Blagoje P.
Dimitrijević
a,
Sunčica Z.
Borozan
b and
Srđan Đ.
Stojanović
*c
aClinic for Ruminants and Pigs, Faculty of Veterinary Medicine, University of Belgrade, Serbia
bDepartment of Chemistry, Faculty of Veterinary Medicine, University of Belgrade, Serbia
cICTM - Department of Chemistry, University of Belgrade, Serbia. E-mail: srdjanst@chem.bg.ac.rs; Fax: +381-11-636-061
First published on 25th October 2012
In this study we have described the π–π and cation–π interactions between the porphyrin ring and the protein part of porphyrin-containing proteins to better understand their stabilizing role. The number of π–π interactions was higher than that of cation–π interactions in the same set of proteins studied. The pyrrole groups of one porphyrin can be involved in π–π interactions with π systems of another porphyrin in the protein. We have found 5.1% cation–π interactions between porphyrin Fe2+ metal cations and π systems of surrounding amino acids as well as the pyrrole rings of other porphyrins. We observed that most of the π–π interactions have an energy in the range −0.5 to −2.0 kcal mol−1, while the cation–π interactions showed an energy in the range −2 to −4 kcal mol−1. Further, an appreciable number of metal/cation–π interaction pairs have an energy in the range −6 to −13 kcal mol−1. The preferred parallel-stacked orientation is found to be more stable than a T-shaped structure for the full set of π–π interaction pairs. In the case of cation–π interactions, it was found that 44% of the cation–π interactions involved planar stacking, 37% of the interactions belonged to the oblique category, and the remaining 19% of the interactions were of the orthogonal type. The separation distance between the cation group and the aromatic ring decreases as the interplanar angle decreases. Furthermore, in the present study we have found that 10.4% of π residues and 3.9% of cationic residues were found to have one or more stabilization centers. Amino acids deployed in the environment of porphyrin rings are deposited in helices and coils. The results from this study might be used for structure-based porphyrin protein prediction and as scaffolds for future porphyrin-containing protein design.
Noncovalent inter- and intramolecular interactions involving aromatic rings are ubiquitous in chemical and biological systems, and span from molecular recognition to self-assembly, and to catalysis and transport.9 Aromatic π–π interactions not only determine biological structures but also modulate the physical properties of residues at enzyme active sites.10 The calculated π–π interaction energies of the parallel, edge–face (T-shaped) and offset stacked are −1.48, −2.46 and −2.48 kcal mol−1, respectively,11 and the major source of attraction is not short-range (such as charge-transfer), but long-range interactions (quadrupole–quadrupole electrostatic and dispersion).12,13 It has been suggested that the perpendicular and parallel-displaced configurations are more common than the sandwich geometry as these, especially the former one, expose three aromatic faces to the outside, offering greater possibility for additional interactions with other groups.14
Cation–π interactions are unique biomolecular binding forces that occur between electron-rich aromatic rings and organic or inorganic (metallic) cations. This type of noncovalent interaction can be very strong, as has been confirmed by solid-state studies of small-molecule crystal structures15,16 and by theoretical and experimental analyses in the gas phase and in aqueous media.9,16 The strength of cation–π interactions ranges between 10 and 150 kcal mol−1,17 which is comparable to the three other major types of molecular interactions: hydrogen bonds, van der Waals interactions, and electrostatic interactions. Cation–π interactions are therefore considered to be an essential force in generating tertiary and quaternary protein structures induced by oligomerization and protein folding.18 While the cation–π interaction is arguably the strongest of the nonbonded interactions, its strength critically depends on the nature of the aromatic system and the charge of the cation.19 Ever since the study of cation–π interactions has gained prominence, the interplay of theory and experiments to see the manifestation of this reasonably strong interaction in chemistry and biology has become inevitable. In a series of pioneering papers Dougherty and co-workers have explored the role, relevance and range of the cation–π interaction between aromatic amino acids and the side chains of arginine, lysine, etc.20 Parallel to this, a series of experimental studies carried out by Rodgers and co-workers to explore the cation–π interaction involving alkali and alkaline metal ions reveal their quantitative strengths.21,22 Computations have played an important role in studying the cation–π interaction in aromatic, heteroaromatic and polycyclic aromatic compounds.23–25 While complexes with a dication have a strength comparable to covalent bonds, monocationic complexes have much smaller binding energies. In contrast, although π–π interactions are much weaker, they have been implicated in a large number of studies where their role has been discussed in eliciting a biological response.26 A strong cooperativity between cation–π interactions involving alkali and alkaline earth metal ions and π–π and hydrogen bonding interactions has also been observed.27,28
Studies of the metal center in heme proteins and model systems have shown that many factors, including noncovalent interactions, may play important roles in the properties of these metalloproteins.29–32 The importance of CH/π interactions of aromatic residues with porphyrin in hemoproteins was shown.33 It was demonstrated that the noncovalent interactions with the propionic groups of porphyrins are very important for the orientations of the imidazoles and that the conformations of the propionic groups have a strong influence on these interactions.34,35 We have previously shown that XH/π interactions are involved in the interactions of every porphyrin ring with surrounding amino acids and that these amino acids are highly conserved, demonstrating that XH/π interactions play important roles in the stability of porphyrin interactions with the protein part.36 Furthermore, we have studied strong hydrogen bonds and hydrophobic interactions,37 and non-canonical interactions38 of porphyrins in porphyrin-containing proteins. In this study, we have extended the analysis of interactions between the porphyrin ring and the protein part of porphyrin-containing proteins on all the π–π and cation–π interactions in order to better understand their stabilizing role. Results from this study might be used for structure-based porphyrin protein prediction and as scaffolds for future porphyrin-containing protein design.
A computer program Discovery Studio Visualizer 3.1 was used for the calculation of various types of π–π and cation–π interactions and their geometrical features with default settings (Fig. 1). π–π interactions are determined following the methodology of McGaughey et al.42 This method finds stacked and staggered π–π interactions by performing the following tests: the distance between the centroid of each pair of π rings is determined to find those which fall within the Center Distance (Rcen) cutoff (default 8 Å). For these, an atom from each ring should be within the Closest Atom Distance (Rclo) cutoff (default 4.5 Å). The angle θ between the normal of one or both rings and the centroid–centroid vector must fall between 0° and ± the Theta Angle cutoff (default 60°), and the angle λ between the normal to each ring must fall between 0° and ± the Lambda Angle cutoff (default 30°). The aromatic rings of the porphyrins and the amino acid side chains of His, Phe, Tyr and Trp are considered to be π systems.
Fig. 1 Parameters for π interactions: π–π interactions (a) and cation–π interactions (b). |
The following tests are performed to find cation–π interactions:43 cations are found; these are considered to be atoms with a formal charge of >0.5 to allow the inclusion of delocalized cationic species such as lysine and arginine side chains. The distance (R) between each cation and the centroid of each π ring is tested to find those within 1 Å more than the Distance Cutoff (i.e., 7 Å by default). In such cases, the distance between the cation and the mean plane of the ring is determined to find those within the Distance Cutoff (default 6 Å). The angle (θ) between the cation–centroid vector and the ring plane should be more than the Minimum Angle (default 45°). The angle α is the dihedral angle between the two extended planes. The positive charge centers of Arg (NE, CZ, CD, NH1, NH2), Lys (CE, NZ) and the metal atoms from porphyrin are considered to be cations in the present study.
We have computed the energetic contribution of all π–π and cation–π interactions for each protein in the dataset. The total interaction energy (Etotal) has been divided into electrostatic (Ees) and van der Waals (Evdw) energies and were computed using the program NAMD 2.8,44 which has implemented a subset of the OPLS force field to calculate the energies.45
Fig. 2 and 3 were prepared using the program Discovery Studio Visualizer 3.1.40
Fig. 2 Details of π–π and cation–π interactions for the HEM1001 of hemoglobin from midge larva Propsilocerus akamusi (PDB ID code 1x3k). The interactions are marked with red lines (π–π interactions: HEM1001:C—Phe46, HEM1001:D—Phe46; cation–π interactions: HEM1001:D—Arg101:NE). |
Fig. 3 Details of the interactions linking the three porphyrins of the cytochrome c-552 from Wolinella succinogenes (PDB ID code 3bnj). The π–π interactions are marked with red lines (HEM513:D—HEM515:D, HEM515:A—HEM516:A). |
Type of interaction | π–π | Cation–π | ||
---|---|---|---|---|
Numbera | %b | Numbera | %b | |
a The number of times a particular amino acid occurs in an appropriate interaction. b Percent of amino acid occurs in an appropriate interaction. c Pyrrole ring. d Chelate metal ion of porphyrin. | ||||
Amino acid | ||||
His | 12 | 5.5 | — | — |
Phe | 86 | 39.6 | — | — |
Trp | 4 | 1.8 | — | — |
Tyr | 24 | 11.1 | — | — |
Arg | — | — | 34 | 15.6 |
Lys | — | — | 10 | 4.6 |
Porphyrinc | 36 | 16.7 | — | — |
Fe2+d | — | — | 11 | 5.1 |
Total | 162 | 74.7 | 55 | 25.3 |
The number of π–π interactions was higher than that of cation–π interactions in the same set of proteins studied. Out of the aromatic residues, Phe is the most common amino acid involved in π–π interactions. This might be because, of the aromatic amino acids, Phe occurs most frequently. It is observed that Arg occurs more frequently than Lys among the cationic residues in the set of porphyrin-containing proteins, suggesting its important role in porphyrin recognition and structural stability. Because the sidechain of arginine is larger and less well water-solvated than that of lysine, it likely benefits from better van der Waals interactions with the aromatic ring. This result is consistent with previous studies of globular proteins,20 sugar-binding proteins,54 single chain “all-alpha” proteins,55 and immunoglobulin proteins.56
We have elucidated the structural motifs containing Phe residues, as the most frequently involved, responsible for the π–π interactions using the PDBeMotif web interface.57 We discovered two types of structural motifs significantly over-represented in porphyrin-containing proteins: beta-turn-ir and niche-4r. It has been shown that these small structural motifs could be involved in protein functions such as determining the conformation of enzyme active sites or binding sites.57–59 Our underlying hypothesis was that Phe residues with a high frequency of π–π interactions are probably linked to structural implications.
It is very interesting that in the proteins that contain more than one porphyrin, the pyrrole groups of one porphyrin can be involved in π–π interactions with the π systems of another porphyrin in the protein. We have found 36 (16.7%) of those interactions. It is likely that these interactions contribute significantly to the overall stability of porphyrin rings. Three iron-porphyrins from the binding pocket of the cytochrome c-552 from Wolinella succinogenes (PDB ID code 3bnj) are shown in Fig. 3. There are two π–π interactions between the porphyrins (HEM513:D—HEM515:D, HEM515:A—HEM516:A).
We have found 11 (5.1%) cation–π interactions between porphyrin Fe2+ metal cations and the π systems of surrounding amino acids as well as the pyrrole rings of other porphyrins. The interactions involve Phe, Tyr, and His amino acid residues as the π systems. Three of them include pyrrole rings. Cation–π interactions between metal cations and π systems recently emerged as one of the fundamental bonding motifs of key importance in molecular recognition.60–64
In our analysis, we have investigated multiple π interactions, an illustrative example is shown in Fig. 2. The analysis shows that about 58% of the total π interactions in the dataset are involved in the formation of multiple π interactions. This conveys that furcation is an inherent characteristic of macromolecular crystal structures.38 The furcation level of π–π interactions is much higher than that of cation–π interactions in porphyrin-containing proteins. The majority of furcated interactions exhibit longer distances than the simple non-furcated interactions as expected (Fig. 5).
We used the ConSurf server to compute the conservation score of residues involved in π–π and cation–π interactions in porphyrin-containing proteins. Among the cationic residues 37.8% of them had a conservation score of ≥6, the cutoff value used to identify the stabilizing residues. None of the cationic residues had the highest score of 9. 8.6% of the aromatic residues had the highest score of 9 and 39.7% of the residues had a conservation score of ≥6 and <9. Analysis of the conservation patterns of π–π and cation–π interactions have shown that the multiple interactions have been conserved more than the single interactions, and it is considered that structurally conserved residues are important in protein stability and folding.65
Fig. 4 Occurrence of residue pairs at different ranges of π–π ( ) and cation–π ( ; ) interaction energy. |
We have calculated the average interaction energy for all possible pairs of porphyrins and aromatic residues and the results are presented in Table 2.
Interacting pair | N (kcal mol−1) | E es (kcal mol−1) | E vdw (kcal mol−1) | E total (kcal mol−1) |
---|---|---|---|---|
a Number of interactions. b Electrostatic energy. c Van der Waals energy. d Total interaction energy. e Pyrrole ring. f Chelate metal ion of porphyrin. The standard deviations are given in parenthesis. | ||||
π–π | ||||
Porphyrin–Hise | 12 | 0.14 (0.4) | −1.05 (0.5) | −0.91 (0.4) |
Porphyrin–Phe | 86 | −0.27 (0.2) | −1.05 (0.4) | −1.32 (0.5) |
Porphyrin–Trp | 4 | −0.35 (1.1) | −1.53 (0.2) | −1.88 (1.2) |
Porphyrin–Tyr | 24 | −0.09 (0.3) | −0.96 (0.4) | −1.05 (0.5) |
Porphyrin–Porphyrin | 36 | 0.03 (0.1) | −0.77 (0.3) | −0.74 (0.3) |
Cation–π | ||||
Porphyrin–Arg | 34 | −2.15 (1.1) | −0.72 (0.5) | −2.87 (1.6) |
Porphyrin–Lys | 10 | −2.09 (0.8) | −0.24 (0.1) | −2.33 (0.9) |
Fe2+–Hisf | 1 | −9.28 | −0.01 | −9.29 |
Fe2+–Phe | 5 | −10.32 (1.4) | −0.01 (0.0) | −10.33 (1.4) |
Fe2+–Tyr | 2 | −6.82 (1.9) | −0.01 (0.0) | −6.83 (1.9) |
Fe2+–Porphyrin | 3 | −7.25 (0.8) | −0.08 (0.0) | −7.33 (0.9) |
The pairwise π–π interaction energy between the porphyrin ring and aromatic residues shows that the porphyrin–Trp (−1.88 kcal mol−1) energy is the strongest and porphyrin–His (−0.91 kcal mol−1) is the lowest among the all possible pairs as shown in Table 2. However, the energy associated with each porphyrin–porphyrin interaction, which is observed in very few proteins, is found to be the least among the pairs (−0.74 kcal mol−1). π–π interactions have a van der Waals energy that is more than five times stronger than the electrostatic energy for the interacting pairs containing porphyrin and aromatic residues. Regardless of the angle, the aromatic side chains orient in a fashion to minimize Rclo between the two rings and thus maximize the van der Waals attraction. The geometrical parameters of porphyrin–porphyrin interactions are not as good, and this might be due to steric reasons.
The average cation–π interaction energy of Arg and Lys is −2.87 and −2.33 kcal mol−1, respectively. The average energy in other globular proteins is −2.9 and −3.3, respectively, for Arg and Lys. In transmembrane helical and strand proteins, the interaction energy for Arg is stronger than Lys.66 The composition of cation–π interaction energy into electrostatic and van der Waals energy terms showed that in our dataset there is a stronger electrostatic energy than van der Waals energy.
The pair-wise cation–π interaction energy between the cationic and aromatic residues shows that the energetic contribution due to pairs including a metal (Fe2+–His and Fe2+–Phe) was the strongest among all possible pairs in porphyrin-containing proteins. From this observation, we are able to infer that the contributions of pairs including a metal are highly significant to the energetic contribution due to cation–π interactions in the dataset studied. The energies associated with cation–π interactions may be important contributors to the overall protein stability.
The native structure is the compromise of a large number of noncovalent interactions that exist in proteins and the geometrical features relating two residue-types are expected to be rather broad. However, based on the distribution of interplanar angles it was suggested that there is non-randomness in the packing of side chains.67 Using two angular parameters to define the relative orientations, each pair of planar residues is indeed found to exhibit a preference for one or two geometries, and similarly avoidance for some.
Fig. 5 Distance distribution of π–π interactions (a), and cation–π interactions (b). |
The angle distributions for the inter-ring orientational angles θ and λ shown in Fig. 6 were generated considering only π interacting pairs between the porphyrin ring and surrounding amino acid side chains. If there was no intrinsic angular energetic preference, the profiles in Fig. 6 would appear flat; instead the θ distribution (Fig. 6a) has a peak near 45° and the λ distributions (Fig. 6b) have peaks around 16 and 26°.
Fig. 6 θ angle distribution of π–π interactions (a), and λ angle distribution of π–π interactions (b). |
The interaction geometries of cation–π interactions were analyzed. In a previous study of interplanar residue contacts, Brocchieri and Karlin18 found it convenient to partition the interplanar angle (α) into three categories, planar (0° < α < 30°), oblique (30° < α < 60°) and orthogonal (60° < α < 90°). When a similar analysis was applied in the present study, it was found that 44% of the cation–π interactions involved planar stacking between the Lys (CE, NZ atoms) and Arg (NE, CZ, CD, NH1, NH2 atoms) and the aromatic ring of porphyrin. Planar interactions had an average α of 19°. Thirty-seven percent (37%) of the interactions belonged to the oblique category with an average α = 41°. The remaining 19% of the interactions were of the orthogonal type with an average α = 75°. In general, the separation distance between the cation group and the aromatic ring decreases as the interplanar angle decreases. This can be expected since coplanar stacking leads to maximal (shortest distance) contact between the interacting groups.
The geometries that are observed in abundance are not necessarily the ones that have the highest interaction energy between the two moieties in a pair, but the ones that can provide the maximum overall stability to the protein structure by the optimum use of all π interactions.
We have computed the stabilization centers for all π–π and cation–π interaction-forming residues in porphyrin-containing proteins. Table 3 shows the percentage contribution of the individual amino acid residue which is part of the stabilizing center involved in π–π and cation–π interactions with porphyrin rings.
Considering the whole data set of 553 stabilizing residues, 44 (8.0%) are involved in building π–π and cation–π interactions with the porphyrin ring. These data suggest that the amino acids that are part of stabilization centers are not involved in the building of interactions with the porphyrin ring, but they are involved in the stabilization of the surrounding structure and overall structure of the protein. We have found that 10.4% of π residues and 3.9% of cationic residues were found to have one or more stabilization centers. Aromatic residues were found to have more stabilization centers than cationic residues in porphyrin-containing proteins. It may be observed that the highest contribution comes from Phe, His, Tyr among aromatic and from Arg among charged amino acids. Interestingly, Trp and Lys were the least abundant amino acids which are stabilizing centers of the whole data set. This could be explained by the fact that these amino acids are involved in various strong and weak hydrogen bonds with the porphyrin ring and is located in the environment.38 From this we infer that these residues might contribute additional stability to the porphyrin-containing proteins in addition to their participation in π–π and cation–π interactions.
Thus, we have analyzed amino acid secondary structure preferences for the whole data set of 74 proteins. In all analyzed proteins, only two types of secondary structures were present: helices and coils. We have observed a significant number of cation–π interactions between the residues in α-helical segments, which is consistent with previous reports that the residues in α-helices have the tendency to form cation–π interactions.20,70 It was found that, among the cationic residues Lys preferred to be coiled and Arg favored coil and helix conformation. From this observation, we infer that Arg residues might stabilize coils and helices by bonding to porphyrin-containing proteins. The π residues (Phe, Tyr, Trp, His) preferred to be in a helix. This analysis indicates that the π–π and cation–π interactions do not occur at random but have a residue-specific preference for a particular secondary structure.
An interesting question concerns the location of cation–π interactions within protein structures. Cationic residues generally prefer to be on the surface of proteins whereas aromatic amino acids prefer to remain in the hydrophobic core. Because a cation–π interaction contains both a cation and an aromatic, it is not clear whether the interacting pairs should prefer to be located on the surfaces of proteins or in the cores. Traditional methods for determining residue surface accessibility rely on calculating the water-exposed surface area for a given amino acid. Because cation–π partners are necessarily in contact with one another, their water-accessible surface is diminished, even though the interacting pair as a unit may be well solvated. We observed that Arg and Lys preferred to be in an exposed region. Among the aromatic residues, it was observed that Phe and Tyr preferred to be in the partially buried region, while Trp preferred to be in the fully buried regions. This observation is quite reasonable in the sense that the aromatic residues are in principle, non polar residues, and tend to be buried. Since Arg and Lys are polar in nature they tend be exposed to the solvent surface. The polar residues might contribute significantly to the stability of the porphyrin-containing proteins as the contribution of the global energy is much greater in solvation.
We observed that the number of π–π interactions was higher than that of cation–π interactions in the same set of proteins studied. Out of the aromatic residues, Phe is the most common amino acid involved in π–π interactions. It is observed that Arg exhibits higher occurrences than Lys among the cationic residues. It is very interesting that in the proteins that contain more than one porphyrin, the pyrrole groups of one porphyrin can be involved in π–π interactions with the π systems of another porphyrin in the protein. We have found 11 (5.1%) cation–π interactions between porphyrin Fe2+ metal cations and the π systems of surrounding amino acids as well as the pyrrole rings of other porphyrins. Analysis of conservation patterns of π–π and cation–π interactions have shown that the multiple interactions have been conserved more than the single interactions. We discovered two types of structural motifs significantly over-represented in porphyrin-containing proteins: beta-turn-ir and niche-4r. It has been shown that these small structural motifs could be involved in protein functions such as determining the conformation of enzyme active sites or binding sites. The analysis of the energetic contribution of the protein interacting residues has revealed that most of the π–π interactions have an energy in the range −0.5 to −2.0 kcal mol−1. In our database it was found that cation–π interactions showed an energy less than −6 kcal mol−1, and most of them have energies in the range −2 to −4 kcal mol−1. The pair-wise cation–π interaction energy between the cationic and aromatic residues shows that the energetic contribution due to pairs including a metal (Fe2+–His and Fe2+–Phe) was the strongest among all possible pairs in porphyrin-containing proteins. The preferred parallel-stacked orientation is found to be more stable than a T-shaped structure by 0.7 kcal mol−1 for the full set of π–π interaction pairs. In the case of cation–π interactions, it was found that 44% of the cation–π interactions involved planar stacking, while 37% of the interactions belonged to the oblique category. The remaining 19% of the interactions were of the orthogonal type. The separation distance between the cation group and the aromatic ring decreases as the interplanar angle decreases. Furthermore, aromatic residues were found to have more stabilization centers than cationic residues in porphyrin-containing proteins. It may be observed that the highest contribution comes from Phe, His, Tyr among aromatic and from Arg among charged amino acids. The secondary structure and solvent accessibility of residues in porphyrin-containing proteins reveal that most of the π–π and cation–π interaction-forming residues prefer the secondary structure of alpha helical segments. It is found that in the porphyrin-containing proteins, the positively charged amino acids, Arg and Lys, involved in π–π and cation–π interactions, prefer to be in the solvent exposed surface but the aromatic amino acids (Phe and Tyr) prefer the partially buried regions of the protein, while Trp prefers to be in the fully buried regions.
In conclusion, observations obtained in this study identify π–π and cation–π interactions and structural motifs that contribute to the stabilization of the porphyrin ring by proteins, relevant to the understanding of structure–function relationships in porphyrin-containing proteins, and useful to the efforts made to design proteins able to incorporate this versatile and ubiquitous prosthetic group.
This journal is © The Royal Society of Chemistry 2012 |