Biomedical subjects

Volkhard Helms

Publications and source records attributed to Volkhard Helms.

At least 19 recordsLinked to original sources

PubMed2007

Protein dynamics tightly connected to the dynamics of surrounding and internal water molecules.

Proteins are key components of biological cells. For example, enzymes catalyze biochemical reactions, membrane transporters are responsible for uptake and release of critical and superfluous components from the cell environment, and structural proteins are responsible for the stability of the cell wall and cytoskeleton. Many of the diverse protein functions involve dynamic transitions ranging from small local atomic displacements up to large allosteric conformational changes. In any conformation, proteins are in contact with the universal solvent medium of cells, water. Water not only surrounds proteins but is often an integral part of proteins and also is involved in key mechanistic steps. This Minireview discusses recent experimental and theoretical results on the role of water for protein dynamics and function.

Hydrogen Bonding↗

PubMed2006

Assembly of transmembrane helices of simple polytopic membrane proteins from sequence conservation patterns.

The transmembrane (TM) domains of most membrane proteins consist of helix bundles. The seemingly simple task of TM helix bundle assembly has turned out to be extremely difficult. This is true even for simple TM helix bundle proteins, i.e., those that have the simple form of compact TM helix bundles. Herein, we present a computational method that is capable of generating native-like structural models for simple TM helix bundle proteins having modest numbers of TM helices based on sequence conservation patterns. Thus, the only requirement for our method is the presence of more than 30 homologous sequences for an accurate extraction of sequence conservation patterns. The prediction method first computes a number of representative well-packed conformations for each pair of contacting TM helices, and then a library of tertiary folds is generated by overlaying overlapping TM helices of the representative conformations. This library is scored using sequence conservation patterns, and a subsequent clustering analysis yields five final models. Assuming that neighboring TM helices in the sequence contact each other (but not that TM helices A and G contact each other), the method produced structural models of Calpha atom root-mean-square deviation (CA RMSD) of 3-5 A from corresponding crystal structures for bacteriorhodopsin, halorhodopsin, sensory rhodopsin II, and rhodopsin. In blind predictions, this type of contact knowledge is not available. Mimicking this, predictions were made for the rotor of the V-type Na(+)-adenosine triphosphatase without such knowledge. The CA RMSD between the best model and its crystal structure is only 3.4 A, and its contact accuracy reaches 55%. Furthermore, the model correctly identifies the binding pocket for sodium ion. These results demonstrate that the method can be readily applied to ab initio structure prediction of simple TM helix bundle proteins having modest numbers of TM helices.

Bacteriorhodopsins↗

PubMed2006

Reconstruction of a kinetic model of the chromatophore vesicles from Rhodobacter sphaeroides.

We present a molecular model of a chromatophore vesicle from Rhodobacter sphaeroides. These vesicles are ideal benchmark systems for molecular and systemic simulations, because they have been well studied, they are small, and they are naturally separated from their cellular environment. To set up a photosynthetic chain working under steady-state conditions, we compiled from the experimental literature the specific activities and geometries that have been determined for their constituents. This data then allowed defining the stoichiometries for all membrane proteins. This article contains the kinetic part of the reconstructed model, while the spatial reconstruction is presented in a companion article. By considering the transport properties of the Cytochrome c(2) and ubiquinone pools, we show that their size and oxidation states allow for an efficient buffering of the statistical fluctuations that arise from the small size of the vesicles. Stoichiometric and kinetic considerations indicate that a typical chromatophore vesicle of Rb. sphaeroides with a diameter of 45 nm should contain approximately five bc(1) monomers.

ATP Synthetase Complexes↗

PubMed2006

A spatial model of the chromatophore vesicles of Rhodobacter sphaeroides and the position of the Cytochrome bc1 complex.

The photosynthetic apparatus of purple bacteria is generally considered a well-studied and understood system. However, recent atomic force microscopy images of flattened chromatophore vesicles from Rhodobacter sphaeroides restarted a debate about the stoichiometry and positions of the membrane proteins, with the interpretations of the observed images only partly being in agreement with earlier models. The most puzzling observation from the recent images is that the Cytochrome bc(1) complex, which is a central part of the photosynthetic apparatus, seems to be missing on the chromatophore vesicles, even when these were extracted from photosynthetically grown bacteria. From the available information on the geometry of the vesicle and of the proteins we reconstructed here a three-dimensional model vesicle at molecular resolution. Its central feature, also determining its diameter of approximately 45 nm, is an equatorial array of LH1 dimers, lined by a region of LH2 rings. This naturally puts the Cytochrome bc(1) complexes and the ATPase at the vesicle's poles. This spatial model may explain why the vesicle's endcaps with the bc(1) complexes are lost during the preparatory steps of the imaging process together with the ATPase and are therefore absent from the available images.

Adenosine Triphosphatases↗

PubMed2006

Tandem repeats in the CpG islands of imprinted genes.

In contrast to most genes in mammalian genomes, imprinted genes are monoallelically expressed depending on the parental origin of the alleles. Imprinted gene expression is regulated by distinct DNA elements that exhibit allele-specific epigenetic modifications, such as DNA methylation. These so-called differentially methylated regions frequently overlap with CpG islands. Thus, CpG islands of imprinted genes may contain special DNA elements that distinguish them from CpG islands of biallelically expressed genes. Here, we present a detailed study of CpG islands of imprinted genes in mouse and in human. Our study shows that imprinted genes more frequently contain tandem repeat arrays in their CpG islands than randomly selected genes in both species. In addition, mouse imprinted genes more frequently possess intragenic CpG islands that may serve as promoters of allele-specific antisense transcripts. This feature is much less pronounced in human, indicating an interspecies variability in the evolution of imprinting control elements.

Animals↗

PubMed2006

How strongly do sequence conservation patterns and empirical scales correlate with exposure patterns of transmembrane helices of membrane proteins?

Given the difficulty in determining high-resolution structures of helical membrane proteins, sequence-based prediction methods can be useful in elucidating diverse physiological processes mediated by this important class of proteins. Predicting the angular orientations of transmembrane (TM) helices about the helix axes, based on the helix parameters from electron microscopy data, is a classical problem in this regard. This problem has triggered the development of a number of different empirical scales. Recently, sequence conservation patterns were also made use of for improved predictions. Empirical scales and sequence conservation patterns (collectively termed as "prediction scales") have also found frequent applications in other research areas of membrane proteins: for example, in structure modeling and in prediction of buried TM helices. This trend is expected to grow in the near future unless there are revolutionary developments in the experimental characterization of membrane proteins. Thus, it is timely and imperative to carry out a comprehensive benchmark test over the prediction scales proposed so far to determine their pros and cons. In the current analysis, we use exposure patterns of TM helices as a golden standard, because if one develops a prediction scale that correlates perfectly with exposure patterns of TM helices, it will enable one to predict buried residues (or buried faces) of TM helices with an accuracy of 100%. Our analysis reveals several important points. (1) It demonstrates that sequence conservation patterns are much more strongly correlated with exposure patterns of TM helices than empirical scales. (2) Scales that were specifically parameterized using structure data (structure-based scales) display stronger correlation than hydrophobicity-based scales, as expected. (3) A nonnegligible difference is observed among the structure-based scales in their correlational property, suggesting that not every learning algorithm is equally effective. (4) A straightforward framework of optimally combining sequence conservation patterns and empirical scales is proposed, which reveals that improvements gained from combining the two sources of information are not dramatic in almost all cases. In turn, this calls for the development of fundamentally different scales that capture the essentials of membrane protein folding for substantial improvements.

Algorithms↗

PubMed2006

Energy barriers of proton transfer reactions between amino acid side chain analogs and water from ab initio calculations.

Proton transfer reactions were studied in all titratable pairs of amino acid side chains where, under physiologically reasonable conditions, one amino acid may function as a donor and the other one as an acceptor. Energy barriers for shifting the proton from donor to acceptor atom were calculated by electronic structure methods at the MP2/6-31++G(d,p) level, and the well-known double-well potentials were characterized. The energy difference between both minima can be expressed by a parabola using as argument the donor-acceptor distance R(DA). In this work, the fit parameters of the quadratic expression are determined for each donor-acceptor pair. Moreover, it was found previously that the energy barriers of the reactions can be expressed by an analytical expression depending on the distance between donor and acceptor and the energy difference between donor and acceptor bound states. The validity of this approach is supported by the extensive new data set. This new parameterization of proton transfer barriers between titratable amino acid side chains allows us to very efficiently estimate proton transfer probabilities in molecular modelling studies or during classical molecular dynamics simulation of biomolecular systems.

Amino Acids↗

PubMed2005

Diffusional encounter of barnase and barstar.

We present an analysis of trajectories from Brownian dynamics simulations of diffusional protein-protein encounter for the well-studied system of barnase and barstar. This analysis reveals details about the optimal association pathways, the regions of the encounter complex, possible differences of the pathways for dissociation and association, the coupling of translational and rotation motion, and the effect of mutations on the trajectories. We found that a small free-energy barrier divides the energetically most favorable region into a region of the encounter complex above the barnase binding interface and a region around a second energy minimum near the RNA binding loop. When entering the region of the encounter complex from the region near the RNA binding loop, barstar has to change its orientation to increase the electrostatic attraction between the proteins. By concentrating the analysis on the successful binding trajectories, we found that the region of the second minimum is not essential for the binding of barstar to barnase. Nevertheless, this region may be helpful to steer barstar into the region of the encounter complex. When applying the same analysis to several barnase mutants, we found that single mutations may drastically change the free-energy landscape and may significantly alter the population of the two minima. Therefore, certain protein-protein pairs may require careful adaptation of the positions of encounter and transition states when interpreting mutation effects on kinetic rates of association and/or dissociation.

Bacterial Proteins↗

PubMed2005

Statistical analysis of predominantly transient protein-protein interfaces.

A non-redundant set of 170 protein-protein interfaces of known structure was statistically analyzed for residue and secondary-structure compositions, pairing preferences and side-chain-backbone interaction frequencies. By focussing mainly on transient protein-protein interfaces, the results underline previous findings for protein-protein interfaces but also show some new interesting aspects of transient interfaces. The residue compositions at interfaces found in this study correlate well with the results of other studies. On average, contacts between pairs of hydrophobic and polar residues were unfavorable, and the charged residues tended to pair subject to charge complementarity. Secondary structure composition analysis shows that neither helices nor beta-sheets are dominantly populated at interfaces. Analyzing the pairing preferences of the secondary structure elements revealed a higher affinity within the same elements and alludes to tight packings. In addition, the results for the side-chain and backbone interaction frequencies, which were measured under more stringent conditions, showed a high occurrence of side-chain-backbone interactions. Taking a closer look at the helix and beta-sheet binding frequencies for a given side-chain and backbone interaction underlined the relevance of tight packings. The polarity of interfaces increased with decreasing interface size. These types of information may be useful for scoring complexes in protein-protein docking studies or for prediction of protein-protein interfaces from the sequences alone.

Amino Acids↗

PubMed2005

Are solvation free energies of homogeneous helical peptides additive?

We investigated the additivity of the solvation free energy of amino acids in homogeneous helices of different length in water and in chloroform. Solvation free energies were computed by multiconfiguration thermodynamic integration involving extended molecular dynamics simulations and by applying the generalized-born surface area solvation model to static helix geometries. The investigation focused on homogeneous peptides composed of uncharged amino acids, where the backbone atoms are kept fixed in an ideal helical conformation. We found nonlinearity especially for short peptides, which does not allow a simple treatment of the interaction of amino acids with their surroundings. For homogeneous peptides longer than five residues, the results from both methods are in quite good agreement and solvation energies are to a good extent additive.

Amino Acids↗

PubMed2005

Dynamical binding of proline-rich peptides to their recognition domains.

Recognition of proline-rich sequences plays an important role for the assembly of multi-protein complexes during the course of eukaryotic signal transduction and is mediated by a set of protein folds that share characteristic features. For many complex systems containing proline-rich sequences, multiple binding modes have been found by theoretical and/or experimental studies. In this review, we discuss the different binding modes as well as the correlated dynamics of the peptides and their recognition domains, and some implications to their biological functions. Furthermore, we give an outlook of the systems in the context of systems biology.

Ligands↗

PubMed2005

Alternative binding modes of proline-rich peptides binding to the GYF domain.

Recognition of proline-rich sequences plays an important role for the assembly of multiprotein complexes during the course of eukaryotic signal transduction and is mediated by a set of protein folds that share characteristic features. The GYF (glycine-tyrosine-phenylalanine) domain is known as a member of the superfamily of recognition domains for proline-rich sequences. Recent studies on the complexation of the CD2BP2-GYF domain with CD2 peptides showed that the peptide adopts an extended conformation and forms a polyproline type-II helix involving residues Pro4-Pro7 [Freund et al. (2002) EMBO J. 21, 5985-5995]. R/K/GxxPPGxR/K is the key signature for the peptides that bind to the GYF domain [Kofler et al. (2004) J. Biol. Chem. 279, 28292-28297]. In our combined theoretical and experimental study, we show that the peptides adopt a polyproline II helical conformation in the unbound form as well as in the complex. From molecular dynamics simulations, we identify a novel binding mode for the G8W mutant and the wild-type peptide (shifted by one proline in register). In contrast, the conformation of the peptide mutant H9M remains close to the experimentally derived wild-type GYF-peptide complex. Possible functional implications of this altered conformation of the bound ligand are discussed in the light of our experimental and theoretical results.

Adaptor Proteins, Signal Transducing↗

PubMed2005

Cyclophilin A binds to linear peptide motifs containing a consensus that is present in many human proteins.

Cyclophilin A (CypA) is a peptidyl-prolyl cis/trans-isomerase that is involved in multiple signaling events of eukaryotic cells. It might either act as a catalyst for prolyl bond isomerization, or it can form stoichiometric complexes with target proteins. We have investigated the linear sequence recognition code for CypA by phage display and found the consensus motif FGPXLp to be selected after five rounds of panning. The peptide FGPDLPAGD showed inhibition of the isomerase reaction and NMR chemical shift mapping experiments highlight the CypA interaction epitope. Ligand docking suggests that the peptide was able to bind to CypA in the cis- and trans-conformation. Protein Data Bank searches reveal that many human proteins contain the consensus motif, and several of these protein motifs are shown to interact with CypA in vitro. These sequences represent putative target sites for binding of CypA to intracellular proteins.

Amino Acid Motifs↗

PubMed2005

Molecular dynamics simulation of truncated bovine adrenodoxin.

The 128 amino acid long soluble protein adrenodoxin (Adx) is a typical member of the ferredoxin protein family that are electron carrier proteins with an iron-sulfur cofactor. Adx carries electrons from adrenodoxin reductase (AdR) to cytochrome P450s. Its binding modes to these proteins were previously characterized by site-directed mutagenesis, by X-ray crystallography for the complex Adx:AdR, and by NMR. However, no clear evidence has been provided for the driving force that promotes Adx detachment from AdR upon reduction. Here, we characterized the conformational dynamics of unbound Adx in the oxidized and reduced forms using 2-20 ns long molecular dynamics simulations. The most noticeable difference between both forms is the enhanced flexibility of the loop (47-51) surrounding the iron-sulfur cluster in the reduced form. Together with several structural displacements at the binding interface, this increased flexibility may be the key factor promoting unbinding of reduced Adx from AdR. This points to an intrinsic property of reduced Adx that drives dissociation.

Adrenodoxin↗

PubMed2005

Titration behavior of residues at the entrance of the D-pathway of cytochrome c oxidase from paracoccus denitrificans investigated by continuum electrostatic calculations.

Continuum electrostatic calculations were employed to investigate the titration curves of the fully oxidized state of wild type and several variants of cytochrome c oxidase from Paracoccus denitrificans (N131D, N131C, N131V, and D124N) for different values of the dielectric constant of the protein. The effects of the mutations at the entrance of the D-proton transfer pathway were found to be quite localized to their immediate surroundings. The results can be well interpreted in the light of the available biochemical and structural data and help understanding the effects of mutations on proton conductivity. The mutations of aspartic acid Asp-I-124 to a neutral residue resulted in a decreased pK(a) value of His-I-28 suggesting that the mutation of His-I-28 may have a significant influence on the coupling of electron and proton transfer in cytochrome c oxidase. We also investigated the effect of the mutations N131D, N131C, and N131V on the residue Glu-I-278 in terms of its pK(a) value and electrostatic interaction energies.

Amino Acid Substitution↗

PubMed2004

Novel scoring function for modeling structures of oligomers of transmembrane alpha-helices.

Specific non-covalent interactions between transmembrane (TM) alpha-helices are important in a variety of biological processes. Experimental and computational studies have shown that van der Waals interactions play an important role in the tight packing between TM alpha-helices, although polar interactions can also be important in some instances. Based on the assumption that van der Waals interaction alone is sufficient for a meso-scale (residue-scale) description of the interaction between TM alpha-helices, we have designed a novel residue-scale scoring function for modeling structures of oligomers of TM alpha-helices. We first calculated atomistic van der Waals interaction energies between two amino acids, X and Y, of a pair of parallel alpha-helices, glycine-X-glycine and glycine-Y-glycine and compiled them according to three variables, the distance between the two C(alpha) atoms and the rotational angles of X and Y about their helical axes. Upon averaging over the rotational angles, we obtained one-dimensional interaction energy profiles that are functions of the distance between C(alpha) atoms only. Each of the interaction energy profiles was fitted with a generic fitting function of the distance between C(alpha) atoms, yielding analytical scoring functions for all possible amino acid pairs. For glycophorin A, neu/erbB-2, and phospholamban, lowest-energy conformations obtained through exhaustive scanning of the entire conformational space using the scoring functions were compatible with available experimental data.

Amino Acid Sequence↗

PubMed2004

Dynamic water networks in cytochrome C oxidase from Paracoccus denitrificans investigated by molecular dynamics simulations.

We present a molecular dynamics study of cytochrome c oxidase from Paracoccus denitrificans in the fully oxidized state, embedded in a fully hydrated dimyristoylphosphatidylcholine lipid bilayer membrane. Parallel simulations with different levels of protein hydration, 1.125 ns each in length, were carried out under conditions of constant temperature and pressure using three-dimensional periodic boundary conditions and full electrostatics to investigate the distribution and dynamics of water molecules and their corresponding hydrogen-bonded networks inside cytochrome c oxidase. The majority of the water molecules had residence times shorter than 100 ps, but a few water molecules are fixed inside the protein for up to 1.125 ns. The hydrogen-bonded network in cytochrome c oxidase is not uniformly distributed, and the degree of water arrangement is variable. The average number of solvent sites in the proton-conducting K- and D-pathways was determined. In contrast to single water files in narrow geometries we observe significant diffusion of individual water molecules along these pathways. The highly fluctuating hydrogen-bonded networks, combined with the significant diffusion of individual water molecules, provide a basis for the transfer of protons in cytochrome c oxidase, therefore leading to a better understanding of the mechanism of proton pumping.

Computer Simulation↗

PubMed2004

A Brownian dynamics study: the effect of a membrane environment on an electron transfer system.

During the past few years, three-dimensional crystal structures of many of the important integral membrane proteins responsible for the bioenergetic processes of photosynthesis and respiration have been determined. Moreover, a few crystal structures of protein-protein complexes have become available that characterize the interaction between those membrane proteins and the electron carrier protein cytochrome c. Here, we address the association kinetics for binding of cytochrome c to cytochrome c oxidase (COX) from Paracoccus denitrificans by Brownian dynamics simulations. The effects of ionic strength and protein mutations were studied for two different cytochrome c species: the positively charged, dipolar horse heart cytochrome c and the negatively charged physiological electron transfer partner cytochrome c(552). We studied association toward "naked" COX and toward membrane-embedded COX where the membrane is represented as an uncharged DPPC bilayer modeled in atomistic detail. For the nonnatural association toward "naked" COX, the association rates are >100 times larger for horse heart cytochrome c than for cytochrome c(552). Interestingly, the presence of the lipid bilayer leads to a dramatic decrease of the association rate of horse heart cytochrome c, but slightly enhances association of cytochrome c(552), leading to very similar association rates of both proteins to membrane-embedded COX. This finding from computational modeling studies may reflect the optimization of surface patches and of the total net charge on electron transfer pairs in nature.

Animals↗