Structure and dynamics of a molten globular enzyme

Conformational fluctuations of UreG, an intrinsically disordered enzyme

UreG proteins are small GTP binding (G) proteins that catalyze the hydrolysis of GTP necessary for the maturation of urease, a virulence factor in bacterial pathogenesis. UreG proteins are the first documented cases of intrinsically disordered enzymes. The comprehension of the dynamics of folding-unfolding events occurring in this protein could shed light on the enzymatic mechanism of UreG. Here, we used the recently developed replica exchange with solute tempering (REST2) computational methodology to explore the conformational space of UreG from Helicobacter pylori (HpUreG) and to identify its structural fluctuations. The same simulation and analysis protocol has been applied to HypB from Methanocaldococcus jannaschii (MjHypB), which is closely related to UreG in both sequence and function, even though it is not intrinsically disordered. A comparison of the two systems reveals that both HpUreG and MjHypB feature a substantial rigidity of the protein regions involved in catalysis, justifying its residual catalytic activity. On the other hand, HpUreG tends to unfold more than MjHypB in portions involved in protein-protein interactions with metallochaperones necessary for the formation of multiprotein complexes known to be involved in urease activation.

Protein conformational disorder and enzyme catalysis

Though lacking a well-defined three-dimensional structure, intrinsically unstructured proteins are ubiquitous in nature. These molecules play crucial roles in many cellular processes, especially signaling and regulation. Surprisingly, even enzyme catalysis can tolerate substantial disorder. This observation contravenes conventional wisdom but is relevant to an understanding of how protein dynamics modulates enzyme function. This chapter reviews properties and characteristics of disordered proteins, emphasizing examples of enzymes that lack defined structures, and considers implications of structural disorder for catalytic efficiency and evolution.

Packing in molten globules and native states

Close packing of hydrophobic residues in the protein interior is an important determinant of protein stability. Cavities introduced by large to small substitutions are known to destabilize proteins. Conversely, native states of proteins and protein fragments can be stabilized by filling in existing cavities. Molten globules (MGs) were initially used to describe a state of protein which has well-defined secondary structure but little or no tertiary packing. Subsequent studies have shown that MGs do have some degree of native-like topology and specific packing. Wet molten globules (WMGs) with hydrated cores and considerably decreased packing relative to the native state have been studied extensively. Recently there has been renewed interest in identification and characterization of dry molten globules (DMGs). These are slightly expanded forms of the native state which show increased conformational flexibility, native-like main-chain hydrogen bonding and dry interiors. The generality of occurrence of DMGs during protein unfolding and the extent and nature of packing in DMGs remain to be elucidated. Packing interactions in native proteins and MGs can be probed through mutations. Next generation sequencing technologies make it possible to determine relative populations of mutants in a large pool. When this is coupled to phenotypic screens or cell-surface display, it becomes possible to rapidly examine large panels of single-site or multi-site mutants. From such studies, residue specific contributions to protein stability and function can be estimated in a highly parallelized fashion. This complements conventional biophysical methods for characterization of packing in native states and molten globules.

Increased phospholipase A2 activity with phosphorylation of peroxiredoxin 6 requires a conformational change in the protein

We have shown previously and confirmed in this study that the phospholipase A(2) (PLA(2)) activity of peroxiredoxin 6 (Prdx6) is markedly increased by phosphorylation. This report evaluates the conformation and thermodynamic stability of Prdx6 protein after phosphorylation to understand the physical basis for increased activity. Phosphorylation resulted in decreased negative far-UV CD, strengthened ANS binding, and a lack of rigid tertiary structure, compatible with a change in conformation to that of a molten globule. The ΔG°(D) was 3.3 ± 0.3 kcal mol(-1) for Prdx6 and 1.7 ± 0.7 kcal mol(-1) for pPrdx6, suggesting that phosphorylation destabilizes the protein. Phosphorylation of Prdx6 changed the conformation of the N-terminal domain exposing Trp 33, as determined by tryptophan fluorescence and NaI fluorescence quenching. The kinetics of interaction of proteins with unilamellar liposomes (50:25:15:10 DPPC:egg PC:cholesterol:PG molar ratio) were evaluated with tryptophan fluorescence. pPrdx6 bound to liposomes with a higher affinity (K(d) = 5.6 ± 1.2 μM) than Prdx6 (K(d) = 24.9 ± 4.5 μM). By isothermal titration calorimetry, pPrdx6 bound to liposomes with a large exothermic heat loss (ΔH = -31.49 ± 0.22 kcal mol(-1)). Correlating our conformational studies with the published crystal structure of oxidized Prdx6 suggests that phosphorylation results in exposure of hydrophobic residues, thereby providing accessibility to the sites for liposome binding. Because binding of the enzyme to the phospholipid substrate interface is a requirement for PLA(2) activity, these results indicate that a change in the conformation of Prdx6 upon its phosphorylation is the basis for enhancement of PLA(2) enzymatic activity.

Structural instability tuning as a regulatory mechanism in protein-protein interactions

Protein-protein interactions mediate a vast number of cellular processes. Here, we present a regulatory mechanism in protein-protein interactions mediated by finely tuned structural instability and coupled with molecular mimicry. We show that a set of type III secretion (TTS) autoinhibited homodimeric chaperones adopt a molten globule-like state that transiently exposes the substrate binding site as a means to become rapidly poised for binding to their cognate protein substrates. Packing defects at the homodimeric interface stimulate binding, whereas correction of these defects results in less labile chaperones that give rise to nonfunctional biological systems. The protein substrates use structural mimicry to offset the weak spots in the chaperones and to counteract their autoinhibitory conformation. This regulatory mechanism of protein activity is evolutionarily conserved among several TSS systems and presents a lucid example of functional advantage conferred upon a biological system by finely tuned structural instability.

Insights in the (un)structural organization of Bacillus pasteurii UreG, an intrinsically disordered GTPase enzyme

In the past, enzymatic activity has always been expected to be dependent on overall protein rigidity, necessary for substrate recognition and optimal orientation. However, increasing evidence is now accumulating, revealing that some proteins characterized by intrinsic disorder are actually able to perform catalysis. Among them, the only known natural intrinsically disordered enzyme is UreG, a GTPase that, in plants and bacteria, is involved in the protein interaction network leading to Ni(2+) ions delivery into the active site of urease. In this paper, we report a detailed analysis of the unfolding behaviour of UreG from Bacillus pasteurii (BpUreG), following its thermal and chemical denaturation with a combination of fluorescence spectroscopy, calorimetry, CD and NMR. The results demonstrate that BpUreG exists as an ensemble of inter-converting conformations, whose degrees of secondary structure depend on temperature and denaturant concentration. In particular, three major types of conformational ensembles with different degrees of residual structure were identified, with major structural characteristics resembling those of a molten globule (low temperature, absence of denaturant), pre-molten globule (high temperature, absence or presence of denaturant) and random coil (low temperature, presence of denaturant). Transitions among these ensembles of conformational states occur non-cooperatively although reversibly, with a gradual loss or acquisition of residual structure depending on the conditions. A possible role of disorder in the biological function of UreG is envisaged and discussed.

An N-terminal protein degradation tag enables robust selection of highly active enzymes

Degradation tags are short peptide sequences that target proteins for destruction by housekeeping proteases. We previously utilized the C-terminal SsrA tag in directed evolution experiments to decrease the intracellular lifetime of a growth-limiting enzyme and thereby facilitate selection of highly active variants. In this study, we examine the N-terminal RepA tag as an alternative degradation signal for laboratory evolution. Although RepA proved to be less effective than SsrA at lowering protein concentrations in the cell, its N-terminal location dramatically reduced the occurrence of truncation and frameshift artifacts in selection experiments. We exploited this improvement to evolve a topologically redesigned chorismate mutase that is intrinsically disordered but already highly active for the conversion of chorismate to prephenate. After three rounds of mutagenesis and high-stringency selection, a robust and more nativelike variant was obtained that exhibited a catalytic efficiency (k(cat)/K(M) = 84000 M(-1) s(-1)) comparable to that of a natural dimeric chorismate mutase. Because of concomitant increases in catalyst yield, the level of intracellular prephenate production increased approximately 30-fold overall over the course of evolution.

Pericyclic reactions catalyzed by chorismate-utilizing enzymes

One of the fundamental questions of enzymology is how catalytic power is derived. This review focuses on recent developments in the structure--function relationships of chorismate-utilizing enzymes involved in siderophore biosynthesis to provide insight into the biocatalysis of pericyclic reactions. Specifically, salicylate synthesis by the two-enzyme pathway in Pseudomonas aeruginosa is examined. The isochorismate-pyruvate lyase is discussed in the context of its homologues, the chorismate mutases, and the isochorismate synthase is compared to its homologues in the MST family (menaquinone, siderophore, or tryptophan biosynthesis) of enzymes. The tentative conclusion is that the activities observed cannot be reconciled by inspection of the active site participants alone. Instead, individual activities must arise from unique dynamic properties of each enzyme that are tuned to promote specific chemistries.

Impact of oxidation on protein therapeutics: conformational dynamics of intact and oxidized acid-β-glucocerebrosidase at near-physiological pH

The solution dynamics of an enzyme acid-β-glucocerebrosidase (GCase) probed at a physiologically relevant (lysosomal) pH by hydrogen/deuterium exchange mass spectrometry (HDX-MS) reveals very uneven distribution of backbone amide protection across the polypeptide chain. Highly mobile segments are observed even within the catalytic cavity alongside highly protective segments, highlighting the importance of the balance between conformational stability and flexibility for enzymatic activity. Forced oxidation of GCase that resulted in a 40-60% reduction in in vitro biological activity affects the stability of some key structural elements within the catalytic site. These changes in dynamics occur on a longer time scale that is irrelevant for catalysis, effectively ruling out loss of structure in the catalytic site as a major factor contributing to the reduction of the catalytic activity. Oxidation also leads to noticeable destabilization of conformation in remote protein segments on a much larger scale, which is likely to increase the aggregation propensity of GCase and affect its bioavailability. Therefore, it appears that oxidation exerts its negative impact on the biological activity of GCase indirectly, primarily through accelerated aggregation and impaired trafficking.

Biology of amyloid: structure, function, and regulation

Amyloids are highly ordered cross-β sheet protein aggregates associated with many diseases including Alzheimer's disease, but also with biological functions such as hormone storage. The cross-β sheet entity comprising an indefinitely repeating intermolecular β sheet motif is unique among protein folds. It grows by recruitment of the corresponding amyloid protein, while its repetitiveness can translate what would be a nonspecific activity as monomer into a potent one through cooperativity. Furthermore, the one-dimensional crystal-like repeat in the amyloid provides a structural framework for polymorphisms. This review summarizes the recent high-resolution structural studies of amyloid fibrils in light of their biological activities. We discuss how the unique properties of amyloids gives rise to many activities and further speculate about currently undocumented biological roles for the amyloid entity. In particular, we propose that amyloids could have existed in a prebiotic world, and may have been the first functional protein fold in living cells.

Test MM-PB/SA on true conformational ensembles of protein-ligand complexes

The molecular mechanics Poisson-Boltzmann surface area (MM-PB/SA) method has been popular for computing protein-ligand binding free energies in recent years. All previous evaluations of the MM-PB/SA method are based upon computer-generated conformational ensembles, which may be affected by the defective computational methods used for preparing these conformational ensembles. In an attempt to reach more convincing conclusions, we have evaluated the MM-PB/SA method on a set of 24 diverse protein-ligand complexes, each of which has a set of conformations derived from NMR spectroscopy. Our results indicate that both MM-PB/SA and molecular mechanics generalized Born surface area (MM-GB/SA) are able to produce a modest correlation between their results and the experimentally measured binding free energies on our test set. In particular, both MM-PB/SA and MM-GB/SA produced better results by using a representative structure (R = 0.72-0.79) rather than averaging over the conformational ensemble of each given complex (R = 0.61-0.74). A head-to-head comparison with four selected scoring functions (X-Score, PLP, ChemScore, and DrugScore) on the same test set reveals that MM-PB/SA and MM-GB/SA results are marginally better than those produced by scoring funcitons, supporting the value of the MM-PB/SA method. Nevertheless, scoring functions are still more cost-effective options, especially for high-throughput tasks.

Enzymatic activity in disordered states of proteins

Although disordered proteins are able to carry out a variety of different functions, particularly those involved in signalling and regulation, they have been observed to perform catalysis only in a small number of cases. The presence of structural disorder is indeed expected to be poorly compatible with enzymatic catalysis, which requires a well-organised environment in the active site of the enzyme in order to facilitate the formation of the transition state of the chemical reaction to be catalysed. Despite this stringent requirement, current evidence suggests that certain partially disordered proteins could be catalytically active by becoming structured in the regions of their active sites, even if their overall states retain a significant degree of conformational heterogeneity. This type of mechanism, however, does not appear to be not very common, perhaps because the time required to the conformational search within a disordered state to establish a catalytic environment in the presence of the substrate should not be longer than the overall turnover time required for optimal function. In addition, the catalytic environment should be maintained for long enough despite the structural fluctuations to enable the catalytic reaction to take place. As some partially unstructured proteins have been reported to be capable of overcoming these severe limitations and act as enzymes, their study can increase our general understanding of the mechanism of enzymatic catalysis, as well as extend our ability to control the range of functions that can be performed by disordered proteins.

Protein functional landscapes, dynamics, allostery: a tortuous path towards a universal theoretical framework

Energy landscape theories have provided a common ground for understanding the protein folding problem, which once seemed to be overwhelmingly complicated. At the same time, the native state was found to be an ensemble of interconverting states with frustration playing a more important role compared to the folding problem. The landscape of the folded protein – the native landscape – is glassier than the folding landscape; hence, a general description analogous to the folding theories is difficult to achieve. On the other hand, the native basin phase volume is much smaller, allowing a protein to fully sample its native energy landscape on the biological timescales. Current computational resources may also be used to perform this sampling for smaller proteins, to build a ‘topographical map’ of the native landscape that can be used for subsequent analysis. Several major approaches to representing this topographical map are highlighted in this review, including the construction of kinetic networks, hierarchical trees and free energy surfaces with subsequent structural and kinetic analyses. In this review, we extensively discuss the important question of choosing proper collective coordinates characterizing functional motions. In many cases, the substates on the native energy landscape, which represent different functional states, can be used to obtain variables that are well suited for building free energy surfaces and analyzing the protein's functional dynamics. Normal mode analysis can provide such variables in cases where functional motions are dictated by the molecule's architecture. Principal component analysis is a more expensive way of inferring the essential variables from the protein's motions, one that requires a long molecular dynamics simulation. Finally, the two popular models for the allosteric switching mechanism, ‘preexisting equilibrium’ and ‘induced fit’, are interpreted within the energy landscape paradigm as extreme points of a continuum of transition mechanisms. Some experimental evidence illustrating each of these two models, as well as intermediate mechanisms, is presented and discussed.

Design, selection, and characterization of a split chorismate mutase

Split proteins are versatile tools for detecting protein-protein interactions and studying protein folding. Here, we report a new, particularly small split enzyme, engineered from a thermostable chorismate mutase (CM). Upon dissecting the helical-bundle CM from Methanococcus jannaschii into a short N-terminal helix and a 3-helix segment and attaching an antiparallel leucine zipper dimerization domain to the individual fragments, we obtained a weakly active heterodimeric mutase. Using combinatorial mutagenesis and in vivo selection, we optimized the short linker sequences connecting the leucine zipper to the enzyme domain. One of the selected CMs was characterized in detail. It spontaneously assembles from the separately inactive fragments and exhibits wild-type like CM activity. Owing to the availability of a well characterized selection system, the simple 4-helix bundle topology, and the small size of the N-terminal helix, the heterodimeric CM could be a valuable scaffold for enzyme engineering efforts and as a split sensor for specifically oriented protein-protein interactions.

Protein dynamics and conformational disorder in molecular recognition

Recognition requires protein flexibility because it facilitates conformational rearrangements and induced-fit mechanisms upon target binding. Intrinsic disorder is an extreme on the continuous spectrum of possible protein dynamics and its role in recognition may seem counterintuitive. However, conformational disorder is widely found in many eukaryotic regulatory proteins involved in processes such as signal transduction and transcription. Disordered protein regions may in fact confer advantages over folded proteins in binding. Rapidly interconverting and diverse conformers may create mean electrostatic fields instead of presenting discrete charges. The resultant "polyelectrostatic" interactions allow for the utilization of post-translational modifications as a means to change the net charge and thereby modify the electrostatic interaction of a disordered region. Plasticity of disordered protein states enables steric advantages over folded proteins and allows for unique binding configurations. Disorder may also have evolutionary advantages, as it facilitates alternative splicing, domain shuffling and protein modularity. As proteins exist in a continuous spectrum of disorder, so do their complexes. Indeed, disordered regions in complexes may control the degree of motion between domains, mask binding sites, be targets of post-translational modifications, permit overlapping binding motifs, and enable transient binding of different binding partners, making them excellent candidates for signal integrators and explaining their prevalence in eukaryotic signaling pathways. "Dynamic" complexes arise if more than two transient protein interfaces are involved in complex formation of two binding partners in a dynamic equilibrium. "Disordered" complexes, in contrast, do not involve significant ordering of interacting protein segments but rely exclusively on transient contacts. The nature of these interactions is not well understood yet but advancements in the structural characterization of disordered states will help us gain insights into their function and their implications for health and disease.

At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis?

Enzymes play a key role in almost all biological processes, accelerating a variety of metabolic reactions as well as controlling energy transduction, the transcription, and translation of genetic information, and signaling. They possess the remarkable capacity to accelerate reactions by many orders of magnitude compared to their uncatalyzed counterparts, making feasible crucial processes that would otherwise not occur on biologically relevant timescales. Thus, there is broad interest in understanding the catalytic power of enzymes on a molecular level. Several proposals have been put forward to try to explain this phenomenon, and one that has rapidly gained momentum in recent years is the idea that enzyme dynamics somehow contributes to catalysis. This review examines the dynamical proposal in a critical way, considering basically all reasonable definitions, including (but not limited to) such proposed effects as "coupling between conformational and chemical motions," "landscape searches" and "entropy funnels." It is shown that none of these proposed effects have been experimentally demonstrated to contribute to catalysis, nor are they supported by consistent theoretical studies. On the other hand, it is clarified that careful simulation studies have excluded most (if not all) dynamical proposals. This review places significant emphasis on clarifying the role of logical definitions of different catalytic proposals, and on the need for a clear formulation in terms of the assumed potential surface and reaction coordinate. Finally, it is pointed out that electrostatic preorganization actually accounts for the observed catalytic effects of enzymes, through the corresponding changes in the activation free energies.

Interactions between PTB RRMs induce slow motions and increase RNA binding affinity

Polypyrimidine tract binding protein (PTB) participates in a variety of functions in eukaryotic cells, including alternative splicing, mRNA stabilization, and internal ribosomal entry site-mediated translation initiation. Its mechanism of RNA recognition is determined in part by the novel geometry of its two C-terminal RNA recognition motifs (RRM3 and RRM4), which interact with each other to form a stable complex (PTB1:34). This complex itself is unusual among RRMs, suggesting that it performs a specific function for the protein. In order to understand the advantage it provides to PTB, the fundamental properties of PTB1:34 are examined here as a comparative study of the complex and its two constituent RRMs. Both RRM3 and RRM4 adopt folded structures that NMR data show to be similar to their structure in PRB1:34. The RNA binding properties of the domains differ dramatically. The affinity of each separate RRM for polypyrimidine tracts is far weaker than that of PTB1:34, and simply mixing the two RRMs does not create an equivalent binding platform. (15)N NMR relaxation experiments show that PTB1:34 has slow, microsecond motions throughout both RRMs including the interdomain linker. This is in contrast to the individual domains, RRM3 and RRM4, where only a few backbone amides are flexible on this time scale. The slow backbone dynamics of PTB1:34, induced by packing of RRM3 and RRM4, could be essential for high-affinity binding to a flexible polypyrimidine tract RNA and also provide entropic compensation for its own formation.

The EVB as a quantitative tool for formulating simulations and analyzing biological and chemical reactions

Recent years have seen dramatic improvements in computer power, allowing ever more challenging problems to be approached. In light of this, it is imperative to have a quantitative model for examining chemical reactivity, both in the condensed phase and in solution, as well as to accurately quantify physical organic chemistry (particularly as experimental approaches can often be inconclusive). Similarly, computational approaches allow for great progress in studying enzyme catalysis, as they allow for the separation of the relevant energy contributions to catalysis. Due to the complexity of the problems that need addressing, there is a need for an approach that can combine reliability with an ability to capture complex systems in order to resolve long-standing controversies in a unique way. Herein, we will demonstrate that the empirical valence bond (EVB) approach provides a powerful way to connect the classical concepts of physical organic chemistry to the actual energies of enzymatic reactions by means of computation. Additionally, we will discuss the proliferation of this approach, as well as attempts to capture its basic chemistry and repackage it under different names. We believe that the EVB approach is the most powerful tool that is currently available for studies of chemical processes in the condensed phase in general and enzymes in particular, particularly when trying to explore the different proposals about the origin of the catalytic power of enzymes.

Peptide-dominated membranes preceding the genetic takeover by RNA: latest thinking on a classic controversy

It is commonly presumed that abiotic membranes were colonized by proteins later on. Yet, hydrophobic peptides could have formed primordial protein-dominated membranes on their own. In a metabolism-first context, "autocatalytically closed" sets of statistical peptides could organize a self-maintaining protometabolism, assisted by an unfolding set of ribotide-related cofactors. Pairwise complementary ribotide cofactors may have formed docking guides for stochastic peptide formation, before replicating RNA emerged from this subset. Tidally recurring wet-drying cycles and an early onset of photosynthetic activities are considered most likely to meet the thermodynamic requirements. Conceivably, the earliest peptide-dominated vesicles were engaged in light harvesting, together with isoprenoid-tethered pigments, rather than providing an external boundary. Early on, the bulk of prebiotic organic matter can have formed a contiguous layer covering the mineral sediment, held in place by colloidal coherence of a hydrogel matrix. This unconventional scenario assumes a late onset of cellular individualization - perhaps from within, resembling endosporogenesis.

Deconstructing the native state: energy landscapes, function, and dynamics of globular proteins

Proteins are highly complex molecules with features exquisitely selected by nature to carry out essential biological functions. Physical chemistry and polymer physics provide us with the tools needed to make sense of this complexity. Upon translation, many proteins fold to a thermodynamically stable form known as the native state. The native state is not static, but consists of a hierarchy of conformations, that are continuously explored through dynamics. In this review we provide a brief introduction to some of the core concepts required in the discussion of the protein native dynamics using energy landscapes ideas. We first discuss recent works which have challenged the structure-function paradigm by demonstrating function in disordered proteins. Next we examine the hierarchical organization in the energy landscapes using atomistic molecular dynamics simulations and principal component analysis. In particular, the role of direct and water-mediated contacts in sculpting the landscape is elaborated. Another approach to studying the native state ensemble is based on choosing high-resolution order parameters for computing one- or two-dimensional free energy surfaces. We demonstrate that 2D free energy surfaces provide rich thermodynamic and kinetic information about the native state ensemble. Brownian dynamics simulations on such a surface indicate that protein conformational dynamics is weakly activated. Finally, we briefly discuss implicit and coarse-grained protein models and emphasize the solvent role in determining native state structure and dynamics.

Structure-function analyses of isochorismate-pyruvate lyase from Pseudomonas aeruginosa suggest differing catalytic mechanisms for the two pericyclic reactions of this bifunctional enzyme

The isochorismate-pyruvate lyase from Pseudomonas aeruginosa (PchB) catalyzes two pericyclic reactions in a single active site. PchB physiologically produces salicylate and pyruvate from isochorismate for ultimate incorporation of the salicylate into the siderophore pyochelin. PchB also produces prephenate from chorismate, most likely due to structural homology to the Escherchia coli chorismate mutase. The molecular basis of catalysis among enzymatic pericyclic reactions is a matter of debate, one view holding that catalysis may be derived from electrostatic transition state stabilization and the opposing view that catalysis is derived from the generation of a reactive substrate conformation. Mutant forms of PchB were generated by site-directed mutagenesis at the site (K42) hypothesized to be key for electrostatic transition state stabilization (K42A, K42Q, K42E, and K42H). The loop containing K42 is mobile, and a mutant to slow loop dynamics was also designed (A43P). Finally, a previously characterized mutation (I87T) was also produced. Circular dichroism was used to assess the overall effect on secondary structure as a result of the mutations, and X-ray crystallographic structures are reported for K42A with salicylate and pyruvate bound and for apo-I87T. The data illustrate that the active site architecture is maintained in K42A-PchB, which indicates that differences in activity are not caused by secondary structural changes or by differences in active site loop conformation but rather by the chemical nature of this key residue. In contrast, the I87T structure demonstrates considerable mobility, suggesting that loop dynamics and conformational plasticity may be important for efficient catalysis. Finally, the mutational effects on k(cat) provide evidence that the two activities of PchB are not covariant and that a single hypothesis may not provide a sufficient explanation for catalysis.

Structure and inhibition of the SARS coronavirus envelope protein ion channel

The envelope (E) protein from coronaviruses is a small polypeptide that contains at least one α-helical transmembrane domain. Absence, or inactivation, of E protein results in attenuated viruses, due to alterations in either virion morphology or tropism. Apart from its morphogenetic properties, protein E has been reported to have membrane permeabilizing activity. Further, the drug hexamethylene amiloride (HMA), but not amiloride, inhibited in vitro ion channel activity of some synthetic coronavirus E proteins, and also viral replication. We have previously shown for the coronavirus species responsible for severe acute respiratory syndrome (SARS-CoV) that the transmembrane domain of E protein (ETM) forms pentameric α-helical bundles that are likely responsible for the observed channel activity. Herein, using solution NMR in dodecylphosphatidylcholine micelles and energy minimization, we have obtained a model of this channel which features regular α-helices that form a pentameric left-handed parallel bundle. The drug HMA was found to bind inside the lumen of the channel, at both the C-terminal and the N-terminal openings, and, in contrast to amiloride, induced additional chemical shifts in ETM. Full length SARS-CoV E displayed channel activity when transiently expressed in human embryonic kidney 293 (HEK-293) cells in a whole-cell patch clamp set-up. This activity was significantly reduced by hexamethylene amiloride (HMA), but not by amiloride. The channel structure presented herein provides a possible rationale for inhibition, and a platform for future structure-based drug design of this potential pharmacological target.

Coronaviruses are viral pathogens that cause a variety of lethal diseases in birds and mammals, and common colds in humans. In 2003, however, an animal coronavirus was able to infect humans and produced severe acute respiratory syndrome (SARS), causing a near pandemic. Such events are likely to reoccur in the future, and new antiviral strategies are necessary. A small coronavirus protein called ‘envelope’ is important for pathogenesis, affecting the formation of the viral envelope and the distribution of the virus in the body. In vitro studies have shown that synthetic coronavirus envelope proteins have channel activity that in some cases has been inhibited by the drug hexamethylene amiloride, but not by amiloride. In the present paper, we have characterized the structure responsible for this channel activity. We have also determined the binding site of the drug hexamethylene amiloride in the channel, and shown that amiloride has only a mild effect on the NMR signals from the protein. The validity of these results is supported using mammalian cells expressing full length SARS-CoV E, where channel activity was inhibited by hexamethylene amiloride, but only mildly by amiloride. The structural model described for this channel provides a valuable insight into coronavirus envelope protein ion channel activity, and could serve as a platform for the development of novel anti-viral drugs.

Structure and function of a complex between chorismate mutase and DAHP synthase: efficiency boost for the junior partner

Chorismate mutase catalyzes a key step in the shikimate biosynthetic pathway towards phenylalanine and tyrosine. Curiously, the intracellular chorismate mutase of Mycobacterium tuberculosis (MtCM; Rv0948c) has poor activity and lacks prominent active-site residues. However, its catalytic efficiency increases >100-fold on addition of DAHP synthase (MtDS; Rv2178c), another shikimate-pathway enzyme. The 2.35 A crystal structure of the MtCM-MtDS complex bound to a transition-state analogue shows a central core formed by four MtDS subunits sandwiched between two MtCM dimers. Structural comparisons imply catalytic activation to be a consequence of the repositioning of MtCM active-site residues on binding to MtDS. The mutagenesis of the C-terminal extrusion of MtCM establishes conserved residues as part of the activation machinery. The chorismate-mutase activity of the complex, but not of MtCM alone, is inhibited synergistically by phenylalanine and tyrosine. The complex formation thus endows the shikimate pathway of M. tuberculosis with an important regulatory feature. Experimental evidence suggests that such non-covalent enzyme complexes comprising an AroQ(delta) subclass chorismate mutase like MtCM are abundant in the bacterial order Actinomycetales.

The empirical valence bond as an effective strategy for computer-aided enzyme design

The ability of the empirical valence bond (EVB) to be used in screening active site residues in enzyme design is explored in a preliminary way. This validation is done by comparing the ability of this approach to evaluate the catalytic contributions of various residues in chorismate mutase. It is demonstrated that the EVB model can serve as an accurate tool in the final stages of computer-aided enzyme design (CAED). The ability of the model to predict quantitatively the catalytic power of enzymes should augment the capacity of current approaches for enzyme design.

Toward accurate screening in computer-aided enzyme design

The ability to design effective enzymes is one of the most fundamental challenges in biotechnology and in some respects in biochemistry. In fact, such ability would be one of the most convincing manifestations of a full understanding of the origin of enzyme catalysis. In this work, we explore the reliability of different simulation approaches, in terms of their ability to rank different possible active site constructs. This validation is done by comparing the ability of different approaches to evaluate the catalytic contributions of various residues in chorismate mutase. It is demonstrated that the empirical valence bond (EVB) model can serve as a practical yet accurate tool in the final stages of computer-aided enzyme design (CAED). Other approaches for fast screening are also examined and found to be less accurate and mainly useful for qualitative screening of ionized residues. It is pointed out that accurate ranking of different options for enzyme design cannot be accomplished by approaches that cannot capture the electrostatic preorganization effect. This is in particular true with regard to current design approaches that use gas phase or small cluster calculations and then estimate the interaction between the enzyme and the transition state (TS) model rather than the TS binding free energy or the relevant activation free energy. The ability of the EVB model to provide a tool for quantitative ranking in the final stage of CAED may help in progressing toward the design of enzymes whose catalytic power is closer to that of native enzymes than to that of the current generation of designer enzymes.

Relative tolerance of an enzymatic molten globule and its thermostable counterpart to point mutation

Enzyme structures reflect the complex interplay between the free energy of unfolding (DeltaG) and catalytic efficiency. Consequently, the effects of point mutations on structure, stability, and function are difficult to predict. It has been proposed that the mutational robustness of homologous enzymes correlates with a higher initial DeltaG. To examine this issue, we compared the tolerance of a natural thermostable chorismate mutase and an engineered molten globular variant to targeted mutation. These mutases possess similar sequence, structure, and catalytic efficiency but dramatically different DeltaG values. We find that analogous point mutations can have widely divergent effects on catalytic activity in these scaffolds. In a set of five rationally designed single-amino acid changes, the thermostable scaffold suffers activity losses ranging from 50-fold smaller, for an aspartate-to-glycine substitution at the active site, to 2-fold greater, for a phenylalanine-to-tryptophan substitution in the hydrophobic core, versus that of the molten globular scaffold. However, biophysical characterization indicates that the variations in catalytic efficiency are not caused by losses of either secondary structural integrity or thermodynamic stability. Rather, the activity differences between variant pairs are very much context-dependent and likely stem from subtle changes in the fine structure of the active site. Thus, in many cases, it may be more productive to focus on changes in local conformation than on global stability when attempting to understand and predict how enzymes respond to point mutations.

Structural reorganization and preorganization in enzyme active sites: comparisons of experimental and theoretically ideal active site geometries in the multistep serine esterase reaction cycle

Many enzymes catalyze reactions with multiple chemical steps, requiring the stabilization of multiple transition states during catalysis. Such enzymes must strike a balance between the conformational reorganization required to stabilize multiple transition states of a reaction and the confines of a preorganized active site in the polypeptide tertiary structure. Here we investigate the compromise between structural reorganization during the catalytic process and preorganization of the active site for a multistep enzyme-catalyzed reaction, the hydrolysis of esters by the Ser-His-Asp/Glu catalytic triad. Quantum mechanical transition states were used to generate ensembles of geometries that can catalyze each individual step in the mechanism. These geometries are compared to each other by superpositions of catalytic atoms to find "consensus" geometries that can catalyze all steps with minimal rearrangement. These consensus geometries are found to be excellent matches for the natural active site. Preorganization is therefore found to be the major defining characteristic of the active site, and reorganizational motions often proposed to promote catalysis have been minimized. The variability of enzyme active sites observed by X-ray crystallography was also investigated empirically. A catalog of geometrical parameters relating active site residues to each other and to bound inhibitors was collected from a set of crystal structures. The crystal-structure-derived values were then compared to the ranges found in quantum mechanically optimized structures along the entire reaction coordinate. The empirical ranges are found to encompass the theoretical ranges when thermal fluctuations are taken into account. Therefore, the active sites are preorganized to a geometry that can be objectively and quantitatively defined as minimizing conformational reorganization while maintaining optimal transition state stabilization for every step during catalysis. The results provide a useful guiding principle for de novo design of enzymes with multistep mechanisms.

On the relationship between folding and chemical landscapes in enzyme catalysis

Elucidating the relationship between the folding landscape of enzymes and their catalytic power has been one of the challenges of modern enzymology. The present work explores this issue by using a simplified folding model to generate the free-energy landscape of an enzyme and then to evaluate the activation barriers for the chemical step in different regions of the landscape. This approach is used to investigate the recent finding that an engineered monomeric chorismate mutase exhibits catalytic efficiency similar to the naturally occurring dimer even though it exhibits the properties of an intrinsically disordered molten globule. It is found that the monomer becomes more confined than its native-like counterpart upon ligand binding but still retains a wider catalytic region. Although the overall rate acceleration is still determined by reduction of the reorganization energy, the detailed contribution of different barriers yields a more complex picture for the chemical process than that of a single path. This work provides insight into the relationship between folding landscapes and catalysis. The computational approach used here may also provide a powerful strategy for modeling single-molecule experiments and designing enzymes.

Kinetics and thermodynamics of ligand binding to a molten globular enzyme and its native counterpart

An engineered monomeric chorismate mutase (mMjCM) has been found to combine high catalytic activity with the characteristics of a molten globule. To gain insight into the dramatic structural changes that accompany binding of a transition-state analog, we examined mMjCM by isothermal calorimetry and compared it with its dimeric parent protein, MjCM (CM from Methanococcus jannaschii), a thermostable and conventionally folded enzyme. As expected for a ligand-induced ordering process, there is a large entropic penalty for binding to the monomer relative to the dimer (-TDeltaDeltaS=5.1+/-0.5 kcal/mol, at 20 degrees C). However, this unfavorable entropy term is largely offset by enthalpic gains (DeltaDeltaH=-3.5+/-0.4 kcal/mol), presumably arising from tightening of non-covalent interactions throughout the monomeric complex. Stopped-flow kinetic measurements further reveal that the catalytic molten globule binds and releases ligands significantly faster than its natural counterpart, demonstrating that partial structural disorder can speed up molecular recognition. These results illustrate how structural plasticity may strongly perturb the thermodynamics and kinetics of transition-state recognition while negligibly affecting catalytic efficiency.

Biological function in a non-native partially folded state of a protein

As structural flexibility is known to be required for enzyme catalysis and pattern recognition and a significant fraction of eukaryotic proteins appear to be unfolded or contain unstructured regions, biological activity of conformational states distinct from fully folded structures could be more common than previously thought. By applying a procedure that allows the recovery of enzymatic activity to be monitored in real time, we show that a non-native state populated transiently during folding of the acylphosphatase from Sulfolobus solfataricus is enzymatically active. The structural characterization of this partially folded state reveals that enzymatic activity is possible even if the catalytic site is structurally heterogeneous, whereas the remainder of the structure acts as a scaffold. These results extend the spectrum of biological functions carried out in the absence of a folded state to include enzyme catalysis.