Structural attributes for the recognition of weak and anomalous regions in coiled-coils of myosins and other motor proteins
© Sunitha et al.; licensee BioMed Central Ltd. 2012
Received: 13 March 2012
Accepted: 28 August 2012
Published: 25 September 2012
Coiled-coils are found in different proteins like transcription factors, myosin tail domain, tropomyosin, leucine zippers and kinesins. Analysis of various structures containing coiled-coils has revealed the importance of electrostatic and hydrophobic interactions. In such domains, regions of different strength of interactions need to be identified since they could be biologically relevant.
We have updated our coiled-coil validation webserver, now called COILCHECK+, where new features were added to efficiently identify the strength of interaction at the interface region and measure the density of charged residues and hydrophobic residues. We have examined charged residues and hydrophobic ladders, using a new algorithm called CHAHO, which is incorporated within COILCHECK + server. CHAHO permits the identification of spatial charged residue patches and the continuity of hydrophobic ladder which stabilizes and destabilizes the coiled-coil structure.
The availability of such computational tools should be useful to understand the importance of spatial clustering of charged residues and the continuity of hydrophobic residues at the interface region of coiled-coil dimers. COILCHECK + is a structure based tool to validate coiled-coil stability; it can be accessed athttp://caps.ncbs.res.in/coilcheckplus.
KeywordsCoiled-coil Charged-patch Hydro ladder Pseudoenergies Charged clusters Heptads Knobs-into-holes packing
Structure and function of proteins are primarily determined from their amino acid sequences. The alpha-helical coiled-coils are simple structural units that consist of repeating blocks of seven residues which are commonly termed as 'heptads'. They are mostly seen as protein-protein interaction domains which mediate many vital functions in the system like oligomerization, cell division, transport of macromolecules, mobility and transcription. Coiled-coils consist of two or more helices wound around as a superhelix - like strand of a rope. On the basis of the number of helices involved in coil formation, they are differentiated into higher order structures like dimers, trimers, tetramers and pentamers[1–4]. The most commonly observed coiled-coil types are dimers with two α-helices and trimers with three α-helices wrapped around each other into a left-handed superhelix.
Coiled-coils usually contain a repeated pattern, “HPPHCPC” of hydrophobic (H), polar (P) and charged (C) amino acid residues referred to as heptad repeat. The position of the heptad repeat is usually labeled as 'abcdefg', where ‘a’ and ‘d’ positions are ideally occupied by the hydrophobic residues and ‘e’ and ‘g’ are occupied by oppositely charged residues. The crucial fact in the structure of coiled-coils is the burial of the hydrophobic residues to form a core region, permitting it to form an amphipathic structure, which provides the thermodynamic driving force for oligomerization[5, 6]. The hallmark of coiled-coils is the distinctive packing of amino acid side chain from one helix (knob) into a space surrounded by four side chains (hole) of the other helix named as ‘knobs-into-holes’. The packing in a coiled-coil interface is exceptionally tight, with almost complete Van der Waals contact between the side chains of residues at ‘a’ and ‘d’ positions. The charged residues at ‘e’ and ‘g’ positions also have an important role in the stability of the coiled-coils. Altogether, the set of interactions observed between ‘a’, ‘d’, ‘e’ and ‘g’ positions of the heptad repeat forms the basis for knobs-into-holes packing in coiled-coils.
The coiled-coil structures can be broadly grouped into short-length and long- length coiled-coils. The short-length coiled-coils mostly act as dimerization domains in transcription factors. The long-length coiled-coils, which consist of several hundreds of amino acid residues, are found in variety of proteins like intermediate filaments, myosins, kinesins and SMC proteins. The non-ideality of the heptad register is a common fact seen in most of the long-length coiled-coil proteins like myosins and this could include either the presence of heptad breaks or the presence of unfavorable amino acids at the specific position of the heptad[9, 10]. These irregularities would lead to the deviation of the perfect knobs-into-holes packing of the coiled-coil domain. In addition to the packing at the core positions of a coiled-coil, the distribution of selected hydrophobic residues at the ’a’ and ‘d’ position of the heptad is also equally important for the stability of two-stranded coiled-coils[11, 12]. Besides the hydrophobic packing, the contributions of long-range electrostatic interactions also play a central role in the overall stability of coiled-coil domains[13, 14]. The strength of interaction between the coiled-coil dimer is crucial, since most of them act as oligomerization units in protein structures. Regions with different strength of interactions need to be identified since they could be biologically relevant. Specificity in coiled-coil heptads and their interactions play an important role in the oligomerization process[15, 16]. Studies on homo- and heterotypic coiled-coils have shown the importance of individual amino acids at particular heptad positions and their contribution to the stability[17, 18]. Specifically positions ‘a’ and ‘d’ are studied extensively by substituting 20 different amino acids in a de novo-designed coiled-coil model to identify the stability and the oligomerization state[19, 20]. Though coiled-coils have a simple geometry forming dimers and trimers predominantly in nature, still they are capable of forming higher order structures. These complex structures are very interesting due to their structural organization, Complexes with an oligomerization state above pentamer are gaining attraction due to their barrel-like appearance and diverse function. These studies show the importance of coiled-coils and the specificity at the interface which is related to the structural organization of these simple domains.
There are different methods currently available that could predict the coiled-coil forming region from sequences, based on amino acid propensities, profile-profile comparisons, pairwise residue correlations and HMM-based approaches[22–30]. Structural analysis programs like SOCKET can determine the extent of knobs-into-holes interaction at the helix interface; thus, regions of tight packing between helices engaged in coiled-coil can be identified efficiently. In addition to the sequence database of coiled-coils from the genome of Arabidopsis thaliana, a relational database CC + of coiled-coil structures processed by the SOCKET program is also available to the public domain[33, 34]. Detailed analysis of long-length coiled-coils, for example tropomyosin structure, would show that many regions are quite deviated from ideality which is still unexplored. In such and other examples, weak/relatively-flexible regions within the coiled-coil are associated with a biological function of the molecule. Such emerging facts provide us an impetus to develop a method to validate a coiled-coil region given the structure based on the strength of interaction between the helices forming coiled-coils. One popular manner of associating strength or stability is to attribute psuedoenergies to the system of interest. The potential energy due to intermolecular interactions at coiled-coils could be viewed as a sum of individual components such as van der Waals and electrostatic intermolecular interactions, as described in our in-house program COILCHECK. The gap between sequence, structure and stability of a coiled-coil molecule is efficiently bridged by the COILCHECK program, which acts as a promising tool to identify and compare weak and strong regions within coiled-coil structures.
The magnitude and nature of the interaction between the coiled-coil dimer is provided as energy per residue in COILCHECK (for details on COILCHECK method and energy range please see). Since then, the program has been updated to include new features and we refer to the new version of webserver as ‘COILCHECK+’. New features include a tool for hydrogen atom fixing and hydrogen bond energy calculation to accurately describe the prevailing stabilizing interaction at the coiled-coil interface (please see details below). An additional option is provided to choose distant-dependent dielectric constant in electrostatic energy calculation. We have also analyzed newer structural entries using the updated method and standardized the energy ranges for acceptable and stable coiled-coil regions. We also report the availability of a new structural analysis algorithm, called CHAHO within COILCHECK + webserver, for the presence of charged residue patches and ‘hydrophobic ladder’ within given coiled-coils. Finally, we illustrate the application of COILCHECK + on examples of long length coiled-coils for which biophysical studies are emerging to show their relative instability.
Benchmarking study of COILCHECK+
COILCHECK + energies before and after SOCKET for the 112 structural entries
COILCHECK + energy before SOCKET (kJ/mol)
COILCHECK + energy after SOCKET (kJ/mol)
COILCHECK + energy before SOCKET (kJ/mol)
COILCHECK + energy after SOCKET (kJ/mol)
COILCHECK + energy before SOCKET (kJ/mol)
COILCHECK + energy after SOCKET (kJ/mol)
1L8D - 1
2CH7 - 1
1C1G - 1
1L8D - 2
2CH7 - 2
1C1G - 2
1CII - 1
1CII - 2
1CII - 3
2GD7 - 1
2GD7 - 2
To illustrate that the COILCHECK + energies are synchronous with energy calculations obtained using molecular mechanics and also to show that COILCHECK energies do reflect the overall quality of the coiled coil structures, we have compared the energies of few structures at different iterations of minimization process against its corresponding COILCHECK + total stabilizing energy. These structures were minimized using Tripos force field (SYBYL 7.1, Tripos Inc) and energies were calculated at ten different minimization steps (10–100 iterations). Minimization was carried out using Powell’s gradient with non-bonded interaction cut-off value of 8 and a distance-dependent dielectric constant equal to 1. Initial optimization was done using Simplex method and minimization was terminated at a convergence of 0.05 kcal mol Å-1. The comparison between SYBYL and COILCHECK + energies are shown inAdditional file 5. There exists a high correlation between both the energy values. This shows that COILCHECK + energies are sensitive enough to reflect the minute structural changes incorporated during the minimization process.
CHAHO algorithm for structural analysis of charged residue patches and hydrophobic ladder
Clusters of charged residues and their role in stability
The charged residues and long-range electrostatic interactions between these residues additionally contribute to the stability of coiled-coil dimers. Thus, using the charged-patch scoring scheme, we were able to identify stabilizing and destabilizing heptads which may have a significant role in the stability of the dimers. For the known crystal structures, where the coiled-coil boundaries are identified using SOCKET, the average charged-patch score was calculated. Most of the structures had stabilizing charged-patch score, while a few had destabilizing score (see Additional file6). Since the charged-patch scoring scheme reflects the variable distance shell and the specific type of interaction between the charged residues, we anticipate that there will be a proper correlation between the charged-patch scores and the electrostatic energy of the proteins. Structures with a net positive score for charged-patch means that they contain larger number of favorable electrostatic interactions, thus we expect them to have good COILCHECK + energy and an inverse correlation between charged patch score and COILCHECK + energy. In order to inspect the level of relationship between these two factors, the charged-patch average score was calculated for all the analyzed structures according to the process explained in the methodology and the electrostatic energy by applying Coulomb’s equation. Ideally, higher the charge-patch score the lower the COILCHECK + energy would be. Thus, we expect negative correlation between these two values. The charge-patch score and COILCHECK + total energy/residue values for the structures analyzed are shown in Additional file7. The correlation between these two values is negative, thereby suggesting that the method is sensitive in indentifying the stabilizing and destabilizing regions. It also emphasizes the high contribution of charged residue (electrostatic) interactions in contributing to the total energy of the system.
Packing, continuity and specificity of core hydrophobic residues
Regions of weak interactions identified by COILCHECK+
Apart from the various programs for predicting the coiled-coil forming region from sequence data, SOCKET is the only program that uses structural information to identify coiled-coils. Despite the apparent sequence level of simplicity in retaining a regular heptad pattern, coiled-coils display considerable degree of structural diversity which leads to the formation of parallel, anti-parallel and higher order coils. These structurally distinct features have been recognized using the straightforward approach of identifying the characteristic knobs-into-holes side chain packing at the coiled-coil interface. COILCHECK package further examines the strength of interactions at coiled-coils and pinpoints the weak and strong regions. Although COILCHECK + is not a coiled-coil locator program, it has wide area of applications in validating the strength of interacting dimers. Regions with relatively low stability of the coiled-coil interface could be vitally required for the function of the protein. For instance, these regions could be protein-protein interaction zones or it could be a nucleation point for conformational changes that are crucial for the biological function. Here, we show the advantage of using COILCHECK + in some examples of the long length coiled-coils; though these proteins have a coiled-coil structure with regular packing, still they have stable and variable regions for a functional role. These stable and unstable regions are not identified by other coiled-coil analysis programs. The identification of ‘unstable’ dimer region in coiled-coils can pave way to explore structural importance of this simple domain which is involved in a wide variety of cellular processes.
The second example is an anti-parallel coiled-coil domain of beclin1 protein (PDB id: 3Q8T). The protein is composed of a series of imperfect a-d pairing at the interface. The coiled-coil dimer domain is shown to be metastable at physiological temperatures using circular dichroism studies and the metastable property of this homodimeric coiled-coil domain is required for its ready transition to the more stable and functional beclin1-Atg14L/UVRAG heterodimer. SOCKET identifies this structure as a proper coiled-coil domain with complete heptad pattern. But the stability of the coiled-coil domain which is related to the function of the protein need to be identified, this is efficiently done by COILCHECK+. The overall COILCHECK + energy for the structure 3Q8T is -1.161 kJ/mol which is quite low showing that the coiled-coil domain is less stable. To gain further details on each heptad and variable regions, the 13-heptad full length structure was divided into eleven parts, where each part has three heptads and one heptad has a overlap with the previous structure. Each of the structure was analyzed using COILCHECK+. The energy values are shown in Figure4b, it can be clearly seen that all the parts have energy above - kJ/mol and certain regions were very variable. These less stable heptads (heptad-3,5 &6) are regions where the imperfect a-d pairing resides and the unfavorable repulsive pair E224a-D221e’ is present in heptad7 which is also found to be less stable. Li and coworkers have shown that these imperfect residue interactions are responsible for the metastable beclin1 homodimer. Additionally, instability is caused by clustered negatively charged residues near heptad- 3 and 10, long range electrostatic interactions between these ionic residues (E184, E188, E189, E190, E240, D242, D243, E244) provide large amount of unfavorable repulsive energy to the system. These unfavorable electrostatic interactions could be another factor that contributes to the lower stability of the coiled-coil domain.
The third example the tropomyosin mid-region (PDB id: 2B9C). Since tropomyosin is wound around the F-actin polymer, it should have flexible zones which would permit the coiled-coil to bend and bind with the actin. The middle region of tropomyosin coiled-coil domain is detected to be flexible which is responsible for the bending of the molecule. Though the structure is identified to be a proper coiled-coil with complete heptad pattern, COILCHECK + energy for this flexible mid-region is very poor (0.13 kJ/mol), showing the unstable nature of this region.
COILCHECK + was able to distinguish the relatively unstable flexible regions that could provide biological importance and value. Such examples also show the importance of bridging the gap between identification and validation of coiled-coils. It is shown that identifying the relatively unstable regions would also throw light in the direction of recognizing functionally important regions which is biologically relevant.
Analysis on specific systems
Study on a well-known structure: tropomyosin
Classification of charged residues of 1C1G tropomyosin structure into four categories: SB (stabilizing buried), SE (stabilizing exposed), DB (destabilizing buried), DE (destabilizing exposed)
Myosin VI medial tail and its importance
Myosins are diverse family of actin-based molecular motors. The general structure of myosins consists of a N-terminal head domain followed by neck domain of variable lengths and the C-terminal domain which has the coiled-coil and the cargo-binding domain. Different biological functions within this large family of molecular motors depend on the differences in their tail domains, commonly rich in the α-helical coiled-coil motif. Myosin VI is a unique class of unconventional myosins that has a reverse directionality on actin movement. In spite of having a short neck region, the molecule takes up ~36 nm step size and acts as a processive molecule travelling long distances along the actin filament. The myosin VI tail is composed of four distinct domains: the proximal tail (PT), medial tail (MT), distal tail (DT) and the cargo binding domain (CBD). The medial tail has gained importance in the recent past, where different views are put forward to address the structure and function of this region. The study of Spink and coworkers had showed that dimerization happens at the CBD and the MT domains, with clusters of charged residues acting as a single α-helix, stable enough to bridge the gap between the two heads in a 36 nm step. It is also noted that the MT is composed of ER/K motif which gives substantial rigidity to the α-helical region. Another study has showed that the PT region adopts a three-helix bundle conformation and the dimerization is initiated by two CBD regions at close proximity. Since the initial part of MT region has properly spaced hydrophobic residues, this region may be engaged in forming a coiled-coil, leading to the unfolding of PT and the lever arm extension is responsible for the large step of the molecule. The fact that MT region is strongly predicted to form coiled-coil by various programs like PAIRCOIL2 and COILS and the region has clusters of charged residues has urged us to study the human myosin VI MT with our developed methods to gain insight on its functional importance.
Charged residues of human myosin VI MT (medial tail) are classified into SB (stabilizing buried), SE (stabilizing exposed), DB (destabilizing buried) and DE (destabilizing exposed) categories along with buried DB (destabilizing buried) residues
COILCHECK + energies of medial tail and 908–940 modelled construct
van der Waals
The entire study on human myosin VI medial tail has shown that the two helices, assumed to form coiled-coil, are very poorly stabilized with each other due to the lack of hydrophobic residues at the core position throughout the complete stretch. Although few hydrophobic residues are present in the initial part of the tail, it appears that the repulsive electrostatic interactions would not permit the helices to form a stable coiled-coil. The dimerization of CBD, which leads to the close proximity of the tail regions, may induce interactions between charged residues permitting them to dimerize. In a physiological state (with water and salt in the protein environment), the exposed charged residues would be shielded by salt and the interaction between the buried residues would be the sole responsible unit for dimerization. However, in the MT region of myosin VI, the presence of buried charged residues suggests that this region is highly destabilized. Even though the MT region of myosin VI is predicted to be a strong coiled-coil by various programs, our analysis and energies show that the interactions between the two helices are weak and unstable. This shows the importance of examining hydrophobic ladder and spatial buried charged clusters to gain a better understanding of the stability of coiled-coils. Likewise, regions in coiled-coil with abundant solvent-exposed spatial charge clusters could be zones of protein-protein interactions. This analysis also provides impetus to carry out detailed examination of the spatial distribution of sequence features on various proteins that contain coiled-coils.
Coiled-coils are important structural motifs playing vital role in various fields of biology. It is important to have a better understanding of the interactions and energetic of the molecule which is the driving force for specific association. Regions of coiled-coils with relatively weak stability and blocks of flexible zones on the long-length coiled-coil proteins, usually coupled with imperfect heptads and/or solvent-exposed spatial charge clusters, are assumed to be potential target binding sites.
In order to understand the contribution of the different interactions to the stability of coiled-coil dimers, we had conducted an extensive structural analysis on known two-stranded coiled-coils using both SOCKET and COILCHECK. The results of the analysis had lead to an understanding that the inter-chain electrostatic interactions (contributed by clusters of charged residues) and the continuity of the core hydrophobic ladder (forming the ideal knobs-into-holes packing) are extremely important to the overall stability of the protein. An explicit analysis of hydrophobic ladder and the presence of unfavorable charged-patches can provide a better understanding of weakly stable or functionally important zones within coiled-coils. This has urged us to develop a dedicated algorithm which will validate the coiled-coil structure based on the energies, evaluate the distribution of spatial clusters of charged residues, to examine the continuity of the hydrophobic residues at the core positions and identify hot spot zones which are stable and destable. Such an algorithm will also act as a promising tool to identify structurally vital residues and to study many important protein families. This can also guide in designing point mutations which could reveal functional importance of pivotal residues, to distinguish particular regions in the coiled-coil motif rich in charged residues and to provide clues about protein-protein interaction sites.
We have updated our coiled-coil validation webserver from COILCHECK to COILCHECK+. The improved version includes a tool for fixing hydrogen atoms to a given PDB structure and calculating the hydrogen bond energy at the interface and an option to choose distant-dependent dielectric in the electrostatic energy calculations. Additionally, we have developed CHAHO a program to identify clustered charged residues and the continuity of core hydrophobic residues. This additional emphasis for the presence of charged residue clusters and detailed examination of the hydrophobic ladder was required in order to detect weakly stabilized regions which may be missed during average pseudoenergy calculations.
Geometric fixing of hydrogen atoms and calculation of hydrogen bond energy in COILCHECK+
The method utilized for fixing hydrogens is as delineated here. The position of hydrogen with respect to the connecting atom has been determined/fixed geometrically using standard bond lengths, angles and torsion angles for all types of relevant atoms to make up hydroxyl, carboxyl, methyl, methylene, tertiary groups and considering sp3, sp2, sp atomic states of hybridization following the published method. The hydrogen bond energy is calculated based on Kabsch and Sander’s equation, where only inter-chain hydrogen bonds are considered to attribute energy values. The hydrogen bonds are categorized into three classes: interchain backbone-backbone (BB), backbone-side chain (BS) and side chain-side chain (SS) interactions. The details of the residues involved in hydrogen bonding at a defined cutoff of 3.2Å (the donor-acceptor distance), type of interaction and the energy contributed are also provided in the output.
Distant-dependent dielectric in electrostatic energy calculation in COILCHECK+
The charges q1 and q2 for each charged residue were taken from CHARRM package, we had used distance dependent dielectric (DDD) constant in the energy calculation where D = 2r and r is the distance between the two charged atoms. The use of DDD in electrostatic energy calculation is able to yield an appropriate energy value depending on the strength of interaction seen between the residues.
Identification of spatial charged-patches
To identify clusters of charged residues which are responsible for stabilizing and destabilizing the intensity of interaction between the coiled-coil dimer, we developed a routine called ‘charged-patch’ within CHAHO. The aim of the charged-patch method is to identify specific regions or charged patches in the coiled-coil dimer which are important for the stability of the protein. Given a PDB structure and the two interacting chain identifiers, the program employs a sliding window of each heptad on both the chains and the residues are first classified into buried and exposed based on their relative solvent accessibility (RSA; identified by the PSA program from JOY package) (please see Figure2 for a complete protocol). For every individual charged residue of an heptad, the number of other interacting charged residues are identified for two distance ranges - a 6Å intra-chain and 12Å inter-chain distant shells. It is well-known that electrostatic interactions (between charged residues) can traverse fairly long distances, but become weak at longer lengths. In order to account for this, the scores are weighed according to the distance between two charged residues, variable distance shells (of 3Å, 6Å, 9Å and 12Å) were considered to add weight to the score to each type of interactions at different distances. Every charged residue is assigned a charged-patch score which ranges from zero to one based upon the electrostatic interactions with its partner at different distant shells. The program adopts a two-way approach of consolidation in which scoring function (SC) takes up the heptad and residue scheme. In the heptad scheme, the score is a cumulative value of the interactions of all charged residue (ACR) in the heptad and in the residue scheme every individual charged residue (ICR) is assigned a score based on its interacting partners and the distance shell. Thus stabilizing/destabilizing heptads and stabilizing/destabilizing residues are identified based on their scores and the set of residues forming the charged-patch for the stability of the system are also recognized.
The presence of charged patches on the surface of coiled coils could be required for interactions with other proteins. However, clusters of like-charged residues within charge patches could contribute to poor stability to a coiled-coil system. Therefore, the solvent burial of charged residues were considered and charged patches were categorized as solvent-exposed (E) or solvent-buried (B). Based on their RSA values and scores, charged residues are categorized into four groups: Stabilizing buried (SB), Stabilizing exposed (SE), Destabilizing buried (DB) and Destabilizing exposed (DE) pairs of residues. Apart from the charged-patch, which states the magnitude of interaction within the coiled-coil dimer interface, the distribution of charged residues in the sequence is also found to be structurally important, because it additionally conveys the information about the probable sites for protein-protein interactions. Hence, at a defined sliding window size, the distribution of positive and negative charged residues and the net charges are identified and displayed. Finally, the identified charged-patch at the residue and heptad levels and the distribution of charged residues are mapped on the structure. The identification of such spatial charged-patches will permit us to recognize hot spots which are crucial for the stability of the structure. This could also enable us to identify key residues that can either stabilize or destabilize the system.
Measurement of hydrophobic content in hydrophobic ladder
The ‘hydrophobic-ladder’ score, part of CHAHO algorithm, is designed to examine the hydrophobic content of ‘a’ and ‘d’ positions of the assigned heptads at both helices that engage to form the coiled-coil. The more hydrophobic these positions are at a heptad, the higher the score will be. To begin with, we assigned a score of 0.25 if any of these four positions in a coiled-coil heptad was occupied by a hydrophobic residue (Ala/Val/Leu/Ile/Phe/Tyr/Trp) and 0 if it is not. A maximum score of 1 can be expected for a heptad. Whereas ideally ‘a’ and ’d’ positions at both helices of a heptad is expected to retain hydrophobic residues for conducive and stabilizing interactions, this may not be the case in real-time sequences. The presence of every hydrophobic residue at any of the four positions in a heptad contributes 0.25 each to the score at a heptad.
Next, we also devised a propensity-based scoring scheme to normalize for the frequency of occurrence of certain hydrophobic residues. The propensity values of the hydrophobic residues from MTIDK matrix are directly incorporated to show the amino acid preferences at the particular positions. Scores were assigned to every heptad and graphs are plotted where the values are smoothened over a five-heptad sliding window. Such smoothened graphs will allow the user to identify regions which are stable and likely to be unstable, thereby providing a clue about the flexibility of the structure at particular regions.
Availability and requirements
Project name: COILCHECK+.
Project home page:http://caps.ncbs.res.in/coilcheckplus.
Operation system(s): Ubuntu Linux 11.10+.
Programming language: Java, C++, Perl, PHP.
Other requirements: None.
Any restrictions to use by non-academics: COILCHECK + can be used free-of-charge by non-academics, provided appropriate citation and credit is given to the authors of this publication.
The desired format of input file would be a PDB file with two chain identifiers. However, the user can select two chains out of a multi-chain PDB file as well. The COILCHECK + server calculates energies between two chains of a PDB file; ideally the parameters have been standardized with dimeric coiled-coils. Parallel/antiparallel, homodimeric, heterodimeric coiled-coils can be validated using COILCHECK + server. Additionally, user can select from different options for type of electrostatics calculation, charged-patch and hydrophobic ladder algorithm, apart from the default options. The outputs can also be selected by the user for particular energy component and visualization of different type of interacting residues are provided. The output page contains energy terms, possible interactions, charged-patch score and hydrophobic ladder score plots. The energy terms provided are hydrogen bond, electrostatics and Van der Waals components, along with the corresponding potential interactions like hydrophobic interactions and short contacts. The hydrophobic ladder score and charge patch scores are averaged over a seven-residue window and can be used best when both the chains start from equivalent corresponding heptads.
Availability of supporting data
Algorithm for the analysis of charged residue clusters and hydrophobic ladders
Distance dependent dielectric constant.
This study was supported by Human Frontier Science Program (HFSP). We thank James A. Spudich for his fruitful discussions during the development of the CHAHO algorithm.
- Mckay DB, Steitz TA: Structure of cataboilite gene activator protein at 2.9Å resolution suggests binding to left-handed B-DNA. Nature. 1981, 290: 744-749. 10.1038/290744a0.PubMedView ArticleGoogle Scholar
- Wilson IA, Skehel JJ, Wiley DC: Structure of the haemagglutinin membrane glycoprotein of influenza virus at 3 Å resolution. Nature. 1981, 289: 366-373. 10.1038/289366a0.PubMedView ArticleGoogle Scholar
- Malashkevich VN, Kammerer RA, Efimov VP, Schulthess T, Engel J: The crystal structure of a five-stranded coiled coil in COMP: a prototype Ion channel. Science. 1996, 274: 761-765. 10.1126/science.274.5288.761.PubMedView ArticleGoogle Scholar
- Liu J, Zheng Q, Deng Y, Cheng CS, Kallenbach NR, Lu M: A seven-helix coiled coil. Proc Natl Acad Sci USA. 2006, 103: 15457-15462. 10.1073/pnas.0604871103.PubMedPubMed CentralView ArticleGoogle Scholar
- Burkhard P, Stetefeld J, Strelkov SV: Coiled coils: a versatile protein folding motif. Trends Cell Biol. 2001, 11: 82-88. 10.1016/S0962-8924(00)01898-5.PubMedView ArticleGoogle Scholar
- Lupas A: Coiled coils:new structures and new functions. Trends Biochem Sci. 1996, 21: 375-382.PubMedView ArticleGoogle Scholar
- Mason JM, Arndt KM: Coiled coil domains: stability, specificity and biological implications. ChemBioChem. 2004, 5: 170-176. 10.1002/cbic.200300781.PubMedView ArticleGoogle Scholar
- Rose A, Meier I: Scaffolds, levers, rods and springs: diverse cellular functions of long coiled-coil proteins. Cell Mol Life Sci. 2004, 61: 1996-2009.PubMedView ArticleGoogle Scholar
- Brown JH: Breaking symmetry in protein dimers: design and functions. Protein Sci. 2006, 15: 1-13. 10.1110/ps.051658406.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown JH, Cohen C, Parry DA: Heptad breaks in alpha-helical coiled coils: stutters and stammers. Protein. 1996, 26: 134-145. 10.1002/(SICI)1097-0134(199610)26:2<134::AID-PROT3>3.0.CO;2-G.View ArticleGoogle Scholar
- Zhu BY, Zhou NE, Kay CM, Hodges RS: Packing and hydrophobicity effects on protein folding and stability: effects of beta-branched amino acids, valine and isoleucine, on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers. Protein Sci. 1993, 2: 383-394.PubMedPubMed CentralView ArticleGoogle Scholar
- Kwok SC, Hodges RS: Stabilizing and destabilizing clusters in the hydrophobic core of long Two-stranded alpha-helical coiled-coils. J Biol Chem. 2004, 279: 21576-21588. 10.1074/jbc.M401074200.PubMedView ArticleGoogle Scholar
- Manson JM, Hagemann UB, Arndt KM: Role of hydrophobic and electrostatic interactions in coiled coil stability and specificity. Biochemistry. 2009, 48: 10380-10388. 10.1021/bi901401e.View ArticleGoogle Scholar
- Sagui C, Darden TA: Molecular dynamics simulations of biomolecules :Long-range electrostatic effects. Annu Rev Biophys Biomol Struct. 1999, 28: 155-179. 10.1146/annurev.biophys.28.1.155.PubMedView ArticleGoogle Scholar
- Grigoryan G, Keating AE: Structural specificity in coiled-coil interactions. Curr Opin Struct Biol. 2008, 18: 477-483. 10.1016/j.sbi.2008.04.008.PubMedPubMed CentralView ArticleGoogle Scholar
- Mason JM, Hagemann UB, Arndt KM: Role of hydrophobic and electrostatic interactions in coiled coil stability and specificity. Biochemistry. 2009, 48: 10380-10388. 10.1021/bi901401e.PubMedView ArticleGoogle Scholar
- Acharya A, Ruvinov SB, Gal J, Moll JR, Vinson C: A heterodimerizing leucine zipper coiled coil system for examining the specificity of a position interactions: amino acids I, V, L, N, A, and K. Biochemistry. 2002, 41: 14122-14131. 10.1021/bi020486r.PubMedView ArticleGoogle Scholar
- Acharya A, Rishi V, Vinson C: Stability of 100 homo and heterotypic coiled-coil a-a' pairs for ten amino acids (A, L, I, V, N, K, S, T, E, and R). Biochemistry. 2006, 45: 11324-11332. 10.1021/bi060822u.PubMedView ArticleGoogle Scholar
- Wagschal K, Tripet B, Lavigne P, Mant C, Hodges RS: The role of position a in determining the stability and oligomerization state of alpha-helical coiled coils: 20 amino acid stability coefficients in the hydrophobic core of proteins. Protein Sci. 1999, 8: 2312-2329.PubMedPubMed CentralView ArticleGoogle Scholar
- Tripet B, Wagschal K, Lavigne P, Mant CT, Hodges RS: Effects of side-chain characteristics on stability and oligomerization state of a de novo-designed model coiled-coil: 20 amino acid substitutions in position "d". J Mol Biol. 2000, 300: 377-402. 10.1006/jmbi.2000.3866.PubMedView ArticleGoogle Scholar
- Woolfson DN, Bartlett GJ, Bruning M, Thomson AR: New currency for old rope: from coiled-coil assemblies to α-helical barrels. Curr Opin Struct Biol. 2012, 22: 1-10. 10.1016/j.sbi.2012.01.001.View ArticleGoogle Scholar
- Lupas A, Van Dyke M, Stock J: Predicting coiled coils from protein sequence. Science. 1991, 252: 1162-1164. 10.1126/science.252.5009.1162.PubMedView ArticleGoogle Scholar
- Lupas A: Prediction and analysis of coiled coil structures. Methods Enzymol. 1996, 266: 513-525.PubMedView ArticleGoogle Scholar
- Berger B, Wilson DB, Wolf E, Tonchev T, Milla M, Kim PS: Predicting coiled coils by use of pairwise residue correlations. Biochemistry. 1995, 92: 8259-8263.Google Scholar
- McDonnell AV, Jiang T, Keating AE, Berger B: Paircoil2:improved prediction of coiled coils from sequence. Bioinformatics. 2006, 22: 356-358. 10.1093/bioinformatics/bti797.PubMedView ArticleGoogle Scholar
- Gruber M, Söding J, Lupas AN: REPPER–repeats and their periodicities in fibrous proteins. Nucleic Acids Res. 2005, 33: W239-W243. 10.1093/nar/gki405.PubMedPubMed CentralView ArticleGoogle Scholar
- Wolf E, Kim PS, Berger B: Multicoil:a program for prediction two- and three- stranded coiled coils. Protein Sci. 1997, 6: 1179-1189. 10.1002/pro.5560060606.PubMedPubMed CentralView ArticleGoogle Scholar
- Trigg J, Gutwin K, Keating AE, Berger B: Multicoil2: predicting coiled coils and their oligomerization states from sequence in the twilight zone. PLoS One. 2011, 6: e23519-10.1371/journal.pone.0023519.PubMedPubMed CentralView ArticleGoogle Scholar
- Delorenzi M, Speed T: An HMM model for coiled coil domains and a comparision with PSSM-based predictions. Bioinformatics. 2002, 18: 617-625. 10.1093/bioinformatics/18.4.617.PubMedView ArticleGoogle Scholar
- Bartoli L, Fariselli P, Krogh A, Casadio R: CCHMM_PROF: a HMM-based coiled-coil predictor with evolutionary information. Bioinformatics. 2009, 25: 2757-2763. 10.1093/bioinformatics/btp539.PubMedView ArticleGoogle Scholar
- Walshaw J, Woolfson DN: SOCKET:A program for identifying and analysing coiled-coil motifs with protein structures. J.Mol.Bio. 2001, 307: 1427-1450. 10.1006/jmbi.2001.4545.View ArticleGoogle Scholar
- Rose A, Manikantan S, Schraegle SJ, Maloy MA, Stahlberg EA, Meier I: Genome-wide identification of Arabidopsis coiled-coil proteins and establishment of the ARABI-COIL database. Plant Physiol. 2004, 134: 927-939. 10.1104/pp.103.035626.PubMedPubMed CentralView ArticleGoogle Scholar
- Moutevelis E, Woolfson DN: A periodic table of coiled-coil protein structures. J. Mol. Bio. 2009, 385: 726-732. 10.1016/j.jmb.2008.11.028.View ArticleGoogle Scholar
- Testa OD, Moutevelis E, Woolfson DN: CC+: a relational database of coiled-coil structures. Nucleic Acids Res. 2009, 37: D315-D322. 10.1093/nar/gkn675.PubMedPubMed CentralView ArticleGoogle Scholar
- Bhaduri A, Pugalenthi G, Gupta N, Sowdhamini R: iMOT: an interactive package for the selection of spatially interacting motifs. Nucleic Acids Res. 2004, 32: W602-W605. 10.1093/nar/gkh375.PubMedPubMed CentralView ArticleGoogle Scholar
- Alva V, Syamala Devi DP, Sowdhamini R: COILCHECK:An interactive server for the analysis of interface regions in coiled coils. Protein Pept Lett. 2008, 15: 33-38. 10.2174/092986608783330314.PubMedView ArticleGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H: The Protein Data Bank. Nucleic Acids Res. 2000, 28: 235-242. 10.1093/nar/28.1.235.PubMedPubMed CentralView ArticleGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.PubMedView ArticleGoogle Scholar
- Glover JN, Harrison SC: Crystal structure of the heterodimeric bZIP transcription factor c-Fos-c-Jun bound to DNA. Nature. 1995, 373: 257-261. 10.1038/373257a0.PubMedView ArticleGoogle Scholar
- Strelkov SV, Herrmann H, Geisler N, Wedig T, Zimbelmann R, Aebi U, Burkhard P: Conserved segments 1A and 2B of the intermediate filament dimer: their atomic structures and role in filament assembly. EMBO J. 2002, 21: 1255-1266. 10.1093/emboj/21.6.1255.PubMedPubMed CentralView ArticleGoogle Scholar
- Dong G, Medkova M, Novick P, Reinisch KM: A catalytic coiled coil: structural insights into the activation of the Rab GTPase Sec4p by Sec2p. Mol Cell. 2007, 25: 455-462. 10.1016/j.molcel.2007.01.013.PubMedPubMed CentralView ArticleGoogle Scholar
- Li X, He L, Che KH, Funderburk SF, Pan L, Pan N, Zhang M, Yue Z, Zhao Y: Imperfect interface of Beclin1 coiled-coil domain regulates homodimer and heterodimer formation with Atg14L and UVRAG. Nat Commun. 2012, 3: 662-PubMedPubMed CentralView ArticleGoogle Scholar
- Brown JH, Zhou Z, Reshetnikova L, Robinson H, Yammani RD, Tobacman LS, Cohen C: Structure of the mid-region of tropomyosin: bending and binding sites for actin. Proc Natl Acad Sci USA. 2005, 102: 18878-18883. 10.1073/pnas.0509269102.PubMedPubMed CentralView ArticleGoogle Scholar
- Li XE, Holmes KC, Lehman W, Jung H, Fischer S: The shape and flexibility of tropomyosin coiled coil: implication for actin filament assembly and regulation. J Mol Bio. 2010, 395: 327-339. 10.1016/j.jmb.2009.10.060.View ArticleGoogle Scholar
- Hitchcock-DeGregori SE, Singh A: What makes tropomyosin an actin binding protein?. A perspective. J. Struct Bio. 2010, 170: 319-324.View ArticleGoogle Scholar
- Stewart M: Structural basis for bending tropomyosin around actin in muscle thin filament. Proc Natl Acad Sci USA. 2001, 98: 8165-8166. 10.1073/pnas.151265198.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown JH: How sequence direct bending in tropomyosin and other two-stranded alpha-helical coiled coils. Protein Sci. 2010, 19: 1366-1375. 10.1002/pro.415.PubMedPubMed CentralView ArticleGoogle Scholar
- Stewart M: Fourteen actin-binding sites on tropomyosin?. Nature. 1975, 257: 331-333. 10.1038/257331a0.PubMedView ArticleGoogle Scholar
- Brown JH, Kim KH, Jun G, Greenfield NJ, Dominguez R, Volkmann N, Hitchcock-DeGregori SE, Cohen C: Deciphering the design of the tropomyosin molecule. Proc Natl Acad Sci USA. 2001, 98: 8496-8501. 10.1073/pnas.131219198.PubMedPubMed CentralView ArticleGoogle Scholar
- Straussman R, Ben-Ya'acov A, Woolfson DN, Ravid S: Kinking the coiled coil — negatively charged residues at the coiled-coil interface. J Mol Biol. 2007, 366: 1232-1242. 10.1016/j.jmb.2006.11.083.PubMedView ArticleGoogle Scholar
- Wells AL, Lin AW, Chen LQ, Safer D, Cain SM, Hasson T, Carragher BO, Milligan RA, Sweeney HL: Myosin VI is an actin-based motor that moves backwards. Nature. 1999, 401: 505-508. 10.1038/46835.PubMedView ArticleGoogle Scholar
- Rock RS, Rice SE, Wells AL, Purcell TJ, Spudich JA, Sweeney HL: Myosin VI is a processive motor with a large step size. Proc Natl Acad Sci USA. 2001, 98: 13655-13659. 10.1073/pnas.191512398.PubMedPubMed CentralView ArticleGoogle Scholar
- Spink BJ, Sivaramakrishnan S, Lipfert J, Doniach S, Spudich JA: Long single alpha-helical tail domains bridge the gap between structure and function of myosin VI. Nat Struct Mol Biol. 2008, 15: 591-597. 10.1038/nsmb.1429.PubMedPubMed CentralView ArticleGoogle Scholar
- Sivaramakrishnan S, Spink BJ, Sim AY, Doniach S, Spudich JA: Dynamic charge interactions create surprising rigidity in the ER/K alpha-helical protein motif. Proc Natl Acad Sci USA. 2008, 105: 13356-13361. 10.1073/pnas.0806256105.PubMedPubMed CentralView ArticleGoogle Scholar
- Mukherjea M, Llinas P, Kim H, Travaglia M, Safer D, Ménétrey J, Franzini-Armstrong C, Selvin PR, Houdusse A, Sweeney HL: Myosin VI dimerization triggers an unfolding of a three-helix bundle in order to extend its reach. Mol Cell. 2009, 35: 305-315. 10.1016/j.molcel.2009.07.010.PubMedPubMed CentralView ArticleGoogle Scholar
- Nardelli M: A calculation program for calculating hydrogen atom coordinates. Comput Chem. 1982, 6: 139-152. 10.1016/0097-8485(82)80008-9.View ArticleGoogle Scholar
- Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983, 22: 2577-2637. 10.1002/bip.360221211.PubMedView ArticleGoogle Scholar
- Brooks BR, Bruccoleri RE, Olafson BD, States DJ, Swaminathan S, Karplus M: CHARMM: a program for macromolecular energy, minimization, and dynamics calculations. J Comp Chem. 1983, 4: 187-217. 10.1002/jcc.540040211.View ArticleGoogle Scholar
- Brooks BR, Brooks CL, Mackerell AD, Nilsson L, Petrella RJ, Roux B, Won Y, Archontis G, Bartels C, Boresch S, Caflisch A, Caves L, Cui Q, Dinner AR, Feig M, Fischer S, Gao J, Hodoscek M, Im W, Kuczera K, Lazaridis T, Ma J, Ovchinnikov V, Paci E, Pastor RW, Post CB, Pu JZ, Schaefer M, Tidor B, Venable RM: CHARMM: the biomolecular simulation program. J Comput Chem. 2009, 30: 1545-1614. 10.1002/jcc.21287.PubMedPubMed CentralView ArticleGoogle Scholar
- Mizuguchi K, Deane CM, Blundell TL, Johnson MS, Overington JP: JOY: protein sequence-structure representation and analysis. Bioinformatics. 1998, 14: 617-623. 10.1093/bioinformatics/14.7.617.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License(http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.