- Open Access
MDcons: Intermolecular contact maps as a tool to analyze the interface of protein complexes from molecular dynamics trajectories
© Abdel-Azeim et al.; licensee BioMed Central Ltd. 2014
- Published: 6 May 2014
Molecular Dynamics (MD) simulations of protein complexes suffer from the lack of specific tools in the analysis step. Analyses of MD trajectories of protein complexes indeed generally rely on classical measures, such as the RMSD, RMSF and gyration radius, conceived and developed for single macromolecules. As a matter of fact, instead, researchers engaged in simulating the dynamics of a protein complex are mainly interested in characterizing the conservation/variation of its biological interface.
On these bases, herein we propose a novel approach to the analysis of MD trajectories or other conformational ensembles of protein complexes, MDcons, which uses the conservation of inter-residue contacts at the interface as a measure of the similarity between different snapshots. A "consensus contact map" is also provided, where the conservation of the different contacts is drawn in a grey scale. Finally, the interface area of the complex is monitored during the simulations. To show its utility, we used this novel approach to study two protein-protein complexes with interfaces of comparable size and both dominated by hydrophilic interactions, but having binding affinities at the extremes of the experimental range. MDcons is demonstrated to be extremely useful to analyse the MD trajectories of the investigated complexes, adding important insight into the dynamic behavior of their biological interface.
MDcons specifically allows the user to highlight and characterize the dynamics of the interface in protein complexes and can thus be used as a complementary tool for the analysis of MD simulations of both experimental and predicted structures of protein complexes.
- Molecular Dynamic Simulation
- Molecular Dynamic Trajectory
- Root Mean Square Fluctuation
- Conservation Rate
- Hydrophilic Residue
The thousands of proteins expressed in cells perform many of their functions through interactions with other proteins. Understanding protein-protein interactions is a crucial step in the investigation of many processes, such as intracellular signaling pathways , antibody-antigen pairing , enzyme-inhibitor interactions . Therefore, the detailed characterization on a structural basis of protein-protein complexes has become an important task for both wet and computational biologists. Thousands of experimental structures of protein-protein complexes are currently available from the wwPDB  and many more can be reliably predicted by computational approaches, as shown in the Critical Assessment of PRedicted Interactions (CAPRI) experiment [5, 6].
The dynamical characterization of such complexes can add valuable information on the recognition process, and Molecular dynamics (MD) has indeed been long established as a useful tool to help understanding biological process at the atomic and molecular levels [7–10]. The recent advances in the computational power and the development of new theoretical methods have moreover made MD simulations of large systems on a large time scale much more affordable than before [11–13]. However, MD simulations of protein complexes suffer from the lack of specific tools in the analysis step. The analysis of MD trajectories of protein complexes indeed generally relies on classical measures, such as root-mean-square deviation (RMSD), root-mean-square-fluctuation (RMSF) and radius of gyration, to list a few, conceived and developed for single macromolecules. As a matter of fact, instead, researchers engaged in simulating the dynamics of a protein complex are mainly interested in characterizing the conservation/variation of its biological interface.
Whereas a large number of specific tools have been developed to analyze single biomolecules, for example to measure the volume of pockets in proteins (Mdpocket  and POVME ), to analyze their hydrogen bonds network (HBonanza ) or to characterize the fine structural details of DNA molecules (MDDNA ) during MD simulations, there is a paucity of methods to specifically analyse MD trajectories of protein complexes. In this context, we have previously shown that intermolecular contact maps (i.e. maps where a black dot is present at the cross-over of two residues belonging to the two interacting molecules if any pair of atoms belonging to the two residues is closer than a chosen cut-off distance), can be successful applied to the analysis of both experimental and predicted structures of protein-protein complexes [18–21]. In particular, we have shown that an intermolecular contact map can identify uniquely and intuitively the surface of interaction, representing a sort of fingerprint of the complex and reporting the crucial information in a ready-to-read form. It therefore allows to easily and intuitively discriminate between similar and different binding solutions . Furthermore, we have shown that intermolecular contact maps, together with the measure of the conservation of inter-molecular contacts, can be used to analyse docking model ensembles  and to reliably extract from them the native-like solutions [20, 21].
Herein we extend this approach to the analysis of MD trajectories of protein complexes. We use the conservation of inter-residue contacts (ICs) at the interface as a measure of the similarity between different snapshots. A "consensus contact map" is also provided, where the conservation of the different contacts is drawn in a grey scale, and the complex interface area along the simulation is monitored. This novel approach is embodied in an algorithm we have named MDcons (Molecular Dynamics CONSensus). In the following, as exemplary cases, we apply MDcons to study two protein-protein complexes with comparable interface areas and binding affinities at the extremes of the experimental range. We show that differences in the dynamic behavior of the analyzed systems, particularly at their biological interface, are efficiently outlined by the MDcons approach. As a final remark, we note that the importance of inter-residue contacts when analyzing 3D structures of protein complexes is well established. In the CAPRI experiment, for instance, the correctness of a prediction, i.e. its similarity to the native structure, is assessed based on a combination of RMSD criteria and of conservation of inter-residue contacts, as compared to the native structure [5, 22]. Furthermore, the fraction of common inter-residue contacts among a set of docking models has been recently shown to successfully apply to their clustering .
Main features of the analysed complexes
Δ G (kcal/mol)
# ICs CRkl = 1
MD Interface area (Å2)
Im7 immunity protein
733 ± 38
641 ± 41
Colicin E7 nuclease - Im7 immunity protein complex (7CEI)
As already mentioned, the interaction between the two partners is very specific and strong, presenting an experimental Gibbs free-energy of binding equal to -19.50 kcal/mol . Both proteins present high percentage of hydrophilic residues at the binding interface. Our COCOMAPS  analysis revealed that 35 of the total 52 ICs indeed occur between two hydrophilic residues, while only one of them involves the two hydrophobic residues significantly buried at the interface, namely Val27/Im7 and Val523/colE7. Importantly, fifteen direct hydrogen bonds are found at the interface, five of them being salt bridges. The Asn26 and Asp31 residues of Im7 contribute significantly to the H-bonds network, the former one participating in three H-bonds, the latter one giving two salt bridges with ColE7 Arg520 and Lys525. Clearly, the charged residues at the interface play an important role in the interaction. It was shown that the triple mutant of Im7 where Asp31, Asp35 and Glu39 are mutated to Asn/Gln completely lacks in vivo inhibitory activity against ColE7, with inhibitory effect of the residues in the order Asp31 > Asp35 > Glu39 . Additionally, two Im7 tyrosines in the middle of the biological interface, Tyr55 and Tyr56, also seem to play a crucial role . They contribute one intermolecular H-bond each, but also give hydrophobic interactions with ColE7. In particular, Tyr55 is 98% and Tyr56 73%buried upon the complex formation, being inserted into ColE7 non-polar pockets (Figure 2).
MD simulations and MDcons analysis. As an indicative measure of the stability and conformational drift of the complex in the simulation, the RMSD of the C-alpha atoms from their initial position was monitored as a function of the simulation time and is reported in the Additional file 1. The RMSD trend quickly reaches a plateau at 1.5 Å after the equilibration steps, indicating that the system remains stable during the 100-ns long simulation. The radius of gyration Rγ, a property linked to the molecular volume and compactness, is also very stable after few ns, showing a very limited variation and indicating that the complex does not undergo significant conformational changes (Additional file 1). The good convergence of the simulation is reflected by the high root mean square inner product (RMSIP) value, 0.79, calculated between the first 10 principal component vectors of the two halves of the trajectory.
The overall conservation of the interface during the dynamics can be visually appreciated by the comparison of the MD consensus map and the inter-molecular contact map of the X-ray structure (reported in Figure 3) and is reflected by the C70 and C90 values. The MDcons analysis results in C70 and C90 values of 0.84 and 0.69, which means that 84% and 69% of the ICs are conserved in at least 70% and 90% of the frames, respectively, thus indicating a limited flexibility of the interface during the dynamics. In particular, of the total 52 ICs present in the X-ray structure, 19 are maintained in all the examined frames (CRkl = 1, Additional file 2). These include contacts between residues involved in 9 out of the total 15 direct inter-molecular hydrogen bonds. The interactions between Tyr55/Im7 and the side chain of Lys528/ColE7, and between Pro57/Im7 and Ser514/ColE7 are also fully conserved. As a remark, these contacts are driven by van der Waals interactions, for example between the aliphatic part of the Lys528 side chain and the aromatic ring of Tyr55, with no evidence of a H-bond or a π-cation interaction between the amino group of Lys528 and the OH group or the aromatic ring of Ty55 in the examined frames. This highlights that tools based on the detection of standard interaction motifs, such as H-bond or π-cation interaction, would completely miss this strongly maintained contact .
The remaining 33 ICs exhibit CRkl in the range 0.99 to 0.07. The least conserved ICs, among those observed in the X-ray structure, correspond to Asn26/Im7-Val523/ColE7 (CRkl = 0.14), Leu34/Im7-Lys525/ColE7 (CRkl = 0.13), Asp35/Im7-Thr539/ColE7 (CRkl = 0.10) and Asp35/Im7-Lys528/ColE7 (CRkl = 0.07). The first three interactions are non-specific and located at the external boundaries of the binding interface. Contributing residues are at a minimum distance above 4 Å in the X-ray structure. Therefore their low conservation is consistent with a small fluctuation of the two proteins around their recognition interface. The last one is instead surprisingly a salt bridge, which is clearly completely lost during the MD simulations. As we contemporarily found a conservation rate of 0.90 between Asp35/Im7 and Lys537/ColE7 giving no H-bond in the X-ray structure (Additional file 2), we decided to monitor the distances between the functional groups of Asp35 of Im7and Lys528/Lys537 of ColE7 in the time, from which it is clear that the salt-bridge between Asp35/Im7 and Lys528/ColE7 is rapidly lost during the simulation and replaced by a Asp35/Im7 salt-bridge with Lys537/ColE7 (Additional file 1). Lys537/ColE7 also establishes a salt-bridge with Glu39/Im7 that is more than 5 Å apart in the X-ray structure (see below), as shown in the inset of Figure 2 for frame 113, corresponding to 5.65 ns of simulation. Therefore, quite surprisingly, some variability in the complex salt-bridge network emerges during the dynamics simulation, around the crucial negatively charged Asp35/Im7 residue.
CD2-CD58 complex (1QA9)
MD simulations and MDcons analysis. RMSD of the C-alpha atoms from their initial position, and the radius of gyration Rγ of the system are reported in the Additional file 1. Both measures indicate that the complex is stable during the dynamics simulation. The RMSD values are pretty stable on values between 1.5 and 2.0 Å, although more fluctuating than the 7CEI ones in the first half of the simulation, where it reaches maximum values of 2.5 Å. The radius of gyration Rγ is also stable. It assumes an absolute value very close to that of the 7CEI system, although showing a slightly higher fluctuation. Also for this simulation, a high RMSIP value, 0.80, was calculated, indicating that it is very well converged.
The C70 and C90 parameters assume the values of 0.78 and 0.49, respectively (Table 1), that means that only 78% of ICs are maintained in at least 70% of the frames and only 49% in at least 90% of them, indicating a certain variability of the ICs during the simulation. This is also apparent from the complex consensus map (Figure 6b), where it is clear that additional spots appear as compared to the contact map of the X-ray structure, i.e. additional ICs are established during the simulation, while most ICs of the X-ray structure are not fully conserved (appear as gray spots in the consensus map).
Of the 38 ICs observed in the X-ray structure, only 9 (24%) are indeed rigidly maintained (CRkl = 1) during the MD simulation. These include only 4 of the 15 intermolecular hydrogen bonds, showing that the remaining 11 hydrogen bonds cannot be established at least in some of the explored frames, as the corresponding residues are more that 5 Å apart. Among the 100% conserved ICs there are also some substantially hydrophobic ones, which cannot be easily characterized with standard analysis tools based on, for example, the characterization of hydrogen-bonds at the interface. For instance, the ICs between Phe46/CD58 and Gly90/CD2 and the one between the aromatic ring of Tyr86/CD2 and the aliphatic portion of Lys34/CD58 have a CRkl = 1.0. The remaining 31 ICs of the starting X-ray structure exhibit CRkl in the range 0.99 to 0.00. The least conserved ICs, among those observed in the X-ray structure, Arg48/CD2-Arg44/CD58 (CRkl = 0.11), Arg48/CD2-Phe46/CD58 (CRkl = 0.00) and Arg48/CD2-Lys34/CD58 (CRkl = 0.00), all involve Arg48 of CD2 in the central region of the CD58 binding interface (see Additional file 2). The disruption of this network of interactions is probably due to repulsive electrostatic interaction between Arg47/CD2 on one side, and Arg44 and Lys34 of CD58 on the other side.
Interestingly, during the dynamics 59 additional ICs are established, which are not present in the starting X-ray structure, as can easily be seen from the corresponding consensus map shown in Figure 6b. Ten of them have a CRkl above 0.50 and 5 of them are very well conserved, exhibiting CRkl values above 0.80. They include ICs made by residues located at the boundary of the binding interface, such as Lys41/CD2 and Lys43/CD2, andPro80/CD58 and Asn81/CD58, but also residues in the middle of the interface, such as Tyr86 in CD2 and Leu27 and Glu39 in CD58 (Additional file 2). From the consensus map (Figure 6b), it is also easy to spot a new contact region around residue 80 of both proteins. In particular, novel contacts emerge between Lys82/CD2 and Glu78/CD58 (CRkl = 0.63) and between Tyr86/CD2 and Glu78/CD58 (CRkl = 0.29). These also are ICs in the middle of the binding interface.
In a small fraction of snapshots (1-3%), we observed ICs between residues, such as Lys91/CD2 and Phe49/CD58, which are more than 10 Å apart in the crystallographic structure, clearly indicating a significant flexibility of the complex interface, not captured by the RMSD or the Rγ. Differently from 7CEI, for 1QA9 the interface area is slightly reduced during the dynamics. The average interface area of the 2000 snapshots is indeed 641Å2 (with a standard deviation of 41), that is 35 Å2 less than the X-ray one, 676Å2 , see Figure 4.
Comparison between the two complexes
MDcons analysis of the two structures highlights a higher flexibility and variability in the interface of CD2-CD58 relative to that of Im7-ColE7. This can easily be appreciated by comparing the consensus maps obtained for the two systems (and by their comparison with the inter-molecular contact maps of the corresponding X-ray structures, see Figures 3,6) and is reflected in the lower values of C50, C70 and C90 parameters obtained for CD2-CD58 (Table 1). Specifically, while the C50 and C70 of CD2-CD58 are marginally smaller than those calculated for Im7-ColE7, its C90 is significantly smaller (0.49 vs. 0.69). Furthermore, only 9 ICs at the interface of the CD2-CD58 complex, corresponding to 24% of the contacts observed in the X-ray structure, are 100% conserved during the dynamics, which compares with the 19 fully conserved ICs at the interface of Im7-ColE7, corresponding to 37% of the contacts of the X-ray structure (Tables S1, S3). Finally, the fraction of X-ray contacts with a CRkl< 0.50 is 12% for 7CEI, indicating that only a few contacts are basically lost during the MD simulation, a value that increases to the clearly higher value of 26% for 1QA9.
Overall, the above analysis clearly indicates that the ICs of 1QA9 are somewhat more flexible than those of 7CEI, highlighting a more dynamic behavior of the interface of 1QA9. Furthermore, the average interface area of 1QA9 is reduced relative to the X-ray structure during the MD simulation, whereas the 7CEI one experiences a small increase (Figure 4).
Both systems acquire novel ICs during the MD simulations, not observed in the X-ray structure (Tables S2,S4). Precisely, 59 and 42 novel contacts are observed for 1QA9 and 7CEI, respectively. However, in CD2-CD58 the novel established ICs do not compensate the ones lost compared to the starting X-ray structure and the interface area is consequently reduced by about 35 Å2. On the contrary, the interface area of the Im7- ColE7 complex is increased on average by 40 Å2 during the MD simulation.
Once they were outlined by MDcons, we could follow and characterize the ICs that exhibited a significant variation compared to the X-ray structure by classical analysis tools, such as distance monitoring along the simulation. Results of this analysis surprisingly showed that the salt-bridge network of Im7-ColE7 is quite dynamic, with Asp35/Im7 losing its H-bond with Lys528/ColE7 during the dynamics, and acquiring one novel H-bond with Lys537/ColE7, which in turns rearranges its side-chain to also give a H-bond with Glu39/Im7.
It is interesting that we found a poor conservation in terms of ICs for the interface of the human CD2-CD58 complex and a decrease in its extension, compared to the static 3D structure, as it has been experimentally shown to have very rapid on and off kinetic rates and to exchange to new partners, after dissociation, rather than rebind to the same CD2-CD58 pair components . The rapid off rate of CD2-CD58 has previously been hypothesized to be related to the virtual absence of hydrophobic contacts at the interface. However, this is a feature that such complex shares with many other complexes, including the Im7-ColE7 one analysed here. We suggest that the CD2-CD58 peculiar binding behavior could also be related to the dynamics of the CD2 and CD58 proteins in the complex, that spontaneously lose a significant fraction of contacts while exploring novel ones.
We performed MD simulations on two protein complexes with interfaces of comparable area and both dominated by hydrophilic residues. To analyse the results of the simulations we used some classical tools and a novel tool, MDcons (freely available for download at https://www.molnac.unisa.it/BioTools/mdcons/index.php/MDcons), we specifically developed for the analysis of the MD trajectories of protein complexes. The two systems under analysis correspond to the high binding affinity Im7-ColE7 structure and to the low affinity CD2-CD58 complex. Analysis through the MDcons tool allowed us to outline interesting dynamic features at the interface of both the analysed complexes, evidencing overall a more rigid interface for Im7-ColE7, and a relatively more flexible and dynamic interface for the CD2-CD58 complex. We have thus demonstrated the utility of MDcons in complementing classical analysis tools, when studying the dynamics of protein complexes. This approach can also be straightforwardly applied to other macromolecular complexes, such as DNA-protein and RNA-protein complexes.
The starting coordinates of the Molecular Dynamics (MD) simulations of the two protein complexes selected for this study were the PDB codes: 7CEI and 1QA9 (resolution 2.30 and 3.20 Å2, respectively). MD simulations were performed using the GROMACS 4.6.2 package with the Amber99SB-ILDN force field . The complexes were slightly relaxed using 50 steps of steepest descent minimizer. The systems were then immersed in an explicit TIP3P water cubic box [45, 46], under periodic boundary conditions, which extended at least 10 Å away in each direction from any atom of the complex and contained 14267 and 12074 water molecules for the 7CEI and 1QA9 systems, respectively. Four sodium ions were added to the 1QA9 system and three chloride ions to the 7CEI system, to neutralize them as needed for the particle mesh Ewald calculation  of the long-range electrostatic interactions, while a cut-off of 10 Å was used for van der Waals and short-range electrostatic interactions. 500 steps of steepest descent minimization were performed to remove bad contacts with the solvent. All bonds involving hydrogen atoms were constrained by the LINCS algorithm . Equilibration of the solvent and ions around the complexes with position constraints of the heavy atoms was performed for 1 ns at constant temperature (300 K) in the NVT ensemble, followed by 1 ns at constant temperature (300 K) and pressure (1 atm) in the NPT ensemble. NVT simulations were carried out using the velocity rescaling thermostat (V-rescale)  and the NPT ones using the Parrinello-Rahman barostat . After equilibration, 100 ns MD simulations were performed in the NPT ensemble. Root mean square deviation from the X-ray coordinates (RMSD), the root mean square fluctuation (RMSF) and the radius of gyration were calculated with standard GROMACS tools. To assess the convergence of the trajectories, we split the trajectories into two halves, and we calculated the root mean square inner product (RMSIP) between the first 10 eigenvectors of PCA analysis on the two halves [51, 52]. The calculated RMSIP values are 0.79 for the 7CEI and 0.80 for the 1QA9 trajectory. 2000 snapshots were then generated for each system from the whole trajectory of 100 ns, by writing the coordinates every 50 ps, for further MDcons analyses. For comparison purposes, one more independent 100 ns-long NPT simulation was run for each system with different initial velocities. For each system, the similarity between the two trajectories was assessed in terms of RMSIP on the first 10 eigenvectors from PCA analysis. The obtained RMSIP value was 0.76 for 7CEI and 0.74 for 1QA9. RMSD and gyration radius values along the simulations were also compared for each system (Fig. S3). As the two trajectories for each system resulted to be very similar, for the sake of simplicity the analysis reported in the text was limited to the first trajectory for each system.
The MDcons package consists of a collection of fortran sources and python scripts. It is freely available for download at the wiki page: https://www.molnac.unisa.it/BioTools/mdcons/index.php/MDcons. It takes as input a set of PDB snapshot files or a single PDB trajectory file representing the protein-protein complex conformations during the MD simulations. MDcons then gives as output: i) a consensus map which represents the frequency of inter-residue contacts for all the MD simulation; ii) a list of the inter-residue contacts with the relative conservation rate during the MD simulations; iii) the C50, C70 and C90 coefficients, which represent respectively the number of inter-residue contacts conserved for at least 50, 70 and 90% of the analysed frames.
Where nckl is the total number of frames where residues k of protein A and l of protein B are in contact, and N is the total number of analysed frames. The conservation rate thus ranges between CRkl = 0, if the contact between residues k and l never happens, to CRkl = 1, if the contact is present in all the frames.
Finally, the interface area for each snapshot is estimated as half of the Accessible Surface Area (ASA) buried upon complex formation as calculated by NACCESS .
LC is grateful to the KAUST for generous baseline research funding. RO was supported by University Parthenope "Potenziamento linee di ricerca dipartimentali" funding.
The publication costs for this article were funded by BRF funds, KAUST.
This article has been published as part of BMC Bioinformatics Volume 15 Supplement 5, 2014: Italian Society of Bioinformatics (BITS): Annual Meeting 2013. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/15/S5
- Pawson T, Nash P: Protein-protein interactions define specificity in signal transduction. Genes Dev. 2000, 14: 1027-1047.PubMedGoogle Scholar
- Wilson IA, Stanfield RL: Antibody-antigen interactions: new structures and new conformational changes. Current opinion in structural biology. 1994, 4: 857-867. 10.1016/0959-440X(94)90267-4.View ArticlePubMedGoogle Scholar
- Jones S, Thornton JM: Principles of protein-protein interactions. Proceedings of the National Academy of Sciences of the United States of America. 1996, 93: 13-20. 10.1073/pnas.93.1.13.PubMed CentralView ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28: 235-242. 10.1093/nar/28.1.235.PubMed CentralView ArticlePubMedGoogle Scholar
- Lensink MF, Mendez R, Wodak SJ: Docking and scoring protein complexes: CAPRI 3rd Edition. Proteins. 2007, 69: 704-718. 10.1002/prot.21804.View ArticlePubMedGoogle Scholar
- Mendez R, Leplae R, Lensink MF, Wodak SJ: Assessment of CAPRI predictions in rounds 3-5 shows progress in docking procedures. Proteins. 2005, 60: 150-169. 10.1002/prot.20551.View ArticlePubMedGoogle Scholar
- Karplus M, McCammon JA: Molecular dynamics simulations of biomolecules. Nature structural biology. 2002, 9: 646-652. 10.1038/nsb0902-646.View ArticlePubMedGoogle Scholar
- Wang W, Donini O, Reyes CM, Kollman PA: Biomolecular simulations: recent developments in force fields, simulations of enzyme catalysis, protein- ligand, protein-protein, and protein-nucleic acid noncovalent interactions. Annual review of biophysics and biomolecular structure. 2001, 30: 211-243. 10.1146/annurev.biophys.30.1.211.View ArticlePubMedGoogle Scholar
- Hansson T, Oostenbrink C, van Gunsteren W: Molecular dynamics simulations. Current opinion in structural biology. 2002, 12: 190-196. 10.1016/S0959-440X(02)00308-1.View ArticlePubMedGoogle Scholar
- Karplus M, Kuriyan J: Molecular dynamics and protein function. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102: 6679-6685. 10.1073/pnas.0408930102.PubMed CentralView ArticlePubMedGoogle Scholar
- Hess B, Kutzner C, van der Spoel D, Lindahl E: GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. Journal of chemical theory and computation. 2008, 4: 435-447. 10.1021/ct700301q.View ArticleGoogle Scholar
- Gotz AW, Williamson MJ, Xu D, Poole D, Le Grand S, Walker RC: Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born. Journal of chemical theory and computation. 2012, 8: 1542-1555. 10.1021/ct200909j.PubMed CentralView ArticlePubMedGoogle Scholar
- Pronk S, Pall S, Schulz R, Larsson P, Bjelkmar P, Apostolov R, Shirts MR, Smith JC, Kasson PM, van der Spoel D: GROMACS 4.5: a high- throughput and highly parallel open source molecular simulation toolkit. Bioinformatics. 2013, 29: 845-854. 10.1093/bioinformatics/btt055.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmidtke P, Bidon-Chanal A, Luque FJ, Barril X: MDpocket: open-source cavity detection and characterization on molecular dynamics trajectories. Bioinformatics. 2011, 27: 3276-3285. 10.1093/bioinformatics/btr550.View ArticlePubMedGoogle Scholar
- Durrant JD, de Oliveira CA, McCammon JA: POVME: an algorithm for measuring binding-pocket volumes. Journal of molecular graphics & modelling. 2011, 29: 773-776. 10.1016/j.jmgm.2010.10.007.View ArticleGoogle Scholar
- Durrant JD, McCammon JA: HBonanza: a computer algorithm for molecular-dynamics-trajectory hydrogen-bond analysis. Journal of molecular graphics & modelling. 2011, 31: 5-9.View ArticleGoogle Scholar
- Dixit SB, Beveridge DL: Structural bioinformatics of DNA: a web-based tool for the analysis of molecular dynamics results and structure prediction. Bioinformatics. 2006, 22: 1007-1009. 10.1093/bioinformatics/btl059.View ArticlePubMedGoogle Scholar
- Vangone A, Spinelli R, Scarano V, Cavallo L, Oliva R: COCOMAPS: a web application to analyze and visualize contacts at the interface of biomolecular complexes. Bioinformatics. 2011, 27: 2915-2916. 10.1093/bioinformatics/btr484.View ArticlePubMedGoogle Scholar
- Vangone A, Oliva R, Cavallo L: CONS-COCOMAPS: a novel tool to measure and visualize the conservation of inter-residue contacts in multiple docking solutions. Bmc Bioinformatics. 2012, 13 (Suppl 4): S19-10.1186/1471-2105-13-S4-S19.PubMed CentralView ArticlePubMedGoogle Scholar
- Oliva R, Vangone A, Cavallo L: Ranking multiple docking solutions based on the conservation of inter-residue contacts. Proteins. 2013, 81: 1571-1584. 10.1002/prot.24314.View ArticlePubMedGoogle Scholar
- Vangone A, Cavallo L, Oliva R: Using a consensus approach based on the conservation of inter-residue contacts to rank CAPRI models. Proteins. 2013, 81: 2210-2220. 10.1002/prot.24423.View ArticlePubMedGoogle Scholar
- Mendez R, Leplae R, De Maria L, Wodak SJ: Assessment of blind predictions of protein-protein interactions: current status of docking methods. Proteins. 2003, 52: 51-67. 10.1002/prot.10393.View ArticlePubMedGoogle Scholar
- Rodrigues JP, Trellet M, Schmitz C, Kastritis P, Karaca E, Melquiond AS, Bonvin AM: Clustering biomolecular complexes by residue contacts similarity. Proteins. 2012, 80: 1810-1817.PubMedGoogle Scholar
- Kastritis PL, Bonvin AMJJ: Are Scoring Functions in Protein-Protein Docking Ready To Predict Interactomes? Clues from a Novel Binding Affinity Benchmark. Journal of Proteome Research. 2010, 9: 2216-2225. 10.1021/pr9009854.View ArticlePubMedGoogle Scholar
- Ko TP, Liao CC, Ku WY, Chak KF, Yuan HS: The crystal structure of the DNase domain of colicin E7 in complex with its inhibitor Im7 protein. Structure with Folding & Design. 1999, 7: 91-102.View ArticleGoogle Scholar
- Wang JH, Smolyar A, Tan K, Liu JH, Kim M, Sun ZY, Wagner G, Reinherz EL: Structure of a heterophilic adhesion complex between the human CD2 and CD58 (LFA-3) counterreceptors. Cell. 1999, 97: 791-803. 10.1016/S0092-8674(00)80790-4.View ArticlePubMedGoogle Scholar
- Schamberger GP, Diez-Gonzalez F: Selection of recently isolated colicinogenic Escherichia coli strains inhibitory to Escherichia coli O157:H7. Journal of food protection. 2002, 65: 1381-1387.PubMedGoogle Scholar
- Murinda SE, Rashid KA, Roberts RF: In vitro assessment of the cytotoxicity of nisin, pediocin, and selected colicins on simian virus 40-transfected human colon and Vero monkey kidney cells with trypan blue staining viability assays. Journal of food protection. 2003, 66: 847-853.PubMedGoogle Scholar
- Kortemme T, Joachimiak LA, Bullock AN, Schuler AD, Stoddard BL, Baker D: Computational redesign of protein-protein interaction specificity. Nature structural & molecular biology. 2004, 11: 371-379. 10.1038/nsmb749.View ArticleGoogle Scholar
- Joachimiak LA, Kortemme T, Stoddard BL, Baker D: Computational design of a new hydrogen bond network and at least a 300-fold specificity switch at a protein-protein interface. Journal of molecular biology. 2006, 361: 195-208. 10.1016/j.jmb.2006.05.022.View ArticlePubMedGoogle Scholar
- Ma JC, Dougherty DA: The Cationminus signpi Interaction. Chem Rev. 1997, 97: 1303-1324. 10.1021/cr9603744.View ArticlePubMedGoogle Scholar
- da Silva AJ, Brickelmaier M, Majeau GR, Li Z, Su L, Hsu YM, Hochman PS: Alefacept, an immunomodulatory recombinant LFA-3/IgG1 fusion protein, induces CD16 signaling and CD2/CD16-dependent apoptosis of CD2(+) cells. J Immunol. 2002, 168: 4462-4471. 10.4049/jimmunol.168.9.4462.View ArticlePubMedGoogle Scholar
- Mrowietz U: Treatment targeted to cell surface epitopes. Clinical and experimental dermatology. 2002, 27: 591-596. 10.1046/j.1365-2230.2002.01171.x.View ArticlePubMedGoogle Scholar
- Spitzer TR, McAfee SL, Dey BR, Colby C, Hope J, Grossberg H, Preffer F, Shaffer J, Alexander SI, Sachs DH: Nonmyeloablative haploidentical stem-cell transplantation using anti-CD2 monoclonal antibody (MEDI-507)- based conditioning for refractory hematologic malignancies. Transplantation. 2003, 75: 1748-1751. 10.1097/01.TP.0000064211.23536.AD.View ArticlePubMedGoogle Scholar
- Braun J, Sieper J: Role of novel biological therapies in psoriatic arthritis: effects on joints and skin. BioDrugs : clinical immunotherapeutics, biopharmaceuticals and gene therapy. 2003, 17: 187-199. 10.2165/00063030-200317030-00005.View ArticleGoogle Scholar
- Weger W: Current status and new developments in the treatment of psoriasis and psoriatic arthritis with biological agents. British journal of pharmacology. 2010, 160: 810-820. 10.1111/j.1476-5381.2010.00702.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Gokhale A, Weldeghiorghis TK, Taneja V, Satyanarayanajois SD: Conformationally constrained peptides from CD2 to modulate protein- protein interactions between CD2 and CD58. Journal of medicinal chemistry. 2011, 54: 5307-5319. 10.1021/jm200004e.PubMed CentralView ArticlePubMedGoogle Scholar
- Arulanandam AR, Moingeon P, Concino MF, Recny MA, Kato K, Yagita H, Koyasu S, Reinherz EL: A soluble multimeric recombinant CD2 protein identifies CD48 as a low affinity ligand for human CD2: divergence of CD2 ligands during the evolution of humans and mice. J Exp Med. 1993, 177: 1439-1450. 10.1084/jem.177.5.1439.View ArticlePubMedGoogle Scholar
- Dustin ML: Adhesive bond dynamics in contacts between T lymphocytes and glass-supported planar bilayers reconstituted with the immunoglobulin- related adhesion molecule CD58. The Journal of biological chemistry. 1997, 272: 15782-15788. 10.1074/jbc.272.25.15782.View ArticlePubMedGoogle Scholar
- Recny MA, Neidhardt EA, Sayre PH, Ciardelli TL, Reinherz EL: Structural and functional characterization of the CD2 immunoadhesion domain. Evidence for inclusion of CD2 in an alpha-beta protein folding class. The Journal of biological chemistry. 1990, 265: 8542-8549.PubMedGoogle Scholar
- Sayre PH, Hussey RE, Chang HC, Ciardelli TL, Reinherz EL: Structural and binding analysis of a two domain extracellular CD2 molecule. J Exp Med. 1989, 169: 995-1009. 10.1084/jem.169.3.995.View ArticlePubMedGoogle Scholar
- Davis SJ, Davies EA, Barclay AN, Daenke S, Bodian DL, Jones EY, Stuart DI, Butters TD, Dwek RA, van der Merwe PA: Ligand binding by the immunoglobulin superfamily recognition molecule CD2 is glycosylation- independent. The Journal of biological chemistry. 1995, 270: 369-375. 10.1074/jbc.270.1.369.View ArticlePubMedGoogle Scholar
- Davis SJ, Ikemizu S, Wild MK, van der Merwe PA: CD2 and the nature of protein interactions mediating cell-cell recognition. Immunol Rev. 1998, 163: 217-236. 10.1111/j.1600-065X.1998.tb01199.x.View ArticlePubMedGoogle Scholar
- Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, Dror RO, Shaw DE: Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins. 2010, 78: 1950-1958.PubMed CentralPubMedGoogle Scholar
- Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: Comparison of Simple Potential Functions for Simulating Liquid Water. Journal of Chemical Physics. 1983, 79: 926-935. 10.1063/1.445869.View ArticleGoogle Scholar
- Jorgensen WL, Duffy EM, Tiradorives J: Computational Investigations of Protein Denaturation - Apomyoglobin and Chaotrope-Arene Interactions. Philos T R Soc A. 1993, 345: 87-96. 10.1098/rsta.1993.0120.View ArticleGoogle Scholar
- Darden T, York D, Pedersen L: Particle Mesh Ewald - an N.Log(N) Method for Ewald Sums in Large Systems. Journal of Chemical Physics. 1993, 98: 10089-10092. 10.1063/1.464397.View ArticleGoogle Scholar
- Hess B, Bekker H, Berendsen HJC, Fraaije JGEM: LINCS: A linear constraint solver for molecular simulations. Journal of computational chemistry. 1997, 18: 1463-1472. 10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-H.View ArticleGoogle Scholar
- Bussi G, Donadio D, Parrinello M: Canonical sampling through velocity rescaling. The Journal of chemical physics. 2007, 126: 014101-10.1063/1.2408420.View ArticlePubMedGoogle Scholar
- Parrinello M, Rahman A: Polymorphic Transitions in Single-Crystals - a New Molecular-Dynamics Method. J Appl Phys. 1981, 52: 7182-7190. 10.1063/1.328693.View ArticleGoogle Scholar
- de Groot BL, van Aalten DM, Amadei A, Berendsen HJ: The consistency of large concerted motions in proteins in molecular dynamics simulations. Biophysical journal. 1996, 71: 1707-1713. 10.1016/S0006-3495(96)79372-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Amadei A, Ceruso MA, Di Nola A: On the convergence of the conformational coordinates basis set obtained by the essential dynamics analysis of proteins' molecular dynamics simulations. Proteins. 1999, 36: 419-424. 10.1002/(SICI)1097-0134(19990901)36:4<419::AID-PROT5>3.0.CO;2-U.View ArticlePubMedGoogle Scholar
- Hubbard SJ, Thornton JM: 'NACCESS' Computer Program. Department of Biochemistry and Molecular Biology. 1993, University College LondonGoogle 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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.