Minimum Energy Pathways for Proton Transfer between Adjacent Sites Exposed to Water

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.

Dynamic water bridging and proton transfer at a surface carboxylate cluster of photosystem II

A hydrogen-bond cluster at a negatively-charged protein interface with a bound protein and long-lived waters might be a proton storage site.

Prediction of perturbed proton transfer networks

The transfer of protons through proton translocating channels is a complex process, for which direct samplings of different protonation states and side chain conformations in a transition network calculation provide an efficient, bias-free description. In principle, a new transition network calculation is required for every unsampled change in the system of interest, e.g. an unsampled protonation state change, which is associated with significant computational costs. Transition networks void of or including an unsampled change are termed unperturbed or perturbed, respectively. Here, we present a prediction method, which is based on an extensive coarse-graining of the underlying transition networks to speed up the calculations. It uses the minimum spanning tree and a corresponding sensitivity analysis of an unperturbed transition network as initial guess and refinement parameter for the determination of an unknown, perturbed transition network. Thereby, the minimum spanning tree defines a sub-network connecting all nodes without cycles and minimal edge weight sum, while the sensitivity analysis analyzes the stability of the minimum spanning tree towards individual edge weight reductions. Using the prediction method, we are able to reduce the calculation costs in a model system by up to 80%, while important network properties are maintained in most predictions.

The mechanism of water/ion exchange at a protein surface: a weakly bound chloride in Helicobacter pylori apoflavodoxin

Binding/unbinding of small ligands, such as ions, to/from proteins influences biochemical processes such as protein folding, enzyme catalysis or protein/ligand recognition. We have investigated the mechanism of chloride/water exchange at a protein surface (that of the apoflavodoxin from Helicobacter pylori) using classical all-atom molecular dynamics simulations. They reveal a variety of chloride exit routes and residence times; the latter is related to specific coordination modes of the anion. The role of solvent molecules in the mechanism of chloride unbinding has been studied in detail. We see no temporary increase in chloride coordination along the release process. Instead, the coordination of new water molecules takes place in most cases after the chloride/protein atom release event has begun. Moreover, the distribution function of water entrance events into the first chloride solvation shell peaks after chloride protein atom dissociation events. All these observations together seem to indicate that water molecules simply fill the vacancies left by the previously coordinating protein residues. We thus propose a step-by-step dissociation pathway in which protein/chloride interactions gradually break down before new water molecules progressively fill the vacant positions left by protein atoms. As observed for other systems, water molecules associated with bound chloride or with protein atoms have longer residence times than those bound to the free anion. The implications of the exchange mechanism proposed for the binding of the FMN (Flavin Mononucleotide) protein cofactor are discussed.

QM/MM investigation of the catalytic mechanism of angiotensin-converting enzyme

Angiotensin-converting enzyme (ACE) converts angiotensin I to angiotensin II and degrades bradykinin and other vasoactive peptides. ACE inhibitors are used to treat diseases such as hypertension and heart failure. It is thus highly desirable to understand the catalytic mechanism of ACE, as this should facilitate the design of more powerful and selective ACE inhibitors. ACE exhibits two different active domains, the C-domain and the N-domain. In this work, we systematically investigated the inhibitor- and substrate-binding patterns in the N-domain of human ACE using a combined quantum mechanical and molecular mechanical approach. The hydrolysis of hippuryl-histidyl-leucine (HHL) as catalyzed by the N-domain of human somatic ACE was explored, and the effects of chloride ion on the overall reaction were also investigated. Two models, one with and one without a chloride ion at the first binding position, were then designed to examine the chloride dependence of inhibitor-substrate binding and the catalytic mechanism. Our calculations indicate that the hydrolysis reaction follows a stepwise general base/general acid catalysis path. The estimated mean free energy barrier height in the two models is about 15.6 kcal/mol, which agrees very well with the experimentally estimated value of 15.8 kcal/mol. Our simulations thus suggest that the N-domain is in a mixed form during ACE-catalyzed hydrolysis, with the single-chloride-ion and the double-chloride-ion forms existing simultaneously. Graphical Abstract Superposition of ACE C- and N- domains.

Effect of sucrose on chemically and thermally induced unfolding of domain-I of human serum albumin: Solvation dynamics and fluorescence anisotropy study

The present study is devoted to understand the effect of sucrose on the hydration dynamics and rotational relaxation dynamics within the domain-I of HSA during chemically as well as thermally induced unfolding. It has been observed that the average solvation time become slower in the presence of sucrose for the lower concentrations of GnHCl, however at higher concentrations of GnHCl the effect of sucrose is almost negligible. From the time resolved fluorescence anisotropy it has been observed that in the lower concentration region of GnHCl the sucrose induced stabilization is small as compared to the higher concentrations of GnHCl. We have concluded that the hydration dynamics plays an important role in the sucrose induced stabilization process at the low concentration region; whereas environmental restriction is responsible at the higher concentration of GnHCl. However, we have observed a negligible stabilizing effect of sucrose towards the temperature induced unfolding.

Energetics and dynamics of proton transfer reactions along short water wires

Proton transfer (pT) reactions in biochemical processes are often mediated by chains of hydrogen-bonded water molecules. We use hybrid density functional calculations to study pT along quasi one-dimensional water arrays that connect an imidazolium-imidazole proton donor-acceptor pair. We characterize the structures of intermediates and transition states, the energetics, and the dynamics of the pT reactions, including vibrational contributions to kinetic isotope effects. In molecular dynamics simulations of pT transition paths, we find that for short water chains with four water molecules, the pT reactions are semi-concerted. The formation of a high-energy hydronium intermediate next to the proton-donating group is avoided by a simultaneous transfer of a proton from the donor to the first water molecule, and from the first water molecule into the water chain. Lowering the dielectric constant of the environment and increasing the water chain length both reduce the barrier for pT. We study the effect of the driving force on the energetics of the pT reaction by changing the proton affinity of the donor and acceptor groups through halogen and methyl substitutions. We find that the barrier of the pT reaction depends linearly on the proton affinity of the donor but is nearly independent of the proton affinity of the acceptor, corresponding to Brønsted slopes of one and zero, respectively.

Water-assisted proton transfer in ferredoxin I

The role of water molecules in assisting proton transfer (PT) is investigated for the proton-pumping protein ferredoxin I (FdI) from Azotobacter vinelandii. It was shown previously that individual water molecules can stabilize between Asp(15) and the buried [3Fe-4S](0) cluster and thus can potentially act as a proton relay in transferring H(+) from the protein to the μ(2) sulfur atom. Here, we generalize molecular mechanics with proton transfer to studying proton transfer reactions in the condensed phase. Both umbrella sampling simulations and electronic structure calculations suggest that the PT Asp(15)-COOH + H(2)O + [3Fe-4S](0) → Asp(15)-COO(-) + H(2)O + [3Fe-4S](0) H(+) is concerted, and no stable intermediate hydronium ion (H(3)O(+)) is expected. The free energy difference of 11.7 kcal/mol for the forward reaction is in good agreement with the experimental value (13.3 kcal/mol). For the reverse reaction (Asp(15)-COO(-) + H(2)O + [3Fe-4S](0)H(+) → Asp(15)-COOH + H(2)O + [3Fe-4S](0)), a larger barrier than for the forward reaction is correctly predicted, but it is quantitatively overestimated (23.1 kcal/mol from simulations versus 14.1 from experiment). Possible reasons for this discrepancy are discussed. Compared with the water-assisted process (ΔE ≈ 10 kcal/mol), water-unassisted proton transfer yields a considerably higher barrier of ΔE ≈ 35 kcal/mol.

The self-organizing fractal theory as a universal discovery method: the phenomenon of life

A universal discovery method potentially applicable to all disciplines studying organizational phenomena has been developed. This method takes advantage of a new form of global symmetry, namely, scale-invariance of self-organizational dynamics of energy/matter at all levels of organizational hierarchy, from elementary particles through cells and organisms to the Universe as a whole. The method is based on an alternative conceptualization of physical reality postulating that the energy/matter comprising the Universe is far from equilibrium, that it exists as a flow, and that it develops via self-organization in accordance with the empirical laws of nonequilibrium thermodynamics. It is postulated that the energy/matter flowing through and comprising the Universe evolves as a multiscale, self-similar structure-process, i.e., as a self-organizing fractal. This means that certain organizational structures and processes are scale-invariant and are reproduced at all levels of the organizational hierarchy. Being a form of symmetry, scale-invariance naturally lends itself to a new discovery method that allows for the deduction of missing information by comparing scale-invariant organizational patterns across different levels of the organizational hierarchy.An application of the new discovery method to life sciences reveals that moving electrons represent a keystone physical force (flux) that powers, animates, informs, and binds all living structures-processes into a planetary-wide, multiscale system of electron flow/circulation, and that all living organisms and their larger-scale organizations emerge to function as electron transport networks that are supported by and, at the same time, support the flow of electrons down the Earth's redox gradient maintained along the core-mantle-crust-ocean-atmosphere axis of the planet. The presented findings lead to a radically new perspective on the nature and origin of life, suggesting that living matter is an organizational state/phase of nonliving matter and a natural consequence of the evolution and self-organization of nonliving matter.The presented paradigm opens doors for explosive advances in many disciplines, by uniting them within a single conceptual framework and providing a discovery method that allows for the systematic generation of knowledge through comparison and complementation of empirical data across different sciences and disciplines.

Probing deuterium isotope effect on structure and solvation dynamics of human serum albumin

The deuterium isotopic effect on the structure and solvation dynamics of the protein, human serum albumin (HSA), has been studied by using circular dichroism (CD), femtosecond up-conversion, FRET, and single-molecule spectroscopy. The CD spectra suggest that D(2)O affects the structure of HSA, leading to a 20% decrease in the helical structure. The FRET study indicates that the distance of C153 from the lone tryptophan residue of HSA is quite similar (≈21 Å) in H(2)O and D(2)O, and hence, the location of the probe in the protein remains the same in the two solvents. The single-molecule study suggests that coumarin 153 (C153) binds almost exclusively (>96%) to one site of HSA. Solvation dynamics of C153 in HSA is found to be markedly retarded in D(2)O compared with H(2)O. In H(2)O, the solvation of C153 bound to HSA is found to be biexponential with one component of 7 ps (30%) and a long component of 350 ps (70%). In D(2)O, we detected a short component of 4 ps (41%) and a long component of 950 ps (59%). Thus, the ultraslow component of the solvation dynamics of C153 bound to HSA in D(2)O (950 ps) is 2.5-fold slower than that in H(2)O (350 ps). The marked deuterium isotope effect has been ascribed to water molecules confined in the protein environment and to a lesser extent to the structural modification of protein by D(2)O.

On the orientation of the catalytic dyad in aspartic proteases

The recent re-refinement of the X-ray structure of apo plasmepsin II from Plasmodium falciparum suggests that the two carboxylate groups in the catalytic dyad are noncoplanar, (Robbins et al., Acta Crystallogr D Biol Crystallogr 2009;65: 294-296) in remarkable contrast with the vast majority of structures of aspartic proteases. Here, evidence for the noncoplanarity of the catalytic aspartates is provided by analysis of multiple explicit water molecular dynamics (MD) simulations of plasmepsin II, human beta-secretase, and HIV-protease. In the MD runs of plasmepsin II, the angle between the planes of the two carboxylates of the catalytic dyad is almost always in the range 60 degrees -120 degrees , in agreement with the perpendicular orientation in the re-refined X-ray structure. The noncoplanar arrangement is prevalent also in the beta-secretase simulations, as well as in the runs with the inhibitor-bound proteases. Quantum-mechanics calculations provide further evidence that before catalysis the noncoplanar arrangement is favored energetically in eukaryotic aspartic proteases. Remarkably, the coplanar orientation of the catalytic dyad is observed in MD simulations of HIV-protease at 100 K but not at 300 K, which indicates that the noncoplanar arrangement is favored by conformational entropy. This finding suggests that the coplanar orientation in the crystal structures of apo aspartic proteases is promoted by the very low temperature used for data collection (usually around 100 K).

Lateral proton transfer between the membrane and a membrane protein

Proton transport across biological membranes is a key step of the energy conservation machinery in living organisms, and it has been proposed that the membrane itself plays an important role in this process. In the present study we have investigated the effect of incorporation of a proton transporter, cytochrome c oxidase, into a membrane on the protonation kinetics of a fluorescent pH-sensitive probe attached at the surface of the protein. The results show that proton transfer to the probe was slightly accelerated upon attachment at the protein surface (approximately 7 x 1010 s(-1) M(-1), compared to the expected value of (1-2) x 10(10) s(-1) M(-1)), which is presumably due to the presence of acidic/His groups in the vicinity. Upon incorporation of the protein into small unilamellar phospholipid vesicles the rate increased by more than a factor of 400 to approximately 3 x 10(13) s(-1) M(-1), which indicates that the protein-attached probe is in rapid protonic contact with the membrane surface. The results indicate that the membrane acts to accelerate proton uptake by the membrane-bound proton transporter.

Phosphate assisted proton transfer in water and sugar glasses: a study using fluorescence of pyrene-1-carboxylate and IR spectroscopy

The role of water's H-bond percolation network in acid-assisted proton transfer was studied in water and glycerol solutions and in sugar glasses. Proton transfer rates were determined by the fluorescence of pyrene-1-carboxylate, a compound with a higher pK in its excited state relative to the ground state. Excitation of pyrene-1-COO- produces fluorescence from pyrene-1-COOH when a proton is accepted during the excited singlet state lifetime of pyrene-1-COO-. The presence of glycerol as an aqueous cosolvent decreases proton transfer rates from phosphoric and acetic acid in a manner that does not follow the Stokes relationship on viscosity. In sugar glass composed of trehalose and sucrose, proton transfer occurs when phosphate is incorporated in the glass. Sugar glass containing phosphate retains water and it is suggested that proton transfer requires this water. The infrared (IR) frequency of water bending mode in sugar glass and in aqueous solution is affected by the presence of phosphate and the IR spectral bands of all phosphate species in water are temperature dependent; both results are consistent with H-bonding between water and phosphate. The fluorescence results, which studied the effect of cosolvent, highlight the role of water in assisting proton transfer in reactions involving biological acids, and the IR results, which give spectroscopic evidence for H-bonding between water and phosphate, are consistent with a mechanism of proton transfer involving H-bonding. The possibility that the phosphate-rich surface of membranes assists in proton equilibration in cells is discussed.

Dynamics of voltage profile in enzymatic ion transporters, demonstrated in electrokinetics of proton pumping rhodopsin

H(+)-pumping rhodopsins mediate a primordial conversion of light to metabolic energy. Bacteriorhodopsin from Halobacterium salinarium is the first identified and (biochemically) best-studied H(+)-pumping rhodopsin. The electrical properties of H(+)-pumping rhodopsins, however, are known in more detail for the homolog Acetabularia rhodopsin, isolated from the eukaryotic green alga Acetabularia acetabulum. Based on data from Acetabularia rhodopsin we present a general reaction kinetic model of H(+)-pumping rhodopsins with only seven independent parameters, which fits the kinetic properties of photocurrents as functions of light, transmembrane voltage, internal and external pH, and time. The model describes fast photoisomerization of retinal with simultaneous H(+) transfer to an H(+) acceptor, reprotonation of retinal from the intracellular face via an H(+) donor, and proton release to the extracellular space via an H(+) release complex. The voltage sensitivities of the individual reaction steps and their temporal changes are treated here by a novel approach, whereby--as in an Ohmic voltage divider--the effective portions of the total transmembrane voltage decrease with the relative velocities of the individual reaction steps. This analysis quantitatively infers dynamic changes of the voltage profile and of the pK values of the H(+)-binding sites involved.

Spatial and temporal control of microwave triggered chemiluminescence: a protein detection platform

We have combined the principles of microwave circuitry and antenna design and our recent work in microwave-triggered metal-enhanced chemiluminescence to now "trigger" chemically and enzyme-catalyzed chemiluminescent reactions with spatial and temporal control. With this technology platform, we achieve spatial and temporal control of enzyme and chemically catalyzed chemiluminescence reactions to achieve more than 500-fold increases in "on-demand" photon flux from chemically catalyzed chemiluminescent reactions. We also report a 6-fold increase in photon flux from HRP-catalyzed assays on disposable coverslips functionalized with HRP and placed proximal to the substrates modified with thin-film aluminum triangle disjointed "bow-tie" structures. In addition, we demonstrate the applicability of this technology to develop multiplexed or high-throughput chemiluminescent assays. We also demonstrate the clinical and biological relevance of this technology platform by affixing aluminum structures in proximity to HRP protein immobilized on nitrocellulose to improve the sensitivity for this model Western blot scheme by 50-fold. We believe analytical applications that rely on enzyme-catalyzed chemiluminescence, such as immunoassays, may greatly benefit from this new platform technology.