Biomedical subjects

Ernst-Walter Knapp

Publications and source records attributed to Ernst-Walter Knapp.

31 records · Page 2Linked to original sources

PubMed2005

Tuning electron transfer by ester-group of chlorophylls in bacterial photosynthetic reaction center.

Accessory chlorophylls (B(A/B)) in bacterial photosynthetic reaction center play a key role in charge-separation. Although light-exposed and dark-adapted bRC crystal structures are virtually identical, the calculated B(A) redox potentials for one-electron reduction differ. This can be traced back to different orientations of the B(A) ester-group. This tuning ability of chlorophyll redox potentials modulates the electron transfer from SP* to B(A).

Chlorophyll↗

PubMed2005

Energetics of proton transfer pathways in reaction centers from Rhodobacter sphaeroides. The Glu-H173 activated mutants.

Electron transfer between the primary and secondary quinones (Q(A), Q(B)) in the bacterial photosynthetic reaction center (bRC) is coupled with proton uptake at Q(B). The protons are conducted from the cytoplasmic side, probably with the participation of two water channels. Mutations of titratable residues like Asp-L213 to Asn (inhibited mutant) or the double mutant Glu-L212 to Ala/Asp-L213 to Ala inhibit these electron transfer-coupled proton uptake events. The inhibition of the proton transfer (PT) process in the single mutant can be restored by a second mutation of Arg-M233 to Cys or Arg-H177 to His (revertant mutant). These revertant mutants shed light on the location of the main proton transfer pathway of wild type bRC. In contrast to the wild type and inhibited mutant bRC, the revertant mutant bRC showed notable proton uptake at Glu-H173 upon formation of the Q(B)- state. In all of these mutants, the pK(a) of Asp-M17 decreased by 1.4-2.4 units with respect to the wild type bRC, whereas a significant pK(a) upshift of up to 5.8 units was observed at Glu-H122, Asp-H170, Glu-H173, and Glu-H230 in the revertant mutants. These residues belonging to the main PT pathway are arranged along water channel P1 localized mainly in subunit H. bRC possesses subunit H, which has no counterpart in photosystem II. Thus, bRC may possess alternative PT pathways involving water channels in subunit H, which becomes active in case the main PT pathway is blocked.

Alanine↗

PubMed2004

A molecular docking study of estrogenically active compounds with 1,2-diarylethane and 1,2-diarylethene pharmacophores.

Numerous selective estrogen receptor modulators (SERMs) have been synthesized and assayed in recent years. The focus of this study is to apply coarse-grain molecular docking procedures coupled with fine-grain all-atom force field optimization strategies to shed light on the binding mechanisms of currently available estrogen receptor-active compounds. Although the mechanics of ligand binding in estrogen receptors is generally well understood, there is room for surprises. In this paper computational evidence corroborating the experimentally observed type I agonistic binding mode for estradiol (E2) and diethylstilbesterol (DES) and the type II antagonistic binding mode for 4-hydroxytamoxifen and raloxifen is presented. Included in this type I agonistic mode are the DES derivatives, transstilbene and 1,2-diaryldiaminoethane. In addition, a novel 'type II agonistic' binding mode for 2,3-diarylimidazolines, 4,5-diarylimidazoles, 2,3-diarylpiperazines is introduced. This mode is stabilized by suggesting alternative hydrogen bond anchor points in the ligand binding domain as potential leads for future drug design.

Computer Simulation↗

PubMed2004

Accurate pKa determination for a heterogeneous group of organic molecules.

Single-molecule studies that allow to compute pKa values, proton affinities (gas-phase acidity/basicity) and the electrostatic energy of solvation have been performed for a heterogeneous set of 26 organic compounds. Quantum mechanical density functional theory (DFT) using the Becke-half&half and B3LYP functionals on optimized molecular geometries have been carried out to investigate the energetics of gas-phase protonation. The electrostatic contribution to the solvation energies of protonated and deprotonated compounds were calculated by solving the Poisson equation using atomic charges generated by fitting the electrostatic potential derived from the molecular wave functions in vacuum. The combination of gas-phase and electrostatic solvation energies by means of the thermodynamic cycle enabled us to compute pKa values for the 26 compounds, which cover six distinct chemical groups (carboxylic acids, benzoic acids, phenols, imides, pyridines and imidazoles). The computational procedure for determining pKa values is accurate and transferable with a root-mean-square deviation of 0.53 and 0.57 pKa units and a maximum error of 1.0 pKa and 1.3 pKa units for Becke-half&half and B3LYP DFT functionals, respectively.

Journal Article↗

PubMed2004

Variation of Ser-L223 hydrogen bonding with the QB redox state in reaction centers from Rhodobacter sphaeroides.

Ser-L223 is close to ubiquinone (Q(B)) in the B-branch of the bacterial photosynthetic reaction center (bRC) from Rhodobacter (Rb) sphaeroides. Therefore, the presence of a hydrogen bond (H bond) between the two was naturally proposed from the crystal structure. The hydrogen bonding pattern of Q(B) from the light-exposed structure was studied by generating hydrogen atom coordinates based on the CHARMM force field. In the Q(B) neutral charge state (Q(B)(0)), no H bond was found between the oxygen of the OH group from Ser-L223 and the carbonyl oxygen of Q(B) that is distal to the non-heme iron. In the reduced state (Q(B)(-)), however, Ser-L213 was found to form an H bond with Q(B) only when Asp-L213 is protonated by more than 0.75 H(+). This indicates the significance of the protonation of Asp-L213 in forming an H bond between Ser-L223 and Q(B). We found that the driving force to form the H bond between Ser-L223 and Q(B) is enhanced by the positively charged Arg-L217. The calculated Q(B) redox potentials with or without this H bond discriminated two ET rates, which are close to the faster and slower time phases observed in UV-Vis and FTIR studies. Together with the calculated redox potential of the quinones, this H-bond formation could play a key role in conformational gating for the ET process from Q(A) to Q(B).

Aspartic Acid↗

PubMed2004

Peptide binding at class I major histocompatibility complex scored with linear functions and support vector machines.

We explore two different methods to predict the binding ability of nonapeptides at the class I major histocompatibility complex using a general linear scoring function that defines a separating hyperplane in the feature space of sequences. In absence of suitable data on non-binding nonapeptides we generated sequences randomly from a selected set of proteins from the protein data bank. The parameters of the scoring function were determined by a generalized least square optimization (LSM) and alternatively by the support vector machine (SVM). With the generalized LSM impaired data for learning with a small set of binding peptides and a large set of non-binding peptides can be treated in a balanced way rendering LSM more successful than SVM, while for symmetric data sets SVM has a slight advantage compared to LSM.

Amino Acid Sequence↗

PubMed2003

Unfolding of the cold shock protein studied with biased molecular dynamics.

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.

Bacterial Proteins↗

PubMed2003

Tuning heme redox potentials in the cytochrome C subunit of photosynthetic reaction centers.

The photosynthetic reaction center (RC) from Rhodopseudomonas viridis contains four cytochrome c hemes. They establish the initial part of the electron transfer (ET) chain through the RC. Despite their chemical identity, their midpoint potentials cover an interval of 440 mV. The individual heme midpoint potentials determine the ET kinetics and are therefore tuned by specific interactions with the protein environment. Here, we use an electrostatic approach based on the solution of the linearized Poisson-Boltzmann equation to evaluate the determinants of individual heme redox potentials. Our calculated redox potentials agree within 25 meV with the experimentally measured values. The heme redox potentials are mainly governed by solvent accessibility of the hemes and propionic acids, by neutralization of the negative charges at the propionates through either protonation or formation of salt bridges, by interactions with other hemes, and to a lesser extent, with other titratable protein side chains. In contrast to earlier computations on this system, we used quantum chemically derived atomic charges, considered an equilibrium-distributed protonation pattern, and accounted for interdependencies of site-site interactions. We provide values for the working potentials of all hemes as a function of the solution redox potential, which are crucial for calculations of ET rates. We identify residues whose site-directed mutation might significantly influence ET processes in the cytochrome c part of the RC. Redox potentials measured on a previously generated mutant could be reproduced by calculations based on a model structure of the mutant generated from the wild type RC.

Binding Sites↗

PubMed2003

Redox potential of quinones in both electron transfer branches of photosystem I.

The redox potentials of the two electron transfer (ET) active quinones in the central part of photosystem I (PSI) were determined by evaluating the electrostatic energies from the solution of the Poisson-Boltzmann equation based on the crystal structure. The calculated redox potentials are -531 mV for A1A and -686 mV for A1B. From these results we conclude the following. (i) Both branches are active with a much faster ET in the B-branch than in the A-branch. (ii) The measured lifetime of 200-290 ns of reduced quinones agrees with the estimate for the A-branch and corroborates with an uphill ET from this quinone to the iron-sulfur cluster as observed in recent kinetic measurements. (iii) The electron paramagnetic resonance spectroscopic data refer to the A-branch quinone where the corresponding ET is uphill in energy. The negative redox potential of A1 in PSI is primarily because of the influence from the negatively charged FX, in contrast to the positive shift on the quinone redox potential in bacterial reaction center and PSII that is attributed to the positively charged non-heme iron atom. The conserved residue Asp-B575 changes its protonation state after quinone reduction. The difference of 155 mV in the quinone redox potentials of the two branches were attributed to the conformation of the backbone with a large contribution from Ser-A692 and Ser-B672 and to the side chain of Asp-B575, whose protonation state couples differently with the formation of the quinone radicals.

Benzoquinones↗

PubMed2003

Redox potential of quinones in photosynthetic reaction centers from Rhodobacter sphaeroides: dependence on protonation of Glu-L212 and Asp-L213.

The absolute values of the one-electron redox potentials of the two quinones (Q(A) and Q(B)) in bacterial photosynthetic reaction centers from Rhodobacter sphaeroides were calculated by evaluating the electrostatic energies from the solution of the linearized Poisson-Boltzmann equation at pH 7.0. The redox potential for Q(A) was calculated to be between -173 and -160 mV, which is close to the lowest measured values that are assumed to refer to nonequilibrated protonation patterns in the redox state Q(A)(-). The redox potential of quinone Q(B) is found to be about 160-220 mV larger for the light-exposed than for the dark-adapted structure. These values of the redox potentials are obtained if Asp-L213 is nearly protonated (probability 0.75-1.0) before and after electron transfer from Q(A) to Q(B), while Glu-L212 is partially protonated (probability 0.6) in the initial state Q(A)(-)Q(B)(0) and fully protonated in the final state Q(A)(0)Q(B)(-). Conversely, if the charge state of the quinones is varied from Q(A)(-)Q(B)(0) to Q(A)(0)Q(B)(-) corresponding to the electron transfer from Q(A) to Q(B), Asp-L213 remains protonated, while Glu-L212 changes its protonation state from 0.15 H(+) to fully protonated. In agreement with results from FTIR spectra, there is proton uptake at Glu-L212 going along with the electron transfer, whereas Asp-L213 does not change its protonation state. However, in our simulations Asp-L213 is considered to be protonated rather than ionized as deduced from FTIR spectra. The calculated redox potential of Q(A) shows little dependence on the charge state of Asp-L213, which is due to a strong coupling with the protonation state of Asp-M17 but increases by 50 mV if Glu-L212 changes from the ionized to the protonated charge state. Both are in agreement with fluorescence measurements observing the decay of SP(+)Q(A)(-) in a wide pH regime. The computed difference in redox potential of Q(B) in the light-exposed and dark-adapted structure was traced back to the hydrogen bond of Q(B) with His-L190 that is lost in the dark-adapted structure and the charge of the non-heme iron atom, which is closer to Q(B) in the light-exposed than in the dark-adapted structure.

Amino Acid Substitution↗

PubMed2002

Energetics of radical transfer in DNA photolyase.

Charge separation and radical transfer in DNA photolyase from Escherichia coli is investigated by computing electrostatic free energies from a solution of the Poisson-Boltzmann equation. For the initial charge separation 450 meV are available. According to recent experiments [Aubert et al. Nature 2000, 405, 586-590] the flavin receives an electron from the proximal tryptophan W382, which consequently forms a cationic radical WH(*)(+)382. The radical state is subsequently transferred along the triad W382-W359-W306 of conserved tryptophans. The radical transfer to the intermediate tryptophan W359 is nearly isoenergetic (58 meV uphill); the radical transfer from the intermediate W359 to the distal W306 is 200 meV downhill in energy, funneling and stabilizing the radical state at W306. The resulting cationic radical WH(*)(+)306 is further stabilized by deprotonation, yielding the neutral radical W(*)306, which is 214 meV below WH(*)(+)306. The time scale of the charge recombination process yielding back the resting enzyme with FADH(*) is governed by reprotonation of W306, with a calculated lifetime of 1.2 ms that correlates well with the measured lifetime of 17 ms. In photolyase from Anacystis nidulans the radical state is partially transferred to a tyrosine [Aubert et al. Proc. Natl. Acad. Sci. U.S.A. 1999, 96, 5423-5427]. In photolyase from Escherichia coli, there is a tyrosine (Y464) close to the distal tryptophan W306 that could play this role. We show that this tyrosine cannot be involved in radical transfer, because the electron transfer from tyrosine to W306 is much too endergonic (750 meV) and a direct hydrogen transfer is likely too slow. Coupling of specific charge states of the tryptophan triad with protonation patterns of titratable residues of photolyase is small.

Deoxyribodipyrimidine Photo-Lyase↗

PubMed2002

Driving forces of protein association: the dimer-octamer equilibrium in arylsulfatase A.

The enzyme arylsulfatase A (ASA) occurs in solution as dimer (alpha(2)) above pH 6 and associates to octamers (alpha(2))(4) below pH 6. The crystal structure of ASA suggests that the (alpha(2))-(alpha(2))(4) equilibrium is regulated by protonation/deprotonation of Glu-424 located at the interface between (alpha(2)) dimers in the octamer. The reason for this assumption is that Glu-424 can be in two different conformers where it forms an intra or intermolecular hydrogen bond, respectively. In the present study we investigate this protein association process theoretically. The electrostatic energies are evaluated by solving the Poisson-Boltzmann equation for the inhomogeneous dielectric of the protein-water system for the dimer and octamer configurations. If a conventional surface energy term is used for the nonelectrostatic interactions, the absolute value of free energy of association fails to agree with experiment. A more detailed treatment that explicitly accounts for hydrophilic and hydrophobic character of the amino acids in the dimer-dimer interface of the octamer can explain this discrepancy qualitatively. The pH dependence of the computed association energy clearly demonstrates that the octamer is more stable at low pH if Glu-424 becomes protonated and forms an intermolecular hydrogen bond. We found a slight preference of Glu-424 to be in a conformation where its acidic group is fully solvent-exposed in the dimer state to form hydrogen bonds with water molecules. Application of the proton linkage model to calculate the association energy from the simulated data yielded results identical to the one obtained from the corresponding direct method.

Binding Sites↗

PubMed2002

Protonation and stability of the globular domain of influenza virus hemagglutinin.

A partial dissociation of the HA1 subunits of influenza virus hemagglutinin (HA) is considered to be the initial step of conformational changes of the HA ectodomain leading to a membrane fusion active conformation (L. Godley, J. Pfeifer, D. Steinhauer, B. Ely, G. Shaw, R. Kaufman, E. Suchanek, C. Pabo, J.J. Skehel, D.C. Wiley, and S. Wharton, 1992, Cell 68:635-645; G.W. Kemble, D.L.Bodian, J. Rose, I.A. Wilson, and J.M. White, 1992, J. Virol. 66:4940-4950). Here, we explore a mechanism that provides an understanding of the physical and chemical basis for such dissociation and relies on two essential observations. First, based on the x-ray structure of HA from X31 (I.A. Wilson, J.J. Skehel, and D.C. Wiley, 1981, Nature 289:366-373), and by employing techniques of molecular modeling, we show that the protonation of the HA1 subunits is enhanced at the conditions known to trigger conformational changes of the HA ectodomain. Second, we found that the dependence of the calculated relative degree of protonation of the HA1 domain on temperature and pH is similar to that observed experimentally for the conformational change of HA assessed by proteinase K sensitivity. We suggest that at the pH-temperature conditions typical for the conformational change of HA and membrane fusion, dissociation of the HA1 subunits is caused by the enhanced protonation of the HA1 subunits leading to an increase in the positive net charge of these subunits and, in turn, to a weakened attraction between them.

Amino Acid Motifs↗