Peptide investigations of pairwise interactions in the collagen triple-helix

Pairwise interactions have been studied for the major secondary structures in proteins. The present work extends the characterization of interactions between side-chains to the context of a collagen triple-helix. In this study, the most frequent Gly-X-Y tripeptide sequences in collagen are characterized in terms of interchain interactions between non-imino acid X and Y residues, through the use of host-guest peptides and statistical frequency analysis. Stabilities predicted on the basis of additivity show good agreement with experimental values for almost half of the peptides, indicating a lack of interaction. A small number of peptides have a stability lower than predicted, while a larger number are more stable than expected. Of all triplets containing residues of opposite charge, only Gly-Lys-Asp and Gly-Arg-Asp exhibit stabilizing electrostatic interactions, and these pairs are found together preferentially in collagens. Repulsion of like charges is observed in Gly-Arg-Lys, Gly-Lys-Arg, and Gly-Glu-Asp sequences, and a small degree of hydrophobic stabilization was observed for the Gly-Leu-Leu guest triplet. The data reported here help clarify basic principles of triple-helix stability. In addition, the experimentally determined stabilities of the tripeptide units found most frequently in collagens constitute a database useful for predicting triple-helix stability in peptides, collagens and other triple-helix-containing proteins.     

Sequence dependence of beta-hairpin structure: comparison of a salt bridge and an aromatic interaction

A comparison of the contributions and position dependence of cross-strand electrostatic and aromatic side-chain interactions to beta-sheet stability has been performed by using nuclear magnetic resonance in a well-folded beta-hairpin peptide of the general sequence XRTVXVdPGOXITQX. Phe-Phe and Glu-Lys pairs were varied at the internal and terminal non-hydrogen-bonded position, and the resulting stability was measured by the effects on alpha-hydrogen and aromatic hydrogen chemical shifts. It was determined that the introduction of a Phe-Phe pair resulted in a more folded peptide, regardless of position, and a more tightly folded core. Substitution of the Glu-Lys pair at the internal position results in a less folded peptide and increased fraying at the terminal residues. Upfield shifting of the aromatic hydrogens provided evidence for an edge-face aromatic interaction, regardless of position of the Phe-Phe pair. In peptides with two Phe-Phe pairs, substitution with Glu-Lys at either position resulted in a weakening of the aromatic interaction and a subsequent decrease in peptide stability. Thermal denaturation of the peptides containing Phe-Phe indicates that the aromatic interaction is enthalpically favored, whereas the folding of hairpins with cross-strand Glu-Lys pairs was less enthalpically favorable but entropically more favorable.     

Electrostatic interactions in the association of proteins: an analysis of the thrombin-hirudin complex.

The role of electrostatic interactions in stabilization of the thrombin-hirudin complex has been investigated by means of two macroscopic approaches: the modified Tanford-Kirkwood model and the finite-difference method for numerical solution of the Poisson-Boltzmann equations. The electrostatic potentials around the thrombin and hirudin molecules were asymmetric and complementary, and it is suggested that these fields influence the initial orientation in the process of the complex formation. The change of the electrostatic binding energy due to mutation of acidic residues in hirudin has been calculated and compared with experimentally determined changes in binding energy. In general, the change in electrostatic binding energy for a particular mutation calculated by the modified Tanford-Kirkwood approach agreed well with the experimentally observed change. The finite-difference approach tended to overestimate changes in binding energy when the mutated residues were involved in short-range electrostatic interactions. Decreases in binding energy caused by mutations of amino acids that do not make any direct ionic interactions (e.g., Glu 61 and Glu 62 of hirudin) can be explained in terms of the interaction of these charges with the positive electrostatic potential of thrombin. Differences between the calculated and observed changes in binding energy are discussed in terms of the crystal structure of the thrombin-hirudin complex.     

Identification and characterization of a bactericidal and proapoptotic peptide from Cycas revoluta seeds with DNA binding properties

Nowadays, novel pharmacies have been screened from plants. Among them are the peptides, which show multiple biotechnological activities. In this report, a small peptide (Ala–Trp–Lys–Leu–Phe–Asp–Asp–Gly–Val) with a molecular mass of 1,050 Da was purified from Cycas revoluta seeds by using reversed‐phase liquid chromatography. This peptide shows clear deleterious effects against human epidermoid cancer (Hep2) and colon carcinoma cells (HCT15). It caused inhibition of cancer cell proliferation and further disruption of nucleosome structures, inducing apoptosis by direct DNA binding. A remarkable antibacterial activity was also observed in this same peptide. Nevertheless, no significant lysis of normal RBC cells was observed in the presence of peptide. Additionally, an acetylation at the N‐termini portion is able to reduce both activities. Bioinformatics tools were also utilized for construction of a three‐dimensional model showing a single amphipathic helix. Since in vitro binding studies show that the target of this peptide seems to be DNA, theoretical docking studies were also performed to better understand the interaction between peptide and nucleic acids and also to shed some light on the acetyl group role. Firstly, binding studies showed that affinity contacts basically occur due to electrostatic attraction. The complex peptide‐ssDNA was clearly oriented by residues Ala¹, Lys³, and Asp⁶, which form several hydrogen bonds that are able to stabilize the complex. When acetyl was added, hydrogen bonds are broken, reducing the peptide affinity. In summary, it seems that information here provided could be used to design a novel derivative of this peptide which a clear therapeutic potential. J. Cell. Biochem. 113: 184–193, 2012. © 2011 Wiley Periodicals, Inc.     

Solvent effects on the secondary structure of the membrane-active zervamicin determined by PELDOR spectroscopy

Zervamicin is a voltage-gated ion-channel-forming peptide. Channels are generally considered to be formed by first insertion of amphipathic molecules into the phospholipid bilayer, followed by self-assembly of a variable number of transmembrane helices. We have studied the length of the peptide structure to address the question whether this peptide is long enough to span the phospholipid bilayer. The pulsed electron-electron double resonance (PELDOR) spectroscopic technique was used to determine the length of the helical molecule in membrane-mimicking solvents. This was achieved from the distance-related dipole-dipole interaction between spin labels, which were located at both ends of the linear peptide chain. The data were obtained by using samples of frozen glassy solutions of MeOH, MeOH/toluene, and MeOH/CHCl(3). Contributions of inter- and intramolecular interactions of spin labels were separated to analyze the intramolecular interaction and the distance distribution function between the labels. It is shown that the main maximum of the distribution functions is located at a distance of ca. 3.3 nm, and this distance appears to be only slightly dependent on the solvent composition. The distribution function was observed to narrow after addition of either CHCl(3) or toluene to MeOH. This effect is rationalized in terms of a decreased mobility of the terminal amino acid residues. By molecular-dynamics simulations, it was shown that the conformation, corresponding with the predominant distance found by PELDOR, agrees well with the mixed alpha/3(10)-helical that was previously determined by NMR. However, in the case toluene was added to the MeOH solution to further increase the hydrophobicity of the environment of the membrane-active peptide, the distribution function gives rise to a minor fraction (7-8%) with a distance of 4.2 nm. This distance corresponds most likely to the more extended 2(7)-helix structure.     

Gain of local structure in an amphipathic peptide does not require a specific tertiary framework

In this work, we studied how an amphipathic peptide of the surface of the globular protein thioredoxin, TRX94‐108, acquires a native‐like structure when it becomes involved in an apolar interaction network. We designed peptide variants where the tendency to form α‐helical conformation is modulated by replacing each of the leucine amino acid residues by an alanine. The induction of structure caused by sodium dodecyl sulfate (SDS) binding was studied by capillary zone electrophoresis, circular dichroism, DOSY‐NMR, and molecular dynamics simulations (MDS). In addition, we analyzed the strength of the interaction between a C18 RP‐HPLC matrix and the peptides. The results presented here reveal that (a) critical elements in the sequence of the wild‐type peptide stabilize a SDS/peptide supramolecular cluster; (b) the hydrophobic nature of the interaction between SDS molecules and the peptide constrains the ensemble of conformations; (c) nonspecific apolar surfaces are sufficient to stabilize peptide secondary structure. Remarkably, MDS shed light on a contact network formed by a limited number of SDS molecules that serves as a structural scaffold preserving the helical conformation of this module. This mechanism might prevail when a peptide with low helical propensity is involved in structure consolidation. We suggest that folding of peptides sharing this feature does not require a preformed tightly‐packed protein core. Thus, the formation of specific tertiary interactions would be the consequence of peptide folding and not its cause. In this scenario, folding might be thought of as a process that includes unspecific rounds of structure stabilization guiding the protein to the native state. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Fluid and condensed ApoA-I/phospholipid monolayers provide insights into ApoA-I membrane insertion

Apolipoprotein A-I (ApoA-I) is a protein implicated in the solubilization of lipids and cholesterol from cellular membranes. The study of ApoA-I in phospholipid (PL) monolayers brings relevant information about ApoA-I/PL interactions. We investigated the influence of PL charge and acyl chain organization on the interaction with ApoA-I using dipalmitoyl-phosphatidylcholine, dioleoyl-phosphatidylcholine and dipalmitoyl-phosphatidylglycerol monolayers coupled to ellipsometric, surface pressure, atomic force microscopy and infrared (polarization modulation infrared reflection-absorption spectroscopy) measurements. We show that monolayer compressibility is the major factor controlling protein insertion into PL monolayers and show evidence of the requirement of a minimal distance between lipid headgroups for insertion to occur, Moreover, we demonstrate that ApoA-I inserts deepest at the highest compressibility of the protein monolayer and that the presence of an anionic headgroup increases the amount of protein inserted in the PL monolayer and prevents the steric constrains imposed by the spacing of the headgroup. We also defined the geometry of protein clusters into the lipid monolayer by atomic force microscopy and show evidence of the geometry dependence upon the lipid charge and the distance between headgroups. Finally, we show that ApoA-I helices have a specific orientation when associated to form clusters and that this is influenced by the character of PL charges. Taken together, our results suggest that the interaction of ApoA-I with the cellular membrane may be driven by a mechanism that resembles that of antimicrobial peptide/lipid interaction.     

Structure of Self-Aggregated Alamethicin in ePC Membranes Detected by Pulsed Electron-Electron Double Resonance and Electron Spin Echo Envelope Modulation Spectroscopies

PELDOR spectroscopy was exploited to study the self-assembled super-structure of the [Glu(OMe)^7,18,19]alamethicin molecules in vesicular membranes at peptide to lipid molar ratios in the range of 1:70-1:200. The peptide molecules were site-specifically labeled with TOAC electron spins. From the magnetic dipole-dipole interaction between the nitroxides of the monolabeled constituents and the PELDOR decay patterns measured at 77 K, intermolecular-distance distribution functions were obtained and the number of aggregated molecules (n ≈ 4) was estimated. The distance distribution functions exhibit a similar maximum at 2.3 nm. In contrast to Alm16, for Alm1 and Alm8 additional maxima were recorded at 3.2 and ∼5.2 nm. From ESEEM experiments and based on the membrane polarity profiles, the penetration depths of the different spin-labeled positions into the membrane were qualitatively estimated. It was found that the water accessibility of the spin-labels follows the order TOAC-1 > TOAC-8 ≈ TOAC-16. The geometric data obtained are discussed in terms of a penknife molecular model. At least two peptide chains are aligned parallel and eight ester groups of the polar Glu(OMe)^18,19 residues are suggested to stabilize the self-aggregate superstructure.     

Distance dependence and salt sensitivity of pairwise, coulombic interactions in a protein

Histidine pK(a) values were measured in charge-reversal (K78E, K97E, K127E, and K97E/K127E) and charge-neutralization (E10A, E101A, and R35A) mutants of staphylococcal nuclease (SNase) by (1)H-NMR spectroscopy. Energies of interaction between pairs of charges (DeltaG(ij)) were obtained from the shifts in pK(a) values relative to wild-type values. The data describe the distance dependence and salt sensitivity of pairwise coulombic interactions. Calculations with a continuum electrostatics method captured the experimental DeltaG(ij) when static structures were used and when the protein interior was treated empirically with a dielectric constant of 20. The DeltaG(ij) when r(ij) < or = 10 A were exaggerated slightly in the calculations. Coulomb's law with a dielectric constant near 80 and a Debye-Hückel term to account for screening by the ionic strength reproduced the salt sensitivity and distance dependence of DeltaG(ij) as well as the structure-based method. In their interactions with each other, surface charges behave as if immersed in water; the Debye length describes realistically the distance where interactions become negligible at a given ionic strength. On average, charges separated by distances (r(ij)) approximately 5 A interacted with DeltaG(ij) approximately 0.6 kcal/mole in 0.01 M KCl, but DeltaG(ij) decayed to < or =0.10 kcal/mole when r(ij) = 20 A. In 0.10 M KCl, DeltaG(ij) approximately 0.10 kcal/mole when r(ij) = 10 A. In 1.5 M KCl, only short-range interactions with r(ij) < or = 5 A persisted. Although at physiological ionic strengths the interactions between charges separated by more than 10 A are extremely weak, in situations where charge imbalance exists many weak interactions can cumulatively produce substantial effects.     

光谱和微量热法分析柑橘苷与BSA分子间作用及构建其理论模型

光谱和微量热法分析柑橘苷(naringin,NAR)与牛血清白蛋白（bovine serum albumin, BSA）分子间作用,构建NAR与BSA分子间作用的理论模型。采用紫外-荧光光谱法解析Fōrster方程求得NAR与BSA分子间作用及分子间作用的临界距离r,等温滴定微量热技术测定NAR与BSA分子间作用的积分量热曲线,获得Δ H并通过Gibbs-Helmholtz方程获取Δ S和Δ G。基于光谱和微量热辅助分析,构建NAR与BSA分子间作用的理论模型。结果表明,光谱法测出NAR与BSA发生分子间作用,NAR与BSA分子间作用的临界距离为2.06 nm,表明NAR与BSA分子间作用为短程分子间作用。微量热法成功测定出NAR与BSA分子间的热力学参数Δ H<0,Δ S>0,Δ G<0,说明NAR与BSA分子间作用是自发进行的放热相互作用。依据Ross理论分析NAR与BSA分子间作用力主要是疏水作用力和静电作用力。模型构建结果说明,NAR与BSA分子间作用主要发生在BSA的domain IIA区域,NAR与BSA分子间作用力主要是静电作用力,兼有范德华作用力和氢键。实验与理论模型构建结果基本一致。本研究工作可为深入了解蛋白质与大分子化合物间的作用以及研究微观药理学机制提供有益的参考。     

Structural insight into the XTACC3/XMAP215 interaction from CD and NMR studies on model peptides

TACC3 is a centrosomal adaptor protein that plays important roles during mitotic spindle assembly. It interacts with chTOG/XMAP215, which catalyzes the addition of tubulin dimers during microtubule growth. A 3D coiled‐coil model for this interaction is available but the structural details are not well described. To characterize this interaction at atomic resolution, we have designed a simplified version of the system based on small peptides. Four different peptides have been studied by circular dichroism and nuclear magnetic resonance both singly and in all possible combinations; namely, five peptide pairs and two trios. In cosolvents, all single peptides tend to adopt helical conformations resembling those of the full‐length protein. However, neither the single peptides nor pairs of peptides form coiled coils. We show that the simultaneous presence of all preformed helices is a prerequisite for binding. The simplest 3D model for the interaction, based on the NMR results, is proposed. Interestingly, the peptide's structure remains unaffected by mutations at essential positions for TACC3 activity. This suggests that the lack of interaction of this TACC3 mutant with XMAP does not correlate with changes in the protein structure and that specific interactions are likely responsible for the interaction and stability of the complex.     

Biological consequences of improving the structural stability of hairpins that have antimicrobial activity

Designing new antimicrobial peptides (AMPs) focuses heavily on the activity of the peptide and less on the elements that stabilize the secondary structure of these peptides. Studies have shown that improving the structure of naturally occurring AMPs can affect activity and so here we explore the relationship between structure and activity of two non‐naturally occurring AMPs. We have used a backbone‐cyclized peptide as a template and designed an uncyclized analogue of this peptide that has antimicrobial activity. We focused on beta‐hairpin‐like structuring features. Improvements to the structure of this peptide reduced the activity of the peptide against gram‐negative, Escherichia coli but improved the activity against gram‐positive, Corynebacterium glutamicum. Distinctions in structuring effects on gram‐negative versus gram‐positive activity were also seen in a second peptide system. Structural improvements resulted in a peptide that was more active than the native against gram‐positive bacterium but less active against gram‐negative bacterium. Our results show that there is not always a correlation between improved hairpin‐structuring and activity. Other factors such as the type of bacteria being targeted as well as net positive charge can play a role in the potency of AMPs. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.     

Optimization of the electrostatic interactions between ionized groups and peptide dipoles in proteins.

The three-dimensional optimization of the electrostatic interactions between the charged amino acid residues and the peptide partial charges was studied by Monte-Carlo simulations on a set of 127 nonhomologous protein structures with known atomic coordinates. It was shown that this type of interaction is very well optimized for all proteins in the data set, which suggests that they are a valuable driving force, at least for the native side-chain conformations. Similar to the optimization of the charge-charge interactions (Spassov VZ, Karshikoff AD, Ladenstein R, 1995, Protein Sci 4:1516-1527), the optimization effect was found more pronounced for enzymes than for proteins without enzymatic function. The asymmetry in the interactions of acidic and basic groups with the peptide dipoles was analyzed and a hypothesis was proposed that the properties of peptide dipoles are a factor contributing to the natural selection of the basic amino acids in the chemical composition of proteins.     

DNA A-tracts bending: polarization effects on electrostatic interactions across their minor groove

Bending by the DNA A-tracts constitutes a contentious issue, suggesting deficiencies in the physics employed so far. Here, we inquire as to the importance in this bending of many-body polarization effects on the electrostatic interactions across their narrow minor groove. We have done this on the basis of the findings of Jarque and Buckingham who developed a procedure based on a Monte Carlo simulation for two charges of the same sign embedded in a polarizable medium. Remarkably, the present analysis reveals that for compact DNA conformations, which result from dynamic effects, an overall attractive interaction operates between the phosphate charges; this interaction is especially strong for the narrow minor groove of the A-tracts, suggesting a tendency for DNA to bend toward this groove. This tendency is in agreement with the conclusions of electrophoretic and NMR solution studies. The present analysis is also consistent with the experimental observations that the minor groove is much more easily compressible than the major groove and the bending propensity of the A-tracts is greatly reduced at “premelting” temperatures. By contrast, the dielectric screening model predicts a repulsion between the phosphate charges and is not consistent with the aforementioned bending tendency or experimental observations.     

Electrostatic interactions in proteins: Calculations of the electrostatic term of free energy and the electrostatic potential field

The pH-dependence of the electrostatic energy of interactions between titratable groups is calculated for some well studied globular proteins: basic pancreatic trypsin inhibitor, sperm whale myoglobin and tuna cytochrome c. The calculations are carried out using a semi-empirical appraach in terms of the macroscopic model based on the Kirkwood-Tanford theory. The results are discussed in the light of their physicochemical and biological properties. It was found that the pH-dependence of the electrostatic energy correlates with the III–IV transition of cytochrome c. The electrostatic field of the cysteine proteinase inhibitor, cystatin, was calculated in two ways. In the first one, the electrostatic field created by the pH dependent charges of the ionizable groups and peptide dipoles was calculated using the approach proposed. In the second one, the finite-difference method was used. The results obtained by the two methods are in overall agreement. The calculated field was discussed in terms of the binding of cystatin to papain.     

Electrodynamics of interactions in electrolyte media. Possible consequences in biological functions

Attention is drawn to the fact that electric charge interactions in aqueous electrolyte solutions, such as physiological media, depend not only on the distance between the interacting charges but also on the frequency at which they move. This fact has important consequences for some biological processes, particularly the kinetic ones. Consideration of charge screening, including the dynamic effects, shows that electrical interactions in quasi-physiological systems are effective even above Debye’s length, provided that charges move at frequencies higher than ～250 MHz. For each electrolyte solution, it is possible to define a threshold frequency, the Maxwell frequency, as the point of crossover between a conducting and a dielectric regime in the solution.     

Computational design of the Fyn SH3 domain with increased stability through optimization of surface charge charge interactions

Computational design of surface charge-charge interactions has been demonstrated to be an effective way to increase both the thermostability and the stability of proteins. To test the robustness of this approach for proteins with predominantly beta-sheet secondary structure, the chicken isoform of the Fyn SH3 domain was used as a model system. Computational analysis of the optimal distribution of surface charges showed that the increase in favorable energy per substitution begins to level off at five substitutions; hence, the designed Fyn sequence contained four charge reversals at existing charged positions and one introduction of a new charge. Three additional variants were also constructed to explore stepwise contributions of these substitutions to Fyn stability. The thermodynamic stabilities of the variants were experimentally characterized using differential scanning calorimetry and far-UV circular dichroism spectroscopy and are in very good agreement with theoretical predictions from the model. The designed sequence was found to have increased the melting temperature, DeltaT (m) = 12.3 +/- 0.2 degrees C, and stability, DeltaDeltaG(25 degrees C) = 7.1 +/- 2.2 kJ/mol, relative to the wild-type protein. The experimental data suggest that a significant increase in stability can be achieved through a very small number of amino acid substitutions. Consistent with a number of recent studies, the presented results clearly argue for a seminal role of surface charge-charge interactions in determining protein stability and suggest that the optimization of surface interactions can be an attractive strategy to complement algorithms optimizing interactions in the protein core to further enhance protein stability.     

Age‐related cleavages of crystallins in human lens cortical fiber cells generate a plethora of endogenous peptides and high molecular weight complexes

Low molecular weight peptides derived from the breakdown of crystallins have been reported in adult human lenses. The proliferation of these LMW peptides coincides with the earliest stages of cataract formation, suggesting that the protein cleavages involved may contribute to the aggregation and insolubilization of crystallins. This study reports the identification of 238 endogenous LMW crystallin peptides from the cortical extracts of four human lenses representing young, middle and old‐age human lenses. Analysis of the peptide terminal amino acids showed that Lys and Arg were situated at the C‐terminus with significantly higher frequency compared to other residues, suggesting that trypsin‐like proteolysis may be active in the lens cortical fiber cells. Selected reaction monitoring analysis of an endogenous αA‐crystallin peptide (αA_57‐65) showed that the concentration of this peptide in the human lens increased gradually to middle age, after which the rate of αA_57‐65 formation escalated significantly. Using 2D gel electrophoresis/nanoLC‐ESI‐MS/MS, 12 protein complexes of 40–150 kDa consisting of multiple crystallin components were characterized from the water soluble cortical extracts of an adult human lens. The detection of these protein complexes suggested the possibility of crystallin cross‐linking, with these complexes potentially acting to stabilize degraded crystallins by sequestration into water soluble complexes. Proteins 2015; 83:1878–1886. © 2015 Wiley Periodicals, Inc.