Propagating conformational changes over long (and short) distances in proteins

Identification of a chemoreceptor for tricarboxylic acid cycle intermediates: differential chemotactic response towards receptor ligands

We report the identification of McpS as the specific chemoreceptor for 6 tricarboxylic acid (TCA) cycle intermediates and butyrate in Pseudomonas putida. The analysis of the bacterial mutant deficient in mcpS and complementation assays demonstrate that McpS is the only chemoreceptor of TCA cycle intermediates in the strain under study. TCA cycle intermediates are abundantly present in root exudates, and taxis toward these compounds is proposed to facilitate the access to carbon sources. McpS has an unusually large ligand-binding domain (LBD) that is un-annotated in InterPro and is predicted to contain 6 helices. The ligand profile of McpS was determined by isothermal titration calorimetry of purified recombinant LBD (McpS-LBD). McpS recognizes TCA cycle intermediates but does not bind very close structural homologues and derivatives like maleate, aspartate, or tricarballylate. This implies that functional similarity of ligands, such as being part of the same pathway, and not structural similarity is the primary element, which has driven the evolution of receptor specificity. The magnitude of chemotactic responses toward these 7 chemoattractants, as determined by qualitative and quantitative chemotaxis assays, differed largely. Ligands that cause a strong chemotactic response (malate, succinate, and fumarate) were found by differential scanning calorimetry to increase significantly the midpoint of protein unfolding (T(m)) and unfolding enthalpy (DeltaH) of McpS-LBD. Equilibrium sedimentation studies show that malate, the chemoattractant that causes the strongest chemotactic response, stabilizes the dimeric state of McpS-LBD. In this respect clear parallels exist to the Tar receptor and other eukaryotic receptors, which are discussed.

Ingestion-controlling network: what's language got to do with it?

The prevalence of obesity has increased dramatically worldwide, whereas the types of treatment and their efficacy have not substantially changed over the last two decades. Additionally, drugs used to control weight gain could occasionally create untoward effects in cardiovascular functions, as well as in behaviors, memory, sleep, and emotions because the molecular machinery responsible for ingestion control is interconnected with or shared by the above domains. How each group of drugs preserves the privacy of its message in the mutual network is not fully understood. In the present essay, the graph theory approach was used to explore some aspects of molecular signaling as though they were a 'language'. Its emphasis is on 'molecular polysemy', a term that refers to the ability of biomolecules to be used like words in natural languages more than one-way. This has physiological and clinical implications, in particular when planning drug designs with "specially engineered shotgun loads" that target a combination of biomolecules that assure a better therapeutic outcome without causing deficits in connected but patho-physiologically irrelevant bystanders.

Unraveling the allosteric mechanism of serine protease inhibition by an antibody

Recent structural studies have outlined the mechanism of protease inhibition by active site-directed antibodies. However, the molecular basis of allosteric inhibition by antibodies has been elusive. Here we report the 2.35 A resolution structure of the trypsin-like serine protease hepatocyte growth factor activator (HGFA) in complex with the allosteric antibody Ab40, a potent inhibitor of HGFA catalytic activity. The antibody binds at the periphery of the substrate binding cleft and imposes a conformational change on the entire 99-loop (chymotrypsinogen numbering). The altered conformation of the 99-loop is incompatible with substrate binding due to the partial collapse of subsite S2 and the reorganization of subsite S4. Remarkably, a single residue deletion of Ab40 abolished inhibition of HGFA activity, commensurate with the reversal of the 99-loop conformation to its "competent" state. The results define an "allosteric switch" mechanism as the basis of protease inhibition by an allosteric antibody.

In situ D-periodic molecular structure of type II collagen

Collagens are essential components of extracellular matrices in multicellular animals. Fibrillar type II collagen is the most prominent component of articular cartilage and other cartilage-like tissues such as notochord. Its in situ macromolecular and packing structures have not been fully characterized, but an understanding of these attributes may help reveal mechanisms of tissue assembly and degradation (as in osteo- and rheumatoid arthritis). In some tissues such as lamprey notochord, the collagen fibrillar organization is naturally crystalline and may be studied by x-ray diffraction. We used diffraction data from native and derivative notochord tissue samples to solve the axial, D-periodic structure of type II collagen via multiple isomorphous replacement. The electron density maps and heavy atom data revealed the conformation of the nonhelical telopeptides and the overall D-periodic structure of collagen type II in native tissues, data that were further supported by structure prediction and transmission electron microscopy. These results help to explain the observed differences in collagen type I and type II fibrillar architecture and indicate the collagen type II cross-link organization, which is crucial for fibrillogenesis. Transmission electron microscopy data show the close relationship between lamprey and mammalian collagen fibrils, even though the respective larger scale tissue architecture differs.

Millisecond dynamics in the allosteric enzyme imidazole glycerol phosphate synthase (IGPS) from Thermotoga maritima

IGPS is a 51 kDa heterodimeric enzyme comprised of two proteins, HisH and HisF, that catalyze the hydrolysis of glutamine to produce NH(3) in the HisH active site and the cyclization of ammonia with N'-[(5'-phosphoribulosyl)formimino]-5-aminoimidazole-4-carboxamide-ribonucleotide (PRFAR) in HisF to produce imidazole glycerol phosphate (IGP) and 5-aminoimidazole-4-carboxamide ribotide (AICAR). Binding of PRFAR and IGP stimulates glutaminase activity in the HisH enzyme over 5,000 and 100-fold, respectively, despite the active sites being >25 A apart. The details of this long-range protein communication process were investigated by solution NMR spectroscopy and CPMG relaxation dispersion experiments. Formation of the heterodimer enzyme results in a reduction in millisecond motions in HisF that extend throughout the protein. Binding of lGP results in an increase in protein-wide millisecond dynamics evidenced as severe NMR line broadening and elevated R (ex) values. Together, these data demonstrate a grouping of flexible residues that link the HisF active site with the protein interface to which HisH binds and provide a model for the path of communication between the IGPS active sites.

Flexibility in the inducer binding region is crucial for allostery in the Escherichia coli lactose repressor

Lactose repressor protein (LacI) utilizes an allosteric mechanism to regulate transcription in Escherichia coli, and the transition between inducer- and operator-bound states has been simulated by targeted molecular dynamics (TMD). The side chains of amino acids 149 and 193 interact and were predicted by TMD simulation to play a critical role in the early stages of the LacI conformational change. D149 contacts IPTG directly, and variations at this site provide the opportunity to dissect its role in inducer binding and signal transduction. Single mutants at D149 or S193 exhibit a minimal change in operator binding, and alterations in inducer binding parallel changes in operator release, indicating normal allosteric response. The observation that the double mutant D149A/S193A exhibits wild-type properties excludes the requirement for inter-residue hydrogen bond formation in the allosteric response. The double mutant D149C/S193C purified from cell extracts shows decreased sensitivity to inducer binding while retaining wild-type binding affinities and kinetic constants for both operator and inducer. By manipulating cysteine oxidation, we show that the more reduced state of D149C/S193C responds to inducer more like the wild-type protein, whereas the more oxidized state displays diminished inducer sensitivity. These features of D149C/S193C indicate that the novel disulfide bond formed in this mutant impedes the allosteric transition, consistent with the role of this region predicted by TMD simulation. Together, these results establish the requirement for flexibility in the spatial relationship between D149 and S193 rather than a specific D149-S193 interaction in the LacI allosteric response to inducer.

Critical role of substrate conformational change in the proton transfer process catalyzed by 4-oxalocrotonate tautomerase

4-Oxalocrotonate tautomerase enzyme (4-OT) catalyzes the isomerization of 2-oxo-4-hexenedioate to 2-oxo-3-hexenedioate. The chemical process involves two proton transfers, one from a carbon of the substrate to the nitrogen of Pro1 and another from this nitrogen atom to a different carbon of the substrate. In this paper the isomerization has been studied using the combined quantum mechanical and molecular mechanical method with a dual-level treatment of the quantum subsystem employing the MPW1BK density functional as the higher level. Exploration of the potential energy surface shows that the process is stepwise, with a stable intermediate state corresponding to the deprotonated substrate and a protonated proline. The rate constant of the overall process has been evaluated using ensemble-averaged variational transition state theory, including the quantized vibrational motion of a primary zone of active-site atoms and a transmission coefficient based on an ensemble of optimized reaction coordinates to account for recrossing trajectories and optimized multidimensional tunneling. The two proton-transfer steps have similar free energy barriers, but the transition state associated with the first proton transfer is found to be higher in energy. The calculations show that reaction progress is coupled to a conformational change of the substrate, so it is important that the simulation allows this flexibility. The coupled conformational change is promoted by changes in the electron distribution of the substrate that take place as the proton transfers occur.

Allosteric communication occurs via networks of tertiary and quaternary motions in proteins

Allosteric proteins bind an effector molecule at one site resulting in a functional change at a second site. We hypothesize that allosteric communication in proteins relies upon networks of quaternary (collective, rigid-body) and tertiary (residue-residue contact) motions. We argue that cyclic topology of these networks is necessary for allosteric communication. An automated algorithm identifies rigid bodies from the displacement between the inactive and the active structures and constructs "quaternary networks" from these rigid bodies and the substrate and effector ligands. We then integrate quaternary networks with a coarse-grained representation of contact rearrangements to form "global communication networks" (GCNs). The GCN reveals allosteric communication among all substrate and effector sites in 15 of 18 multidomain and multimeric proteins, while tertiary and quaternary networks exhibit such communication in only 4 and 3 of these proteins, respectively. Furthermore, in 7 of the 15 proteins connected by the GCN, 50% or more of the substrate-effector paths via the GCN are "interdependent" paths that do not exist via either the tertiary or the quaternary network. Substrate-effector "pathways" typically are not linear but rather consist of polycyclic networks of rigid bodies and clusters of rearranging residue contacts. These results argue for broad applicability of allosteric communication based on structural changes and demonstrate the utility of the GCN. Global communication networks may inform a variety of experiments on allosteric proteins as well as the design of allostery into non-allosteric proteins.

Allosteric regulation of proteases

Allostery is a basic principle of control of enzymatic activities based on the interaction of a protein or small molecule at a site distinct from an enzyme's active center. Allosteric modulators represent an alternative approach to the design and synthesis of small-molecule activators or inhibitors of proteases and are therefore of wide interest for medicinal chemistry. The structural bases of some proteinaceous and small-molecule allosteric protease regulators have already been elucidated, indicating a general mechanism that might be exploitable for future rational design of small-molecule effectors.

The 'Tyranny of choices' in the ingestion-controlling network

BACKGROUND

Currently used antiobesity remedies offer only a modest weight reduction, and have untoward effects that can complicate treatment efforts. Motivated by the needs of the pharmacotherapy of obesity, the study explored the role of neuropeptide Y, leptin, and corticotrophin-releasing hormone.

METHOD

The study used Ingenuity Pathway Analysis which is a tool for automated discovery and visualization of molecular interactions.

RESULTS

In ingestion-controlling networks, neuropeptide Y, leptin, and corticotrophin-releasing hormone molecules are commonly combined into the units designated as 'maximal motifs'. The analysis of this triad allowed suggesting that maximal motifs are not more than a compendium of admission rules and transmission alternatives of their nodes catalogued in the dataset. Nonetheless, these options seem to endow them with the flexibility needed to respond dynamically as a functional unit to changing internal (metabolic) conditions or environmental challenges.

CONCLUSION

Thus far, each peptide represents a separate target for pharmaceutical interventions (as judged by US patents scanned). The study concludes with predictions regarding designs of 'multitargeted' antiobesity agents since only by hitting a combination of targets can an appropriate therapeutic effect be achieved.

A coupled equilibrium shift mechanism in calmodulin-mediated signal transduction

We used nuclear magnetic resonance data to determine ensembles of conformations representing the structure and dynamics of calmodulin (CaM) in the calcium-bound state (Ca(2+)-CaM) and in the state bound to myosin light chain kinase (CaM-MLCK). These ensembles reveal that the Ca(2+)-CaM state includes a range of structures similar to those present when CaM is bound to MLCK. Detailed analysis of the ensembles demonstrates that correlated motions within the Ca(2+)-CaM state direct the structural fluctuations toward complex-like substates. This phenomenon enables initial ligation of MLCK at the C-terminal domain of CaM and induces a population shift among the substates accessible to the N-terminal domain, thus giving rise to the cooperativity associated with binding. Based on these results and the combination of modern free energy landscape theory with classical allostery models, we suggest that a coupled equilibrium shift mechanism controls the efficient binding of CaM to a wide range of ligands.

Conformational flexibility of a model protein upon immobilization on self-assembled monolayers

The present study reports on the retention of conformational flexibility of a model allosteric protein upon immobilization on self-assembled monolayers (SAMs) on gold. Organothiolated SAMs of different compositions were utilized for adsorptive and covalent attachment of bovine liver glutamate dehydrogenase (GDH), a well-characterized allosteric enzyme. Sensitive fluorimetric assays were developed to determine immobilization capacity, specific activity, and allosteric properties of the immobilized preparations as well as the potential for repeated use and continuous catalytic transformations. The allosteric response of the free and immobilized forms towards ADP, L-leucine and high concentrations of NAD(+), some of the well-known activators for this enzyme, were determined and compared. The enzyme immobilized by adsorption or chemical binding responded similarly to the activators with a greater degree of activation, as compared to the free form. Also loss of activity involving the two immobilization procedures were similar, suggesting that residues essential for catalytic activity or allosteric properties of GDH remained unchanged in the course of chemical modification. A recently established method was used to predict GDH orientation upon immobilization, which was found to explain some of the experimental results presented. The general significance of these observations in connection with retention of native properties of protein structures upon immobilization on SAMs is discussed.

Contact rearrangements form coupled networks from local motions in allosteric proteins

Allosteric proteins bind an effector molecule at one site resulting in a functional change at a second site. We hypothesize that networks of contacts altered, formed, or broken are a significant contributor to allosteric communication in proteins. In this work, we identify which interactions change significantly between the residue-residue contact networks of two allosteric structures, and then organize these changes into graphs. We perform the analysis on 15 pairs of allosteric structures with effector and substrate each present in at least one of the two structures. Most proteins exhibit large, dense regions of contact rearrangement, and the graphs form connected paths between allosteric effector and substrate sites in five of these proteins. In the remaining 10 proteins, large-scale conformational changes such as rigid-body motions are likely required in addition to contact rearrangement networks to account for substrate-effector communication. On average, clusters which contain at least one substrate or effector molecule comprise 20% of the protein. These allosteric graphs are small worlds; that is, they typically have mean shortest path lengths comparable to those of corresponding random graphs and average clustering coefficients enhanced relative to those of random graphs. The networks capture 60-80% of known allostery-perturbing mutants in three proteins, and the metrics degree and closeness are statistically good discriminators of mutant residues from nonmutant residues within the networks in two of these three proteins. For two proteins, coevolving clusters of residues which have been hypothesized to be allosterically important differ from the regions with the most contact rearrangement. Residues and contacts which modulate normal mode fluctuations also often participate in the contact rearrangement networks. In summary, residue-residue contact rearrangement networks provide useful representations of the portions of allosteric pathways resulting from coupled local motions.

Activation of G protein-coupled receptors: beyond two-state models and tertiary conformational changes

Transformation of G protein-coupled receptors (GPCRs) from a quiescent to an active state initiates signal transduction. All GPCRs share a common architecture comprising seven transmembrane-spanning alpha-helices, which accommodates signal propagation from a diverse repertoire of external stimuli across biological membranes to a heterotrimeric G protein. Signal propagation through the transmembrane helices likely involves mechanistic features common to all GPCRs. The structure of the light receptor rhodopsin may serve as a prototype for the transmembrane architecture of GPCRs. Early biochemical, biophysical, and pharmacological studies led to the conceptualization of receptor activation based on the context of two-state equilibrium models and conformational changes in protein structure. More recent studies indicate a need to move beyond these classical paradigms and to consider additional aspects of the molecular character of GPCRs, such as the oligomerization and dynamics of the receptor.

Structural identification of the pathway of long-range communication in an allosteric enzyme

Allostery is a common mechanism of regulation of enzyme activity and specificity, and its signatures are readily identified from functional studies. For many allosteric systems, structural evidence exists of long-range communication among protein domains, but rarely has this communication been traced to a detailed pathway. The thrombin mutant D102N is stabilized in a self-inhibited conformation where access to the active site is occluded by a collapse of the entire 215-219 beta-strand. Binding of a fragment of the protease activated receptor PAR1 to exosite I, 30-A away from the active site region, causes a large conformational change that corrects the position of the 215-219 beta-strand and restores access to the active site. The crystal structure of the thrombin-PAR1 complex, solved at 2.2-A resolution, reveals the details of this long-range allosteric communication in terms of a network of polar interactions.

The intrinsic dynamics and function of nickel-binding regulatory protein: insights from elastic network analysis

Nickel-responsive protein NikR regulates the nickel uptake in nickel-dependent bacteria by interacting with the operator of nikABCDE and subsequently repressing the transcription of NikABCDE, an ABC-type nickel transporter system. The function of NikR and its affinity for the operator DNA is highly conformation-dependent, which has been confirmed by three independent crystallographic studies on NikR proteins from different bacteria. Depending on the intracellular nickel concentration, NikR is able to adopt either the open form or one of the two closed forms (cis and trans) that differ in the domain-domain arrangement. Only the closed cis form is optimal for DNA binding. We examined the low-resolution vibrational spectrum of NikR in each conformational form using the elastic network model and observed large-scale domain-domain vibrations that are closely related to the conformational transitions required for function, particularly the symmetric bending mode and the asymmetric twisting mode. This analysis on the intrinsic dynamics coded in the three-dimensional molecular construct allows us to examine the proposed mechanisms of NikR regulation from the standpoint of protein collective motions. Our findings further support the three-state equilibrium hypothesis proposed by others, and imply that an isolated closed cis form may be dynamically unstable but can be stabilized by DNA binding. However, we also found that the simple C(alpha)-model used in the current analysis is insufficient to capture the impact of nickel binding on the protein dynamics, for which an all-atom model with detailed atom typing is more appropriate.

A potential molecular switch in an alpha-helical coiled coil

The yeast DNA-binding protein GCN4 forms a homo-dimer through a self-complementary coiled-coil interface. In this article, we describe how such coiled-coils might be bistable and, through Molecular Dynamics computations on the GCN4 coiled coil, we show that the coiled coil can indeed switch between the two states by a pathway in which there is a progressive "flipping" of consecutive steps along the interface. We discuss the general implications of potentially bistable coiled-coil interfaces for allosteric signal-transmission mechanisms along homo-dimeric coiled coils and for the packing of helices in globular proteins.

Myosin dynamics on the millisecond time scale

Myosin is a motor protein associating with actin and ATP. It translates along actin filaments against a force by transduction of free energy liberated with ATP hydrolysis. Various myosin crystal structures define time points during ATPase showing the protein undergoes large conformation change during transduction over a cycle with approximately 10 ms periodicity. The protein conformation trajectory between two intermediates in the cycle is surmised by non-equilibrium Monte Carlo simulation utilizing free-energy minimization. The trajectory shows myosin transduction of free energy to mechanical work giving evidence for: (i) a causal relationship between product release and work production in the native isoform that is correctly disrupted in a chemically modified protein, (ii) the molecular basis of ATP-sensitive tryptophan fluorescence enhancement and acrylamide quenching, (iii) an actin-binding site peptide containing the free-energy barrier to ATPase product release defining the rate limiting step and, (iv) a scenario for actin-activation of myosin ATPase.

Using molecular dynamics to map interaction networks in an aminoacyl-tRNA synthetase

Long-range functional communication is a hallmark of many enzymes that display allostery, or action-at-a-distance. Many aminoacyl-tRNA synthetases can be considered allosteric, in that their trinucleotide anticodons bind the enzyme at a site removed from their catalytic domains. Such is the case with E. coli methionyl-tRNA synthase (MetRS), which recognizes its cognate anticodon using a conserved tryptophan residue 50 A away from the site of tRNA aminoacylation. The lack of details regarding how MetRS and tRNA(Met) interact has limited efforts to deconvolute the long-range communication that occurs in this system. We have used molecular dynamics simulations to evaluate the mobility of wild-type MetRS and a Trp-461 variant shown previously by experiment to be deficient in tRNA aminoacylation. The simulations reveal that MetRS has significant mobility, particularly at structural motifs known to be involved in catalysis. Correlated motions are observed between residues in distant structural motifs, including the active site, zinc binding motif, and anticodon binding domain. Both mobility and correlated motions decrease significantly but not uniformly upon substitution at Trp-461. Mobility of some residues is essentially abolished upon removal of Trp-461, despite being tens of Angstroms away from the site of mutation and solvent exposed. This conserved residue does not simply participate in anticodon binding, as demonstrated experimentally, but appears to mediate the protein's distribution of structural ensembles. Finally, simulations of MetRS indicate that the ligand-free protein samples conformations similar to those observed in crystal structures with substrates and substrate analogs bound. Thus, there are low energetic barriers for MetRS to achieve the substrate-bound conformations previously determined by structural methods.

Ligand binding and allostery can emerge simultaneously

A heterotropic allosteric effect involves an effector molecule that is distinct from the substrate or ligand of the protein. How heterotropic allostery originates is an unanswered question. We have previously created several heterotropic allosteric enzymes by recombining the genes for TEM1 beta-lactamase (BLA) and maltose binding protein (MBP) to create BLAs that are positively or negatively regulated by maltose. We show here that one of these engineered enzymes has approximately 10(6) M(-1) affinity for Zn(2+), a property that neither of the parental proteins possesses. Furthermore, Zn(2+) is a negative effector that noncompetitively switches off beta-lactam hydrolysis activity. Mutagenesis experiments indicate that the Zn(2+)-binding site does not involve a histidine or a cysteine, which is atypical of natural Zn(2+)-binding sites. These studies also implicate helices 1 and 12 of the BLA domain in allosteric signal propagation. These results support a model for the evolution of heterotropic allostery in which effector affinity and allosteric signaling emerge simultaneously.

Local motions in a benchmark of allosteric proteins

Allosteric proteins have been studied extensively in the last 40 years, but so far, no systematic analysis of conformational changes between allosteric structures has been carried out. Here, we compile a set of 51 pairs of known inactive and active allosteric protein structures from the Protein Data Bank. We calculate local conformational differences between the two structures of each protein using simple metrics, such as backbone and side-chain Cartesian displacement, and torsion angle change and rearrangement in residue-residue contacts. Thresholds for each metric arise from distributions of motions in two control sets of pairs of protein structures in the same biochemical state. Statistical analysis of motions in allosteric proteins quantifies the magnitude of allosteric effects and reveals simple structural principles about allostery. For example, allosteric proteins exhibit substantial conformational changes comprising about 20% of the residues. In addition, motions in allosteric proteins show strong bias toward weakly constrained regions such as loops and the protein surface. Correlation functions show that motions communicate through protein structures over distances averaging 10-20 residues in sequence space and 10-20 A in Cartesian space. Comparison of motions in the allosteric set and a set of 21 nonallosteric ligand-binding proteins shows that nonallosteric proteins also exhibit bias of motion toward weakly constrained regions and local correlation of motion. However, allosteric proteins exhibit twice as much percent motion on average as nonallosteric proteins with ligand-induced motion. These observations may guide efforts to design flexibility and allostery into proteins.

Reconciling the "old" and "new" views of protein allostery: a molecular simulation study of chemotaxis Y protein (CheY)

A combination of thirty-two 10-ns-scale molecular dynamics simulations were used to explore the coupling between conformational transition and phosphorylation in the bacteria chemotaxis Y protein (CheY), as a simple but representative example of protein allostery. Results from these simulations support an activation mechanism in which the beta4-alpha4 loop, at least partially, gates the isomerization of Tyr106. The roles of phosphorylation and the conserved Thr87 are deemed indirect in that they stabilize the active configuration of the beta4-alpha4 loop. The indirect role of the activation event (phosphorylation) and/or conserved residues in stabilizing, rather than causing, specific conformational transition is likely a feature in many signaling systems. The current analysis of CheY also helps to make clear that neither the "old" (induced fit) nor the "new" (population shift) views for protein allostery are complete, because they emphasize the kinetic (mechanistic) and thermodynamic aspects of allosteric transitions, respectively. In this regard, an issue that warrants further analysis concerns the interplay of concerted collective motion and sequential local structural changes in modulating cooperativity between distant sites in biomolecules.

Synthetic multivalent ligands as probes of signal transduction

Cell-surface receptors acquire information from the extracellular environment and coordinate intracellular responses. Many receptors do not operate as individual entities, but rather as part of dimeric or oligomeric complexes. Coupling the functions of multiple receptors may endow signaling pathways with the sensitivity and malleability required to govern cellular responses. Moreover, multireceptor signaling complexes may provide a means of spatially segregating otherwise degenerate signaling cascades. Understanding the mechanisms, extent, and consequences of receptor co-localization and interreceptor communication is critical; chemical synthesis can provide compounds to address the role of receptor assembly in signal transduction. Multivalent ligands can be generated that possess a variety of sizes, shapes, valencies, orientations, and densities of binding elements. This Review focuses on the use of synthetic multivalent ligands to characterize receptor function.

Is allostery an intrinsic property of all dynamic proteins?

Allostery involves coupling of conformational changes between two widely separated binding sites. The common view holds that allosteric proteins are symmetric oligomers, with each subunit existing in "at least" two conformational states with a different affinity for ligands. Recent observations such as the allosteric behavior of myoglobin, a classical example of a nonallosteric protein, call into question the existing allosteric dogma. Here we argue that all (nonfibrous) proteins are potentially allosteric. Allostery is a consequence of re-distributions of protein conformational ensembles. In a nonallosteric protein, the binding site shape may not show a concerted second-site change and enzyme kinetics may not reflect an allosteric transition. Nevertheless, appropriate ligands, point mutations, or external conditions may facilitate a population shift, leading a presumably nonallosteric protein to behave allosterically. In principle, practically any potential drug binding to the protein surface can alter the conformational redistribution. The question is its effectiveness in the redistribution of the ensemble, affecting the protein binding sites and its function. Here, we review experimental observations validating this view of protein allostery.

Thrombin: a paradigm for enzymes allosterically activated by monovalent cations

Enzymes activated by monovalent cations are abundantly represented in plants and in the animal world. The mechanism, of activation involves formation of a ternary intermediate with the enzyme-substrate complex, or binding of the cation to an allosteric site in the protein. Thrombin is a Na+-activated enzyme with procoagulant, anticoagulant and signaling roles. The binding of Na+ influences allosterically thrombin function and offers a paradigm for regulatory control of protease activity and specificity. Here we review the molecular basis of thrombin allostery as recently emerged from mutagenesis and structural studies. The role of Na+ in blood coagulation and the evolution of serine proteases are also discussed.

Directed evolution of protein switches and their application to the creation of ligand-binding proteins

We describe an iterative approach for creating protein switches involving the in vitro recombination of two nonhomologous genes. We demonstrate this approach by recombining the genes coding for TEM1 beta-lactamase (BLA) and the Escherichia coli maltose binding protein (MBP) to create a family of MBP-BLA hybrids in which maltose is a positive or negative effector of beta-lactam hydrolysis. Some of these MBP-BLA switches were effectively "on-off" in nature, with maltose altering catalytic activity by as much as 600-fold. The ability of these switches to confer an effector-dependent growth/no growth phenotype to E. coli cells was exploited to rapidly identify, from a library of 4 x 10(6) variants, MBP-BLA switch variants that respond to sucrose as the effector. The transplantation of these mutations into wild-type MBP converted MBP into a "sucrose-binding protein," illustrating the switches potential as a tool to rapidly identify ligand-binding proteins.

Engineering allosteric protein switches by domain insertion

Domain insertion is proving to be an effective way to construct hybrid proteins exhibiting switch-like behavior. In this strategy, two existing domains, the first exhibiting a signal recognition function and the second containing the function to be modulated, are fused such that the recognition of the signal by the first domain is transmitted to the second domain, thereby modulating its activity. Recent directed evolution experiments indicate that the structural space comprised of the recombination of unrelated protein domains may be rich in switching behavior, particularly when the circular permutation of domains is also employed. This bodes well for potential basic science, sensing and therapeutic applications of molecular switches.

Biophysical and kinetic characterization of HemAT, an aerotaxis receptor from Bacillus subtilis

HemAT from Bacillus subtilis is a new type of heme protein responsible for sensing oxygen. The structural and functional properties of the full-length HemAT protein, the sensor domain (1-178), and Tyr-70 mutants have been characterized. Kinetic and equilibrium measurements reveal that both full-length HemAT and the sensor domain show two distinct O(2) binding components. The high-affinity component has a K(dissociation) approximately 1-2 microM and a normal O(2) dissociation rate constant, k(O2) = 50-80 s(-1). The low-affinity component has a K(dissociation) approximately 50-100 microM and a large O(2) dissociation rate constant equal to approximately 2000 s(-1). The low n-value and biphasic character of the equilibrium curve indicate that O(2) binding to HemAT involves either independent binding to high- and low-affinity subunits in the dimer or negative cooperativity. Replacement of Tyr-70(B10) with Phe, Leu, or Trp in the sensor domain causes dramatic increases in k(O2) for both the high- and low-affinity components. In contrast, the rates and affinity for CO binding are little affected by loss of the Tyr-70 hydroxyl group. These results suggest highly dynamic behavior for the Tyr-70 side chain and the fraction of the "up" versus "down" conformation is strongly influenced by the nature of the iron-ligand complex. As a result of having both high- and low-affinity components, HemAT can respond to oxygen concentration gradients under both hypoxic (0-10 microM) and aerobic (50-250 microM) conditions, a property which could, in principle, be important for a robust sensing system. The unusual ligand-binding properties of HemAT suggest that asymmetry and apparent negative cooperativity play an important role in the signal transduction pathway.

The importance of being dimeric

Why are there so many dimeric proteins and enzymes? While for heterodimers a functional explanation seems quite reasonable, the case of homodimers is more puzzling. The number of homodimers found in all living organisms is rapidly increasing. A thorough inspection of the structural data from the available literature and stability (measured from denaturation-renaturation experiments) allows one to suggest that homodimers can be divided into three main types according to their mass and the presence of a (relatively) stable monomeric intermediate in the folding-unfolding pathway. Among other explanations, we propose that an essential advantage for a protein being dimeric may be the proper and rapid assembly in the cellular milieu.

Molecular dissection of Na+ binding to thrombin

Na(+) binding near the primary specificity pocket of thrombin promotes the procoagulant, prothrombotic, and signaling functions of the enzyme. The effect is mediated allosterically by a communication between the Na(+) site and regions involved in substrate recognition. Using a panel of 78 Ala mutants of thrombin, we have mapped the allosteric core of residues that are energetically linked to Na(+) binding. These residues are Asp-189, Glu-217, Asp-222, and Tyr-225, all in close proximity to the bound Na(+). Among these residues, Asp-189 shares with Asp-221 the important function of transducing Na(+) binding into enhanced catalytic activity. None of the residues of exosite I, exosite II, or the 60-loop plays a significant role in Na(+) binding and allosteric transduction. X-ray crystal structures of the Na(+)-free (slow) and Na(+)-bound (fast) forms of thrombin, free or bound to the active site inhibitor H-d-Phe-Pro-Arg-chloromethyl-ketone, document the conformational changes induced by Na(+) binding. The slow --> fast transition results in formation of the Arg-187:Asp-222 ion pair, optimal orientation of Asp-189 and Ser-195 for substrate binding, and a significant shift of the side chain of Glu-192 linked to a rearrangement of the network of water molecules that connect the bound Na(+) to Ser-195 in the active site. The changes in the water network and the allosteric core explain the thermodynamic signatures linked to Na(+) binding and the mechanism of thrombin activation by Na(+). The role of the water network uncovered in this study establishes a new paradigm for the allosteric regulation of thrombin and other Na(+)-activated enzymes involved in blood coagulation and the immune response.

Reaction-path energetics and kinetics of the hydride transfer reaction catalyzed by dihydrofolate reductase

We have studied the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) and the coenzyme nicotinamide adenine dinucleotide phosphate (NADPH); the substrate is 5-protonated 7,8-dihydrofolate, and the product is tetrahydrofolate. The potential energy surface is modeled by a combined quantum mechanical-molecular mechanical (QM/MM) method employing Austin model 1 (AM1) and a simple valence bond potential for 69 QM atoms and employing the CHARMM22 and TIP3P molecular mechanics force fields for the other 21 399 atoms; the QM and MM regions are joined by two boundary atoms treated by the generalized hybrid orbital (GHO) method. All simulations are carried out using periodic boundary conditions at neutral pH and 298 K. In stage 1, a reaction coordinate is defined as the difference between the breaking and forming bond distances to the hydride ion, and a quasithermodynamic free energy profile is calculated along this reaction coordinate. This calculation includes quantization effects on bound vibrations but not on the reaction coordinate, and it is used to locate the variational transition state that defines a transition state ensemble. Then, the key interactions at the reactant, variational transition state, and product are analyzed in terms of both bond distances and electrostatic energies. The results of both analyses support the conclusion derived from previous mutational studies that the M20 loop of DHFR makes an important contribution to the electrostatic stabilization of the hydride transfer transition state. Third, transmission coefficients (including recrossing factors and multidimensional tunneling) are calculated and averaged over the transition state ensemble. These averaged transmission coefficients, combined with the quasithermodynamic free energy profile determined in stage 1, allow us to calculate rate constants, phenomenological free energies of activation, and primary and secondary kinetic isotope effects. A primary kinetic isotope effect (KIE) of 2.8 has been obtained, in good agreement with the experimentally determined value of 3.0 and with the value 3.2 calculated previously. The primary KIE is mainly a consequence of the quantization of bound vibrations. In contrast, the secondary KIE, with a value of 1.13, is almost entirely due to dynamical effects on the reaction coordinate, especially tunneling.

Identification of substrate binding sites in enzymes by computational solvent mapping

Enzyme structures determined in organic solvents show that most organic molecules cluster in the active site, delineating the binding pocket. We have developed algorithms to perform solvent mapping computationally, rather than experimentally, by placing molecular probes (small molecules or functional groups) on a protein surface, and finding the regions with the most favorable binding free energy. The method then finds the consensus site that binds the highest number of different probes. The probe-protein interactions at this site are compared to the intermolecular interactions seen in the known complexes of the enzyme with various ligands (substrate analogs, products, and inhibitors). We have mapped thermolysin, for which experimental mapping results are also available, and six further enzymes that have no experimental mapping data, but whose binding sites are well characterized. With the exception of haloalkane dehalogenase, which binds very small substrates in a narrow channel, the consensus site found by the mapping is always a major subsite of the substrate-binding site. Furthermore, the probes at this location form hydrogen bonds and non-bonded interactions with the same residues that interact with the specific ligands of the enzyme. Thus, once the structure of an enzyme is known, computational solvent mapping can provide detailed and reliable information on its substrate-binding site. Calculations on ligand-bound and apo structures of enzymes show that the mapping results are not very sensitive to moderate variations in the protein coordinates.

G541R within the 4070A TM protein regulates fusion in murine leukemia viruses

The mutation G541R within the ectodomain of TM was isolated in three independent chimeric enveloped murine leukemia virus (MuLV) viral populations originally impaired in viral passage and in wild-type 4070A. Isolation of G541R in multiple populations suggested it played a critical role in viral envelope function. Using a viral vector system, the observed effects of the G541R mutation within MuLV envelope proteins were pleiotropic and included effects on the regulation of SU-TM interactions and membrane fusion. G541R suppresses enhanced cell-cell fusion events attributable to the absence of the R-peptide yet does not adversely affect virus titers. The ability to suppress cell-cell fusion is dependent on the presence of the C terminus of the amphotropic 4070A SU protein. Within the wild-type 4070A envelope background, the mutation results in a decreased level of Env at the cell surface that is mirrored in the virion. The TM mutation alters recognition of the SU C terminus by a monoclonal antibody, suggestive of an altered conformation. The presence of G541R allowed the virus to achieve a balance between cytopathogenicity and replication and restored productive viral entry.

Structure of the oxygen sensor in Bacillus subtilis: signal transduction of chemotaxis by control of symmetry

Much is now known about chemotaxis signaling transduction for Escherichia coli and Salmonella typhimurium. The mechanism of chemotaxis of Bacillus subtilis is, in a sense, reversed. Attractant binding strengthens the activity of histidine kinase in B. subtilis, instead of an inhibition reaction. The HemAT from B. subtilis can detect oxygen and transmit the signal to regulatory proteins that control the direction of flagella rotation. We have determined the crystal structures of the HemAT sensor domain in liganded and unliganded forms at 2.15 A and 2.7 A resolution, respectively. The liganded structure reveals a highly symmetrical organization. Tyrosine70 shows distinct conformational changes on one subunit when ligands are removed. Our study suggests that disruption of the symmetry of HemAT plays an important role in initiating the chemotaxis signaling transduction cascade.

Crystal structure of the bromide-bound D85S mutant of bacteriorhodopsin: principles of ion pumping

We report the crystal structure of a bromide-bound form of the D85S mutant of bacteriorhodopsin, bR(D85S), a protein that uses light energy rather than ATP to pump halide ions across the cell membrane. Comparison of the structure of the halide-bound and halide-free states reveals that both displacements of individual side-chain positions and concerted helical movements occur on the extracellular side of the protein. Analysis of these structural changes reveals how this ion pump first facilitates ion uptake deep within the cell membrane and then prevents the backward escape of ions later in the pumping cycle. Together with the information provided by structures of intermediate states in the bacteriorhodopsin photocycle, this study also suggests the overall design principles that are necessary for ion pumping.

Structural basis of multiple drug-binding capacity of the AcrB multidrug efflux pump

Multidrug efflux pumps cause serious problems in cancer chemotherapy and treatment of bacterial infections. Yet high-resolution structures of ligand transporter complexes have previously been unavailable. We obtained x-ray crystallographic structures of the trimeric AcrB pump from Escherichia coli with four structurally diverse ligands. The structures show that three molecules of ligands bind simultaneously to the extremely large central cavity of 5000 cubic angstroms, primarily by hydrophobic, aromatic stacking and van der Waals interactions. Each ligand uses a slightly different subset of AcrB residues for binding. The bound ligand molecules often interact with each other, stabilizing the binding.

Engineering and design of ligand-induced conformational change in proteins

The ability to manipulate ligand-induced conformational change, although representing a major challenge to the protein engineer, is an essential end point in efforts to produce novel functional proteins for biotechnology and therapeutic applications. Progress towards this goal requires determining not only what factors control the fold and stability of a protein, but also how ligand binding alters the complex conformational/energetic landscape. Important strides are being made on several fronts, including understanding the origin of long-range effects and allosteric structural mechanisms, using both experimental and theoretical approaches.

Drug interactions at GABA(A) receptors

Neurotransmitter receptor systems have been the focus of intensive pharmacological research for more than 20 years for basic and applied scientific reasons, but only recently has there been a better understanding of their key features. One of these systems includes the type A receptor for the gamma-aminobutyric acid (GABA), which forms an integral anion channel from a pentameric subunit assembly and mediates most of the fast inhibitory neurotransmission in the adult vertebrate central nervous system. Up to now, depending on the definition, 16-19 mammalian subunits have been cloned and localized on different genes. Their assembly into proteins in a poorly defined stoichiometry forms the basis of functional and pharmacological GABA(A) receptor diversity, i.e. the receptor subtypes. The latter has been well documented in autoradiographic studies using ligands that label some of the receptors' various binding sites, corroborated by recombinant expression studies using the same tools. Significantly less heterogeneity has been found at the physiological level in native receptors, where the subunit combinations have been difficult to dissect. This review focuses on the characteristics, use and usefulness of various ligands and their binding sites to probe GABA(A) receptor properties and to gain insight into the biological function from fish to man and into evolutionary conserved GABA(A) receptor heterogeneity. We also summarize the properties of the novel mouse models created for the study of various brain functions and review the state-of-the-art imaging of brain GABA(A) receptors in various human neuropsychiatric conditions. The data indicate that the present ligands are only partly satisfactory tools and further ligands with subtype-selective properties are needed for imaging purposes and for confirming the behavioral and functional results of the studies presently carried out in gene-targeted mice with other species, including man.