A new method is presented for simulating the simultaneous binding equilibrium of electrons and protons on protein molecules, which makes it possible to study the full equilibrium thermodynamics of redox and protonation processes, including electron-proton coupling. The simulations using this method reflect directly the pH and electrostatic potential of the environment, thus providing a much closer and realistic connection with experimental parameters than do usual methods. By ignoring the full binding equilibrium, calculations usually overlook the twofold effect that binding fluctuations have on the behavior of redox proteins: first, they affect the energy of the system by creating partially occupied sites; second, they affect its entropy by introducing an additional empty/occupied site disorder (here named occupational entropy). The proposed method is applied to cytochrome c3 of Desulfovibrio vulgaris Hildenborough to study its redox properties and electron-proton coupling (redox-Bohr effect), using a continuum electrostatic method based on the linear Poisson-Boltzmann equation. Unlike previous studies using other methods, the full reduction order of the four hemes at physiological pH is successfully predicted. The sites more strongly involved in the redox-Bohr effect are identified by analysis of their titration curves/surfaces and the shifts of their midpoint redox potentials and pKa values. Site-site couplings are analyzed using statistical correlations, a method much more realistic than the usual analysis based on direct interactions. The site found to be more strongly involved in the redox-Bohr effect is propionate D of heme I, in agreement with previous studies; other likely candidates are His67, the N-terminus, and propionate D of heme IV. Even though the present study is limited to equilibrium conditions, the possible role of binding fluctuations in the concerted transfer of protons and electrons under nonequilibrium conditions is also discussed. The occupational entropy contributions to midpoint redox potentials and pKa values are computed and shown to be significant.

Solution pH is a determinant parameter on protein function and stability, and its inclusion in molecular dynamics simulations is attractive for studies at the molecular level. Current molecular dynamics simulations can consider pH only in a very limited way, through a somewhat arbitrary choice of a set of fixed charges on the titrable sites. Conversely, continuum electrostatic methods that explicitly treat pH effects assume a single protein conformation whose choice is not clearly defined. In this paper we describe a general method that combines both titration and conformational freedom. The method is based on a potential of mean force for implicit titration and combines both usual molecular dynamics and pH-dependent calculations based on continuum methods. A simple implementation of the method, using a mean field approximation, is presented and applied to the bovine pancreatic trypsin inhibitor. We believe that this constant-pH molecular dynamics method, by correctly sampling both charges and conformation, can become a valuable help in the understanding of the dependence of protein function and stability on pH.

Natural populations of widely distributed organisms often exhibit genetic clinal variation over their geographical ranges. The European anchovy, Engraulis encrasicolus, illustrates this by displaying a two-clade mitochondrial structure clinally arranged along the eastern Atlantic. One clade has low frequencies at higher latitudes, whereas the other has an anti-tropical distribution, with frequencies decreasing towards the tropics. The distribution pattern of these clades has been explained as a consequence of secondary contact after an ancient geographical isolation. However, it is not unlikely that selection acts on mitochondria whose genes are involved in relevant oxidative phosphorylation processes. In this study, we performed selection tests on a fragment of 1044 bp of the mitochondrial cytochrome b gene using 455 individuals from 18 locations. We also tested correlations of six environmental features: temperature, salinity, apparent oxygen utilization and nutrient concentrations of phosphate, nitrate and silicate, on a compilation of mitochondrial clade frequencies from 66 sampling sites comprising 2776 specimens from previously published studies. Positive selection in a single codon was detected predominantly (99%) in the anti-tropical clade and temperature was the most relevant environmental predictor, contributing with 59% of the variance in the geographical distribution of clade frequencies. These findings strongly suggest that temperature is shaping the contemporary distribution of mitochondrial DNA clade frequencies in the European anchovy.

Solution pH is a determinant parameter on protein function and stability, and its inclusion in molecular dynamics simulations is attractive for studies at the molecular level. Current molecular dynamics simulations can consider pH only in a very limited way, through a somewhat arbitrary choice of a set of fixed charges on the titrable sites. Conversely, continuum electrostatic methods that explicitly treat pH effects assume a single protein conformation whose choice is not clearly defined. In this paper we describe a general method that combines both titration and conformational freedom. The method is based on a potential of mean force for implicit titration and combines both usual molecular dynamics and pH-dependent calculations based on continuum methods. A simple implementation of the method, using a mean field approximation, is presented and applied to the bovine pancreatic trypsin inhibitor. We believe that this constant-pH molecular dynamics method, by correctly sampling both charges and conformation, can become a valuable help in the understanding of the dependence of protein function and stability on pH.

A comparative study of the pH-dependent redox mechanisms of several members of the cytochrome c3 family has been carried out. In a previous work, the molecular determinants of this dependency (the so-called redox-Bohr effect) were investigated for one species using continuum electrostatic methods to find groups with a titrating range and strength of interaction compatible with a mediating role in the redox-Bohr effect. Here we clarify these aspects in the light of new and improved pKa calculations, our findings supporting the hypothesis of propionate D from heme I being the main effector in the pH-dependent modulation of the cytochrome c3 redox potentials in all the c3 molecules studied here. However, the weaker (but significant) role of other titrating groups cannot be excluded, their importance and identity changing with the particular molecule under study. We also calculate the relative redox potentials of the four heme centers among the selected members of the c3 family, using a continuum electrostatic method that takes into account both solvation and interaction effects. Comparison of the calculated values with available data for the microscopic redox potentials was undertaken, the quality of the agreement being dependent upon the choice of the dielectric constant for the protein interior. We find that high dielectric constants give best correlations, while low values result in better magnitudes for the calculated potentials. The possibility that the crystallographic calcium ion in c3 from Desulfovibrio gigas may be present in the solution structure was tested, and found to be likely.

Decavanadate, a vanadate oligomer, is known to interact with myosin and to inhibit the ATPase activity, but the putative binding sites and the mechanism of inhibition are still to be clarified. We have previously proposed that the decavanadate (V 10 O 28 6− ) inhibition of the actin-stimulated myosin ATPase activity is non-competitive towards both actin and ATP. A likely explanation for these results is that V 10 binds to the so-called back-door at the end of the Pi-tube opposite to the nucleotide-binding site. In order to further investigate this possibility, we have carried out molecular docking simulations of the V 10 oligomer on three different structures of the myosin motor domain of Dictyostelium discoideum, representing distinct states of the ATPase cycle. The results indicate a clear preference of V 10 to bind at the back-door, but only on the "open" structures where there is access to the phosphate binding-loop. It is suggested that V 10 acts as a "back-door stop" blocking the closure of the 50-kDa cleft necessary to carry out ATP-γ-phosphate hydrolysis. This provides a simple explanation to the non-competitive behavior of V 10 and spurs the use of the oligomer as a tool to elucidate myosin back-door conformational changes in the process of muscle contraction.

No abstract

No abstract