Determination of a unique solution to parallel proton transfer reactions using the genetic algorithm

Reaction within the coulomb-cage; science in retrospect

The Coulomb-cage is defined as the space where the electrostatic interaction between two bodies is more intensive than the thermal energy (k_BT). For small molecule, the Coulomb-cage is a small sphere, extending only few water molecules towards the bulk and its radius is sensitive to the ionic strength of the solution. For charged proteins or membranal structures, the Coulomb-cage can engulf large fraction of the surface and provides a preferred pathway for ion propagation along the surface. Similarly, electrostatic potential at the inner space of a channel can form preferential trajectories passage for ions. The dynamics of ions inside the Coulomb-cage of ions was formulated by the studies of proton-anion recombination of excited photoacids. In the present article, we recount the study of intra- Coulomb-cage reaction taking place on the surface of macro-molecular bodies like micelles, membranes, proteins and intra-protein cavities. The study progressed stepwise, tracing the dynamics of a proton ejected from a photo-acid molecule located at defined sites (on membrane, inter-membrane space, active site of enzyme, inside Large Pore Channels etc.). Accumulation of experimental observations encouraged us to study of the reaction mechanism by molecular dynamics simulations of ions within the Coulomb-cage of proteins surface or inside large pores. The intra-Coulomb-cage proton transfer events follows closely the fine structure of the electrostatic field inside the cage and reflects the shape of nearby dielectric boundaries, the temporal ordering of the solvent molecules and the structural fluctuations of the charged side chains. The article sums some 40 years of research, which in retrospect clarifies the intra-Coulomb-cage reaction mechanism.

A simple yeast-based strategy to identify host cellular processes targeted by bacterial effector proteins

Bacterial effector proteins, which are delivered into the host cell via the type III secretion system, play a key role in the pathogenicity of Gram-negative bacteria by modulating various host cellular processes to the benefit of the pathogen. To identify cellular processes targeted by bacterial effectors, we developed a simple strategy that uses an array of yeast deletion strains fitted into a single 96-well plate. The array is unique in that it was optimized computationally such that despite the small number of deletion strains, it covers the majority of genes in the yeast synthetic lethal interaction network. The deletion strains in the array are screened for hypersensitivity to the expression of a bacterial effector of interest. The hypersensitive deletion strains are then analyzed for their synthetic lethal interactions to identify potential targets of the bacterial effector. We describe the identification, using this approach, of a cellular process targeted by the Xanthomonas campestris type III effector XopE2. Interestingly, we discover that XopE2 affects the yeast cell wall and the endoplasmic reticulum stress response. More generally, the use of a single 96-well plate makes the screening process accessible to any laboratory and facilitates the analysis of a large number of bacterial effectors in a short period of time. It therefore provides a promising platform for studying the functions and cellular targets of bacterial effectors and other virulence proteins.

Evaluation of the heterogeneous reactivity of the syntaxin molecules on the inner leaflet of the plasma membrane

The soluble N-ethylmaleimide-sensitive fusion (NSF) attachment protein (SNAP) receptor (SNARE) protein syntaxin 1A forms nano-sized clusters (membrane rafts) on the plasma membrane (PM) that are in equilibrium with freely diffusing syntaxin molecules. SNARE-complex formation between syntaxin 1A and SNAP-25 (synaptosome-associated protein of 25 kDa) on the PM and synaptobrevin 2 on the vesicles (trans-SNAREs) is crucial for vesicle priming and fusion. This process might be impeded by the spontaneous accumulation of non-fusogenic cis-SNARE complexes formed when all three SNARE proteins reside on the PM. We investigated the kinetics of cis-SNARE complex assembly and disassembly and both exhibited biphasic behavior. The experimental measurements were analyzed through integration of differential rate equations pertinent to the reaction mechanism and through the application of a heuristic search for time constants and concentrations using a genetic algorithm. Reconstruction of the measurements necessitated the partitioning of syntaxin into two phases that might represent the syntaxin clusters and free syntaxin outside the clusters. The analysis suggests that most of the syntaxin in the clusters is concentrated in a nonreactive form. Consequently, cis-SNARE complex assembly in the clusters is substantially slower than outside the rafts. Interestingly, the clusters also mediate efficient disassembly of cis-SNARE complexes possibly attributable to the high local concentration of complexes in the clusters area that allows efficient disassembly by the enzymatic reaction of NSF.

Endoplasmic reticulum glucosidase II is inhibited by its end products

The calnexin/calreticulin cycle is a quality control system responsible for promoting the folding of newly synthesized glycoproteins entering the endoplasmic reticulum (ER). The association of calnexin and calreticulin with the glycoproteins is regulated by ER glucosidase II, which hydrolyzes Glc 2Man X GlcNAc 2 glycans to Glc 1Man X GlcNAc 2 and further to Glc 0Man X GlcNAc 2 ( X represents any number between 5 and 9). To gain new insights into the reaction mechanism of glucosidase II, we developed a kinetic model that describes the interactions between glucosidase II, calnexin/calreticulin, and the glycans. Our model accurately reconstructed the hydrolysis of glycans with nine mannose residues and glycans with seven mannose residues, as measured by Totani et al. [Totani, K., Ihara, Y., Matsuo, I., and Ito, Y. (2006) J. Biol. Chem. 281, 31502-31508]. Intriguingly, our model predicted that glucosidase II was inhibited by its nonglucosylated end products, where the inhibitory effect of Glc 0Man 7GlcNAc 2 was much stronger than that of Glc 0Man 9GlcNAc 2. These predictions were confirmed experimentally. Moreover, our model suggested that glycans with a different number of mannose residues can be equivalent substrates of glucosidase II, in contrast to what had been previously thought. We discuss the possibility that nonglucosylated glycans, existing in the ER, might regulate the entry of newly synthesized glycoproteins into the calnexin/calreticulin cycle. Our model also shows that glucosidase II does not interact with monoglucosylated glycans while they are bound to calnexin or calreticulin.

Systematic search for the rate constants that control the exocytotic process from chromaffin cells by a Genetic Algorithm

We have recently created a kinetic model that reproduces the dynamics of exocytosis with high accuracy. The reconstruction necessitated a search, in a multi-dimensional parameter space, for 37 parameters that described the system, with no assurance that the parameters, which reconstructed the observations, are a unique set. In the present study, a Genetic Algorithm (GA) was used for a thorough search in the unknown parameter space, using a strategy of gradual increase of the complexity of the analyzed input data. Upon systematic incorporation of one to four measurable parameters, used as input signals for the analysis, the constraint set on the GA search imposed the convergence of the free parameters into a single narrow range. The mean values for each adjustable parameter represent a minimum for the fitness function in the multi-dimensional parameter space. The GA search demonstrates that the parameters that control the kinetics of exocytosis are the rate constants of the steps downstream to synaptotagmin binding, and that the equilibrium constant of the binding of calcium to Munc13 controls the calcium-dependent priming process. Thus, the systematic use of the GA creates a link between specific reactions in the process of exocytosis and experimental phenotypes.

Photophysics of a Series of Efficient Fluorescent pH Probes for Dual-Emission-Wavelength Measurements in Aqueous Solutions

This paper evaluates the 5-aryl-2-pyridyloxazole backbone to engineer donor-acceptor fluorescent pH probes after one- or two-photon absorption. Parent fluorophores, as well as derivatives that can be used to label biomolecules, can be easily obtained in good yields. These molecules exhibit a large one-photon absorption in the near-UV range, and a strong fluorescence emission that covers the whole visible domain. The 5-aryl-2-pyridyloxazole derivatives also possess significant cross sections for two-photon absorption. Upon pyridine protonation, large shifts were observed in the absorption spectra after one- and two-photon excitation, as well as in the emission spectra. This feature was used to measure the pK(a) of the investigated compounds that range between 2 and 8. In most of the investigated derivatives, the pK(a) increased upon light excitation and protonation exchanges took place during the lifetime of the excited state, as shown by phase-modulation fluorometry analysis. Several 5-aryl-2-pyridyloxazole derivatives are suggested as efficient probes to reliably measure the pH of aqueous solutions by means of ratiometric methods that are dependent on fluorescence emission.

Reactant Concentrations from Fluorescence Correlation Spectroscopy with Tailored Fluorescent Probes. An Example of Local Calibration-Free pH Measurement

The present account is concerned with the measurement of local reactant concentrations by observing specific fluorescent probes in fluorescence correlation spectroscopy (FCS). The Theoretical Analysis section revisits the photophysical, thermodynamic, and kinetic information that is contained in the corresponding FCS correlation curves. In particular, we examine the conditions under which FCS is revealed as a superior tool to measure concentrations of reactive species. Careful molecular engineering of the specific fluorescent probes that simultaneously integrates photophysical, thermodynamic, and kinetic constraints will be required to benefit most from FCS. We illustrate the FCS titration approach with a series of fluorescent probes that we tailored to measure pH at around 4-6 by FCS after two-photon excitation. We show that an optimal design allows one to access pH without any preliminary calibrations such as the determination of the protonation constant or the photophysical properties of the fluorescent probe.

A new platform to study the molecular mechanisms of exocytosis

The exocytotic process in neurons and neuroendocrine cells consists of a sequence of reactions between well defined proteins. In the present study, we have created for the first time a comprehensive kinetic model that demonstrates the dynamics of interactions between key synaptic proteins that are associated with exocytosis. The interactions between the synaptic proteins were transformed into differential rate equations that, after their integration over time, reconstructed the experimental signal. The model can perfectly reconstruct the kinetics of exocytosis, the calcium-dependent priming and fusion processes, and the effects of genetic manipulation of synaptic proteins. The model suggests that fusion occurs from two parallel pathways and assigns precise, non-identical synaptic protein complexes to the two pathways. In addition, it provides a unique opportunity to study the dynamics of intermediate protein complexes during the fusion process, a possibility that is hidden in most experimental systems. We thus developed a novel approach that allows detailed characterization of the temporal relationship between synaptic protein complexes. This model provides an excellent platform for prediction and quantification of the effects of protein manipulations on exocytosis and opens new avenues for experimental investigation of exocytosis.