Kinetic proofreading models for cell signaling predict ways to escape kinetic proofreading

In the context of cell signaling, kinetic proofreading was introduced to explain how cells can discriminate among ligands based on a kinetic parameter, the ligand-receptor dissociation rate constant. In the kinetic proofreading model of cell signaling, responses occur only when a bound receptor undergoes a complete series of modifications. If the ligand dissociates prematurely, the receptor returns to its basal state and signaling is frustrated. We extend the model to deal with systems where aggregation of receptors is essential to signal transduction, and present a version of the model for systems where signaling depends on an extrinsic kinase. We also investigate the kinetics of signaling molecules, “messengers,” that are generated by aggregated receptors but do not remain associated with the receptor complex. We show that the extended model predicts modes of signaling that exhibit kinetic discrimination for some range of parameters but for other parameter values show little or no discrimination and thus escape kinetic proofreading. We compare model predictions with experimental data.

Excited-state dynamics in photosystem II: Insights from the x-ray crystal structure

The heart of oxygenic photosynthesis is photosystem II (PSII), a multisubunit protein complex that uses solar energy to drive the splitting of water and production of molecular oxygen. The effectiveness of the photochemical reaction center of PSII depends on the efficient transfer of excitation energy from the surrounding antenna chlorophylls. A kinetic model for PSII, based on the x-ray crystal structure coordinates of 37 antenna and reaction center pigment molecules, allows us to map the major energy transfer routes from the antenna chlorophylls to the reaction center chromophores. The model shows that energy transfer to the reaction center is slow compared with the rate of primary electron transport and depends on a few bridging chlorophyll molecules. This unexpected energetic isolation of the reaction center in PSII is similar to that found in the bacterial photosystem, conflicts with the established view of the photophysics of PSII, and may be a functional requirement for primary photochemistry in photosynthesis. In addition, the model predicts a value for the intrinsic photochemical rate constant that is 4 times that found in bacterial reaction centers.

A thermodynamic coupling mechanism for GroEL-mediated unfolding.

Chaperonins prevent the aggregation of partially folded or misfolded forms of a protein and, thus, keep it competent for productive folding. It was suggested that GroEL, the chaperonin of Escherichia coli, exerts this function 1 unfolding such intermediates, presumably in a catalytic fashion. We investigated the kinetic mechanism of GroEL-induced protein unfolding by using a reduced and carbamidomethylated variant of RNase T1, RCAM-T1, as a substrate. RCAM-T1 cannot fold to completion, because the two disulfide bonds are missing, and it is, thus, a good model for long-lived folding intermediates. RCAM-T1 unfolds when GroEL is added, but GroEL does not change the microscopic rate constant of unfolding, ruling out that it catalyzes unfolding. GroEL unfolds RCAM-T1 because it binds with high affinity to the unfolded form of the protein and thereby shifts the overall equilibrium toward the unfolded state. GroEL can unfold a partially folded or misfolded intermediate by this thermodynamic coupling mechanism when the Gibbs free energy of the binding to GroEL is larger than the conformational stability of the intermediate and when the rate of its unfolding is high.

Conformational dependence of electron transfer across de novo designed metalloproteins.

Flash photolysis and pulse radiolysis measurements demonstrate a conformational dependence of electron transfer rates across a 16-mer helical bundle (three-helix metalloprotein) modified with a capping CoIII(bipyridine)3 electron acceptor at the N terminus and a 1-ethyl-1'-ethyl-4,4'- bipyridinium donor at the C terminus. For the CoIII(peptide)3-1-ethyl-1'-ethyl-4,4'-bipyridinium maquettes, the observed transfer is a first order, intramolecular process, independent of peptide concentration or laser pulse energy. In the presence of 6 M urea, the random coil bundle (approximately 0% helicity) has an observed electron transfer rate constant of kobs = 900 +/- 100 s-1. In the presence of 25% trifluoroethanol (TFE), the helicity of the peptide is 80% and the kobs increases to 2000 +/- 200 s-1. Moreover, the increase in the rate constant in TFE is consistent with the observed decrease in donor-acceptor distance in this solvent. Such bifunctional systems provide a class of molecules for testing the effects of conformation on electron transfer in proteins and peptides.

Kinetics of photo-induced electron transfer from high-potential iron-sulfur protein to the photosynthetic reaction center of the purple phototroph Rhodoferax fermentans.

The kinetics of photo-induced electrontransfer from high-potential iron-sulfur protein (HiPIP) to the photosynthetic reaction center (RC) of the purple phototroph Rhodoferarfermentans were studied. The rapid photooxidation of heme c-556 belonging to RC is followed, in the presence of HiPIP, by a slower reduction having a second-order rate constant of 4.8 x 10(7) M(-1) x s(-1). The limiting value of kobs at high HiPIP concentration is 95 s(-1). The amplitude of this slow process decreases with increasing HiPIP concentration. The amplitude of a faster phase, observed at 556 and 425 nm and involving heme c-556 reduction, increases proportionately. The rate constant of this fast phase, determined at 425 and 556 nm, is approximately 3 x 10(5) s(-1). This value is not dependent on HiPIP concentration, indicating that it is related to a first-order process. These observations are interpreted as evidence for the formation of a HiPIP-RC complex prior to the excitation flash, having a dissociation constant of -2.5 microM. The fast phase is absent at high ionic strength, indicating that the complex involves mainly electrostatic interactions. The ionic strength dependence of kobs for the slow phase yields a second-order rate constant at infinite ionic strength of 5.4 x 10(6) M(-1) x s(-1) and an electrostatic interaction energy of -2.1 kcal/mol (1 cal = 4.184 J). We conclude that Rhodoferar fermentans HiPIP is a very effective electron donor to the photosynthetic RC.

The nucleation of receptor-mediated endocytosis.

A theory of the mechanical origins of receptor-mediated endocytosis shows that a spontaneous membrane complex formation can provide the stimulus for a local membrane motion toward the cytosol. This motion is identified with a nucleation stage of receptor-mediated endocytosis. When membrane complexes cluster, membrane deformation is predicted to be most rapid. The rate of growth of membrane depressions depends upon the relative rates of approach of aqueous cytosolic and extracellular fluids toward the cell membrane. With cytosolic and extracellular media characterized by apparent viscosities, the rate of growth of membrane depressions is predicted to increase as the extracellular viscosity nears the apparent viscosity of the cytosol and then to decrease when the extracellular viscosity exceeds that of the cytosol. To determine whether these trends would be apparent in the overall endocytosis rate constant, an experimental study of transferrin-mediated endocytosis in two different cell lines was conducted. The experimental results reveal the same dependence of internalization rate on extracellular viscosity as predicted by the theory. These and other comparisons with experimental data suggest that the nucleation stage of receptor-mediated endocytosis is important in the overall endocytosis process.

Phosphinate analogs of D-, D-dipeptides: slow-binding inhibition and proteolysis protection of VanX, a D-, D-dipeptidase required for vancomycin resistance in Enterococcus faecium.

VanX is a D-Ala-D-Ala dipeptidase that is essential for vancomycin resistance in Enterococcus faecium. Contrary to most proteases and peptidases, it prefers to hydrolyze the amino substrate but not the related kinetically and thermodynamically more favorable ester substrate D-Ala-D-lactate. The enzymatic activity of VanX was previously found to be inhibited by the phosphinate analogs of the proposed tetrahedral intermediate for hydrolysis of D-Ala-D-Ala. Here we report that such phosphinates are slow-binding inhibitors. D-3-[(1-Aminoethyl)phosphinyl]-D-2-methylpropionic acid I showed a time-dependent onset of inhibition of VanX and a time-dependent return to uninhibited steady-state rates upon dilution of the enzyme/inhibitor mixture. The initial inhibition constant Ki after immediate addition of VanX to phosphinate I to form the E-I complex is 1.5 microM but is then lowered by a relatively slow isomerization step to a second complex, E-I*, with a final K*i of 0.47 microM. This slow-binding inhibition reflects a Km/K*i ratio of 2900:1. The rate constant for the slow dissociation of complex E-I* is 0.24 min-1. A phosphinate analog with an ethyl group replacing what would be the side chain of the second D-alanyl residue in the normal tetrahedral adduct gives a K*i value of 90 nM. Partial proteolysis of VanX reveals two protease-sensitive loop regions that are protected by the intermediate analog phosphinate, indicating that they may be part of the VanX active site.

3T3 fibroblasts transfected with a cDNA for mitochondrial aspartate aminotransferase express plasma membrane fatty acid-binding protein and saturable fatty acid uptake.

To explore the relationship between mitochondrial aspartate aminotransferase (mAspAT; EC 2.6.1.1) and plasma membrane fatty acid-binding protein (FABPpm) and their role in cellular fatty acid uptake, 3T3 fibroblasts were cotransfected with plasmid pMAAT2, containing a full-length mAspAT cDNA downstream of a Zn(2+)-inducible metallothionein promoter, and pFR400, which conveys methotrexate resistance. Transfectants were selected in methotrexate, cloned, and exposed to increasing methotrexate concentrations to induce gene amplification. Stably transfected clones were characterized by Southern blotting; those with highest copy numbers of pFR400 alone (pFR400) or pFR400 and pMAAT2 (pFR400/pMAAT2) were expanded for further study. [3H]Oleate uptake was measured in medium containing 500 microM bovine serum albumin and 125-1000 microM total oleate (unbound oleate, 18-420 nM) and consisted of saturable and nonsaturable components. pFR400/pMAAT2 cells exhibited no increase in the rate constant for nonsaturable oleate uptake or in the uptake rate of [14C]octanoate under any conditions. By contrast, Vmax (fmol/sec per 50,000 cells) of the saturable oleate uptake component increased 3.5-fold in pFR400/pMAAT2 cells compared to pFR400, with a further 3.2-fold increase in the presence of Zn2+. Zn2+ had no effect in pFR400 controls (P > 0.5). The overall increase in Vmax between pFR400 and pFR400/pMAAT2 in the presence of Zn2+ was 10.4-fold (P < 0.01) and was highly correlated (r = 0.99) with expression of FABPpm in plasma membranes as determined by Western blotting. Neither untransfected 3T3 nor pFR400 cells expressed cell surface FABPpm detectable by immunofluorescence. By contrast, plasma membrane immunofluorescence was detected in pFR400/pMAAT2 cells, especially if cultured in 100 microM Zn2+. The data support the dual hypotheses that mAspAT and FABPpm are identical and mediate saturable long-chain free fatty acid uptake.

Negative activation enthalpies in the kinetics of protein folding.

Although the rates of chemical reactions become faster with increasing temperature, the converse may be observed with protein-folding reactions. The rate constant for folding initially increases with temperature, goes through a maximum, and then decreases. The activation enthalpy is thus highly temperature dependent because of a large change in specific heat (delta Cp). Such a delta Cp term is usually presumed to be a consequence of a large decrease in exposure of hydrophobic surfaces to water as the reaction proceeds from the denatured state to the transition state for folding: the hydrophobic side chains are surrounded by "icebergs" of water that melt with increasing temperature, thus making a large contribution to the Cp of the denatured state and a smaller one to the more compact transition state. The rate could also be affected by temperature-induced changes in the conformational population of the ground state: the heat required for the progressive melting of residual structure in the denatured state will contribute to delta Cp. By examining two proteins with different refolding mechanisms, we are able to find both of these two processes; barley chymotrypsin inhibitor 2, which refolds from a highly unfolded state, fits well to a hydrophobic interaction model with a constant delta Cp of activation, whereas barnase, which refolds from a more structured denatured state, deviates from this ideal behavior.

Submillisecond folding of monomeric lambda repressor.

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.

The folding of GroEL-bound barnase as a model for chaperonin-mediated protein folding.

We have analyzed the pathway of folding of barnase bound to GroEL to resolve the controversy of whether proteins can fold while bound to chaperonins (GroEL or Cpn60) or fold only after their release into solution. Four phases in the folding were detected by rapid-reaction kinetic measurements of the intrinsic fluorescence of both wild type and barnase mutants. The phases were assigned from their rate laws, sensitivity to mutations, and correspondence to regain of catalytic activity. At high ratios of denatured barnase to GroEL, 4 mol of barnase rapidly bind per 14-mer of GroEL. At high ratios of GroEL to barnase, 1 mol of barnase binds with a rate constant of 3.5 x 10(7) s-1.M-1. This molecule then refolds with a low rate constant that changes on mutation in parallel with the rate constant for the folding in solution. This rate constant corresponds to the regain of the overall catalytic activity of barnase and increases 15-fold on the addition of ATP to a physiologically relevant value of approximately 0.4 s-1. The multiply bound molecules of barnase that are present at high ratios of GroEL to barnase fold with a rate constant that is also sensitive to mutation but is 10 times higher. If the 110-residue barnase can fold when bound to GroEL and many moles can bind simultaneously, then smaller parts of large proteins should be able to fold while bound.

Basis of guanylate cyclase activation by carbon monoxide.

Kinetics of CO association with guanylate cyclase [GTP pyrophosphate-lyase (cyclizing), EC 4.6.1.2] and dissociation from carboxy guanylate cyclase have been studied at pH 7.5 by flash photolysis, yielding rate constants at 23 degrees C of 1.2 +/- 0.1 x 10(5) M-1.sec-1 and 28 +/- 2 sec-1, respectively. While the CO combination rate constant is the same as for the T state of hemoglobin, the CO dissociation rate constant is much higher than expected for a six-coordinate carboxyheme protein; yet the absorption spectrum is indicative of a six-coordinate heme. The two observations are reconciled by a reaction mechanism in which CO dissociation proceeds via a five-coordinate intermediate. This intermediate is structurally very similar to the five-coordinate nitrosyl heme derivative of guanylate cyclase and is presumably responsible for the observed 4-fold activation of guanylate cyclase by CO. Thus, we provide a model that explains enzyme activities of the nitrosyl and carboxy forms of the enzyme on the basis of a common mechanism.

Kinetics and mechanism of polyamide ("peptide") nucleic acid binding to duplex DNA.

To elucidate the mechanism of recognition of double-stranded DNA (dsDNA) by homopyrimidine polyamide ("peptide") nucleic acid (PNA) leading to the strand-displacement, the kinetics of the sequence-specific PNA/DNA binding have been studied. The binding was monitored with time by the gel retardation and nuclease S1 cleavage assays. The experimental kinetic curves obey pseudo-first-order kinetics and the dependence of the pseudo-first-order rate constant, kps, on PNA concentration, P, obeys a power law kps approximately P gamma with 2 < gamma < 3. The kps values for binding of decamer PNA to dsDNA target sites with one mismatch are hundreds of times slower than for the correct site. A detailed kinetic scheme for PNA/DNA binding is proposed that includes two major steps of the reaction of strand invasion: (i) a transient partial opening of the PNA binding site on dsDNA and incorporation of one PNA molecule with the formation of an intermediate PNA/DNA duplex and (ii) formation of a very stable PNA2/DNA triplex. A simple theoretical treatment of the proposed kinetic scheme is performed. The interpretation of our experimental data in the framework of the proposed kinetic scheme leads to the following conclusions. The sequence specificity of the recognition is essentially provided at the "search" step of the process, which consists in the highly reversible transient formation of duplex between one PNA molecule and the complementary strand of duplex DNA while the other DNA strand is displaced. This search step is followed by virtually irreversible "locking" step via PNA2/DNA triplex formation. The proposed mechanism explains how the binding of homopyrimidine PNA to dsDNA meets two apparently mutually contradictory features: high sequence specificity of binding and remarkable stability of both correct and mismatched PNA/DNA complexes.

Investigação do mecanismo cinético da reação de redução de oxigênio em solventes não aquosos

O aumento no consumo energético e a crescente preocupação ambiental frente à emissão de gases poluentes criam um apelo mundial favorável para pesquisas de novas tecnologias não poluentes de fontes de energia. Baterias recarregáveis de lítio-ar em solventes não aquosos possuem uma alta densidade de energia teórica (5200 Wh kg-1), o que as tornam promissoras para aplicação em dispositivos estacionários e em veículos elétricos. Entretanto, muitos problemas relacionados ao cátodo necessitam ser contornados para permitir a aplicação desta tecnologia, por exemplo, a baixa reversibilidade das reações, baixa potência e instabilidades dos materiais empregados nos eletrodos e dos solventes eletrolíticos. Assim, neste trabalho um modelo cinético foi empregado para os dados experimentais de espectroscopia de impedância eletroquímica, para a obtenção das constantes cinéticas das etapas elementares do mecanismo da reação de redução de oxigênio (RRO), o que permitiu investigar a influência de parâmetros como o tipo e tamanho de partícula do eletrocatalisador, o papel do solvente utilizado na RRO e compreender melhor as reações ocorridas no cátodo dessa bateria. A investigação inicial se deu com a utilização de sistemas menos complexos como uma folha de platina ou eletrodo de carbono vítreo como eletrodos de trabalho em 1,2-dimetoxietano (DME)/perclorato de lítio (LiClO4). A seguir, sistemas complexos com a presença de nanopartículas de carbono favoreceu o processo de adsorção das moléculas de oxigênio e aumentou ligeiramente (uma ordem de magnitude) a etapa de formação de superóxido de lítio (etapa determinante de reação) quando comparada com os eletrodos de platina e carbono vítreo, atribuída à presença dos grupos laterais mediando à transferência eletrônica para as moléculas de oxigênio. No entanto, foi observada uma rápida passivação da superfície eletrocatalítica através da formação de filmes finos de Li2O2 e Li2CO3 aumentando o sobrepotencial da bateria durante a carga (diferença de potencial entre a carga e descarga > 1 V). Adicionalmente, a incorporação das nanopartículas de platina (Ptnp), ao invés da folha de platina, resultou no aumento da constante cinética da etapa determinante da reação em duas ordens de magnitude, o qual pode ser atribuído a uma mudança das propriedades eletrônicas na banda d metálica em função do tamanho nanométrico das partículas, e estas modificações contribuíram para uma melhor eficiência energética quando comparado ao sistema sem a presença de eletrocatalisador. Entretanto, as Ptnp se mostraram não específicas para a RRO, catalisando as reações de degradação do solvente eletrolítico e diminuindo rapidamente a eficiência energética do dispositivo prático, devido ao acúmulo de material no eletrodo. O emprego de líquido iônico como solvente eletrolítico, ao invés de DME, promoveu uma maior estabilização do intermediário superóxido formado na primeira etapa de transferência eletrônica, devido à interação com os cátions do líquido iônico em solução, o qual resultou em um valor de constante cinética da formação do superóxido de três ordens de magnitude maior que o obtido com o mesmo eletrodo de carbono vítreo em DME, além de diminuir as reações de degradação do solvente. Estes fatores podem contribuir para uma maior potência e ciclabilidade da bateria de lítio-ar operando com líquidos iônicos.

Estimation Models for Nonuniformly Sampled Signals–Application to Multiple Delay Estimation

This letter presents a method to model propagation channels for estimation, in which the sampling scheme can be arbitrary. Additionally, the method yields accurate models, with a size that converges to the channel duration, measured in Nyquist periods. It can be viewed as an improvement on the usual discretization based on regular sampling at the Nyquist rate. The method is introduced in the context of multiple delay estimation using the MUSIC estimator, and is assessed through a numerical example.