G. Matthias Ullmann

Visualization of the Encounter Ensemble of the Transient Electron Transfer Complex of Cytochrome c and Cytochrome c Peroxidase

Recent studies have provided experimental evidence for the existence of an encounter complex, a transient intermediate in the formation of protein complexes. We use paramagnetic relaxation enhancement NMR spectroscopy in combination with Monte Carlo simulations to characterize and visualize the ensemble of encounter orientations in the short-lived electron transfer complex of yeast cytochrome c (Cc) and cytochrome c peroxidase (CcP). The complete conformational space sampled by the protein molecules during the dynamic part of the interaction was mapped experimentally. The encounter complex was described by an electrostatic ensemble of orientations based on Monte Carlo calculations, considering the protein structures in atomic detail. We demonstrate that this visualization of the encounter complex, in combination with the specific complex, is in excellent agreement with the experimental data. Our results indicate that Cc samples only about 15% of the surface area of CcP, surrounding the specific binding interface. The encounter complex is populated for 30% of the time, representing a mere 0.5 kcal/mol difference in the free energies between the two states. This delicate balance is interpreted to be a consequence of the conflicting requirements of fast electron transfer and high turnover of the complex.

Structure of the non-redox-active tungsten/[4Fe:4S] enzyme acetylene hydratase

The tungsten–iron–sulfur enzyme acetylene hydratase stands out from its class because it catalyzes a nonredox reaction, the hydration of acetylene to acetaldehyde. Sequence comparisons group the protein into the dimethyl sulfoxide reductase family, and it contains a bis-molybdopterin guanine dinucleotide-ligated tungsten atom and a cubane-type [4Fe:4S] cluster. The crystal structure of acetylene hydratase at 1.26 Å now shows that the tungsten center binds a water molecule that is activated by an adjacent aspartate residue, enabling it to attack acetylene bound in a distinct, hydrophobic pocket. This mechanism requires a strong shift of pKa of the aspartate, caused by a nearby low-potential [4Fe:4S] cluster. To access this previously unrecognized W–Asp active site, the protein evolved a new substrate channel distant from where it is found in other molybdenum and tungsten enzymes.

McVol - A program for calculating protein volumes and identifying cavities by a Monte Carlo algorithm

In this paper, we describe a Monte Carlo method for determining the volume of a molecule. A molecule is considered to consist of hard, overlapping spheres. The surface of the molecule is defined by rolling a probe sphere over the surface of the spheres. To determine the volume of the molecule, random points are placed in a three-dimensional box, which encloses the whole molecule. The volume of the molecule in relation to the volume of the box is estimated by calculating the ratio of the random points placed inside the molecule and the total number of random points that were placed. For computational efficiency, we use a grid-cell based neighbor list to determine whether a random point is placed inside the molecule or not. This method in combination with a graph-theoretical algorithm is used to detect internal cavities and surface clefts of molecules. Since cavities and clefts are potential water binding sites, we place water molecules in the cavities. The potential water positions can be used in molecular dynamics calculations as well as in other molecular calculations. We apply this method to several proteins and demonstrate the usefulness of the program. The described methods are all implemented in the program McVol, which is available free of charge from our website at http://www.bisb.uni-bayreuth.de/software.html.

pH-Dependent pKa Values in Proteins—A Theoretical Analysis of Protonation Energies with Practical Consequences for Enzymatic Reactions

Because of their central importance for understanding enzymatic mechanisms, pK(a) values are of great interest in biochemical research. It is common practice to determine pK(a) values of amino acid residues in proteins from NMR or FTIR titration curves by determining the pH at which the protonation probability is 50%. The pH dependence of the free energy required to protonate this residue is then determined from the linear relationship DeltaG(prot) = RT ln 10 (pH-pK(a)), where R is the gas constant and T the absolute temperature. However, this approach neglects that there can be important electrostatic interactions in the proteins that can shift the protonation energy. Even if the titration curves seem to have a standard sigmoidal shape, the protonation energy of an individual site in a protein may depend nonlinearly on pH. To account for this nonlinear dependence, we show that it is required to introduce pK(a) values for individual sites in proteins that depend on pH. Two different definitions are discussed. One definition is based on a rearranged Henderson-Hasselbalch equation, and the other definition is based on an equation that was used by Tanford and Roxby to approximate titration curves of proteins. In the limiting case of weak interactions, the two definitions lead to nearly the same pK(a) value. We discuss how these two differently defined pK(a) values are related to the free energy change required to protonate a site. Using individual site pK(a) values, we demonstrate on simple model systems that the interactions between protonatable residues in proteins can help to maintain the energy required to protonate a site in the protein nearly constant over a wide pH range. We show with the example of RNase T1 that such a mechanism to keep the protonation energy constant is used in enzymes. The pH dependence of pK(a) values may be an important concept in enzyme catalysis. Neglecting this concept, important features of enzymes may be missed, and the enzymatic mechanism may not be fully understood.

Structural Basis for a Kolbe-Type Decarboxylation Catalyzed by a Glycyl Radical Enzyme.

4-Hydroxyphenylacetate decarboxylase is a [4Fe-4S] cluster containing glycyl radical enzyme proposed to use a glycyl/thiyl radical dyad to catalyze the last step of tyrosine fermentation in clostridia. The decarboxylation product p-cresol (4-methylphenol) is a virulence factor of the human pathogen Clostridium difficile . Here we describe the crystal structures at 1.75 and 1.81 Å resolution of substrate-free and substrate-bound 4-hydroxyphenylacetate decarboxylase from the related Clostridium scatologenes . The structures show a (βγ)(4) tetramer of heterodimers composed of a catalytic β-subunit harboring the putative glycyl/thiyl dyad and a distinct small γ-subunit with two [4Fe-4S] clusters at 40 Å distance from the active site. The γ-subunit comprises two domains displaying pseudo-2-fold symmetry that are structurally related to the [4Fe-4S] cluster-binding scaffold of high-potential iron-sulfur proteins. The N-terminal domain coordinates one cluster with one histidine and three cysteines, and the C-terminal domain coordinates the second cluster with four cysteines. Whereas the C-terminal cluster is buried in the βγ heterodimer interface, the N-terminal cluster is not part of the interface. The previously postulated decarboxylation mechanism required the substrates hydroxyl group in the vicinity of the active cysteine residue. In contrast to expectation, the substrate-bound state shows a direct interaction between the substrates carboxyl group and the active site Cys503, while His536 and Glu637 at the opposite side of the active site pocket anchor the hydroxyl group. This state captures a possible catalytically competent complex and suggests a Kolbe-type decarboxylation for p-cresol formation.

Shifting the Equilibrium between the Encounter State and the Specific Form of a Protein Complex by Interfacial Point Mutations

Recent experimental studies have confirmed a long-held view that protein complex formation proceeds via a short-lived encounter state. The population of this transient intermediate, stabilized mainly by long-range electrostatic interactions, varies among different complexes. Here we show that the occupancy of the encounter state can be modulated across a broad range by single point mutations of interfacial residues. Using a combination of Monte Carlo simulations and paramagnetic relaxation enhancement NMR spectroscopy, we illustrate that it is possible to both enhance and diminish the binding specificity in an electron transfer complex of yeast cytochrome c (Cc) and cytochrome c peroxidase. The Cc T12A mutation decreases the population of the encounter to 10% as compared with 30% in the wild-type complex. More dramatically, the Cc R13A substitution reverses the relative occupancies of the stereospecific and the encounter forms, with the latter now being the dominant species with the population of 80%. This finding indicates that the encounter state can make a large contribution to the stability of a protein complex. Also, it appears that by adjusting the amount of the encounter through a judicious choice of point mutations, we can remodel the energy landscape of a protein complex and tune its binding specificity.

A molecular mechanics force field for biologically important sterols.

A parameterization has been performed of the biologically important sterols cholesterol, ergosterol, and lanosterol for the CHARMM27 all‐atom molecular mechanics force field. An automated parameterization method was used that involves fitting the potential to vibrational frequencies and eigenvectors derived from quantum‐chemical calculations. The partial charges were derived by fitting point charges to quantum‐chemically calculated electrostatic potentials. To model the dynamics of the hydroxyl groups of the sterols correctly, the parameter set was refined to reproduce the energy barrier for the rotation of the hydroxyl group around the carbon connected to the hydroxyl of each sterol. The frequency‐matching plots show good agreement between the CHARMM and quantum chemical normal modes. The parameters are tested in a molecular dynamics simulation of the cholesterol crystal structure. The experimental geometry and cell dimensions are well reproduced. The force field derived here is also useful for simulating other sterols such as the phytosterols sigmasterol, and campesterol, and a variety of steroids.

GMCT : A Monte Carlo simulation package for macromolecular receptors

Generalized Monte Carlo titration (GMCT) is a versatile suite of computer programs for the efficient simulation of complex macromolecular receptor systems as for example proteins. The computational model of the system is based on a microstate description of the receptor and an average description of its surroundings in terms of chemical potentials. The receptor can be modeled in great detail including conformational flexibility and many binding sites with multiple different forms that can bind different ligand types. Membrane embedded systems can be modeled including electrochemical potential gradients. Overall properties of the receptor as well as properties of individual sites can be studied with a variety of different Monte Carlo (MC) simulation methods. Metropolis MC, Wang‐Landau MC and efficient free energy calculation methods are included. GMCT is distributed as free open source software at www.bisb.uni‐bayreuth.de under the terms of the GNU Affero General Public License.

Understanding the energetics of helical peptide orientation in membranes.

Understanding the energetic factors determining the positioning and orientation of single‐helical peptides in membranes is of fundamental interest in structural biology. Here, a simple 5‐slab continuum dielectric model for the membrane is examined that distinguishes between the solvent, headgroup, and core regions. An analytical solution for the electrostatic solvation of a single dipole and an all‐atom model of N‐methylacetamide are used to demonstrate the effect of the dielectric boundaries in the system on peptide dipole orientation. The dipole orientation energy is shown to dominate the electrostatic solvation energy of a polyalanine helix in the membrane. With an additional surface‐area‐dependent term to account for the cavity formation in the aqueous region, the continuum electrostatics description is used to examine several helical peptides, the atoms of which are explicitly represented with a molecular mechanics force field. The experimentally determined tilt angles of a number of peptides of alternating alanine and leucine residues, and of glycophorin and melittin, are accurately reproduced by the model. The factors determining the tilt angles and their fluctuations are analyzed. The tilt angles of the simpler peptides are found to increase approximately linearly with peptide length; this effect is also rationalized. The analysis and model presented here provide a step toward the prediction of helical membrane protein structure. Proteins 2005.

Electrostatic potential at the retinal of three archaeal rhodopsins: Implications for their different absorption spectra

The color tuning mechanism of the rhodopsin protein family has been in the focus of research for decades. However, the structural basis of the tuning mechanism in general and of the absorption shift between rhodopsins in particular remains under discussion. It is clear that a major determinant for spectral shifts between different rhodopsins are electrostatic interactions between the chromophore retinal and the protein. Based on the Poisson‐Boltzmann equation, we computed and compared the electrostatic potential at the retinal of three archaeal rhodopsins: bacteriorhodopsin (BR), halorhodopsin (HR), and sensory rhodopsin II (SRII) for which high‐resolution structures are available. These proteins are an excellent test case for understanding the spectral tuning of retinal. The absorption maxima of BR and HR are very similar, whereas the spectrum of SRII is considerably blue shifted—despite the structural similarity between these three proteins. In agreement with their absorption maxima, we find that the electrostatic potential is similar in BR and HR, whereas significant differences are seen for SRII. The decomposition of the electrostatic potential into contributions of individual residues, allowed us to identify seven residues that are responsible for the differences in electrostatic potential between the proteins. Three of these residues are located in the retinal binding pocket and have in fact been shown to account for part of the absorption shift between BR and SRII by mutational studies. One residue is located close to the β‐ionone ring of retinal and the remaining three residues are more than 8 Å away from the retinal. These residues have not been discussed before, because they are, partly because of their location, no obvious candidates for the spectral shift among BR, HR, and SRII. However, their contribution to the differences in electrostatic potential is evident. The counterion of the Schiff base, which is frequently discussed to be involved in the spectral tuning, does not contribute to the dissimilarities between the electrostatic potentials. Proteins 2005.