Biomedical subjects

Günther H Peters

Publications and source records attributed to Günther H Peters.

17 recordsLinked to original sources

PubMed2007

Simulation of the coupling between nucleotide binding and transmembrane domains in the ATP binding cassette transporter BtuCD.

The nucleotide-induced structural rearrangements in ATP binding cassette (ABC) transporters, leading to substrate translocation, are largely unknown. We have modeled nucleotide binding and release in the vitamin B(12) importer BtuCD using perturbed elastic network calculations and biased molecular dynamics simulations. Both models predict that nucleotide release decreases the tilt between the two transmembrane domains and opens the cytoplasmic gate. Nucleotide binding has the opposite effect. The observed coupling may be relevant for all ABC transporters because of the conservation of nucleotide binding domains and the shared role of ATP in ABC transporters. The rearrangements in the cytoplasmic gate region do not provide enough space for B(12) to diffuse from the transporter pore into the cytoplasm, which could suggest that peristaltic forces are needed to exclude B(12) from the transporter pore.

ATP-Binding Cassette Transporters↗

PubMed2006

Ammonium recruitment and ammonia transport by E. coli ammonia channel AmtB.

To investigate substrate recruitment and transport across the Escherichia coli Ammonia transporter B (AmtB) protein, we performed molecular dynamics simulations of the AmtB trimer. We have identified residues important in recruitment of ammonium and intraluminal binding sites selective of ammonium, which provide a means of cation selectivity. Our results indicate that A162 guides translocation of an extraluminal ammonium into the pore lumen. We propose a mechanism for transporting the intraluminally recruited proton back to periplasm. Our mechanism conforms to net transport of ammonia and can explain why ammonia conduction is lost upon mutation of the conserved residue D160. We unify previous suggestions of D160 having either a structural or an ammonium binding function. Finally, our simulations show that the channel lumen is hydrated from the cytoplasmic side via the formation of single file water, while the F107/F215 stack at the inner-most part of the periplasmic vestibule constitutes a hydrophobic filter preventing AmtB from conducting water.

Ammonia↗

PubMed2006

Molecular packing in 1-hexanol-DMPC bilayers studied by molecular dynamics simulation.

The structure and molecular packing density of a "mismatched" solute, 1-hexanol, in lipid membranes of dimyristoyl phosphatidylcholine (DMPC) was studied by molecular dynamics simulations. We found that the average location and orientation of the hexanol molecules matched earlier experimental data on comparable systems. The local density or molecular packing in DMPC-hexanol was elucidated through the average Voronoi volumes of all heavy (non-hydrogen) atoms. Analogous analysis was conducted on trajectories from simulations of pure 1-hexanol and pure (hydrated) DMPC bilayers. The results suggested a positive volume change, DeltaV(m), of 4 cm(3) mol(-1) hexanol partitioned at 310 K in good accordance with experimental values. Analysis of the apparent volumes of each component in the pure and mixed states further showed that DeltaV(m) reflects a balance between a substantial increase in the packing density of the alcohol upon partitioning and an even stronger loosening in the packing of the lipid. Furthermore, analysis of Voronoi volumes along the membrane normal identifies a distinctive depth dependence of the changes in molecular packing. The outer (interfacial) part of the lipid acyl chains (up to C8) is stretched by about 4%. Concomitantly, the average lateral area per chain decreases and these two effects compensate so that the overall packing density in the outer region, where the hexanol molecules are located, remains practically constant. The core of the bilayer (C9-C13) is slightly thinned. The average lateral area per chain in this region expands, resulting in a looser packing density. The net effect in the core is a 2-3% decrease in density corresponding to a total volume increase of approximately 14 cm(3) mol(-1) hexanol partitioned.

Computer Simulation↗

PubMed2006

The effect of calcium on the properties of charged phospholipid bilayers.

We have performed molecular dynamics simulations to investigate the structure and dynamics of charged bilayers as well as the distribution of counterions at the bilayer interface. For this, we have considered the negatively charged di-myristoyl-phosphatidyl-glycerol (DMPG) and di-myristoyl-phosphatidyl-serine (DMPS) bilayers as well as a protonated di-myristoyl-phosphatidyl-serine (DMPSH) bilayer. We were particularly interested in calcium ions due to their important role in biological systems. Simulations performed in the presence of calcium ions (DMPG, DMPS) or sodium ions (DMPS) were run for 45-60 ns. Simulation results for DMPG are compared with fluorescence measurements. The average areas per molecule were 47.4+/-0.5 A2 (DMPG with calcium), 47.3+/-0.5 A2 (DMPS with calcium), 51.3+/-1.0 A2 (DMPS with sodium) and 45.3+/-0.5 A2 (DMPSH). The structure of the negatively charged lipids is significantly affected by the counterions, where calcium ions have a more pronounced effect than sodium ions. Calcium ions were found to be tightly bound to the anionic groups of the lipid molecules and as such appear to constitute an integral part of the membrane interface on nanoseconds time scales. In contrast to sodium ions, calcium ions are localised in a narrow (approximately 10 A) band around the phosphate group. The interaction of calcium with the lipid molecules enhances the molecular packing of the PG and PS lipids. This observation is in good agreement with emission spectra of the membrane partitioning probe Laurdan in DMPG multilamellar vesicles that indicate an increase in the ordering of the DMPG bilayer due to the presence of calcium. Our results indicate that calcium ions, which often function as a second messengers in living cells have a pronounced effect on membrane structures, which may have implications during signal transduction events.

Calcium↗

PubMed2006

Domain-induced activation of human phospholipase A2 type IIA: local versus global lipid composition.

Secretory human phospholipase A2 type IIA (PLA2-IIA) catalyzes the hydrolysis of the sn-2 ester bond in glycerolipids to produce fatty acids and lysolipids. The enzyme is coupled to the inflammatory response, and its specificity toward anionic membrane interfaces suggests a role as a bactericidal agent. PLA2-IIA may also target perturbed native cell membranes that expose anionic lipids to the extracellular face. However, anionic lipid contents in native cells appear lower than the threshold levels necessary for activation. By using phosphatidylcholine/phosphatidylglycerol model systems, we show that local enrichment of anionic lipids into fluid domains triggers PLA2-IIA activity. In addition, the compositional range of enzyme activity is shown to be related to the underlying lipid phase diagram. A comparison is done between PLA2-IIA and snake venom PLA2, which in contrast to PLA2-IIA hydrolyzes both anionic and zwitterionic membranes. In general, this work shows that PLA2-IIA activation can be accomplished through local enrichment of anionic lipids into domains, indicating a mechanism for PLA2-IIA to target perturbed native membranes with low global anionic lipid contents. The results also show that the underlying lipid phase diagram, which determines the lipid composition at a local level, can be used to predict PLA2-IIA activity.

Animals↗

PubMed2005

Methodological problems in pressure profile calculations for lipid bilayers.

From molecular dynamics simulations of a dipalmitoyl-phosphatidyl-choline (DPPC) lipid bilayer in the liquid crystalline phase, pressure profiles through the bilayer are calculated by different methods. These profiles allow us to address two central and unresolved problems in pressure profile calculations: The first problem is that the pressure profile is not uniquely defined since the expression for the local pressure involves an arbitrary choice of an integration contour. We have investigated two different choices leading to the Irving-Kirkwood (IK) and Harasima (H) expressions for the local pressure tensor. For these choices we find that the pressure profile is almost independent of the contour used, which indicates that the local pressure is well defined for a DPPC bilayer in the liquid crystalline phase. This may not be the case for other systems and we therefore suggest that both the IK and H profiles are calculated in order to test the uniqueness of the profile. The second problem is how to include electrostatic interactions in pressure profile calculations when the simulations are conducted without truncating the electrostatic potential, i.e., using the Ewald summation technique. Based on the H expression for the local pressure, we present a method for calculating the contribution to the lateral components of the local pressure tensor from electrostatic interactions evaluated by the Ewald summation technique. Pressure profiles calculated with an electrostatic potential truncation (cutoff) from simulations conducted with Ewald summation are shown to depend on the cutoff in a subtle manner which is attributed to the existence of long-ranged charge ordering in the system. However, the pressure profiles calculated with relatively long cutoffs are qualitatively similar to the Ewald profile for the DPPC bilayer studied here.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2004

Residue 259 in protein-tyrosine phosphatase PTP1B and PTPalpha determines the flexibility of glutamine 262.

To study the flexibility of the substrate-binding site and in particular of Gln262, we have performed adiabatic conformational search and molecular dynamics simulations on the crystal structure of the catalytic domain of wild-type protein-tyrosine phosphatase (PTP) 1B, a mutant PTP1B(R47V,D48N,M258C,G259Q), and a model of the catalytically active form of PTPalpha. For each molecule two cases were modeled: the Michaelis-Menten complex with the substrate analogue p-nitrophenyl phosphate (p-PNPP) bound to the active site and the cysteine-phosphor complex, each corresponding to the first and second step of the phosphate hydrolysis. Analyses of the trajectories revealed that in the cysteine-phosphor complex of PTP1B, Gln262 oscillates freely between the bound phosphate group and Gly259 frequently forming, as observed in the crystal structure, a hydrogen bond with the backbone oxygen of Gly259. In contrast, the movement of Gln262 is restricted in PTPalpha and the mutant due to interactions with Gln259 reducing the frequency of the oscillation of Gln262 and thereby delaying the positioning of this residue for the second step in the catalysis, as reflected experimentally by a reduction in k(cat). Additionally, in the simulation with the Michaelis-Menten complexes, we found that a glutamine in position 259 induces steric hindrance by pushing the Gln262 side chain further toward the substrate and thereby negatively affecting K(m) as indicated by kinetic studies. Detailed analysis of the water structure around Gln262 and the active site Cys215 reveals that the probability of finding a water molecule correctly positioned for catalysis is much larger in PTP1B than in PTP1B(R47V,D48N,M258C,G259Q) and PTPalpha, in accordance with experiments.

Animals↗

PubMed2004

The hydrophobic effect: molecular dynamics simulations of water confined between extended hydrophobic and hydrophilic surfaces.

Structural and dynamic properties of water confined between two parallel, extended, either hydrophobic or hydrophilic crystalline surfaces of n-alkane C(36)H(74) or n-alcohol C(35)H(71)OH, are studied by molecular dynamics simulations. Electron density profiles, directly compared with corresponding experimental data from x-ray reflectivity measurements, reveal a uniform weak de-wetting characteristic for the extended hydrophobic surface, while the hydrophilic surface is weakly wetted. These microscopic data are consistent with macroscopic contact angle measurements. Specific water orientation is present at both surfaces. The ordering is characteristically different between the surfaces and of longer range at the hydrophilic surface. Furthermore, the dynamic properties of water are different at the two surfaces and different from the bulk behavior. In particular, at the hydrophobic surface, time-correlation functions reveal that water molecules have characteristic diffusive behavior and orientational ordering due to the lack of hydrogen bonding interactions with the surface. These observations suggest that the altered dynamical properties of water in contact with extended hydrophobic surfaces together with a partial drying of the surfaces are more indicative of the hydrophobic effect than structural ordering, which we suggest to be independent of surface topology.

Journal Article↗

PubMed2004

Mapping of epitopes for autoantibodies to the type 1 diabetes autoantigen IA-2 by peptide phage display and molecular modeling: overlap of antibody and T cell determinants.

IA-2 is a major target of autoimmunity in type 1 diabetes. IA-2 responsive T cells recognize determinants within regions represented by amino acids 787-817 and 841-869 of the molecule. Epitopes for IA-2 autoantibodies are largely conformational and not well defined. In this study, we used peptide phage display and homology modeling to characterize the epitope of a monoclonal IA-2 Ab (96/3) from a human type 1 diabetic patient. This Ab competes for IA-2 binding with Abs from the majority of patients with type 1 diabetes and therefore binds a region close to common autoantibody epitopes. Alignment of peptides obtained after screening phage-displayed peptide libraries with purified 96/3 identified a consensus binding sequence of Asn-x-Glu-x-x-(aromatic)-x-x-Gly. The predicted surface on a three-dimensional homology model of the tyrosine phosphatase domain of IA-2 was analyzed for clusters of Asn, Glu, and aromatic residues and amino acids contributing to the epitope investigated using site-directed mutagenesis. Mutation of each of amino acids Asn(858), Glu(836), and Trp(799) reduced 96/3 Ab binding by >45%. Mutations of these residues also inhibited binding of serum autoantibodies from IA-2 Ab-positive type 1 diabetic patients. This study identifies a region commonly recognized by autoantibodies in type 1 diabetes that overlaps with dominant T cell determinants.

Amino Acid Sequence↗

PubMed2004

Residue 182 influences the second step of protein-tyrosine phosphatase-mediated catalysis.

Previous enzyme kinetic and structural studies have revealed a critical role for Asp181 (PTP1B numbering) in PTP (protein-tyrosine phosphatase)-mediated catalysis. In the E-P (phosphoenzyme) formation step, Asp181 functions as a general acid, while in the E-P hydrolysis step it acts as a general base. Most of our understanding of the role of Asp181 is derived from studies with the Yersinia PTP and the mammalian PTP1B, and to some extent also TC (T-cell)-PTP and the related PTPa and PTPe. The neighbouring residue 182 is a phenylalanine in these four mammalian enzymes and a glutamine in Yersinia PTP. Surprisingly, little attention has been paid to the fact that this residue is a histidine in most other mammalian PTPs. Using a reciprocal single-point mutational approach with introduction of His182 in PTP1B and Phe182 in PTPH1, we demonstrate here that His182-PTPs, in comparison with Phe182-PTPs, have significantly decreased kcat values, and to a lesser degree, decreased kcat/Km values. Combined enzyme kinetic, X-ray crystallographic and molecular dynamics studies indicate that the effect of His182 is due to interactions with Asp181 and with Gln262. We conclude that residue 182 can modulate the functionality of both Asp181 and Gln262 and therefore affect the E-P hydrolysis step of PTP-mediated catalysis.

Amino Acid Sequence↗

PubMed2004

Evolution of a rippled membrane during phospholipase A2 hydrolysis studied by time-resolved AFM.

The sensitivity of phospholipase A(2) (PLA(2)) for lipid membrane curvature is explored by monitoring, through time-resolved atomic force microscopy, the hydrolysis of supported double bilayers in the ripple phase. The ripple phase presents a corrugated morphology. PLA(2) is shown to have higher activity toward the ripple phase compared to the gel phase in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes, indicating its preference for this highly curved membrane morphology. Hydrolysis of the stable and metastable ripple structures is monitored for equimolar DMPC/1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC)-supported double bilayers. As shown by high-performance liquid chromatography results, DSPC is resistant to hydrolysis at this temperature, resulting in a more gradual hydrolysis of the surface that leads to a change in membrane morphology without loss of membrane integrity. This is reflected in an increase in ripple spacing, followed by a sudden flattening of the lipid membrane during hydrolysis. Hydrolysis of the ripple phase results in anisotropic holes running parallel to the ripples, suggesting that the ripple phase has strip regions of higher sensitivity to enzymatic attack. Bulk high-performance liquid chromatography measurements indicate that PLA(2) preferentially hydrolyzes DMPC in the DMPC/DSPC ripples. We suggest that this leads to the formation of a flat gel-phase lipid membrane due to enrichment in DSPC. The results point to the ability of PLA(2) for inducing a compositional phase transition in multicomponent membranes through preferential hydrolysis while preserving membrane integrity.

Anisotropy↗

PubMed2003

A novel strategy for the development of selective active-site inhibitors of the protein tyrosine phosphatase-like proteins islet-cell antigen 512 (IA-2) and phogrin (IA-2beta).

Islet-cell antigen 512 (IA-2) and phogrin (IA-2beta) are atypical members of the receptor protein tyrosine phosphatase (PTP) family that are characterized by a lack of activity against conventional PTP substrates. The physiological role(s) of these proteins remain poorly defined, although recent studies indicate that IA-2 may be involved in granule trafficking and exocytosis. To further understand their function, we have embarked upon developing low-molecular-mass inhibitors of IA-2 and IA-2beta. Previously, we have shown that a general PTP inhibitor, 2-(oxalylamino)benzoic acid (OBA), can be developed into highly selective and potent inhibitors of PTP1B. However, since wild-type IA-2 and IA-2beta lack conventional PTP activity, a novel strategy was designed whereby catalytically active species were generated by 'back-mutating' key non-consensus catalytic region residues to those of PTP1B. These mutants were then used as tools with which to test the potency and selectivity of OBA and a variety of its derivatives. Catalytically competent IA-2 and IA-2beta species were generated by 'back-mutation' of only three key residues (equivalent to Tyr(46), Asp(181) and Ala(217) using the human PTP1B numbering) to those of PTP1B. Importantly, enzyme kinetic analyses indicated that the overall fold of both mutant and wild-type IA-2 and IA-2beta was similar to that of classic PTPs. In particular, one derivative of OBA, namely 7-(1,1-dioxo-1 H -benzo[ d ]isothiazol-3-yloxymethyl)-2-(oxalylamino)-4,7-dihydro-5 H -thieno[2,3- c ]pyran-3 -carboxylic acid ('Compound 6 ' shown in the main paper), which inhibited IA-2beta((S762Y/Y898P/D933A)) (IA-2beta in which Ser(762) has been mutated to tyrosine, Tyr(898) to proline, and Asp(933) to alanine) with a K (i) value of approximately 8 microM, appeared ideal for future lead optimization. Thus molecular modelling of this classical, competitive inhibitor in the catalytic site of wild-type IA-2beta identified two residues (Ser(762) and Asp(933)) that offer the possibility for unique interaction with an appropriately modified 'Compound 6 '. Such a compound has the potential to be a highly selective and potent active-site inhibitor of wild-type IA-2beta.

Amino Acid Sequence↗

PubMed2003

Water in contact with extended hydrophobic surfaces: direct evidence of weak dewetting.

X-ray reflectivity measurements reveal a significant dewetting of a large hydrophobic paraffin surface floating on water. The dewetting phenomenon extends less than 15 A into the bulk water phase and results in an integrated density deficit of about one water molecule per 25-30 A(2) of water in contact with the paraffin surface. The results are supported by molecular dynamics simulations and related to the hydrophobic effect.

Hydrophobic and Hydrophilic Interactions↗

PubMed2003

Ligand-induced conformational changes: improved predictions of ligand binding conformations and affinities.

A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.

Binding Sites↗

PubMed2003

Enzyme kinetic characterization of protein tyrosine phosphatases.

Protein tyrosine phosphatases (PTPs) play a central role in cellular signaling processes, resulting in an increased interest in modulating the activities of PTPs. We therefore decided to undertake a detailed enzyme kinetic evaluation of various transmembrane and cytosolic PTPs (PTPalpha, PTPbeta, PTPepsilon, CD45, LAR, PTP1B and SHP-1), using pNPP as substrate. Most noticeable is the increase in the turnover number for PTPbeta with increasing pH and the weak pH-dependence of the turnover number of CD45. The kinetic data for PTPalpha-D1 and PTPalpha-D1D2 suggest that D2 affects the catalysis of pNPP. PTPepsilon and the closely homologous PTPalpha behave differently. The K(m) data were lower for PTPepsilon than those for PTPalpha, while the inverse was observed for the catalytic efficiencies.

Animals↗

PubMed2002

Orientation and conformation of a lipase at an interface studied by molecular dynamics simulations.

Electron density profiles calculated from molecular dynamics trajectories are used to deduce the orientation and conformation of Thermomyces lanuginosa lipase and a mutant adsorbed at an air-water interface. It is demonstrated that the profiles display distinct fine structures, which uniquely characterize enzyme orientation and conformation. The density profiles are, on the nanosecond timescale, determined by the average enzyme conformation. We outline a computational scheme that from a single molecular dynamics trajectory allows for extraction of electron density profiles referring to different orientations of the lipase relative to an implicit interface. Profiles calculated for the inactive and active conformations of the lipase are compared with experimental electron density profiles measured by x-ray reflectivity for the lipase adsorbed at an air-water interface. The experimental profiles contain less fine structural information than the calculated profiles because the resolution of the experiment is limited by the intrinsic surface roughness of water. Least squares fits of the calculated profiles to the experimental profiles provide areas per adsorbed enzyme and suggest that Thermomyces lanuginosa lipase adsorbs to the air-water interface in a semiopen conformation with the lid oriented away from the interface.

Ascomycota↗

PubMedDate not supplied

Essential motions in a fungal lipase with bound substrate, covalently attached inhibitor and product.

As an aid to understanding the influence of dynamic fluctuations during esterolytic catalysis, we follow protein flexibility at three different steps along the catalytic pathway from substrate binding to product clearance via a covalently attached inhibitor, which represents a transition-state mimic. We have applied a classical approach, using molecular dynamics simulations to monitor protein dynamics in the nanosecond regime. We filter out small amplitude fluctuations and focus on the anharmonic contributions to the overall dynamics. This 'essential dynamics' analysis reveals different modes of response along the pathway suggesting that binding, catalysis and product clearance occur along different energy surfaces. Motions in the enzyme with a covalently attached ligand are more complex and occur along several eigenvectors. The magnitudes of the fluctuations in these individual subspaces are significantly smaller than those observed for the substrate and product molecules, indicating that the energy surface is shallow and that a relatively large number of conformational substates are accessible. On the other hand, substrate binding and product release occur at distinct modes of the protein flexibility suggesting that these processes occur along rough energy surfaces with only a few minima. Detailed energetic analyses along the trajectories indicated that in all cases binding is dominated by van der Waals interactions. The carboxylate form of the product is stabilized by a tight hydrogen bond network involving in particular Ser82, which may be a potential cause of product inhibition. Considerations such as these should aid the understanding of mechanisms of substrate, inhibitor or product recognition and could become of importance in the design of new substrates or inhibitors for enzymes.

Enzyme Inhibitors↗