Design of fast enzymes by optimizing interaction potential in active site

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Rate Constants and Mechanisms of Protein-Ligand Binding

Whereas protein-ligand binding affinities have long-established prominence, binding rate constants and binding mechanisms have gained increasing attention in recent years. Both new computational methods and new experimental techniques have been developed to characterize the latter properties. It is now realized that binding mechanisms, like binding rate constants, can and should be quantitatively determined. In this review, we summarize studies and synthesize ideas on several topics in the hope of providing a coherent picture of and physical insight into binding kinetics. The topics include microscopic formulation of the kinetic problem and its reduction to simple rate equations; computation of binding rate constants; quantitative determination of binding mechanisms; and elucidation of physical factors that control binding rate constants and mechanisms.

Design of a protein tag and fluorogenic probe with modular structure for live-cell imaging of intracellular proteins

Conditional fluorescence imaging is a powerful technique for precise spatiotemporal analysis of proteins in live cells upon administration of a synthetic probe. To be applicable to various biological phenomena, probes must be membrane-permeable, have a high signal-to-noise ratio, and work quickly. To date, few probes meet all of these requirements. Here, we designed a fluorogenic probe (AcFCANB) that could label intracellular proteins fused to the photoactive yellow protein (PYP) tag in live cells within 30 min and used it to image heterochromatin protein 1 localization in nuclei. AcFCANB is based on a modular platform consisting of fluorophore, ligand and quencher. We accelerated the labeling reaction by strategic mutations of charged residues on the surface of PYP. A simple model based on molecular dynamics simulations quantitatively reproduced the cooperative effect of multiple mutations on labeling rate.

Molecular cloning, characterization and positively selected sites of the glutathione S-transferase family from Locusta migratoria

Glutathione S-transferases (GSTs) are multifunctional enzymes that are involved in the metabolism of endogenous and exogenous compounds and are related to insecticide resistance. The purpose of this study was to provide new information on the molecular characteristics and the positive selection of locust GSTs. Based on the transcriptome database, we sequenced 28 cytosolic GSTs and 4 microsomal GSTs from the migratory locust (Locusta migratoria). We assigned the 28 cytosolic GSTs into 6 classes—sigma, epsilon, delta, theta, omega and zeta, and the 4 microsomal GSTs into 2 subclasses—insect and MGST3. The tissue- and stage-expression patterns of the GSTs differed at the mRNA level. Further, the substrate specificities and kinetic constants of the cytosolic GSTs differed markedly at the protein level. The results of likelihood ratio tests provided strong evidence for positive selection in the delta class. The result of Bayes Empirical Bayes analysis identified 4 amino acid sites in the delta class as positive selection sites. These sites were located on the protein surface. Our findings will facilitate the elucidation of the molecular characteristics and evolutionary aspects of insect GST superfamily.

A method for computing association rate constants of atomistically represented proteins under macromolecular crowding

In cellular environments, two protein molecules on their way to form a specific complex encounter many bystander macromolecules. The latter molecules, or crowders, affect both the energetics of the interaction between the test molecules and the dynamics of their relative motion. In earlier work (Zhou and Szabo 1991 J. Chem. Phys. 95 5948-52), it has been shown that, in modeling the association kinetics of the test molecules, the presence of crowders can be accounted for by their energetic and dynamic effects. The recent development of the transient-complex theory for protein association in dilute solutions makes it possible to easily incorporate the energetic and dynamic effects of crowders. The transient complex refers to a late on-pathway intermediate, in which the two protein molecules have near-native relative separation and orientation, but have yet to form the many short-range specific interactions of the native complex. The transient-complex theory predicts the association rate constant as k(a) = k(a0)exp(-ΔG*(el)/k(B)T), where k(a0) is the 'basal' rate constant for reaching the transient complex by unbiased diffusion, and the Boltzmann factors captures the influence of long-range electrostatic interactions between the protein molecules. Crowders slow down the diffusion, therefore reducing the basal rate constant (to k(ac0)), and induce an effective interaction energy ΔG(c). We show that the latter interaction energy for atomistic proteins in the presence of spherical crowders is 'long'-ranged, allowing the association rate constant under crowding to be computed as k(ac) = k(ac0)exp[-(ΔG*(el) + ΔG*(c))/k(B)T]. Applications demonstrate that this computational method allows for realistic modeling of protein association kinetics under crowding.

Design, mechanism of action, bioavailability and therapeutic effects of mn porphyrin-based redox modulators

Based on aqueous redox chemistry and simple in vivo models of oxidative stress, Escherichia coli and Saccharomyces cerevisiae, the cationic Mn(III) N-substituted pyridylporphyrins (MnPs) have been identified as the most potent cellular redox modulators within the porphyrin class of drugs; their efficacy in animal models of diseases that have oxidative stress in common is based on their high ability to catalytically remove superoxide, peroxynitrite, carbonate anion radical, hypochlorite, nitric oxide, lipid peroxyl and alkoxyl radicals, thus suppressing the primary oxidative event. While doing so MnPs could couple with cellular reductants and redox-active proteins. Reactive species are widely accepted as regulators of cellular transcriptional activity: minute, nanomolar levels are essential for normal cell function, while submicromolar or micromolar levels impose oxidative stress, which is evidenced in increased inflammatory and immune responses. By removing reactive species, MnPs affect redox-based cellular transcriptional activity and consequently secondary oxidative stress, and in turn inflammatory processes. The equal ability to reduce and oxidize superoxide during the dismutation process and recently accumulated results suggest that pro-oxidative actions of MnPs may also contribute to their therapeutic effects. All our data identify the superoxide dismutase-like activity, estimated by log k(cat)O2-*), as a good measure for the therapeutic efficacy of MnPs. Their accumulation in mitochondria and their ability to cross the blood-brain barrier contribute to their remarkable efficacy. We summarize herein the therapeutic effects of MnPs in cancer, central nervous system injuries, diabetes, their radioprotective action and potential for imaging. Few of the most potent modulators of cellular redox-based pathways, MnTE2-PyP5+, MnTDE-2-ImP5+, MnTnHex-2-PyP5+ and MnTnBuOE-2-PyP5+, are under preclinical and clinical development.

Cellular level models as tools for cytokine design

Cytokines and growth factors are critical regulators that connect intracellular and extracellular environments through binding to specific cell-surface receptors. They regulate a wide variety of immunological, growth, and inflammatory response processes. The overall signal initiated by a population of cytokine molecules over long time periods is controlled by the subtle interplay of binding, signaling, and trafficking kinetics. Building on the work of others, we abstract a simple kinetic model that captures relevant features from cytokine systems as well as related growth factor systems. We explore a large range of potential biochemical behaviors, through systematic examination of the model's parameter space. Different rates for the same reaction topology lead to a dramatic range of biochemical network properties and outcomes. Evolution might productively explore varied and different portions of parameter space to create beneficial behaviors, and effective human therapeutic intervention might be achieved through altering network kinetic properties. Quantitative analysis of the results reveals the basis for tensions among a number of different network characteristics. For example, strong binding of cytokine to receptor can increase short-term receptor activation and signal initiation but decrease long-term signaling due to internalization and degradation. Further analysis reveals the role of specific biochemical processes in modulating such tensions. For instance, the kinetics of cytokine binding and receptor activation modulate whether ligand-receptor dissociation can generally occur before signal initiation or receptor internalization. Beyond analysis, the same models and model behaviors provide an important basis for the design of more potent cytokine therapeutics by providing insight into how binding kinetics affect ligand potency.

Rate theories for biologists

Some of the rate theories that are most useful for modeling biological processes are reviewed. By delving into some of the details and subtleties in the development of the theories, the review will hopefully help the reader gain a more than superficial perspective. Examples are presented to illustrate how rate theories can be used to generate insight at the microscopic level into biomolecular behaviors. An attempt is made to clear up a number of misconceptions in the literature regarding popular rate theories, including the appearance of Planck's constant in the transition-state theory and the Smoluchowski result as an upper limit for protein-protein and protein-DNA association rate constants. Future work in combining the implementation of rate theories through computer simulations with experimental probes of rate processes, and in modeling effects of intracellular environments so that theories can be used for generating rate constants for systems biology studies is particularly exciting.

Improved binding of raf to Ras.GDP is correlated with biological activity

The GTP-binding protein Ras plays a central role in the regulation of various cellular processes, acting as a molecular switch that triggers signaling cascades. Only Ras bound to GTP is able to interact strongly with effector proteins like Raf kinase, phosphatidylinositol 3-kinase, and RalGDS, whereas in the GDP-bound state, the stability of the complex is strongly decreased, and signaling is interrupted. To determine whether this process is only controlled by the stability of the complex, we used computer-aided protein design to improve the interaction between Ras and effector. We challenged the Ras.Raf complex in this study because Raf among all effectors shows the highest Ras affinity and the fastest association kinetics. The proposed mutations were characterized as to their changes in dynamics and binding strength. We demonstrate that Ras-Raf interaction can only be improved at the cost of a loss in specificity of Ras.GTP versus Ras.GDP. As shown by NMR spectroscopy, the Raf mutation A85K leads to a shift of Ras switch I in the GTP-bound as well as in the GDP-bound state, thereby increasing the complex stability. In a luciferase-based reporter gene assay, Raf A85K is associated with higher signaling activity, which appears to be a mere matter of Ras-Raf affinity.

Cross-species analysis of the glycolytic pathway by comparison of molecular interaction fields

The electrostatic potential of an enzyme is a key determinant of its substrate interactions and catalytic turnover. Here we invoke comparative analysis of protein electrostatic potentials, along with sequence and structural analysis, to classify and characterize all the enzymes in an entire pathway across a set of different organisms. The electrostatic potentials of the enzymes from the glycolytic pathway of 11 eukaryotes were analyzed by qPIPSA (quantitative protein interaction property similarity analysis). The comparison allows the functional assignment of neuron-specific isoforms of triosephosphate isomerase from zebrafish, the identification of unusual protein surface interaction properties of the mosquito glucose-6-phosphate isomerase and the functional annotation of ATP-dependent phosphofructokinases and cofactor-dependent phosphoglycerate mutases from plants. We here show that plants possess two parallel pathways to convert glucose. One is similar to glycolysis in humans, the other is specialized to let plants adapt to their environmental conditions. We use differences in electrostatic potentials to estimate kinetic parameters for the triosephosphate isomerases from nine species for which published parameters are not available. Along the core glycolytic pathway, phosphoglycerate mutase displays the most conserved electrostatic potential. The largest cross-species variations are found for glucose-6-phosphate isomerase, enolase and fructose-1,6-bisphosphate aldolase. The extent of conservation of electrostatic potentials along the pathway is consistent with the absence of a single rate-limiting step in glycolysis.

Kinetic and thermodynamic properties of the folding and assembly of formate dehydrogenase

The folding mechanism and stability of dimeric formate dehydrogenase from Candida methylica was analysed by exposure to denaturing agents and to heat. Equilibrium denaturation data yielded a dissociation constant of about 10(-13)M for assembly of the protein from unfolded chains and the kinetics of refolding and unfolding revealed that the overall process comprises two steps. In the first step a marginally stable folded monomeric state is formed at a rate (k(1)) of about 2x10(-3)s(-1) (by deduction k(-1) is about 10(-4)s(-1)) and assembles into the active dimeric state with a bimolecular rate constant (k(2)) of about 2x10(4)M(-1)s(-1). The rate of dissociation of the dimeric state in physiological conditions is extremely slow (k(-2) approximately 3x10(-7)s(-1)).

Protein-ion binding process on finite macromolecular concentration. A Poisson-Boltzmann and Monte Carlo study

Electrostatic interactions are one of the key driving forces for protein-ligands complexation. Different levels for the theoretical modeling of such processes are available on the literature. Most of the studies on the Molecular Biology field are performed within numerical solutions of the Poisson-Boltzmann Equation and the dielectric continuum models framework. In such dielectric continuum models, there are two pivotal questions: (a) how the protein dielectric medium should be modeled, and (b) what protocol should be used when solving this effective Hamiltonian. By means of Monte Carlo (MC) and Poisson-Boltzmann (PB) calculations, we define the applicability of the PB approach with linear and nonlinear responses for macromolecular electrostatic interactions in electrolyte solution, revealing some physical mechanisms and limitations behind it especially due the raise of both macromolecular charge and concentration out of the strong coupling regime. A discrepancy between PB and MC for binding constant shifts is shown and explained in terms of the manner PB approximates the excess chemical potentials of the ligand, and not as a consequence of the nonlinear thermal treatment and/or explicit ion-ion interactions as it could be argued. Our findings also show that the nonlinear PB predictions with a low dielectric response well reproduce the pK shifts calculations carried out with an uniform dielectric model. This confirms and completes previous results obtained by both MC and linear PB calculations.

Impact of electrostatics in redox modulation of oxidative stress by Mn porphyrins: protection of SOD-deficient Escherichia coli via alternative mechanism where Mn porphyrin acts as a Mn carrier

Understanding the factors that determine the ability of Mn porphyrins to scavenge reactive species is essential for tuning their in vivo efficacy. We present herein the revised structure-activity relationships accounting for the critical importance of electrostatics in the Mn porphyrin-based redox modulation systems and show that the design of effective SOD mimics (per se) based on anionic porphyrins is greatly hindered by inappropriate electrostatics. A new strategy for the beta-octabromination of the prototypical anionic Mn porphyrins Mn(III) meso-tetrakis(p-carboxylatophenyl)porphyrin ([Mn(III)TCPP](3-) or MnTBAP(3-)) and Mn(III) meso-tetrakis(p-sulfonatophenyl)porphyrin ([Mn(III)TSPP](3-)), to yield the corresponding anionic analogues [Mn(III)Br(8)TCPP](3-) and [Mn(III)Br(8)TSPP](3-), respectively, is described along with characterization data, stability studies, and their ability to substitute for SOD in SOD-deficient Escherichia coli. Despite the Mn(III)/Mn(II) reduction potential of [Mn(III)Br(8)TCPP](3-) and [Mn(III)Br(8)TSPP](3-) being close to the SOD-enzyme optimum and nearly identical to that of the cationic Mn(III) meso-tetrakis(N-methylpyridinium-2-yl)porphyrin (Mn(III)TM-2-PyP(5+)), the SOD activity of both anionic brominated porphyrins ([Mn(III)Br(8)TCPP](3-), E(1/2)=+213 mV vs NHE, log k(cat)=5.07; [Mn(III)Br(8)TSPP](3-), E(1/2)=+209 mV, log k(cat)=5.56) is considerably lower than that of Mn(III)TM-2-PyP(5+) (E(1/2)=+220 mV, log k(cat)=7.79). This illustrates the impact of electrostatic guidance of O(2)(-) toward the metal center of the mimic. With low k(cat), the [Mn(III)TCPP](3-), [Mn(III)TSPP](3-), and [Mn(III)Br(8)TCPP](3-) did not rescue SOD-deficient E. coli. The striking ability of [Mn(III)Br(8)TSPP](3-) to substitute for the SOD enzymes in the E. coli model does not correlate with its log k(cat). In fact, the protectiveness of [Mn(III)Br(8)TSPP](3-) is comparable to or better than that of the potent SOD mimic Mn(III)TM-2-PyP(5+), even though the dismutation rate constant of the anionic complex is 170-fold smaller. Analyses of the medium and E. coli cell extract revealed that the major species in the [Mn(III)Br(8)TSPP](3-) system is not the Mn complex, but the free-base porphyrin [H(2)Br(8)TSPP](4-) instead. Control experiments with extracellular MnCl(2) showed the lack of E. coli protection, indicating that "free" Mn(2+) cannot enter the cell to a significant extent. We proposed herein the alternative mechanism where a labile Mn porphyrin [Mn(III)Br(8)TSPP](3-) is not an SOD mimic per se but carries Mn into the E. coli cell.

Electrostatic rate enhancement and transient complex of protein-protein association

The association of two proteins is bounded by the rate at which they, via diffusion, find each other while in appropriate relative orientations. Orientational constraints restrict this rate to approximately 10(5)-10(6) M(-1) s(-1). Proteins with higher association rates generally have complementary electrostatic surfaces; proteins with lower association rates generally are slowed down by conformational changes upon complex formation. Previous studies (Zhou, Biophys J 1997;73:2441-2445) have shown that electrostatic enhancement of the diffusion-limited association rate can be accurately modeled by $k_{\bf D}$ = $k_{D}0\ {exp} ( - \langle U_{el} \rangle;{\star}/k_{B} T),$ where k(D) and k(D0) are the rates in the presence and absence of electrostatic interactions, respectively, U(el) is the average electrostatic interaction energy in a "transient-complex" ensemble, and k(B)T is the thermal energy. The transient-complex ensemble separates the bound state from the unbound state. Predictions of the transient-complex theory on four protein complexes were found to agree well with the experiment when the electrostatic interaction energy was calculated with the linearized Poisson-Boltzmann (PB) equation (Alsallaq and Zhou, Structure 2007;15:215-224). Here we show that the agreement is further improved when the nonlinear PB equation is used. These predictions are obtained with the dielectric boundary defined as the protein van der Waals surface. When the dielectric boundary is instead specified as the molecular surface, electrostatic interactions in the transient complex become repulsive and are thus predicted to retard association. Together these results demonstrate that the transient-complex theory is predictive of electrostatic rate enhancement and can help parameterize PB calculations.

Redox modulation of oxidative stress by Mn porphyrin-based therapeutics: the effect of charge distribution

We evaluate herein the impact of positive charge distribution on the in vitro and in vivo properties of Mn porphyrins as redox modulators possessing the same overall 5+ charge and of minimal stericity demand: Mn(III) meso-tetrakis(trimethylanilinium-4-yl)porphyrin (MnTTriMAP(5+)), Mn(III) meso-tetrakis(N,N'-dimethylpyrazolium-4-yl)porphyrin (MnTDM-4-PzP(5+)), Mn(III) meso-tetrakis(N,N'-dimethylimidazolium-2-yl)porphyrin (MnTDM-2-ImP(5+)), and the ortho and para methylpyridinium complexes Mn(III) meso-tetrakis(N-methylpyridinium-4-yl)porphyrin (MnTM-4-PyP(5+)) and Mn(III) meso-tetrakis(N-methylpyridinium-2-yl)porphyrin (MnTM-2-PyP(5+)). Both Mn(III)/Mn(II) reduction potential and SOD activity within the series follow the order: MnTTriMAP(5+)<MnTDM-4-PzP(5+)<MnTM-4-PyP(5+)<MnTM-2-PyP(5+)<MnTDM-2-ImP(5+). The kinetic salt effect (KSE) on the catalytic rate constant for superoxide dismutation (k(cat)) indicates that the electrostatic contribution to the O(2)*(-) dismutation is the greatest with MnTM-2-PyP(5+) and follows the order: MnTM-4-PyP(5+)<MnTDM-4-PzP(5+) approximately MnTDM-2-ImP(5+)<MnTM-2-PyP(5+). The KSE observed on k(cat) suggests that the charges are relatively confined within specific regions of the aryl rings. Whereas the charges in imidazolium, pyrazolium, and MnTM-4-PyP(5+) compounds are distributed in-plane with the porphyrin ring, the charges of MnTM-2-PyP(5+) are either above or below the plane, which channels the negatively-charged superoxide toward the axial positions of the Mn porphyrin more efficiently, and leads to the highest KSE. This mimics the tunneling effect observed in the SOD enzymes themselves. The modulation of the reactivity of the Mn center by the electronic perturbations caused by the meso-aryl substituent could be explained by DFT calculation, whereby a correlation between the Mn(III)/Mn(II) reduction potential (and/or SOD activity) and meso-aryl fragment softness descriptors for nucleophilic (s(f)(+)) and radical (s(f)(o)) attacks was observed. MnTDM-4-PzP(5+) and MnTM-4-PyP(5+) did not protect SOD-deficient E. coli grown aerobically, which is in agreement with their low k(cat). MnTM-2-PyP(5+) and MnTDM-2-ImP(5+) have similar high k(cat), but MnTDM-2-ImP(5+) was significantly less protective to E. coli, probably due to its bulkier size, decreased cellular uptake, and/or observed toxicity. The placement of charges closer to the metal center and spatial charge localization increases both the in vitro and the in vivo SOD activity of the compound.

Prediction of salt and mutational effects on the association rate of U1A protein and U1 small nuclear RNA stem/loop II

We have developed a computational approach for predicting protein-protein association rates (Alsallaq and Zhou, Structure 2007, 15, 215). Here we expand the range of applicability of this approach to protein-RNA binding and report the first results for protein-RNA binding rates predicted from atomistic modeling. The system studied is the U1A protein and stem/loop II of the U1 small nuclear RNA. Experimentally it was observed that the binding rate is significantly reduced by increasing salt concentration while the dissociation changes little with salt concentration, and charges distant from the binding site make marginal contribution to the binding rate. These observations are rationalized. Moreover, predicted effects of salt and charge mutations are found to be in quantitative agreement with experimental results.

qPIPSA: relating enzymatic kinetic parameters and interaction fields

Background

The simulation of metabolic networks in quantitative systems biology requires the assignment of enzymatic kinetic parameters. Experimentally determined values are often not available and therefore computational methods to estimate these parameters are needed. It is possible to use the three-dimensional structure of an enzyme to perform simulations of a reaction and derive kinetic parameters. However, this is computationally demanding and requires detailed knowledge of the enzyme mechanism. We have therefore sought to develop a general, simple and computationally efficient procedure to relate protein structural information to enzymatic kinetic parameters that allows consistency between the kinetic and structural information to be checked and estimation of kinetic constants for structurally and mechanistically similar enzymes.

Results

We describe qPIPSA: quantitative Protein Interaction Property Similarity Analysis. In this analysis, molecular interaction fields, for example, electrostatic potentials, are computed from the enzyme structures. Differences in molecular interaction fields between enzymes are then related to the ratios of their kinetic parameters. This procedure can be used to estimate unknown kinetic parameters when enzyme structural information is available and kinetic parameters have been measured for related enzymes or were obtained under different conditions. The detailed interaction of the enzyme with substrate or cofactors is not modeled and is assumed to be similar for all the proteins compared. The protein structure modeling protocol employed ensures that differences between models reflect genuine differences between the protein sequences, rather than random fluctuations in protein structure.

Conclusion

Provided that the experimental conditions and the protein structural models refer to the same protein state or conformation, correlations between interaction fields and kinetic parameters can be established for sets of related enzymes. Outliers may arise due to variation in the importance of different contributions to the kinetic parameters, such as protein stability and conformational changes. The qPIPSA approach can assist in the validation as well as estimation of kinetic parameters, and provide insights into enzyme mechanism.

Prediction of protein-protein association rates from a transition-state theory

We recently developed a theory for the rates of protein-protein association. The theory is based on the concept of a transition state, which separates the bound state, with numerous short-range interactions but restricted translational and rotational freedom, and the unbound state, with, at most, a small number of interactions but expanded configurational freedom. When not accompanied by large-scale conformational changes, protein-protein association becomes diffusion limited. The association rate is then predicted as k(a)=k(a)(0)exp(-DeltaG(el)(double dagger)/k(B)T), where DeltaG(el)(double dagger) is the electrostatic interaction free energy in the transition state, k(a)(0) is the rate in the absence of electrostatic interactions, and k(B)T is thermal energy. Here, this transition-state theory is used to predict the association rates of four protein complexes. The predictions for the wild-type complexes and 23 mutants are found to agree closely with experimental data over wide ranges of ionic strength.

Electrostatic contribution in the catalysis of O2*- dismutation by superoxide dismutase mimics. MnIIITE-2-PyP5+ versus MnIIIBr8T-2-PyP+

The Mn(III) meso-tetrakis(N-ethylpyridinium-2-yl)porphyrin (Mn(III)TE-2-PyP(5+)) is a potent superoxide dismutase (SOD) mimic in vitro and was beneficial in rodent models of oxidative stress pathologies. Its high activity has been ascribed to both the favorable redox potential of its metal center and to the electrostatic facilitation assured by the four positive charges encircling the metal center. Its comparison with the non-alkylated, singly charged analogue Mn(III) beta-octabromo meso-tetrakis(2-pyridyl)porphyrin (Mn(III)Br(8)T-2-PyP(+)) enabled us to evaluate the electrostatic contribution to the catalysis of O(2)() dismutation. Both compounds exhibit nearly identical metal-centered redox potential for Mn(III)/Mn(II) redox couple: +228 mV for Mn(III)TE-2-PyP(5+) and +219 mV versus NHE for Mn(III)Br(8)T-2-PyP(+). The eight electron-withdrawing beta pyrrolic bromines contribute equally to the redox properties of the parent Mn(III)T-2-PyP(+) as do four quaternized cationic meso ortho pyridyl nitrogens. However, the SOD-like activity of the highly charged Mn(III)TE-2-PyP(5+) is >100-fold higher (log k(cat) = 7.76) than that of the singly charged Mn(III)Br(8)T-2-PyP(+) (log k(cat) = 5.63). The kinetic salt effect showed that the catalytic rate constants of the Mn(III)TE-2-PyP(5+) and of its methyl analogue, Mn(III)TM-2-PyP(5+), are exactly 5-fold more sensitive to ionic strength than is the k(cat) of Mn(III)Br(8)T-2-PyP(+), which parallels the charge ratio of these compounds. Interestingly, only a small effect of ionic strength on the rate constant was found in the case of penta-charged para (Mn(III)TM-4-PyP(5+)) and meta isomers (Mn(III)TM-3-PyP(5+)), indicating that the placement of the positive charges in the close proximity of the metal center (ortho position) is essential for the electrostatic facilitation of O(2)() dismutation.

Simulating enzyme reactions: challenges and perspectives

Elucidating how enzymes enhance the rates of the reactions that they catalyze is a major goal of contemporary biochemistry, and it is an area in which computational and theoretical techniques can make a major contribution. This article outlines some of the processes that need to be investigated if enzyme catalysis is to be understood, reviews the current state-of-the-art in enzyme simulation work, and highlights challenges for the future.

New insights into the mechanism of protein-protein association

Association of a protein complex follows a two-step mechanism, with the first step being the formation of an encounter complex that evolves into the final complex. Here, we analyze recent experimental data of the association of TEM1-beta-lactamase with BLIP using theoretical calculations and simulation. We show that the calculated Debye-Hückel energy of interaction for a pair of proteins during association resembles an energy funnel, with the final complex at the minima. All attraction is lost at inter-protein distances of 20 A, or rotation angles of >60 degrees from the orientation of the final complex. For faster-associating protein complexes, the energy funnel deepens and its volume increases. Mutations with the largest impact on association (hotspots for association) have the largest effect on the size and depth of the energy funnel. Analyzing existing evidence, we suggest that the transition state along the association pathway is the formation of the final complex from the encounter complex. Consequently, pairs of proteins forming an encounter complex will tend to dissociate more readily than to evolve into the final complex. Increasing directional diffusion by increasing favorable electrostatic attraction results in a faster forming and slower dissociating encounter complex. The possible applicability of electrostatic calculations for protein-protein docking is discussed.

Disparate ionic-strength dependencies of on and off rates in protein-protein association

Electrostatic interactions have been observed to play important roles in the kinetics of protein-protein association. Ionic strength, by its ability to modulate the magnitude of electrostatic interactions, has often been conveniently used to test their presence. From experiments on a wide range of associating proteins, a common feature has emerged: the on rates show strong dependence on ionic strength whereas the off rates are relatively insensitive. Here this feature is explained by an explicit description of a transition state for the association process and the suggestion that this transition is near the final bound state of two proteins. The molecular basis of the transition state in the bimolecular process lies in the fact that the bound state is characterized by local specific (e.g., van der Waals, hydrophobic, and electrostatic) interactions, whereas the unbound state is characterized by translational and rotational freedom. In the transition state the protein-protein pair encounters a free-energy maximum since its translational-rotational entropy is reduced while the specific interactions are not yet attained. In this formalism of the protein-protein association process, the enhancement of on rates by long-range electrostatic interactions can be written (analogous to an ordinary transition-state theory) in the form k(on) = k(0)(on)exp(-G(el)/k(B)T), where G(el) is the electrostatic free energy of the transition state.

Kinetics of diffusion-assisted reactions in microheterogeneous systems

This review is focused on the basic theory of diffusion-assisted reactions in microheterogeneous systems, from porous solids to self-organized colloids and biomolecules. Rich kinetic behaviors observed experimentally are explained in a unified fashion using simple concepts of competing distance and time scales of the reaction and the embedding structure. We mainly consider pseudo-first-order reactions, such as luminescence quenching, described by the Smoluchowski type of equation for the reactant pair distribution function with a sink term defined by the reaction mechanism. Microheterogeneity can affect the microscopic rate constant. It also enters the evolution equation through various spatial constraints leading to complicated boundary conditions and, possibly, to the reduction of dimensionality of the diffusion space. The reaction coordinate and diffusive motion along this coordinate are understood in a general way, depending on the problem at hand. Thus, the evolution operator can describe translational and rotational diffusion of molecules in a usual sense, it can be a discrete random walk operator when dealing with hopping of adsorbates in solids, or it can correspond to conformational fluctuations in proteins. Mathematical formulation is universal but physical consequences can be different. Understanding the principal features of reaction kinetics in microheterogeneous systems enables one to extract important structural and dynamical information about the host environments by analyzing suitably designed experiments, it helps building effective strategies for computer simulations, and ultimately opens possibilities for designing systems with controllable reactivity properties.

Predicting the rate enhancement of protein complex formation from the electrostatic energy of interaction

The rate of association of proteins is dictated by diffusion, but can be enhanced by favorable electrostatic forces. Here the relationship between the electrostatic energy of interaction, and the kinetics of protein-complex formation was analyzed for the protein pairs of: hirudin-thrombin, acetylcholinesterase-fasciculin and barnase-barstar, and for a panel of point mutants of these proteins. Electrostatic energies of interaction were calculated as the difference between the electrostatic energy of the complex and the sum of the energies of the two individual proteins, using the computer simulation package DelPhi. Calculated electrostatic energies of interaction were compared to experimentally determined rates of association. One kcal/mol of Coulombic interaction energy increased the rate of association by a factor of 2.8, independent of the protein-complex or mutant analyzed. Electrostatic energies of interaction were also determined from the salt dependence of the association rate constant, using the same basic equation as for the theoretical calculation. A Brönsted analysis of the electrostatic energies of interactions plotted versus experimentally determined ln(rate)s of association shows a linear relation between the two, with a beta value close to 1. This is interpreted as the energy of the transition state varies according to the electrostatic interaction energy, fitting a two state model for the association reaction. Calculating electrostatic rate enhancement from the electrostatic interaction energy can be used as a powerful tool to design protein complexes with altered rates of association and affinities.

Evolutionary constraints for dimer formation in prokaryotic Cu,Zn superoxide dismutase

Prokaryotic Cu,Zn superoxide dismutases are characterized by a distinct quaternary structure, as compared to that of the homologous eukaryotic enzymes. Here we report a newly determined crystal structure of the dimeric Cu,Zn superoxide dismutase from Photobacterium leiognathi (crystallized in space group R32, refined at 2.5 A resolution, R-factor 0.19) and analyse it in comparison with that of the monomeric enzyme from Escherichia coli. The dimeric assembly, observed also in a previously studied monoclinic crystal form of P. leiognathi Cu,Zn superoxide dismutase, is based on a ring-shaped subunit contact region, defining a solvated interface cavity. Three clusters of neighbouring residues play a direct role in the stabilization of the quaternary assembly. The present analysis, extended to the amino acid sequences of the other 11 known prokaryotic Cu,Zn superoxide dismutases, shows that at least in five other prokaryotic enzymes the interface residue clusters are under strong evolutionary constraint, suggesting the attainment of a quaternary structure coincident with that of P. leiognathi Cu,Zn superoxide dismutase. Calculation of electrostatic fields for both the enzymes from E. coli and P. leiognathi shows that the monomeric/dimeric association behaviour displayed by prokaryotic Cu, Zn superoxide dismutases is related to the distribution of surface charged residues. Moreover, Brownian dynamics simulations reproduce closely the observed enzyme:substrate association rates, highlighting the role of the active site neighbouring residues in determining the dismutase catalytic properties.

Role of the electrostatic loop charged residues in Cu,Zn superoxide dismutase

We have expressed and characterized a mutant of Xenopus laevis Cu,Zn superoxide dismutase in which four highly conserved charged residues belonging to the electrostatic loop have been replaced by neutral side chains: Lys120 --> Leu, Asp130 --> Gln, Glu131 --> Gln, and Lys134 --> Thr. At low ionic strength, the mutant enzyme is one of the fastest superoxide dismutases ever assayed (k = 6.7 x 10(9) M(-1) s(-1), at pH 7 and mu = 0.02 M). Brownian dynamics simulations give rise to identical enzyme-substrate association rates for both wild-type and mutant enzymes, ruling out the possibility that enhancement of the activity is due to pure electrostatic factors. Comparative analysis of the experimental catalytic rate of the quadruple and single mutants reveals the nonadditivity of the mutation effects, indicating that the hyperefficiency of the mutant is due to a decrease of the energy barrier and/or to an alternative pathway for the diffusion of superoxide within the active site channel. At physiological ionic strength the catalytic rate of the mutant at neutral pH is similar to that of the wild-type enzyme as it is to the catalytic rate pH dependence. Moreover, mutation effects are additive. These results show that, at physiological salt conditions, electrostatic loop charged residues do not influence the diffusion pathway of the substrate and, if concomitantly neutralized, are not essential for high catalytic efficiency of the enzyme, pointing out the role of the metal cluster and of the invariant Arg141 in determining the local electrostatic forces facilitating the diffusion of the substrate towards the active site.

Role of P225 and the C136-C201 disulfide bond in tissue plasminogen activator

The protease domain of tissue plasminogen activator (tPA), a key fibrinolytic enzyme, was expressed in Escherichia coli with a yield of 1 mg per liter of media. The recombinant protein was titrated with the Erythrina caraffa trypsin inhibitor (ETI) and characterized in its interaction with plasminogen and the natural inhibitor plasminogen activator inhibitor-1 (PAI-1). Analysis of the catalytic properties of tPA using a library of chromogenic substrates carrying substitutions at P1, P2, and P3 reveals a strong preference for Arg over Lys at P1, unmatched by other serine proteases like thrombin or trypsin. In contrast to these proteases and plasmin, tPA shows little or no preference for Pro over Gly at P2. A specific inhibition of tPA by Cu2+ was discovered. The divalent cation presumably binds to H188 near D189 in the primary specificity pocket and inhibits substrate binding in a competitive manner with a Kd = 19 microM. In an attempt to engineer Na+ binding and enhanced catalytic activity in tPA, P225 was replaced with Tyr, the residue present in Na+-dependent allosteric serine proteases. The P225Y mutation did not result in cation binding, but caused a significant loss of specificity (up to 100-fold) toward chromogenic substrates and plasminogen and considerably reduced the inhibition by PAI-1 and ETI. Interestingly, the P225Y substitution enhanced the ability of Cu2+ to inhibit the enzyme. Elimination of the C136-C201 disulfide bond, that is absent in all Na+-dependent allosteric serine proteases, significantly enhanced the yield (5 mg per liter of media) of expression in E. coli, but caused no changes in the properties of the enzyme whether residue 225 was Pro or Tyr. These findings point out an unanticipated crucial role for residue 225 in controlling the catalytic activity of tPA, and suggest that engineering of a Na+-dependent allosteric enhancement of catalytic activity in this enzyme, must involve substantial changes in the region homologous to the Na+ binding site of allosteric serine proteases.

Electrostatic steering and ionic tethering in enzyme-ligand binding: insights from simulations

To bind at an enzyme's active site, a ligand must diffuse or be transported to the enzyme's surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and beta-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as "ionic tethering." We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme's surroundings even when the substrate is nonpolar.

Electrostatic enhancement of diffusion-controlled protein-protein association: comparison of theory and experiment on barnase and barstar

The electrostatic enhancement of the association rate of barnase and barstar is calculated using a transition-state theory like expression and atomic-detail modeling of the protein molecules. This expression predicts that the rate enhancement is simply the average Boltzmann factor in the region of configurational space where association occurs instantaneously in the diffusion-controlled limit. Based on experimental evidence, this "transition state" is defined by configurations in which, relative to the stereospecifically bound complex, the two proteins are shifted apart by approximately 8 A (so a layer of water can be accommodated in the interface) and the two binding surfaces are rotated away by 0 degrees to 3 degrees. The values of the average Boltzmann factor, calculated by solving the Poisson-Boltzmann equation, for the wild-type complex and 16 complexes with single mutations are found to correlate well with experimental results for the electrostatic rate enhancement. The predicted rate enhancement is found to be somewhat insensitive to the precise definition of the transition state, due to the long-range nature of electrostatic interactions. The experimental ionic strength dependence of the rate enhancement is also reasonably reproduced.

Correlation between rate of enzyme-substrate diffusional encounter and average Boltzmann factor around active site

The utility of the average Boltzmann factor around the active site of an enzyme as the predictor of the electrostatic enhancement of the substrate binding rate is tested on a set of data on wild-type acetylcholinesterase and 18 charge mutants recently obtained by Brownian dynamics simulations. A good correlation between the average Boltzmann factors and the substrate binding rate constants is found. The effects of single charge mutations on both the Boltzmann factor and the substrate binding rate constant are modest, i.e., < 5 fold increase or decrease. This is consistent with the experimental results of Shafferman et al. but does not support their suggestion that the overall rate of the catalytic reaction is not limited by the diffusional encounter of acetylcholinesterase and its substrate.

Theory of biomolecular recognition

Specific, noncovalent binding of biomolecules can only be understood by considering structural, thermodynamic, and kinetic issues. The theoretical foundations for such analyses have been clarified in the past year. Computational techniques for both particle-based and continuum models continue to improve and to yield useful insights into an ever wider range of biomolecular systems.

Structural adaptations of the cold-active citrate synthase from an Antarctic bacterium

BACKGROUND

The structural basis of adaptation of enzymes to low temperature is poorly understood. Dimeric citrate synthase has been used as a model enzyme to study the structural basis of thermostability, the structure of the enzyme from organisms living in habitats at 55 degrees C and 100 degrees C having previously been determined. Here the study is extended to include a citrate synthase from an Antarctic bacterium, allowing us to explore the structural basis of cold activity and thermostability across the whole temperature range over which life is known to exit.

RESULTS

We report here the first crystal structure of a cold-active enzyme, citrate synthase, isolated from an Antarctic bacterium, at a resolution of 2.09 A. In comparison with the same enzyme from a hyperthermophilic host, the cold-active enzyme has a much more accessible active site, an unusual electrostatic potential distribution and an increased relative flexibility of the small domain compared to the large domain. Several other features of the cold-active enzyme were also identified: reduced subunit interface interactions with no intersubunit ion-pair networks; loops of increased length carrying more charge and fewer proline residues; an increase in solvent-exposed hydrophobic residues; and an increase in intramolecular ion pairs.

CONCLUSIONS

Enzymes from organisms living at the temperature extremes of life need to avoid hot or cold denaturation yet maintain sufficient structural integrity to allow catalytic efficiency. For hyperthermophiles, thermal denaturation of the citrate synthase dimer appears to be resisted by complex networks of ion pairs at the dimer interface, a feature common to other hyperthermophilic proteins. For the cold-active citrate synthase, cold denaturation appears to be resisted by an increase in intramolecular ion pairs compared to the hyperthermophilic enzyme. Catalytic efficiency of the cold-active enzyme appears to be achieved by a more accessible active site and by an increase in the relative flexibility of the small domain compared to the large domain.