Relating Protein Motion to Catalysis

Normal-mode-analysis-guided investigation of crucial intersubunit contacts in the cAMP-dependent gating in HCN channels

Protein structures define a complex network of atomic interactions in three dimensions. Direct visualization of the structure and analysis of the interaction potential energy are not straightforward approaches to pinpoint the atomic contacts that are crucial for protein function. We used the tetrameric hyperpolarization-activated cAMP-regulated (HCN) channel as a model system to study the intersubunit contacts in cAMP-dependent gating. To obtain a systematic survey of the contacts between each pair of residues, we used normal-mode analysis, a computational approach for studying protein dynamics, and constructed the covariance matrix for C-α atoms. The significant contacts revealed by covariance analysis were further investigated by means of mutagenesis and functional assays. Among the mutant channels that show phenotypes different from those of the wild-type, we focused on two mutant channels that express opposite changes in cAMP-dependent gating. Subsequent biochemical assays on isolated C-terminal fragments, including the cAMP binding domain, revealed only minimal effects on cAMP binding, suggesting the necessity of interpreting the cAMP-dependent allosteric regulation at the whole-channel level. For this purpose, we applied the patch-clamp fluorometry technique and observed correlated changes in the dynamic, state-dependent cAMP binding in the mutant channels. This study not only provides further understanding of the intersubunit contacts in allosteric coupling in the HCN channel, it also illustrates an effective strategy for delineating important atomic contacts within a structure.

Mass Modulation of Protein Dynamics Associated with Barrier Crossing in Purine Nucleoside Phosphorylase

The role of protein dynamics on different time scales in enzyme catalysis remains an area of active debate. The connection between enzyme dynamics on the femtosecond time scale and transition state formation has been demonstrated in human purine nucleoside phosphorylase (PNP) through the study of a mass-altered enzyme. Isotopic substitution in human PNP (heavy PNP) decreased the rate of on-enzyme chemistry but did not alter either the transition state structure or steady-state kinetic parameters. Here we investigate the underlying atomic motions associated with altered barrier crossing probability for heavy PNP. Transition path sampling was employed to illuminate the molecular differences between barrier crossing in light and heavy enzymes. The mass effect is apparent in promoting vibrations that polarize the N-ribosidic bond, and that promote the stability of the purine leaving group. These motions facilitate barrier crossing.

Protein dynamics studied with ultrafast two-dimensional infrared vibrational echo spectroscopy

Proteins, enzymes, and other biological molecules undergo structural dynamics as an intrinsic part of their biological functions. While many biological processes occur on the millisecond, second, and even longer time scales, the fundamental structural dynamics that eventually give rise to such processes occur on much faster time scales. Many decades ago, chemical kineticists focused on the inverse of the reaction rate constant as the important time scale for a chemical reaction. However, through transition state theory and a vast amount of experimental evidence, we now know that the key events in a chemical reaction can involve structural fluctuations that take a system of reactants to its transition state, the crossing of a barrier, and the eventual relaxation to product states. Such dynamics occur on very fast time scales. Today researchers would like to investigate the fast structural fluctuations of biological molecules to gain an understanding of how biological processes proceed from simple structural changes in biomolecules to the final, complex biological function. The study of the fast structural dynamics of biological molecules requires experiments that operate on the appropriate time scales, and in this Account, we discuss the application of ultrafast two-dimensional infrared (2D IR) vibrational echo spectroscopy to the study of protein dynamics. The 2D IR vibrational echo experiment is akin to 2D NMR, but it operates on time scales many orders of magnitude faster. In the experiments, a particular vibrational oscillator serves as a vibrational dynamics probe. As the structure of the protein evolves in time, the structural changes are manifested as time-dependent changes in the frequency of the vibrational dynamics probe. The 2D IR vibrational echo experiments can track the vibrational frequency evolution, which we then relate to the time evolution of the protein structure. In particular, we measured protein substate interconversion for mutants of myoglobin using 2D IR chemical exchange spectroscopy and observed well-defined substate interconversion on a sub-100 ps time scale. In another study, we investigated the influence of binding five different substrates to the enzyme cytochrome P450(cam). The various substrates affect the enzyme dynamics differently, and the observed dynamics are correlated with the enzyme's selectivity of hydroxylation of the substrates and with the substrate binding affinity.

Comparing proteins by their internal dynamics: exploring structure-function relationships beyond static structural alignments

The growing interest for comparing protein internal dynamics owes much to the realisation that protein function can be accompanied or assisted by structural fluctuations and conformational changes. Analogously to the case of functional structural elements, those aspects of protein flexibility and dynamics that are functionally oriented should be subject to evolutionary conservation. Accordingly, dynamics-based protein comparisons or alignments could be used to detect protein relationships that are more elusive to sequence and structural alignments. Here we provide an account of the progress that has been made in recent years towards developing and applying general methods for comparing proteins in terms of their internal dynamics and advance the understanding of the structure-function relationship.

A remote mutation affects the hydride transfer by disrupting concerted protein motions in thymidylate synthase

The role of protein flexibility in enzyme-catalyzed activation of chemical bonds is an evolving perspective in enzymology. Here we examine the role of protein motions in the hydride transfer reaction catalyzed by thymidylate synthase (TSase). Being remote from the chemical reaction site, the Y209W mutation of Escherichia coli TSase significantly reduces the protein activity, despite the remarkable similarity between the crystal structures of the wild-type and mutant enzymes with ligands representing their Michaelis complexes. The most conspicuous difference between these two crystal structures is in the anisotropic B-factors, which indicate disruption of the correlated atomic vibrations of protein residues in the mutant. This dynamically altered mutant allows a variety of small thiols to compete for the reaction intermediate that precedes the hydride transfer, indicating disruption of motions that preorganize the protein environment for this chemical step. Although the mutation causes higher enthalpy of activation of the hydride transfer, it only shows a small effect on the temperature dependence of the intrinsic KIE, suggesting marginal changes in the geometry and dynamics of the H-donor and -acceptor at the tunneling ready state. These observations suggest that the mutation disrupts the concerted motions that bring the H-donor and -acceptor together during the pre- and re-organization of the protein environment. The integrated structural and kinetic data allow us to probe the impact of protein motions on different time scales of the hydride transfer reaction within a complex enzymatic mechanism.

Protein dynamics and enzyme catalysis: the ghost in the machine?

One of the most controversial questions in enzymology today is whether protein dynamics are significant in enzyme catalysis. A particular issue in these debates is the unusual temperature-dependence of some kinetic isotope effects for enzyme-catalysed reactions. In the present paper, we review our recent model [Glowacki, Harvey and Mulholland (2012) Nat. Chem. 4, 169-176] that is capable of reproducing intriguing temperature-dependences of enzyme reactions involving significant quantum tunnelling. This model relies on treating multiple conformations of the enzyme-substrate complex. The results show that direct 'driving' motions of proteins are not necessary to explain experimental observations, and show that enzyme reactivity can be understood and accounted for in the framework of transition state theory.

Possible origin of life between mica sheets: does life imitate mica?

The mica hypothesis for the origin of life proposes that life originated between the sheets of muscovite mica. This paper elaborates on two ways that life resembles what might have originated between mica sheets. First, enzymes: The configurations and dynamics of enzymes, with their substrates, cofactors, and sometimes transition metal ions, often resemble mica sheets, with their open-and-shut motions, acting on small molecules between them, sometimes assisted by transition metal ions. Second, organisms: Mica world had the potential to be a community or ecosystem of prebiotic organisms in a way unlike other models for the origin of life.

Spatial distribution of radiation damage to crystalline proteins at 25-300 K

The spatial distribution of radiation damage (assayed by increases in atomic B factors) to thaumatin and urease crystals at temperatures ranging from 25 to 300 K is reported. The nature of the damage changes dramatically at approximately 180 K. Above this temperature the role of solvent diffusion is apparent in thaumatin crystals, as solvent-exposed turns and loops are especially sensitive. In urease, a flap covering the active site is the most sensitive part of the molecule and nearby loops show enhanced sensitivity. Below 180 K sensitivity is correlated with poor local packing, especially in thaumatin. At all temperatures, the component of the damage that is spatially uniform within the unit cell accounts for more than half of the total increase in the atomic B factors and correlates with changes in mosaicity. This component may arise from lattice-level, rather than local, disorder. The effects of primary structure on radiation sensitivity are small compared with those of tertiary structure, local packing, solvent accessibility and crystal contacts.

Gauging the flexibility of the active site in soybean lipoxygenase-1 (SLO-1) through an atom-centered density matrix propagation (ADMP) treatment that facilitates the sampling of rare events

We present a computational methodology to sample rare events in large biological enzymes that may involve electronically polarizing, reactive processes. The approach includes simultaneous dynamical treatment of electronic and nuclear degrees of freedom, where contributions from the electronic portion are computed using hybrid density functional theory and the computational costs are reduced through a hybrid quantum mechanics/molecular mechanics (QM/MM) treatment. Thus, the paper involves a QM/MM dynamical treatment of rare events. The method is applied to probe the effect of the active site elements on the critical hydrogen transfer step in the soybean lipoxygenase-1 (SLO-1) catalyzed oxidation of linoleic acid. It is found that the dynamical fluctuations and associated flexibility of the active site are critical toward maintaining the electrostatics in the regime where the reactive process can occur smoothly. Physical constraints enforced to limit the active site flexibility are akin to mutations and, in the cases studied, have a detrimental effect on the electrostatic fluctuations, thus adversely affecting the hydrogen transfer process.

Conformational dynamics and stability of HP35 studied with 2D IR vibrational echoes

Two-dimensional infrared (2D IR) vibrational echo spectroscopy was used to measure the fast dynamics of two variants of chicken villin headpiece 35 (HP35). The CN of cyanophenylalanine residues inserted in the hydrophobic core were used as a vibrational probe. Experiments were performed on both singly (HP35-P) and doubly CN-labeled peptide (HP35-P(2)) within the wild-type sequence, as well as on HP-35 containing a singly labeled cyanophenylalanine and two norleucine mutations (HP35-P NleNle). There is a remarkable similarity between the dynamics measured in singly and doubly CN-labeled HP35, demonstrating that the presence of an additional CN vibrational probe does not significantly alter the dynamics of the small peptide. The substitution of two lysine residues by norleucines markedly improves the stability of HP35 by replacing charged with nonpolar residues, stabilizing the hydrophobic core. The results of the 2D IR experiments reveal that the dynamics of HP35-P are significantly faster than those of HP35-P NleNle. These observations suggest that the slower structural fluctuations in the hydrophobic core, indicating a more tightly structured core, may be an important contributing factor to HP35-P NleNle's increased stability.

Joint toxicity of alkoxyethanol mixtures: contribution of in silico applications

Exposure to chemicals occurs often as mixtures. Presented in this paper is information on alkoxyethanols and the impact they might have on human health in combination with some commonly found aliphatic and aromatic compounds. Our studies to evaluate the joint toxicity of these chemicals among themselves and in combination with other chemicals reveal a variety of possible outcomes depending on the exposure scenario. The interactions are predominantly based on metabolic pathways and are common among several solvents and organic compounds. Quantitative structure activity relationship (QSAR) analysis can be used with high confidence to identify chemicals that will interact to influence overall joint toxicity. Potential human exposure to a combination of alkoxyethanol, toluene and substituted benzenes may increase reproductive and developmental disease conditions. Inheritable gene alterations result in changes in the enzyme function in different subpopulations causing variations in quantity and/or quality of particular isoenzymes. These changes are responsible for differential metabolism of chemicals in species, genders, and life stages and are often the basis of a population's susceptibility. Unique genotypes introduced as a function of migration can alter the genetic makeup of any given population. Hence special consideration should be given to susceptible populations while conducting chemical health risk assessments.

Global and local molecular dynamics of a bacterial carboxylesterase provide insight into its catalytic mechanism

Carboxylesterases (CEs) are ubiquitous enzymes responsible for the detoxification of xenobiotics. In humans, substrates for these enzymes are far-ranging, and include the street drug heroin and the anticancer agent irinotecan. Hence, their ability to bind and metabolize substrates is of broad interest to biomedical science. In this study, we focused our attention on dynamic motions of a CE from B. subtilis (pnbCE), with emphasis on the question of what individual domains of the enzyme might contribute to its catalytic activity. We used a 10 ns all-atom molecular dynamics simulation, normal mode calculations, and enzyme kinetics to understand catalytic consequences of structural changes within this enzyme. Our results shed light on how molecular motions are coupled with catalysis. During molecular dynamics, we observed a distinct C-C bond rotation between two conformations of Glu310. Such a bond rotation would alternately facilitate and impede protonation of the active site His399 and act as a mechanism by which the enzyme alternates between its active and inactive conformation. Our normal mode results demonstrate that the distinct low-frequency motions of two loops in pnbCE, coil_5 and coil_21, are important in substrate conversion and seal the active site. Mutant CEs lacking these external loops show significantly reduced rates of substrate conversion, suggesting this sealing motion prevents escape of substrate. Overall, the results of our studies give new insight into the structure-function relationship of CEs and have implications for the entire family of α/β fold family of hydrolases, of which this CE is a member.

Mechanistic insights into the retaining glucosyl-3-phosphoglycerate synthase from mycobacteria

Considerable progress has been made in recent years in our understanding of the structural basis of glycosyl transfer. Yet the nature and relevance of the conformational changes associated with substrate recognition and catalysis remain poorly understood. We have focused on the glucosyl-3-phosphoglycerate synthase (GpgS), a "retaining" enzyme, that initiates the biosynthetic pathway of methylglucose lipopolysaccharides in mycobacteria. Evidence is provided that GpgS displays an unusually broad metal ion specificity for a GT-A enzyme, with Mg(2+), Mn(2+), Ca(2+), Co(2+), and Fe(2+) assisting catalysis. In the crystal structure of the apo-form of GpgS, we have observed that a flexible loop adopts a double conformation L(A) and L(I) in the active site of both monomers of the protein dimer. Notably, the L(A) loop geometry corresponds to an active conformation and is conserved in two other relevant states of the enzyme, namely the GpgS·metal·nucleotide sugar donor and the GpgS·metal·nucleotide·acceptor-bound complexes, indicating that GpgS is intrinsically in a catalytically active conformation. The crystal structure of GpgS in the presence of Mn(2+)·UDP·phosphoglyceric acid revealed an alternate conformation for the nucleotide sugar β-phosphate, which likely occurs upon sugar transfer. Structural, biochemical, and biophysical data point to a crucial role of the β-phosphate in donor and acceptor substrate binding and catalysis. Altogether, our experimental data suggest a model wherein the catalytic site is essentially preformed, with a few conformational changes of lateral chain residues as the protein proceeds along the catalytic cycle. This model of action may be applicable to a broad range of GT-A glycosyltransferases.

Active site hydrophobic residues impact hydrogen tunneling differently in a thermophilic alcohol dehydrogenase at optimal versus nonoptimal temperatures

A growing body of data suggests that protein motion plays an important role in enzyme catalysis. Two highly conserved hydrophobic active site residues in the cofactor-binding pocket of ht-ADH (Leu176 and V260) have been mutated to a series of hydrophobic side chains of smaller size, as well as one deletion mutant, L176Δ. Mutations decrease k(cat) and increase K(M)(NAD(+)). Most of the observed decreases in effects on k(cat) at pH 7.0 are due to an upward shift in the optimal pH for catalysis; a simple electrostatic model is invoked that relates the change in pK(a) to the distance between the positively charged nicotinamide ring and bound substrate. Structural modeling of the L176Δ and V260A variants indicates the development of a cavity behind the nicotinamide ring without any significant perturbation of the secondary structure of the enzyme relative to that of the wild type. Primary kinetic isotope effects (KIEs) are modestly increased for all mutants. Above the dynamical transition at 30 °C for ht-ADH [Kohen, A., et al. (1999) Nature 399, 496], the temperature dependence of the KIE is seen to increase with a decrease in side chain volume at positions 176 and 260. Additionally, the relative trends in the temperature dependence of the KIE above and below 30 °C appear to be reversed for the cofactor-binding pocket mutants in relation to wild-type protein. The aggregate results are interpreted in the context of a full tunneling model of enzymatic hydride transfer that incorporates both protein conformational sampling (preorganization) and active site optimization of tunneling (reorganization). The reduced temperature dependence of the KIE in the mutants below 30 °C indicates that at low temperatures, the enzyme adopts conformations refractory to donor-acceptor distance sampling.

Role of active site rigidity in activity: MD simulation and fluorescence study on a lipase mutant

Relationship between stability and activity of enzymes is maintained by underlying conformational flexibility. In thermophilic enzymes, a decrease in flexibility causes low enzyme activity while in less stable proteins such as mesophiles and psychrophiles, an increase in flexibility is associated with enhanced enzyme activity. Recently, we identified a mutant of a lipase whose stability and activity were enhanced simultaneously. In this work, we probed the conformational dynamics of the mutant and the wild type lipase, particularly flexibility of their active site using molecular dynamic simulations and time-resolved fluorescence techniques. In contrast to the earlier observations, our data show that active site of the mutant is more rigid than wild type enzyme. Further investigation suggests that this lipase needs minimal reorganization/flexibility of active site residues during its catalytic cycle. Molecular dynamic simulations suggest that catalytically competent active site geometry of the mutant is relatively more preserved than wild type lipase, which might have led to its higher enzyme activity. Our study implies that widely accepted positive correlation between conformation flexibility and enzyme activity need not be stringent and draws attention to the possibility that high enzyme activity can still be accomplished in a rigid active site and stable protein structures. This finding has a significant implication towards better understanding of involvement of dynamic motions in enzyme catalysis and enzyme engineering through mutations in active site.

Ribonuclease S dynamics measured using a nitrile label with 2D IR vibrational echo spectroscopy

A nitrile-labeled amino acid, p-cyanophenylalanine, is introduced near the active site of the semisynthetic enzyme ribonuclease S to serve as a probe of protein dynamics and fluctuations. Ribonuclease S is the limited proteolysis product of subtilisin acting on ribonuclease A, and consists of a small fragment including amino acids 1-20, the S-peptide, and a larger fragment including residues 21-124, the S-protein. A series of two-dimensional vibrational echo experiments performed on the nitrile-labeled S-peptide and the RNase S are described. The time-dependent changes in the two-dimensional infrared vibrational echo line shapes are analyzed using the center line slope method to obtain the frequency-frequency correlation function (FFCF). The observations show that the nitrile probe in the S-peptide has dynamics that are similar to, but faster than, those of the single amino acid p-cyanophenylalanine in water. In contrast, the dynamics of the nitrile label when the peptide is bound to form ribonuclease S are dominated by homogeneous dephasing (motionally narrowed) contributions with only a small contribution from very fast inhomogeneous structural dynamics. The results provide insights into the nature of the structural dynamics of the ribonuclease S complex. The equilibrium dynamics of the nitrile labeled S-peptide and the ribonuclease S complex are also investigated by molecular dynamics simulations. The experimentally determined FFCFs are compared to the FFCFs obtained from the molecular dynamics simulations, thereby testing the capacity of simulations to determine the amplitudes and time scales of protein structural fluctuations on fast time scales under thermal equilibrium conditions.

Resolving the complex role of enzyme conformational dynamics in catalytic function

Despite growing evidence suggesting the importance of enzyme conformational dynamics (ECD) in catalysis, a consensus on how precisely ECD influences the chemical step and reaction rates is yet to be reached. Here, we characterize ECD in Cyclophilin A, a well-studied peptidyl-prolyl cis-trans isomerase, using normal and accelerated, atomistic molecular dynamics simulations. Kinetics and free energy landscape of the isomerization reaction in solution and enzyme are explored in unconstrained simulations by allowing significantly lower torsional barriers, but in no way compromising the atomistic description of the system or the explicit solvent. We reveal that the reaction dynamics is intricately coupled to enzymatic motions that span multiple timescales and the enzyme modes are selected based on the energy barrier of the chemical step. We show that Kramers' rate theory can be used to present a clear rationale of how ECD affects the reaction dynamics and catalytic rates. The effects of ECD can be incorporated into the effective diffusion coefficient, which we estimate to be about ten times slower in enzyme than in solution. ECD thereby alters the preexponential factor, effectively impeding the rate enhancement. From our analyses, the trend observed for lower torsional barriers can be extrapolated to actual isomerization barriers, allowing successful prediction of the speedup in rates in the presence of CypA, which is in notable agreement with experimental estimates. Our results further reaffirm transition state stabilization as the main effect in enhancing chemical rates and provide a unified view of ECD's role in catalysis from an atomistic perspective.

Role of coupled dynamics in the catalytic activity of prokaryotic-like prolyl-tRNA synthetases

Prolyl-tRNA synthetases (ProRSs) have been shown to activate both cognate and some noncognate amino acids and attach them to specific tRNA(Pro) substrates. For example, alanine, which is smaller than cognate proline, is misactivated by Escherichia coli ProRS. Mischarged Ala-tRNA(Pro) is hydrolyzed by an editing domain (INS) that is distinct from the activation domain. It was previously shown that deletion of the INS greatly reduced cognate proline activation efficiency. In this study, experimental and computational approaches were used to test the hypothesis that deletion of the INS alters the internal protein dynamics leading to reduced catalytic function. Kinetic studies with two ProRS variants, G217A and E218A, revealed decreased amino acid activation efficiency. Molecular dynamics studies showed motional coupling between the INS and protein segments containing the catalytically important proline-binding loop (PBL, residues 199-206). In particular, the complete deletion of INS, as well as mutation of G217 or E218 to alanine, exhibited significant effects on the motion of the PBL. The presence of coupled dynamics between neighboring protein segments was also observed through in silico mutations and essential dynamics analysis. Altogether, this study demonstrates that structural elements at the editing domain-activation domain interface participate in coupled motions that facilitate amino acid binding and catalysis by bacterial ProRSs, which may explain why truncated or defunct editing domains have been maintained in some systems, despite the lack of catalytic activity.

Influence of myeloperoxidase-catalyzing reaction on the binding between myeloperoxidase and anti-myeloperoxidase antibodies

In the current study, whether myeloperoxidase (MPO)-catalyzing reaction could influence the antigenicity of MPO was investigated. Hypochlorite acid, the main product of the catalytic reaction, could lower the binding between MPO-antineutrophil cytoplasmic antibodies (ANCA) and MPO when the available chlorine was higher than 0.031×10(-3) g/l. After MPO-catalyzing reaction with H(2)O(2) lower than 0.469 g/l, the binding level between MPO-ANCAcontaining plasma and MPO increased slightly. The peak binding level was 1.135 ± 0.205 (expressed by the absorbance value at 405 nm). However, with the existence of hydrogen donor (o-phenylenediamine) in the reaction system, the peak binding level between MPO-ANCA-containing plasma and post-catalyzing MPO was significantly higher (1.367 ± 0.321 vs 1.135 ± 0.205, p = 0.023). Moreover, at the approximately physical concentration of H(2)O(2) (0.02 g/l), MPO-ANCA exhibited higher titer to post-catalyzing MPO than to pre-catalyzing MPO (3.91 ± 0.84 vs 3.57 ± 0.84, p < 0.001, expressed as the lgT). These data demonstrated that MPO-catalyzing reaction could potentially increase the antigenicity of MPO.

Taking Ockham's razor to enzyme dynamics and catalysis

The role of protein dynamics in enzyme catalysis is a matter of intense current debate. Enzyme-catalysed reactions that involve significant quantum tunnelling can give rise to experimental kinetic isotope effects with complex temperature dependences, and it has been suggested that standard statistical rate theories, such as transition-state theory, are inadequate for their explanation. Here we introduce aspects of transition-state theory relevant to the study of enzyme reactivity, taking cues from chemical kinetics and dynamics studies of small molecules in the gas phase and in solution--where breakdowns of statistical theories have received significant attention and their origins are relatively better understood. We discuss recent theoretical approaches to understanding enzyme activity and then show how experimental observations for a number of enzymes may be reproduced using a transition-state-theory framework with physically reasonable parameters. Essential to this simple model is the inclusion of multiple conformations with different reactivity.

Primary alcohols activate human TRPA1 channel in a carbon chain length-dependent manner

Transient receptor potential ankyrin 1 (TRPA1) is a calcium-permeable non-selective cation channel that is mainly expressed in primary nociceptive neurons. TRPA1 is activated by a variety of noxious stimuli, including cold temperatures, pungent compounds such as mustard oil and cinnamaldehyde, and intracellular alkalization. Here, we show that primary alcohols, which have been reported to cause skin, eye or nasal irritation, activate human TRPA1 (hTRPA1). We measured intracellular Ca(2+) changes in HEK293 cells expressing hTRPA1 induced by 1 mM primary alcohols. Higher alcohols (1-butanol to 1-octanol) showed Ca(2+) increases proportional to the carbon chain length. In whole-cell patch-clamp recordings, higher alcohols (1-hexanol to 1-octanol) activated hTRPA1 and the potency increased with the carbon chain length. Higher alcohols evoked single-channel opening of hTRPA1 in an inside-out configuration. In addition, cysteine at 665 in the N terminus and histidine at 983 in the C terminus were important for hTRPA1 activation by primary alcohols. Furthermore, straight-chain secondary alcohols increased intracellular Ca(2+) concentrations in HEK293 cells expressing hTRPA1, and both primary and secondary alcohols showed hTRPA1 activation activities that correlated highly with their octanol/water partition coefficients. On the other hand, mouse TRPA1 did not show a strong response to 1-hexanol or 1-octanol, nor did these alcohols evoke significant pain in mice. We conclude that primary and secondary alcohols activate hTRPA1 in a carbon chain length-dependent manner. TRPA1 could be a sensor of alcohols inducing skin, eye and nasal irritation in human.

Dynamic mechanism of proton transfer in mannitol 2-dehydrogenase from Pseudomonas fluorescens: mobile GLU292 controls proton relay through a water channel that connects the active site with bulk solvent

The active site of mannitol 2-dehydrogenase from Pseudomonas fluorescens (PfM2DH) is connected with bulk solvent through a narrow protein channel that shows structural resemblance to proton channels utilized by redox-driven proton pumps. A key element of the PfM2DH channel is the "mobile" Glu(292), which was seen crystallographically to adopt distinct positions up and down the channel. It was suggested that the "down → up" conformational change of Glu(292) could play a proton relay function in enzymatic catalysis, through direct proton shuttling by the Glu or because the channel is opened for water molecules forming a chain along which the protons flow. We report evidence from site-directed mutagenesis (Glu(292) → Ala) substantiated by data from molecular dynamics simulations that support a role for Glu(292) as a gate in a water chain (von Grotthuss-type) mechanism of proton translocation. Occupancy of the up and down position of Glu(292) is influenced by the bonding and charge state of the catalytic acid base Lys(295), suggesting that channel opening/closing motions of the Glu are synchronized to the reaction progress. Removal of gatekeeper control in the E292A mutant resulted in a selective, up to 120-fold slowing down of microscopic steps immediately preceding catalytic oxidation of mannitol, consistent with the notion that formation of the productive enzyme-NAD(+)-mannitol complex is promoted by a corresponding position change of Glu(292), which at physiological pH is associated with obligatory deprotonation of Lys(295) to solvent. These results underscore the important role of conformational dynamics in the proton transfer steps of alcohol dehydrogenase catalysis.

3-picolyl azide adenine dinucleotide as a probe of femtosecond to picosecond enzyme dynamics

Functionally relevant femtosecond to picosecond dynamics in enzyme active sites can be difficult to measure because of a lack of spectroscopic probes that can be located in the active site without altering the behavior of the enzyme. We have developed a new NAD(+) analog 3-Picolyl Azide Adenine Dinucleotide (PAAD(+)), which has the potential to be a general spectroscopic probe for NAD-dependent enzymes. This analog is stable and binds in the active site of a typical NAD-dependent enzyme formate dehydrogenase (FDH) with characteristics similar to those of natural NAD(+). It has an isolated infrared transition with high molar absorptivity that makes it suitable for observing enzyme dynamics using 2D IR spectroscopy. 2D IR experiments show that in aqueous solution, the analog undergoes complete spectral diffusion within hundreds of femtoseconds consistent with the water hydrogen bonding dynamics that would be expected. When bound to FDH in a binary complex, it shows picosecond fluctuations and a large static offset, consistent with previous studies of the binary complexes of this enzyme. These results show that PAAD(+) is an excellent probe of local dynamics and that it should be a general tool for probing the dynamics of a wide range of NAD-dependent enzymes.

Structural position correlation analysis (SPCA) for protein family

Background

The proteins in a family, which perform the similar biological functions, may have very different amino acid composition, but they must share the similar 3D structures, and keep a stable central region. In the conservative structure region similar biological functions are performed by two or three catalytic residues with the collaboration of several functional residues at key positions. Communication signals are conducted in a position network, adjusting the biological functions in the protein family.

Methodology

A computational approach, namely structural position correlation analysis (SPCA), is developed to analyze the correlation relationship between structural segments (or positions). The basic hypothesis of SPCA is that in a protein family the structural conservation is more important than the sequence conservation, and the local structural changes may contain information of biology functional evolution. A standard protein P⁽⁰⁾ is defined in a protein family, which consists of the most-frequent amino acids and takes the average structure of the protein family. The foundational variables of SPCA is the structural position displacements between the standard protein P⁽⁰⁾ and individual proteins P_i of the family. The structural positions are organized as segments, which are the stable units in structural displacements of the protein family. The biological function differences of protein members are determined by the position structural displacements of individual protein P_i to the standard protein P⁽⁰⁾. Correlation analysis is used to analyze the communication network among segments.

Conclusions

The structural position correlation analysis (SPCA) is able to find the correlation relationship among the structural segments (or positions) in a protein family, which cannot be detected by the amino acid sequence and frequency-based methods. The functional communication network among the structural segments (or positions) in protein family, revealed by SPCA approach, well illustrate the distantly allosteric interactions, and contains valuable information for protein engineering study.

Flexibility, diversity, and cooperativity: pillars of enzyme catalysis

This brief review discusses our current understanding of the molecular basis of enzyme catalysis. A historical development is presented, beginning with steady state kinetics and progressing through modern fast reaction methods, nuclear magnetic resonance, and single-molecule fluorescence techniques. Experimental results are summarized for ribonuclease, aspartate aminotransferase, and especially dihydrofolate reductase (DHFR). Multiple intermediates, multiple conformations, and cooperative conformational changes are shown to be an essential part of virtually all enzyme mechanisms. In the case of DHFR, theoretical investigations have provided detailed information about the movement of atoms within the enzyme-substrate complex as the reaction proceeds along the collective reaction coordinate for hydride transfer. A general mechanism is presented for enzyme catalysis that includes multiple intermediates and a complex, multidimensional standard free energy surface. Protein flexibility, diverse protein conformations, and cooperative conformational changes are important features of this model.

Conformational variability of organophosphorus hydrolase upon soman and paraoxon binding

The bacterial enzyme organophosphorus hydrolase (OPH) exhibits both catalytic and substrate promiscuity. It hydrolyzes bonds in a variety of phosphotriester (P-O), phosphonothioate (P-S), phosphofluoridate (P-F), and phosphonocyanate (F-CN) compounds. However, its catalytic efficiency varies markedly for different substrates, limiting the broad-range application of OPH as catalyst in the bioremediation of pesticides and chemical war agents. In the present study, pK(a) calculations and multiple explicit-solvent molecular dynamics (MD) simulations were performed to characterize and contrast the structural dynamics of OPH bound to two substrates hydrolyzed with very distinct catalytic efficiencies: the nerve agent soman (O-pinacolylmethylphosphonofluoridate) and the pesticide paraoxon (diethyl p-nitrophenyl phosphate). pK(a) calculations for the substrate-bound and unbound enzyme showed a significant pK(a) shift from standard values (ΔpK(a) = ±3 units) for residues His254 and Arg275. MD simulations of protonated His254 revealed a dynamic hydrogen bond network connecting the catalytic residue Asp301 via His254 to Asp232, Asp233, Arg275, and Asp235, and is consistent with a previously postulated proton relay mechanism to ferry protons away from the active site with substrates that do not require activation of the leaving group. Hydrogen bonds between Asp301 and His254 were persistent in the OPH-paraoxon complex but not in the OPH-soman one, suggesting a potential role for such interaction in the more efficient hydrolysis of paraoxon over soman by OPH. These results are in line with previous mutational studies of residue His254, which led to an increase of the catalytic efficiency of OPH over soman yet decreased its efficiency for paraoxon. In addition, comparative analysis of the molecular trajectories for OPH bound to soman and paraoxon suggests that binding of the latter facilitates the conformational transition of OPH from the open to the closed substate promoting a tighter binding of paraoxon.

Link between allosteric signal transduction and functional dynamics in a multisubunit enzyme: S-adenosylhomocysteine hydrolase

S-adenosylhomocysteine hydrolase (SAHH), a cellular enzyme that plays a key role in methylation reactions including those required for maturation of viral mRNA, is an important drug target in the discovery of antiviral agents. While targeting the active site is a straightforward strategy of enzyme inhibition, evidence of allosteric modulation of active site in many enzymes underscores the molecular origin of signal transduction. Information of co-evolving sequences in SAHH family and the key residues for functional dynamics that can be identified using native topology of the enzyme provide glimpses into how the allosteric signaling network, dispersed over the molecular structure, coordinates intra- and intersubunit conformational dynamics. To study the link between the allosteric communication and functional dynamics of SAHHs, we performed Brownian dynamics simulations by building a coarse-grained model based on the holo and ligand-bound structures. The simulations of ligand-induced transition revealed that the signal of intrasubunit closure dynamics is transmitted to form intersubunit contacts, which in turn invoke a precise alignment of active site, followed by the dimer-dimer rotation that compacts the whole tetrameric structure. Further analyses of SAHH dynamics associated with ligand binding provided evidence of both induced fit and population shift mechanisms and also showed that the transition-state ensemble is akin to the ligand-bound state. Besides the formation of enzyme-ligand contacts at the active site, the allosteric couplings from the residues distal to the active site are vital to the enzymatic function.

Mass-dependent bond vibrational dynamics influence catalysis by HIV-1 protease

Protein motions that occur on the microsecond to millisecond time scale have been linked to enzymatic rates observed for catalytic turnovers, but not to transition-state barrier crossing. It has been hypothesized that enzyme motions on the femtosecond time scale of bond vibrations play a role in transition state formation. Here, we perturb femtosecond motion by substituting all nonexchangeable carbon, nitrogen, and hydrogen atoms with (13)C, (15)N, and (2)H and observe the catalytic effects in HIV-1 protease. According to the Born-Oppenheimer approximation, isotopic substitution alters vibrational frequency with unchanged electrostatic properties. With the use of a fluorescent peptide to report on multiple steps in the reaction, we observe significantly reduced rates in the heavy enzyme relative to the light enzyme. A possible interpretation of our results is that there exists a dynamic link between mass-dependent bond vibrations of the enzyme and events in the reaction coordinate.

Evolutionarily conserved linkage between enzyme fold, flexibility, and catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein's designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design.

Femtosecond dynamics coupled to chemical barrier crossing in a Born-Oppenheimer enzyme

Contributions of fast (femtosecond) dynamic motion to barrier crossing at enzyme catalytic sites is in dispute. Human purine nucleoside phosphorylase (PNP) forms a ribocation-like transition state in the phosphorolysis of purine nucleosides and fast protein motions have been proposed to participate in barrier crossing. In the present study, (13)C-, (15)N-, (2)H-labeled human PNP (heavy PNP) was expressed, purified to homogeneity, and shown to exhibit a 9.9% increase in molecular mass relative to its unlabeled counterpart (light PNP). Kinetic isotope effects and steady-state kinetic parameters were indistinguishable for both enzymes, indicating that transition-state structure, equilibrium binding steps, and the rate of product release were not affected by increased protein mass. Single-turnover rate constants were slowed for heavy PNP, demonstrating reduced probability of chemical barrier crossing from enzyme-bound substrates to enzyme-bound products. In a second, independent method to probe barrier crossing, heavy PNP exhibited decreased forward commitment factors, also revealing mass-dependent decreased probability for barrier crossing. Increased atomic mass in human PNP alters bond vibrational modes on the femtosecond time scale and reduces on-enzyme chemical barrier crossing. This study demonstrates coupling of enzymatic bond vibrations on the femtosecond time scale to barrier crossing.

Enzymatic transition states, transition-state analogs, dynamics, thermodynamics, and lifetimes

Experimental analysis of enzymatic transition-state structures uses kinetic isotope effects (KIEs) to report on bonding and geometry differences between reactants and the transition state. Computational correlation of experimental values with chemical models permits three-dimensional geometric and electrostatic assignment of transition states formed at enzymatic catalytic sites. The combination of experimental and computational access to transition-state information permits (a) the design of transition-state analogs as powerful enzymatic inhibitors, (b) exploration of protein features linked to transition-state structure, (c) analysis of ensemble atomic motions involved in achieving the transition state, (d) transition-state lifetimes, and (e) separation of ground-state (Michaelis complexes) from transition-state effects. Transition-state analogs with picomolar dissociation constants have been achieved for several enzymatic targets. Transition states of closely related isozymes indicate that the protein's dynamic architecture is linked to transition-state structure. Fast dynamic motions in catalytic sites are linked to transition-state generation. Enzymatic transition states have lifetimes of femtoseconds, the lifetime of bond vibrations. Binding isotope effects (BIEs) reveal relative reactant and transition-state analog binding distortion for comparison with actual transition states.

Characterization of a novel flooding stress-responsive alcohol dehydrogenase expressed in soybean roots

Alcohol dehydrogenase (Adh) is the key enzyme in alcohol fermentation. We analyzed Adh expression in order to clarify the role of Adh of soybeans (Glycine max) to flooding stress. Proteome analysis confirmed that expression of Adh is significantly upregulated in 4-day-old soybean seedlings subjected to 2 days of flooding. Southern hybridization analysis and soybean genome database search revealed that soybean has at least 6 Adh genes. The GmAdh2 gene that responded to flooding was isolated from soybean cultivar Enrei. Adh2 expression was markedly increased 6 h after flooding and decreased 24 h after floodwater drainage. In situ hybridization and Western blot indicated that flooding strongly induces Adh2 expression in RNA and protein levels in the root apical meristem. Osmotic, cold, or drought stress did not induce expression of Adh2. These results indicate that Adh2 is a flooding-response specific soybean gene expressed in root tissue.

Dark progression reveals slow timescales for radiation damage between T = 180 and 240 K

Can radiation damage to protein crystals be `outrun' by collecting a structural data set before damage is manifested? Recent experiments using ultra-intense pulses from a free-electron laser show that the answer is yes. Here, evidence is presented that significant reductions in global damage at temperatures above 200 K may be possible using conventional X-ray sources and current or soon-to-be available detectors. Specifically, `dark progression' (an increase in damage with time after the X-rays have been turned off) was observed at temperatures between 180 and 240 K and on timescales from 200 to 1200 s. This allowed estimation of the temperature-dependent timescale for damage. The rate of dark progression is consistent with an Arrhenius law with an activation energy of 14 kJ mol(-1). This is comparable to the activation energy for the solvent-coupled diffusive damage processes responsible for the rapid increase in radiation sensitivity as crystals are warmed above the glass transition near 200 K. Analysis suggests that at T = 300 K data-collection times of the order of 1 s (and longer at lower temperatures) may allow significant reductions in global radiation damage, facilitating structure solution on crystals with liquid solvent. No dark progression was observed below T = 180 K, indicating that no important damage process is slowed through this timescale window in this temperature range.

Self-assembled enzymatic monolayer directly bound to a gold surface: activity and molecular recognition force spectroscopy studies

Escherichia coli dihydrofolate reductase (ecDHFR) has one surface cysteine, C152, located opposite and distal to the active site. Here, we show that the enzyme spontaneously assembles on an ultraflat gold surface as a homogeneous, covalently bound monolayer. Surprisingly, the activity of the gold-immobilized ecDHFR as measured by radiographic analysis was found to be similar to that of the free enzyme in solution. Molecular recognition force spectroscopy was used to study the dissociation forces involved in the rupture of AFM probe-tethered methotrexate (MTX, a tight-binding inhibitor of DHFR) from the gold-immobilized enzyme. Treatment of the ecDHFR monolayer with free MTX diminished the interaction of the functionalized tip with the surface, suggesting that the interaction was indeed active-site specific. These findings demonstrate the viability of a simple and direct enzymatic surface-functionalization without the use of spacers, thus, opening the door to further applications in the area of biomacromolecular force spectroscopy.