Joanna Makowska

Mechanism of formation of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. I. Importance of hydrophobic interactions in stabilization of beta-hairpin structure.

We previously studied a 16‐amino acid‐residue fragment of the C‐terminal β‐hairpin of the B3 domain (residues 46–61), [IG(46–61)] of the immunoglobulin binding protein G from Streptoccocus, and found that hydrophobic interactions and the turn region play an important role in stabilizing the structure. Based on these results, we carried out systematic structural studies of peptides derived from the sequence of IG (46–61) by systematically shortening the peptide by one residue at a time from both the C‐ and the N‐terminus. To determine the structure and stability of two resulting 12‐ and 14‐amino acid‐residue peptides, IG(48–59) and IG(47–60), respectively, we carried out circular dichroism, NMR, and calorimetric studies of these peptides in pure water. Our results show that IG(48–59) possesses organized three‐dimensional structure stabilized by hydrophobic interactions (Tyr50–Phe57 and Trp48–Val59) at T = 283 and 305 K. At T = 313 K, the structure breaks down because of increased chain entropy, but the turn region is preserved in the same position observed for the structure of the whole protein. The breakdown of structure occurs near the melting temperature of this peptide (Tm = 310 K) measured by differential scanning calorimetry (DSC). The melting temperature of IG(47–60) determined by DSC is Tm = 330 K and its structure is similar to that of the native β‐hairpin at all (lower) temperatures examined (283–313 K). Both of these truncated sequences are conserved in all known amino acid sequences of the B domains of the immunoglobulin binding protein G from bacteria. Thus, this study contributes to an understanding of the mechanism of folding of this whole family of proteins, and provides information about the mechanism of formation and stabilization of a β‐hairpin structural element. Proteins 2009.

Influence of charge and size of terminal amino‐acid residues on local conformational states and shape of alanine‐based peptides

We present results of conformational studies by Circular dichroism and NMR spectroscopy, differential scanning calorimetry, and molecular dynamics, of three alanine‐based peptides: Ac‐KK‐(A)7‐KK‐NH2 (KAK), Ac‐OO‐(A)7‐DD‐NH2 (OAD), and Ac‐KK‐(A)7‐EE‐NH2 (KAE), where A, K, O, D, and E, denote alanine, lysine, ornithine, aspartic acid, and glutamic acid residues, respectively. For OAD and KAE, canonical MD simulations with time‐averaged NMR‐derived restraints demonstrate the presence of an ensemble of structures with a variety of conformational states (polyproline II, α‐helical, α′, and extended, turn); for KAK the conformational states are predominantly polyproline II and extended. The OAD peptide exhibits a bent shape with its ends close to each other, whereas KAK and KAE are more extended. The bent shape was also observed in our earlier study of the Ac‐XX‐(A)7‐OO‐NH2 (XAO) peptide, where X denotes the diaminobutyric acid residue; therefore, the shape seems to depend on the size of the charged side chains at the ends of the alanine sequence and not on their kind. This suggests that the bent shape of the alanine sequence is formed to enable screening of this nonpolar sequence from the solvent by sufficiently short charged side chains. As in our previous study of the XAO peptide, no long polyproline II segments were observed.

Acidic‐basic properties of three alanine‐based peptides containing acidic and basic side chains: Comparison between theory and experiment

The purpose of this work was to evaluate the effect of the nature of the ionizable end groups, and the solvent, on their acid‐base properties in alanine‐based peptides. Hence, the acid‐base properties of three alanine‐based peptides: Ac‐KK‐(A)7‐KK‐NH2 (KAK), Ac‐OO‐(A)7‐DD‐NH2 (OAD), Ac‐KK‐(A)7‐EE‐NH2 (KAE), where A, D, E, K, and O denote alanine, aspartic acid, glutamic acid, lysine, and ornithine, respectively, were determined in water and in methanol by potentiometry. With the availability of these data, the ability of two theoretical methods to simulate pH‐metric titration of those peptides was assessed: (i) the electrostatically driven Monte Carlo method with the ECEPP/3 force field and the Poisson‐Boltzmann approach to compute solvation energy (EDMC/PB/pH), and (ii) the molecular dynamics method with the AMBER force field and the Generalized Born model (MD/GB/pH). For OAD and KAE, pKa1 and pKa2 correspond to the acidic side chains. For all three compounds in both solvents, the pKa1 value is remarkably lower than the pKa of a compound modeling the respective isolated side chain, which can be explained by the influence of the electrostatic field from positively charged ornithine or lysine side chains. The experimental titration curves are reproduced well by the MD/GB/pH approach, the agreement being better if restraints derived from NMR measurements are incorporated in the conformational search. Poorer agreement is achieved by the EDMC/PB/pH method.

Preliminary studies on trigonelline as potential anti-Alzheimer disease agent: Determination by hydrophilic interaction liquid chromatography and modeling of interactions with beta-amyloid☆

For trigonelline, a quaternary-base pyridine alkaloid of presumed Alzheimers disease-preventing activity, a method of determination has been proposed, based on the hydrophilic interaction chromatography (HILIC). That method might be applied to study the agents bioavailability, in particular its permeation through blood-brain barrier, which is an inevitable property for the potential central nervous system affecting drugs. Providing that trigonelline possesses the requested pharmacokinetic properties, once attaining pharmacodynamic phase it must interact effectively with the relevant molecular site in the brain, which is characteristic to neurodegenerative diseases, namely the beta-amyloid peptide. Here it was demonstrated by molecular modeling that affinity of trigonelline to the Aβ(1-42) peptide is high and similar to that of an anti-Alzheimers disease drug candidate - cotinine.

Binding of Cu(II) ions to peptides studied by fluorescence spectroscopy and isothermal titration calorimetry

Steady-state and time-resolved fluorescence quenching measurements supported by Isothermal Titration Calorimetry (ITC) were used to study the interactions of Cu(2+) with four peptides. Two of them were taken from the N-terminal part of the FBP28 protein (formin binding protein) WW domain: Tyr-Lys-Thr-Ala-Asp-Gly-Lys-Thr-Tyr-NH2 (D9) and its mutant Tyr-Lys-Thr-Ala-Asn-Gly-Lys-Thr-Tyr-NH2 (D9_M) as well as two mutated peptides from the B3 domain of the immunoglobulin binding protein G derived from Streptococcus: Asp-Val-Ala-Thr-Tyr-Thr-NH2 (J1) and Glu-Val-Ala-Thr-Tyr-Thr-NH2 (J2). The measurements were carried out at 298.15K in 20mM 2-(N-morpholino)ethanesulfonic acid (MES) buffer solution with a pH of 6. The fluorescence of all peptides was quenched by Cu(2+) ions. The stoichiometry, conditional stability constants and thermodynamic parameters for the interactions of the Cu(2+) ions with D9 and D9_M were determined from the calorimetric data. The values of the conditional stability constants were additionally determined from fluorescence quenching measurements and compared with those obtained from calorimetric studies. There was a good correlation between data obtained from the two techniques. On the other hand, the studies revealed that J1 and J2 do not exhibit an affinity towards metal ions. The obtained results prove that fluorescence quenching experiments may be successfully used in order to determine stability constants of complexes with fluorescent ligands. Finally, based on the obtained results, the coordinating properties of the peptides towards the Cu(2+) ions are discussed.

A Study of the Influence of Charged Residues on β-Hairpin Formation by Nuclear Magnetic Resonance and Molecular Dynamics

Chain reversals are often nucleation sites in protein folding. The β-hairpins of FBP28 WW domain and IgG are stable and have been proved to initiate the folding and are, therefore, suitable for studying the influence of charged residues on β-hairpin conformation. In this paper, we carried out NMR examination of the conformations in solution of two fragments from the FPB28 protein (PDB code: 1E0L) (N-terminal part) namely KTADGKT-NH2 (1E0L 12–18, D7) and YKTADGKTY-NH2 (1E0L 11–19, D9), one from the B3 domain of the protein G (PDB code: 1IGD), namely DDATKT-NH2 (1IGD 51–56) (Dag1), and three variants of Dag1 peptide: DVATKT-NH2 (Dag2), OVATKT-NH2 (Dag3) and KVATKT-NH2 (Dag4), respectively, in which the original charged residue were replaced with non-polar residues or modified charged residues. It was found that both the D7 and D9 peptides form a large fraction bent conformations. However, no hydrophobic contacts between the terminal Tyr residues of D9 occur, which suggests that the presence of a pair of like-charged residues stabilizes chain reversal. Conversely, only the Dag1 and Dag2 peptides exhibit some chain reversal; replacing the second aspartic-acid residue with a valine and the first one with a basic residue results in a nearly extended conformation. These results suggest that basic residues farther away in sequence can result in stabilization of chain reversal owing to screening of the non-polar core. Conversely, smaller distance in sequence prohibits this screening, while the presence oppositely-charged residues can stabilize a turn because of salt-bridge formation.

Like-charged residues at the ends of oligoalanine sequences might induce a chain reversal

We have examined the effect of like-charged residues on the conformation of an oligoalanine sequence. This was facilitated by circular dichroism (CD) and NMR spectroscopic and differential scanning calorimetric (DSC) measurements, and molecular dynamics calculations of the following three alanine-based peptides: Ac-K-(A)(5) -K-NH(2) (KAK5), Ac-K-(A)(4) -K-NH(2) (KAK4), Ac-K-(A)(3) -K-NH(2) (KAK3), where A and K denote alanine and lysine residues, respectively. Our earlier studies suggested that the presence of like-charged residues at the end of a short polypeptide chain composed of nonpolar residues can induce a chain reversal. For all three peptides, canonical molecular dynamics simulations with NMR-derived restraints demonstrate the presence of ensembles of structures with a tendency to form a chain reversal. The KAK3 peptide exhibits a bent shape with its ends close to each other, while KAK4 and KAK5 are more extended. In the KAK5 peptide, the lysine residues do not have any influence on each other and are very mobile. Nevertheless, the tendency to form a more or less pronounced chain reversal is observed and it seems to be stable in all three peptides. This chain reversal seems to be caused by screening of the nonpolar core from the solvent by the hydrated charged residues.

Biochemical Characterization and Validation of a Catalytic Site of a Highly Thermostable Ts2631 Endolysin from the Thermus scotoductus Phage vB_Tsc2631.

Phage vB_Tsc2631 infects the extremophilic bacterium Thermus scotoductus MAT2631 and uses the Ts2631 endolysin for the release of its progeny. The Ts2631 endolysin is the first endolysin from thermophilic bacteriophage with an experimentally validated catalytic site. In silico analysis and computational modelling of the Ts2631 endolysin structure revealed a conserved Zn2+ binding site (His30, Tyr58, His131 and Cys139) similar to Zn2+ binding site of eukaryotic peptidoglycan recognition proteins (PGRPs). We have shown that the Ts2631 endolysin lytic activity is dependent on divalent metal ions (Zn2+ and Ca2+). The Ts2631 endolysin substitution variants H30N, Y58F, H131N and C139S dramatically lost their antimicrobial activity, providing evidence for the role of the aforementioned residues in the lytic activity of the enzyme. The enzyme has proven to be not only thermoresistant, retaining 64.8% of its initial activity after 2 h at 95°C, but also highly thermodynamically stable (Tm = 99.82°C, ΔHcal = 4.58 × 104 cal mol-1). Substitutions of histidine residues (H30N and H131N) and a cysteine residue (C139S) resulted in variants aggregating at temperatures ≥75°C, indicating a significant role of these residues in enzyme thermostability. The substrate spectrum of the Ts2631 endolysin included extremophiles of the genus Thermus but also Gram-negative mesophiles, such as Escherichia coli, Salmonella panama, Pseudomonas fluorescens and Serratia marcescens. The broad substrate spectrum and high thermostability of this endolysin makes it a good candidate for use as an antimicrobial agent to combat Gram-negative pathogens.

Theoretical calculations of homoconjugation equilibrium constants in systems modeling acid-base interactions in side chains of biomolecules using the potential of mean force.

The potentials of mean force (PMFs) were determined for systems forming cationic and anionic homocomplexes composed of acetic acid, phenol, isopropylamine, n‐butylamine, imidazole, and 4(5)‐methylimidazole, and their conjugated bases or acids, respectively, in three solvents with different polarity and hydrogen‐bonding propensity: acetonitrile (AN), dimethyl sulfoxide (DMSO), and water (H2O). For each pair and each solvent a series of umbrella‐sampling molecular dynamics simulations with the AMBER force field, explicit solvent, and counterions added to maintain a zero net charge of a system were carried out and the PMF was calculated by using the Weighted Histogram Analysis Method (WHAM). Subsequently, homoconjugation‐equilibrium constants were calculated by numerical integration of the respective PMF profiles. In all cases but imidazole stable homocomplexes were found to form in solution, which was manifested as the presence of contact minima corresponding to hydrogen‐bonded species in the PMF curves. The calculated homoconjugation constants were found to be greater for complexes with the OHO bridge (acetic acid and phenol) than with the NHN bridge and they were found to decrease with increasing polarity and hydrogen‐bonding propensity of the solvent (i.e., in the series AN > DMSO > H2O), both facts being in agreement with the available experimental data. It was also found that interactions with counterions are manifested as the broadening of the contact minimum or appearance of additional minima in the PMF profiles of the acetic acid‐acetate, phenol/phenolate system in acetonitrile, and the 4(5)‐methylimidazole/4(5)‐methylimidzole cation conjugated base system in dimethyl sulfoxide.

Interplay of charge distribution and conformation in peptides: Comparison of theory and experiment†

We assessed the correlation between charge distribution and conformation of flexible peptides by comparing the theoretically calculated potentiometric‐titration curves of two model peptides, Ac–Lys5–NHMe (a model of poly‐L‐lysine) and Ac–Lys–Ala11–Lys–Gly2–Tyr–NH2 (P1) in water and methanol, with the experimental curves. The calculation procedure consisted of three steps: (i) global conformational search of the peptide under study using the electrostatically driven Monte Carlo (EDMC) method with the empirical conformational energy program for peptides (ECEPP)/3 force field plus the surface‐hydration (SRFOPT) or the generalized Born surface area (GBSA) solvation model as well as a molecular dynamics method with the assisted model building and energy refinement (AMBER)99/GBSA force field; (ii) reevaluation of the energy in the pH range considered by using the modified Poisson–Boltzmann approach and taking into account all possible protonation microstates of each conformation, and (iii) calculation of the average degree of protonation of the peptide at a given pH value by Boltzmann averaging over conformations. For Ac–Lys5–NHMe, the computed titration curve agrees qualitatively with the experimental curve of poly‐L‐lysine in 95% methanol. The experimental titration curves of peptide P1 in water and methanol indicate a remarkable downshift of the first pKa value compared to the values for reference compounds (n‐butylamine and phenol, respectively), suggesting the presence of a hydrogen bond between the tyrosine hydroxyl oxygen and the Hϵ proton of a protonated lysine side chain. The theoretical titration curves agree well with the experimental curves, if conformations with such hydrogen bonds constitute a significant part of the ensemble; otherwise, the theory predicts too small a downward pH shift.