PubMed Health⌕ Search

Biomedical subjects

Menachem Gutman

Publications and source records attributed to Menachem Gutman.

At least 19 recordsLinked to original sources
PubMed2006

Proton antenna effect of the gamma-cyclodextrin outer surface, measured by excited state proton transfer.

The reversible proton dissociation and geminate recombination of the common photoacid, 8-hydroxypyrene-1,3,6-trisulfonate (pyranine), either in dilute aqueous solution or when forming a complex with gamma-cyclodextrin (gamma-CD), has been studied by time-resolved fluorescence spectroscopy and supplemented by molecular modeling and dynamics simulations. We find that the dissociation rate of the proton from the excited molecule was decreased to about approximately 50% of its value in water, while the rate of recombination was doubled. These observations were evaluated by molecular modeling of the reactants at atomic resolution. The combination of the two methodologies indicates that the pyranine in the complex can assume more than one level of interaction with the solvent. The polysugar torus surrounding the pyranine perturbs the hydrogen bond in the dye's immediate vicinity and deforms the electrostatic potential inside the Coulomb cage, causing major deviations from a simple spheric symmetry. These observations can account for the special kinetic features measured for the complex. We suggest that this system can be used as a basic model for evaluating the mechanism of proton transfer in non-homogeneous systems, such as the surface of proteins or biomembranes.

Arylsulfonates↗
PubMed2006

A molecular dynamics study and free energy analysis of complexes between the Mlc1p protein and two IQ motif peptides.

The Mlc1p protein from the budding yeast Saccharomyces cerevisiae is a Calmodulin-like protein, which interacts with IQ-motif peptides located at the yeast's myosin neck. In this study, we report a molecular dynamics study of the Mlc1p-IQ2 protein-peptide complex, starting with its crystal structure, and investigate its dynamics in an aqueous solution. The results are compared with those obtained by a previous study, where we followed the solution structure of the Mlc1p-IQ4 protein-peptide complex by molecular dynamics simulations. After the simulations, we performed an interaction free-energy analysis using the molecular mechanics Poisson-Boltzmann surface area approach. Based on the dynamics of the Mlc1p-IQ protein-peptide complexes, the structure of the light-chain-binding domain of myosin V from the yeast S. cerevisiae is discussed.

Amino Acid Motifs↗
PubMed2006

Molecular dynamics study of a calmodulin-like protein with an IQ peptide: spontaneous refolding of the protein around the peptide.

The Calmodulin (CaM) is a small (16.7 kDa), highly acidic protein that is crucial to all eukaryotes by serving as a prototypical calcium sensor. In the present study, we investigated, through molecular dynamics simulations, the dynamics of a complex between the Mlc1p protein, which is a CaM-like protein, and the IQ4 peptide. This protein-peptide interaction is of high importance because IQ motifs are widely distributed among different kinds of CaM-binding proteins. The Mlc1p-IQ4 complex, which had been resolved by crystallography to 2.1 A, confers to a Ca(+2)-independent stable structure. During the simulations, the complex undergoes a complicated modulation process, which involves bending of the angles between the alpha-helices of the protein, breaking of the alpha-helical structure of the IQ4 peptide into two sections, and formation of new contact points between the protein and the peptide. The dynamics of the process consist of fast sub picosecond events and much slower ones that take a few nanoseconds to completion. Our study expands the information embedded in the crystal structure of the Mlc1p-IQ4 complex by describing its dynamic behavior as it evolves from the crystal structure to a form stable in solution. The article shows that careful application of molecular dynamics simulations can be used for extending the structural information presented by the crystal structure, thereby revealing the dynamic configuration of the protein in its physiological environment.

Algorithms↗
PubMed2006

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.

Algorithms↗
PubMed2006

A molecular dynamics study of the effect of Ca2+ removal on calmodulin structure.

Calmodulin is a small (148 residues), ubiquitous, highly-conserved Ca(2+) binding protein serving as a modulator of many calcium-dependent processes. In this study, we followed, by means of molecular dynamics, the structural stability of the protein when one of its four bound Ca(2+) ions is removed, and compared it to a simulation of the fully Ca(2+) bound protein. We found that the removal of a single Ca(2+) ion from the N-lobe of the protein, which has a lower affinity for the ion, is sufficient to initiate a considerable structural rearrangement. Although the overall structure of the fully 4 Ca(2+) bound protein remained intact in the extended conformation, the Ca(2+)-removed protein changed its conformation into a compact state. The observation that the 3 Ca(2+) loaded protein assumes a compacted solution state is in accord with experimental observation that the NSCP protein, which binds only three Ca(2+) ions, is natively in a compact state. Examination of the folding dynamics reveals a cooperation between the C-lobe, N-lobe, and the interdomain helix that enable the conformation change. The forces driving this conformational change are discussed.

Animals↗
PubMed2006

The mechanism of proton transfer between adjacent sites on the molecular surface.

The surface of a protein, or a membrane, is spotted with a multitude of proton binding sites, some of which are only few A apart. When a proton is released from one site, it propagates through the water by a random walk under the bias of the local electrostatic potential determined by the distribution of the charges on the protein. Eventually, the released protons are dispersed in the bulk, but during the first few nanoseconds after the dissociation, the protons can be trapped by encounter with nearby acceptor sites. While the study of this reaction on the surface of a protein suffers from experimental and theoretical difficulties, it can be investigated with simple model compounds like derivatives of fluorescein. In the present study, we evaluate the mechanism of proton transfer reactions that proceed, preferentially, inside the Coulomb cage of the dye molecules. Kinetic analysis of the measured dynamics reveals the role of the dimension of the Coulomb cage on the efficiency of the reaction and how the ordering of the water molecules by the dye affects the kinetic isotope effect.

Binding Sites↗
PubMed2005

Fatty acid binding proteins: same structure but different binding mechanisms? Molecular dynamics simulations of intestinal fatty acid binding protein.

Fatty acid binding proteins (FABPs) carry fatty acids (FAs) and other lipids in the cellular environment, and are thus involved in processes such as FA uptake, transport, and oxidation. These proteins bind either one or two ligands in a binding site, which appears to be inaccessible from the bulk. Thus, the entry of the substrate necessitates a conformational change, whose nature is still unknown. A possible description of the ligand binding process is given by the portal hypothesis, which suggests that the FA enters the protein through a dynamic area known as the portal region. On the other hand, recent simulations of the adipocyte lipid binding protein (ALBP) suggested a different entry site (the alternative portal). In this article, we discuss molecular dynamics simulations of the apo-intestinal-FABP (I-FABP) in the presence of palmitate molecule(s) in the simulation box. The simulations were carried out to study whether the FA can enter the protein during the simulations (as in the ALBP) and where the ligand entry site is (the portal region, the alternative portal or a different domain). The analysis of the simulations revealed a clear difference between the ALBP and the I-FABP. In the latter case, the palmitate preferentially adsorbed to the portal region, which was more mobile than the rest of the protein. However, no ligand entry was observed in the multi-nanosecond-long simulations, in contrast to ALBP. These findings suggest that, although the main structural motif of the FABPs is common, the fine details of each individual protein structure grossly modulate its reactivity.

Adsorption↗
PubMed2005

Application of classical molecular dynamics for evaluation of proton transfer mechanism on a protein.

Proton transfer reactions on surfaces are prevalent in biology, chemistry and physics. In the present study, we employed classical Molecular Dynamics simulations to search for the presence of transient configurations that enable proton transfer, or proton sharing, between adjacent carboxylate groups on the protein surface. The results demonstrate that, during random fluctuations of the residues on the surface, there are repeated situations in which nearby carboxylates either share a common proton through a hydrogen bond, or are connected by a few water molecules that form conducting networks. These networks do not extend out of the common Coulomb cage of the participating residues and the lifetimes of the bridged structures are sufficiently long to allow passage of a proton between the carboxylates. The detection of domains capable of supporting a rapid proton transfer on a protein supports the notion that clusters of carboxylates are the operative elements of proton collecting antennae, as in bacteriorhodopsin, cytochrome c oxidase or the photosynthetic reaction center.

Aspartic Acid↗
PubMed2005

The mechanism of proton transfer between adjacent sites exposed to water.

The surface of a protein, or a membrane, is spotted with a multitude of proton-binding sites, some of which are only a few angstroms apart. When a proton is released from one site, it propagates through the water by a random walk under the bias of the local electrostatic potential determined by the distribution of the charges on the protein. Some of the released protons disperse into the bulk, but during the first few nanoseconds, the released protons can be trapped by encounter with nearby acceptor sites. This process resembles a scenario which corresponds with the time-dependent Debye-Smoluchowski equation. In the present study, we investigated the mechanism of proton transfer between sites that are only a few angstroms apart, using as a model the proton exchange between sites on a small molecule, fluorescein, having two, spectrally distinguishable, proton-binding sites. The first site is the oxyanion on the chromophore ring structure. The second site is the carboxylate moiety on the benzene ring of the molecule. Through our experiments, we were able to reconstruct the state of protonation at each site and the velocity of proton transfer between them. The fluorescein was protonated by a few nanosecond long proton pulse under specific conditions that ensured that the dye molecules would be protonated only by a single proton. The dynamics of the protonation of the chromophore were measured under varying initial conditions (temperature, ionic strength, and different solvents (H(2)O or D(2)O)), and the velocity of the proton transfer between the two sites was extracted from the overall global analysis of the signals. The dynamics of the proton transfer between the two proton-binding sites of the fluorescein indicated that the efficiency of the site-to-site proton transfer is very sensitive to the presence of the screening electrolyte and has a very high kinetic isotope effect (KIE = 55). These two parameters clearly distinguish the mechanism from proton diffusion in bulk water. The activation energy of the reaction (E(a) = 11 kcal mol(-1)) is also significantly higher than the activation energy for proton dissociation in bulk water (E(a) approximately 2.5 kcal mol(-1)). These observations are discussed with respect to the effect of the solute on the water molecules located within the solvation layer.

Anions↗
PubMed2005

Molecular dynamics of a protein surface: ion-residues interactions.

Time-resolved measurements indicated that protons could propagate on the surface of a protein or a membrane by a special mechanism that enhanced the shuttle of the proton toward a specific site. It was proposed that a suitable location of residues on the surface contributes to the proton shuttling function. In this study, this notion was further investigated by the use of molecular dynamics simulations, where Na(+) and Cl(-) are the ions under study, thus avoiding the necessity for quantum mechanical calculations. Molecular dynamics simulations were carried out using as a model a few Na(+) and Cl(-) ions enclosed in a fully hydrated simulation box with a small globular protein (the S6 of the bacterial ribosome). Three independent 10-ns-long simulations indicated that the ions and the protein's surface were in equilibrium, with rapid passage of the ions between the protein's surface and the bulk. However, it was noted that close to some domains the ions extended their duration near the surface, thus suggesting that the local electrostatic potential hindered their diffusion to the bulk. During the time frame in which the ions were detained next to the surface, they could rapidly shuttle between various attractor sites located under the electrostatic umbrella. Statistical analysis of the molecular dynamics and electrostatic potential/entropy consideration indicated that the detainment state is an energetic compromise between attractive forces and entropy of dilution. The similarity between the motion of free ions next to a protein and the proton transfer on the protein's surface are discussed.

Amino Acids↗
PubMed2005

Molecular dynamics simulations of the adipocyte lipid binding protein reveal a novel entry site for the ligand.

The adipocyte lipid binding protein (ALBP) binds fatty acids (FA) in a cavity that is inaccessible from the bulk. Therefore, the penetration of the FA necessitates conformational changes whose nature is still unknown. It was suggested that the lipid first enters through a "portal region" which consists of the alphaII helix and the adjacent tight turns. The initial event in the ligand binding must be the interaction of the lipid with the protein surface. To analyze this interaction, we have carried out three molecular dynamics simulations of the apo-ALBP, with a palmitate ion initially located at different positions near the protein surface. The simulation indicated that the ligand could adsorb to the protein in more than one location. Yet, in one case, the ligand managed to penetrate the protein by entering a newly formed cavity some 10 A deep. The entry site is located near the N-terminus, at the junction between the loops connecting the beta-strands. Further progression of the penetration seems to be arrested by the need for desolvation of the COOH end of the headgroup. Evolutionary analysis showed that amino acids in this entry site are well conserved. On the basis of these observations, we suggest that the ligand may enter the protein from a site other than the portal region. Furthermore, the rate-limiting step is proposed to be the desolvation of the FA polar headgroup.

Adipocytes↗
PubMed2004

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.

Animals↗
PubMed2003

The role of small intraprotein cavities in the catalytic cycle of bacteriorhodopsin.

The last phase of the proton transfer cycle of bacteriorhodopsin calls for a passage of a proton from D38 to D96. This reaction utilizes a narrow shaft approximately 10-A long that connects the two carboxylates that cross through a very hydrophobic domain. As the shaft is too narrow to be permanently hydrated, there are two alternatives for the proton propagation into the channel. The proton may propagate through the shaft without solvation at the expense of a high electrostatic barrier; alternatively, the shaft will expand to accommodate some water molecules, thus lowering the Born energy for the insertion of the charge into the protein (B. Schätzler, N. A. Dencher, J. Tittor, D. Oesterhelt, S. Yaniv-Checover, E. Nachliel, and G. Gutman, 2003, BIOPHYS: J. 84:671-686). A comparative study of nine published crystal-structures of bacteriorhodopsin identified, next to the shaft, microcavities in the protein whose position and surrounding atoms are common to the reported structures. Some of the cavities either shrink or expand during the photocycle. It is argued that the plasticity of the cavities provides a working space needed for the transient solvation of the shaft, thus reducing the activation energy necessary for the solvation of the shaft. This suggestion is corroborated by the recent observations of Klink et al. (B. U. Klink, R. Winter, M. Engelhard, and I. Chizhov, 2002, BIOPHYS: J. 83:3490-3498) that the late phases of the photocycle (tau >/=1 ms) are strongly inhibited by external pressure.

Bacteriorhodopsins↗
PubMed2003

Subsecond proton-hole propagation in bacteriorhodopsin.

The dynamics of proton transfer between the surface of purple membrane and the aqueous bulk have recently been investigated by the Laser Induced Proton Pulse Method. Following a Delta-function release of protons to the bulk, the system was seen to regain its state of equilibrium within a few hundreds of microseconds. These measurements set the time frame for the relaxation of any state of acid-base disequilibrium between the bacteriorhodopsin's surface and the bulk. It was also deduced that the released protons react with the various proton binding within less than 10 micro s. In the present study, we monitored the photocycle and the proton-cycle of photo-excited bacteriorhodopsin, in the absence of added buffer, and calculated the proton balance between the Schiff base and the bulk phase in a time-resolved mode. It was noticed that the late phase of the M decay (beyond 1 ms) is characterized by a slow (subsecond) relaxation of disequilibrium, where the Schiff base is already reprotonated but the pyranine still retains protons. Thus, it appears that the protonation of D96 is a slow rate-limiting process that generates a "proton hole" in the cytoplasmic section of the protein. The velocity of the hole propagation is modulated by the ionic strength of the solution and by selective replacements of charged residues on the interhelical loops of the protein, at domains that seems to be remote from the intraprotein proton conduction trajectory.

Arylsulfonates↗
PubMed2002

Excited state proton transfer in reverse micelles.

The aqueous phase of water/AOT reversed micelles having varying diameters was probed by a single free diffusing proton that was released form a hydrophilic photoacid molecule (2-naphthol-6,8-disulfonate). The fluorescence decay signals were reconstructed through the geminate recombination algorithm, accounting for the reversible nature of the proton-transfer reactions at the surface of the excited molecule and at the water/detergent interface. The radial diffusion of the proton inside the aqueous phase was calculated accounting for both the entropy of dilution and the total electrostatic energy of the ion pair, consisting of the pair-energy and self-energy of the ions. The analysis implied that micellar surface must be modeled with atomic resolution, assuming that the sulfono residue protrudes above the water/hydrocarbon interface by approximately 2 A. The analysis of the fluorescence decay curves implies that the molecule is located in a solvent with physical-chemical properties very similar to bulk water, except for the dielectric constant. For reversed micelles with r(max) > or = 16 A, the dielectric constant of the aqueous phase was approximately 70 and for smaller micelles, where approximately 60% of the water molecule is in contact with the van der Waals surface of the micelle, it is as low as 60. This reduction is a reflection of the increased fraction of water molecule that is in close interaction with the micelle surface.

Journal Article↗
PubMed2002

Stability of an ion channel in lipid bilayers: implicit solvent model calculations with gramicidin.

Gramicidin is a helical peptide, 15 residues in length, which dimerizes to form ion-conducting channels in lipid bilayers. Here we report calculations of its free energy of transfer from the aqueous phase into bilayers of different widths. The electrostatic and nonpolar contributions to the desolvation free energy were calculated using implicit solvent models, in which gramicidin was described in atomic detail and the hydrocarbon region of the membrane was described as a slab of hydrophobic medium embedded in water. The free energy penalties from the lipid perturbation and membrane deformation effects, and the entropy loss associated with gramicidin immobilization in the bilayer, were estimated from a statistical thermodynamic model of the bilayer. The calculations were carried out using two classes of experimentally observed conformations: a head-to-head dimer of two single-stranded (SS) beta-helices and a double-stranded (DS) intertwined double helix. The calculations showed that gramicidin is likely to partition into the bilayer in all of these conformations. However, the SS conformation was found to be significantly more stable than the DS in the bilayer, in agreement with most of the experimental data. We tested numerous transmembrane and surface orientations of gramicidin in bilayers of various widths. Our calculations indicate that the most favorable orientation is transmembrane, which is indeed to be expected from a channel-forming peptide. The calculations demonstrate that gramicidin insertion into the membrane is likely to involve a significant deformation of the bilayer to match the hydrophobic width of the peptide (22 A), again in good agreement with experimental data. Interestingly, deformation of the bilayer was induced by all of the gramicidin conformations.

Gramicidin↗
PubMed2002

Proton transfer dynamics on the surface of the late M state of bacteriorhodopsin.

The cytoplasmic surface of the BR (initial) state of bacteriorhodopsin is characterized by a cluster of three carboxylates that function as a proton-collecting antenna. Systematic replacement of most of the surface carboxylates indicated that the cluster is made of D104, E161, and E234 (Checover, S., Y. Marantz, E. Nachliel, M. Gutman, M. Pfeiffer, J. Tittor, D. Oesterhelt, and N. Dencher. 2001. Biochemistry. 40:4281-4292), yet the BR state is a resting configuration; thus, its proton-collecting antenna can only indicate the presence of its role in the photo-intermediates where the protein is re-protonated by protons coming from the cytoplasmic matrix. In the present study we used the D96N and the triple (D96G/F171C/F219L) mutant for monitoring the proton-collecting properties of the protein in its late M state. The protein was maintained in a steady M state by continuous illumination and subjected to reversible pulse protonation caused by repeated excitation of pyranine present in the reaction mixture. The re-protonation dynamics of the pyranine anion was subjected to kinetic analysis, and the rate constants of the reaction of free protons with the surface groups and the proton exchange reactions between them were calculated. The reconstruction of the experimental signal indicated that the late M state of bacteriorhodopsin exhibits an efficient mechanism of proton delivery to the unoccupied-most basic-residue on its cytoplasmic surface (D38), which exceeds that of the BR configuration of the protein. The kinetic analysis was carried out in conjunction with the published structure of the M state (Sass, H., G. Büldt, R. Gessenich, D. Hehn, D. Neff, R. Schlesinger, J. Berendzen, and P. Ormos. 2000. Nature. 406:649-653), the model that resolves most of the cytoplasmic surface. The combination of the kinetic analysis and the structural information led to identification of two proton-conducting tracks on the protein's surface that are funneling protons to D38. One track is made of the carboxylate moieties of residues D36 and E237, while the other is made of D102 and E232. In the late M state the carboxylates of both tracks are closer to D38 than in the BR (initial) state, accounting for a more efficient proton equilibration between the bulk and the protein's proton entrance channel. The triple mutant resembles in the kinetic properties of its proton conducting surface more the BR-M state than the initial state confirming structural similarities with the BR-M state and differences to the BR initial state.

Bacteriorhodopsins↗