Conformational and functional analysis of molecular dynamics trajectories by Self-Organising Maps
- Domenico Fraccalvieri^{1},
- Alessandro Pandini^{2},
- Fabio Stella^{3} and
- Laura Bonati^{1}Email author
https://doi.org/10.1186/1471-2105-12-158
© Fraccalvieri et al; licensee BioMed Central Ltd. 2011
Received: 30 August 2010
Accepted: 14 May 2011
Published: 14 May 2011
Abstract
Background
Molecular dynamics (MD) simulations are powerful tools to investigate the conformational dynamics of proteins that is often a critical element of their function. Identification of functionally relevant conformations is generally done clustering the large ensemble of structures that are generated. Recently, Self-Organising Maps (SOMs) were reported performing more accurately and providing more consistent results than traditional clustering algorithms in various data mining problems. We present a novel strategy to analyse and compare conformational ensembles of protein domains using a two-level approach that combines SOMs and hierarchical clustering.
Results
The conformational dynamics of the α-spectrin SH3 protein domain and six single mutants were analysed by MD simulations. The Cα's Cartesian coordinates of conformations sampled in the essential space were used as input data vectors for SOM training, then complete linkage clustering was performed on the SOM prototype vectors. A specific protocol to optimize a SOM for structural ensembles was proposed: the optimal SOM was selected by means of a Taguchi experimental design plan applied to different data sets, and the optimal sampling rate of the MD trajectory was selected. The proposed two-level approach was applied to single trajectories of the SH3 domain independently as well as to groups of them at the same time. The results demonstrated the potential of this approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space in a simple 2D visualisation, as well as of effectively highlighting differences in the conformational dynamics directly related to biological functions.
Conclusions
The use of a two-level approach combining SOMs and hierarchical clustering for conformational analysis of structural ensembles of proteins was proposed. It can easily be extended to other study cases and to conformational ensembles from other sources.
Keywords
Background
Protein dynamics plays a central role in cell life. In many cases biological function involves molecular motion [1] and it was recently suggested that intrinsic dynamics also defines the ability of proteins to adapt and evolve new functions [2]. Therefore, a full understanding of protein function and evolution will require a deeper insight into biomolecular atomistic dynamics.
Significant contributions in this direction have come from computational methods, in particular from Molecular dynamics (MD) simulations [3, 4], by which a large ensemble of molecular structures can be generated to sample the accessible conformational space of a protein. Analysis of this ensemble can provide information about average physico-chemical and geometrical properties, as well as allowing identification of recurring conformations and transitions between them.
As recurrence of and transition between conformations is difficult to extract from the raw ensemble data [5], grouping the conformations becomes a necessity. The most desirable strategy would be to use kinetic clustering[6–11] where conformations are grouped according to their transition probabilities during the simulation and identified clusters are directly related to the free energy landscape. A limitation to this approach arises from the need of an exhaustive sampling with convergence of all pair-wise transition probabilities [11]. A more affordable solution is to use geometrical clustering, because only a representative sampling of the accessible conformations is required. The underlying assumption is that structurally similar conformations lie in the same basin of the free energy surface. While often this is an acceptable approximation, a recent study has suggested caution in interpreting the clustering results [11].
Geometrical clustering for conformational analysis was introduced when simulation time increased up to nanoseconds generating tens of thousands of structures [12–14] and has been extensively used since then [5, 11]. Several data-mining algorithms have been adopted but according to a recent survey [5] no general strategy is available: clustering results are often influenced by the type of algorithm and the choice of optimal parameters is mostly left to user experience and the specific case. Originally adapted to analyse protein folding simulations [12], these algorithms have mainly been implemented for multiple trajectories of the same system. However, recently a great interest has emerged in the comparison of protein flexibility of functionally related proteins as well as in studying the evolutionary conservation and specialization of protein dynamics across distant homologous proteins [2, 15–21]. This new interest emphasises the need for more advanced tools to compare conformational ensembles of different protein domains especially when derived from extensive MD simulations [15, 16, 18, 20].
Among the data-mining algorithms, Self-Organising Maps (SOMs) [22] were recently applied to conformational analysis of lipid molecules [23, 24] and automatic clustering of protein-ligand docking poses [25]. In the context of MD analysis, a recent study [5] compared several methods, including SOMs and traditional clustering algorithms. While the SOMs were identified among the best performing methods, no algorithm emerged as the optimal solution.
In this contribution, we present the use of a two-level approach [26] combining SOMs and hierarchical clustering [27] for the analysis and comparison of multiple MD trajectories of a protein domain. First we illustrate the protocol we developed to characterize and optimize the parameters of the SOM learning process for structural ensembles. Then we present the application of the proposed two-level approach, with the optimized parameters, to the conformational and functional analysis of a test case composed of the α-spectrin SH3 (Spc-SH3) protein domain and a group of its single-site mutants. This is an interesting study case of a small intra-cellular signaling domain [28] where ligand binding activity is modulated by single-mutations greatly affecting the conformational dynamics [29, 30]. This test demonstrated the potential of the proposed approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space embedded in a simple 2D visualisation, given by the SOM stage, as well as the capability of effectively highlighting differences in the conformational dynamics directly related to biological functions, given by the combination of SOMs and the post-clustering stage.
Methods
Conformational sampling and ensemble generation
Molecular Dynamics Simulations
The atomistic dynamics of the set of representative proteins was simulated using the GROMACS package (version 3.3.3) [31–33] with the GROMOS96 43a2 force field. All structures were inserted into an octahedral box with explicit solvent and simulated with periodic boundary conditions. Water molecules were described by a simple point charge (SPC) model [34] and the box size was set to ensure a distance of at least 1.2 nm between the protein and the box boundaries. The solvent was relaxed with a 5 ps MD simulation, then the systems were neutralized by insertion of counter ions, and a short minimization with steepest descent was performed up to convergence on maximum force lower than 1000 kJ/(mol*nm). The resulting systems were simulated for 40 ns in the NPT ensemble (constant number of atoms N, pressure P, and temperature T). Long-range electrostatic interactions were calculated with the particle mesh Ewald (PME) summation method [35, 36], with a 9 Å cut-off for the direct space s ums, a 1.2 Å FFT grid spacing and a 4-order interpolation polynomial for the reciprocal space sums. A thermal bath was independently coupled with protein and solvent using a Berendsen thermostat at 300 K with coupling period of 0.1 ps. The internal degrees of freedom of water were constrained by the Settle algorithm [37], while all bond distances in the protein were constrained by the LINCS algorithm [38]. The integration step was set to 2 fs.
where tr is the trace of the matrix. When s(A,B) is equal to 1, the two spanned subspaces are identical, whereas a value of 0 indicates complete orthogonality. In this work, the overlap between the covariance matrix of the atomic coordinates in each half of a simulation and in the overall trajectory [40] was evaluated for each protein.
Analysis of conformational flexibility
where r_{ i }is the coordinate vector for atom i and the < > indicates a time-averaging. To increase the signal/noise ratio, the RMSF calculation was performed after Essential Dynamics (ED) analysis. ED is a widely applied technique based on principal component analysis (PCA) of conformational ensembles [41], aimed to extract informative directions of motion. Only Cα atoms were included in the analysis [41].
where d is the distance value, i and j are the indices of the selected atoms and N is the total number of computed distances.
Optimization of the Self-Organising Maps for structural ensemble analysis
In the proposed approach the SOM input data are large ensembles of conformations extracted from the MD trajectories of one or more domains. Each conformation is described by the Cartesian coordinates of the Cα atoms of the protein, after projection of the original data on the eigenvectors that define the Essential Space. At the end of the SOM training process (see "Brief introduction to SOMs" for details) the original conformational ensemble ends up projected on a bi-dimensional feature map where similar elements are associated to neighbour neurons. Therefore the whole input space is described by a limited number of vectors (the prototype vectors of the neurons). These are then submitted to linkage clustering and the representative conformations are extracted for conformational and functional analysis.
A optimization of the SOM parameters for the specific type of data is required. The optimization we designed and performed consists of two steps: first, the SOM parameters were characterized and optimized by experimental design to obtain the optimal SOM and then these were used to select the optimal sampling rate of the MD trajectory.
Brief introduction to SOMs
A SOM is a powerful tool for the grouping and visualization of high-dimensional data. It is a specific architecture of artificial neural networks consisting of a low-dimensional (usually bi-dimensional) grid of so-called neurons, called map [42]. The number of neurons may vary from few units up to thousands and each neuron is represented by a d-dimensional model vector (the prototype vector), where d is the dimension of the input data vectors (in this application d is the number of Cα coordinates). The neurons are connected to adjacent neurons by a neighborhood relation, which dictates the topology of the SOM map.
The SOM learning algorithm could be divided into five steps: a) initial randomization of each prototype vector; b) random selection of an input vector; c) neurons "competition" to "win" the input vector; d) update of the prototype vector of the "winning" neuron and of its topological neighbors by pulling them closer to the "won" input vector; e) repeat from b) for each input vector for a given number of epochs sufficient to ensure a convergent training of the map. As a result, the prototype vector of each output neuron represents a particular feature (conformation) drawn from the input space [42].
A useful tool to analyse the results is a pictorial representation of the trained map. In this work we choose a topographic map where each neuron is represented by an hexagon. To visualize the number of input data won by each neuron (hits), black hexagons with size proportional to the neuron population were superimposed to the map. Visual inspection of the map allows the analysis of the final distribution of the original data: neighbor neurons (hexagons) are associated to similar groups of input data (conformations); map regions with highly-populated neurons represent conformations highly sampled during the trajectory.
The main interpretational advantage over hierarchical clustering is that maps are more easily understood and interpreted by humans than dendrograms, which become already very complex in the case of a few hundred data points.
In this work, the analysis was performed using the SOM toolbox in MATLAB [43], but in principle it can be reproduced using other statistical packages, such as R [44]. The calculations were performed on an AMD Opteron 254 1.8 GHz.
Experimental design and map optimization
The result of a SOM learning process depends on a set of parameters and the problem of finding the optimal SOM, i.e. the SOM model which "best characterizes" the underlying input space structure, can be extremely complex. The goal is first to discover the parameters that significantly affect the performance of the SOM model, and then to find the values of such parameters that define the optimal SOM. Unfortunately this is an ill-posed optimization problem [45], so several efforts have been made to provide an efficient and effective solution [46, 47].
where d represents Euclidean distance, j and C_{j} are the j^{ th } cluster of neurons and the corresponding set of neurons, w_{ ij }is the weight vector associated with the i^{ th } neuron of the j^{ th } cluster, μ_{ j }and μ are the centroid of the j^{ th } cluster and of the overall map. The performance measure (E4) is a modification of the clustering efficiency index described in [48] where measures the average cluster tightness while measures the separation of centroids. However, the term scales with the number of output neurons and thus it was normalized through R_{ tot } to account for different SOM sizes. The performance measure (E4) has a unique minimum and was successful against expert-based assessment of clustering [48].
The problem of finding the optimal value of the performance measure is computationally expensive and Artificial Neural Networks are generally affected by the curse of dimensionality. This problem can be solved following different strategies [49], among which Design of Experiments or Experimental Design [50, 51] was selected in this study. In this technique the experiments are designed to have rational relationships to the purposes, needs and physical limitations of the problem. Variables and their values are tested in combinations that ensure to answer the questions of interest as clearly and efficiently as possible. This offers advantages in the economy of the experimentation and provides straightforward estimates of experimental effects and variance [52]. In the present study a set of SOM models with different values of the parameters were selected according to an experimental design plan and the output was used to identify the SOM's most influential parameters and thus to select the optimal SOM model (see Figure 1).
SOM parameters together with allowed ranges.
SOM parameter | Range |
---|---|
Map size | [100,400] |
Lattice type | {hexagonal, rectangular} |
Shape | {sheet, cylinder, toroid} |
Learning algorithm | {batch, sequential} |
Neighbour function | {gaussian, bubble, ep} |
Alpha type | {inverse, linear, power} |
Radius | {1, 2, 3} |
Training length | [1000,5000] |
Starting alpha | [0.01, 0.09] |
Taguchi design plan.
Map size | Lattice type | Shape | Learning algorithm | Neighbour function | Alpha type | Radius | Training length | Starting alpha | Time (s) |
---|---|---|---|---|---|---|---|---|---|
100 | hexagonal | sheet | batch | Gaussian | inverse | 2 | 1000 | 0.01 | 20 |
225 | hexagonal | cylinder | batch | Bubble | power | 3 | 3000 | 0.05 | 125 |
400 | hexagonal | toroid | batch | Ep | linear | 1 | 5000 | 0.09 | 450 |
100 | hexagonal | sheet | sequential | Gaussian | linear | 3 | 1000 | 0.09 | 270 |
225 | hexagonal | cylinder | sequential | Bubble | inverse | 1 | 3000 | 0.01 | 820 |
400 | hexagonal | toroid | sequential | Ep | power | 2 | 5000 | 0.05 | 12000 |
100 | rectangular | cylinder | sequential | Gaussian | linear | 1 | 5000 | 0.05 | 1350 |
225 | rectangular | toroid | sequential | Bubble | inverse | 2 | 1000 | 0.09 | 600 |
400 | rectangular | sheet | sequential | Ep | power | 3 | 3000 | 0.01 | 7700 |
100 | rectangular | toroid | batch | Gaussian | power | 1 | 3000 | 0.01 | 55 |
225 | rectangular | sheet | batch | Bubble | linear | 2 | 5000 | 0.05 | 205 |
400 | rectangular | cylinder | batch | Ep | inverse | 3 | 1000 | 0.09 | 95 |
100 | hexagonal | toroid | sequential | Bubble | power | 1 | 1000 | 0.05 | 290 |
225 | hexagonal | sheet | sequential | Ep | linear | 2 | 3000 | 0.09 | 1800 |
400 | hexagonal | cylinder | sequential | Gaussian | inverse | 3 | 5000 | 0.01 | 12800 |
100 | rectangular | toroid | sequential | Bubble | inverse | 3 | 3000 | 0.09 | 600 |
225 | rectangular | sheet | sequential | Ep | power | 1 | 5000 | 0.01 | 2500 |
400 | rectangular | cylinder | sequential | Gaussian | linear | 2 | 1000 | 0.05 | 2800 |
100 | rectangular | sheet | sequential | Bubble | linear | 3 | 5000 | 0.01 | 1400 |
225 | rectangular | cylinder | sequential | Ep | inv | 1 | 1000 | 0.05 | 650 |
400 | rectangular | toroid | sequential | Gaussian | power | 2 | 3000 | 0.09 | 8500 |
100 | rectangular | cylinder | batch | Bubble | power | 2 | 5000 | 0.09 | 85 |
225 | rectangular | toroid | batch | Ep | linear | 3 | 1000 | 0.01 | 40 |
400 | rectangular | sheet | batch | Gaussian | inverse | 1 | 3000 | 0.05 | 300 |
100 | hexagonal | cylinder | batch | Ep | power | 2 | 1000 | 0.01 | 17 |
225 | hexagonal | toroid | batch | Gaussian | linear | 3 | 3000 | 0.05 | 140 |
400 | hexagonal | sheet | batch | Bubble | inverse | 1 | 5000 | 0.09 | 450 |
100 | hexagonal | cylinder | batch | Ep | linear | 1 | 3000 | 0.09 | 50 |
225 | hexagonal | toroid | batch | Gaussian | inverse | 2 | 5000 | 0.01 | 230 |
400 | hexagonal | sheet | batch | Bubble | power | 3 | 1000 | 0.05 | 90 |
100 | hexagonal | toroid | sequential | Ep | inverse | 3 | 5000 | 0.05 | 1400 |
225 | hexagonal | sheet | sequential | Gaussian | power | 1 | 1000 | 0.09 | 630 |
400 | hexagonal | cylinder | sequential | Bubble | linear | 2 | 3000 | 0.01 | 7600 |
100 | rectangular | sheet | batch | Ep | inverse | 2 | 3000 | 0.05 | 49 |
225 | rectangular | cylinder | batch | Gaussian | power | 3 | 5000 | 0.09 | 236 |
400 | rectangular | toroid | batch | Bubble | linear | 1 | 1000 | 0.01 | 80 |
The results of the numerical experiments associated with the design plan are used to find a linear regression model of the performance measure (E4) as function of the SOM's statistically significant parameters. The linear regression is performed in two stages: a stepwise regression [54] to select the statistically significant parameters and a linear regression fitting only on these to actually model the unknown mapping between them and the performance measure (E4). Stepwise regression parameters were: Prob to Enter = 0.05, Prob to Leave = 0.05 and Direction = Mixed and Rules = No Rules. The linear regression model is validated by comparing actual and predicted performance. JMP software was used for data analysis and linear regression fitting [55]. On the basis of this model the optimal values of the SOM parameters were derived.
Optimal sampling rate of MD trajectories
When dealing with a nanoseconds MD trajectory, it is desirable to train the SOM with a minimum number of selected structures while maintaining a reliable picture of the protein dynamics. Therefore, the effect of the sampling rate on the SOM learning process was assessed and the minimum number of frames to extract from a trajectory, without a statistically significant loss of information, was identified.
In this study such a test was performed on three replicas by randomly selecting n conformations of the trajectory, each belonging to one of n equal width trajectory intervals. To better clarify the statistical procedure, let Q_{ i } ^{ (n) } be the i^{ th } random sample extracted from Q and SOM_{ i } ^{ (n) } be the corresponding SOM trained on it by using the optimal parameters determined by experimental design. Furthermore, let h_{ i } ^{ (n) } and H_{ i } ^{ (n) } be the SOM hit vectors in the case where Q_{ i } ^{ (n) } and Q are presented to SOM_{ i } ^{ (n) } . Then, the χ_{ i } ^{ (n) } Chi-square statistic is computed on non-null cells of the SOM hit vectors H_{ i } ^{ (n) } and h_{ i } ^{ (n) } while the null hypothesis is tested by the n Cressie-Read Chi-square test where the significance level has been set to α = 0.05 [56]. A diagram of the overall procedure is shown in Figure 2.
Cluster analysis of the SOM prototype vectors
The optimal SOM resulting from the previous optimization steps summarizes the conformational variability in the input data (structural ensembles generated by MD) through a set of prototype vectors. The next step of our approach consists on combining the prototype vectors (that can be interpreted as "protoclusters" [26]) in a reduced number of final clusters, where each original conformation belongs to the same cluster as its nearest prototype. To this end, two hierarchical agglomerative clustering algorithms, the complete linkage and the average linkage, were tested, following what described in [26].
where sd_{ d } is the standard deviation, is the mean, and z is a specified constant.
A range of values for z were tested: 1.25, 1.50, 1.75, 2.00, 2.25, 2.50, and 2.75 for both complete and average linkage. Different values of the optimal number of clusters were obtained: from 5 to 8 for the complete linkage and from 4 to 10 for the average linkage. The most informative cluster separation was obtained with complete linkage and z = 2.75, and with average linkage and z = 2.50, which agrees with that suggested in [58].
The stability and better performance of the proposed two-level approach, i.e. SOM and linkage, with respect to other clustering methods, as well as the relative performances of the complete vs. average linkage algorithms were verified by the analysis reported in the Methods and Results sections on "Comparison of cluster analysis methods".
Hierarchical clustering was performed using the MATLAB Statistics Toolbox [59] and the Mojena's rule was applied in its MATLAB implementation, made available by courtesy of Prof. Fernández http://ima.udg.edu/~jamf/MARTIN_MF.htm.
Comparison of cluster analysis methods
The proposed two-level approach for clustering structural ensembles (combining SOM with the average or the complete linkage methods) was compared with other clustering methods applied to the original input data vectors: GROMOS nearest neighbour [60], average and complete linkage methods [27]. GROMOS algorithm is implemented in the GROMACS package [31–33]: neighbours of each data point are defined according to a cutoff distance c; the point with largest neighbourhood defines the first cluster medoid; this point and its neighbours are removed and the algorithm is iterated until all data have been assigned to a cluster.
The quality of the clustering was assessed by the Silhouette (S) index [61] and the Davies-Bouldin (DB) index [62]. The two indices score the compactness and separations of the clusters, and were previously used in similar studies [5, 26]. The S index can also be calculated per data point or cluster. Its value is in the range [-1, 1] and is generally interpreted as evidence in support of the cluster structure [63]: strong (S > 0.7), reasonable (0.5 < S ≤ 0.7), weak (0.25 < S ≤ 0.5) or no significant evidence (S ≤ 0.25). The DB is a measure of relative scattering of the data points. A lower value of DB indicates more compact and separated clusters.
The performance of the methods was compared for the number of clusters ranging from 2 to 15. In the case of GROMOS algorithm, solutions were generated for c ranging from 0.01 to 0.40 with increments of 0.005. Then the best solution for each number of clusters was selected according to the relevant quality index.
Results
Molecular Dynamics of the SH3 domain
The six mutants of the Spc-SH3 domain analysed here (and their PDB codes) are: R21A, R21G, A56G and A56S (2F2W, 2F2X, 2F2V and 2CDT) [30], N47A and N47G (1QKX, 1QKW) [67]. These mutants had been designed to explore how the local perturbations produced by single-point mutations affect both the stability and the p41 affinity of the domain [30].
In particular, it was observed that mutation of R21, at the tip of the RT loop that flanks the binding specificity pocket (Figure 3A), to alanine had favourable effects on the p41 affinity, probably due to the replacement of the bulky arginine side chain by a small side chain. On the contrary, a slight reduction in affinity was observed for the R21G mutant. Both mutations of N47, in the distal loop (Figure 3A), produced significant changes in the stability of the domain and a reduction of p41 affinity. Finally, while mutation of A56, in the 3_{10} helix belonging to one of the hydrophobic binding grooves (Figure 3A), to Ser did not alter significantly the p41 affinity, the mutation to Gly produced a great reduction in affinity.
The conformational dynamics of the wild type structure and the six mutated Spc-SH3 domains were analysed by 40 ns MD simulations (see the Methods section for details). Analysis of the root mean square distance (RMSD) to the starting structure confirmed a general stability of all trajectories, with an equilibration time around 1 ns (see Additional File 1) and a temperature around 300 K during the simulation. The covariance matrix overlap between the conformational spaces spanned by the two halves of the simulation and the overall trajectory for all the SH3 systems (Additional File 2) shows values in the range 0.70 - 0.87, confirming that the conformational subspaces explored in each half are representative of the complete space. Similar values were previously reported in cases of satisfactory sampling [15, 39, 68].
The most informative directions of motion were extracted by Essential Dynamics analysis. Analysis of the fraction of total motion described by different subspaces (Additional File 3) indicates that in all cases more than 79% is described by the first 20 eigendirections (corresponding to 12% of the total space), while the first 30 directions (18% of total space) can explain more than 85% (in most cases more than 90%) of the conformational flexibility of all domains. Therefore, the 30 dimensional essential space of each system was selected for the following analyses. In this way, the noise given by not informative high-frequency motion was removed.
The choice of projecting the trajectories in the essential space led to an increased signal to noise ratio in the plots of the RMSF on the positions of the Cα atoms (E2). A set of comparative RMSF plots is shown in Additional File 4: the WT SH3 structure is rather constrained, with only a region of flexibility in the n-src and the distal loops (Additional File 4a), while the effects of single-point mutations in each of the three positions, R21, N47 and A56 (see Figure 3A), can be clearly observed. While R21A mutation does not affect the overall domain flexibility, the R21G mutation causes a relevant increase of flexibility in the long RT loop and in the region including β3, the distal loop, β4 and the helix 3_{10} (Additional File 4b). These differences suggest that flexibility may have a role in the different binding affinities of the two mutants, that could not be explained only on the basis of the very similar stereo-electronic characteristics of their binding sites [30]. The mutation to Ala at position 47 has little effect on the flexibility, while a change to Gly in the same position causes an increased flexibility in the surrounding region as well as in the whole RT loop (Additional File 4c). This confirms the hypothesis of a long-range propagation of the local perturbation to the binding site, that may have a role in the N47G binding to p41 [69]. Different from previous cases, a mutation to Gly in position 56, in the 3_{10} helix, does not alter either the local or the global flexibility of the domain; an increase of flexibility in the distal loop region is observed for the A56S mutant (Additional File 4d). From these results, the conformational freedom of these two mutants does not seem to be related to the p41 binding affinity [30].
SOM optimization to analyse the SH3 domain dynamics
Following the protocol described in the Methods section, a SOM was optimized for the analysis of the MD trajectories of the WT SH3 and its six mutants. Each sampled conformation was described by the Cartesian coordinates of the Cα atoms projected on the Essential Space. The 55 Cα atoms corresponding to structurally equivalent positions in all the domains were included in the analysis; therefore, the input data presented to the SOM were vectors of 165 elements.
The optimal values for the SOM parameters were found by experimental design using a Taguchi design plan (see Table 2) which consists of 36 runs. Each run was replicated three times on four data sets:
○ WT
○ R21G
○ WT and R21G
○ ALL (WT and the six mutants),
thus resulting in a total of 432 experiments. The choice of multiple and combined data sets aimed at optimizing the map both for the analysis of a single trajectory and for multiple comparisons of different trajectories.
Linear regression, summary of fit.
Parameter | Value |
---|---|
R^{2} | 0.937 |
R^{2}_{adj} | 0.936 |
Mean of Response | 395.01 |
RMSE (Root Mean Square Error) | 27.61 |
Observations | 432 |
Linear regression, effect tests [55].
Source | Nparm | DF | Sum of Squares | F Ratio | Prob > F |
---|---|---|---|---|---|
Map size | 1 | 1 | 1022810.9 | 1341.853 | <0.0001 |
Radius | 1 | 1 | 39779.1 | 52.187 | <0.0001 |
Training length | 1 | 1 | 17056.8 | 22.377 | <0.0001 |
Neighbour function | 2 | 2 | 2935552.0 | 1925.615 | <0.0001 |
TYPEMOL | 3 | 3 | 805604.4 | 352.298 | <0.0001 |
Optimal SOM parameter values.
Case | Map size | Radius | Training length | Neighbour function | Actual | Predicted | Time (s) |
---|---|---|---|---|---|---|---|
WT | 100 | 3 | 5000 | gaussian | 158 | 156 | 50 |
R21G | 100 | 3 | 5000 | gaussian | 265 | 270 | 50 |
WT+R21G | 100 | 3 | 5000 | gaussian | 259 | 238 | 105 |
ALL | 100 | 3 | 5000 | gaussian | 254 | 246 | 340 |
Test on the effects of different sampling rates for the WT SH3 and the R21G and N47A mutants.
statistic | p-value | result | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
number of conformations | sample | WT | R21G | N47A | WT | R21G | N47A | WT | R21G | N47A |
20000 | 1 | 25 | 27 | 16 | 1.00 | 1.00 | 1.00 | - | - | - |
2 | 17 | 28 | 17 | 1.00 | 1.00 | 1.00 | - | - | - | |
3 | 19 | 29 | 22 | 1.00 | 1.00 | 1.00 | - | - | - | |
10000 | 1 | 29 | 35 | 34 | 1.00 | 1.00 | 0.98 | - | - | - |
2 | 28 | 37 | 36 | 1.00 | 1.00 | 0.98 | - | - | - | |
3 | 33 | 36 | 21 | 0.99 | 1.00 | 1.00 | - | - | - | |
5000 | 1 | 30 | 33 | 30 | 1.00 | 1.00 | 0.98 | - | - | - |
2 | 49 | 49 | 28 | 0.91 | 0.98 | 0.98 | - | - | - | |
3 | 34 | 37 | 30 | 0.99 | 1.00 | 0.99 | - | - | - | |
2500 | 1 | 63 | 53 | 23 | 0.26 | 0.92 | 0.99 | - | - | - |
2 | 51 | 48 | 56 | 0.76 | 0.99 | 0.22 | - | - | - | |
3 | 39 | 44 | 34 | 0.92 | 0.99 | 0.98 | - | - | - | |
1600 | 1 | 43 | 49 | 31 | 0.82 | 0.93 | 0.96 | - | - | - |
2 | 42 | 63 | 28 | 0.86 | 0.74 | 0.95 | - | - | - | |
3 | 45 | 48 | 41 | 0.90 | 0.97 | 0.77 | - | - | - | |
800 | 1 | 40 | 57 | 73 | 0.92 | 0.89 | 0.06 | - | - | - |
2 | 57 | 43 | 49 | 0.52 | 0.99 | 0.47 | - | - | - | |
3 | 54 | 44 | 36 | 0.48 | 0.98 | 0.80 | - | - | - | |
400 | 1 | 48 | 73 | 39 | 0.63 | 0.11 | 0.49 | - | - | - |
2 | 50 | 49 | 55 | 0.55 | 0.93 | 0.16 | - | - | - | |
3 | 34 | 39 | 46 | 0.96 | 0.99 | 0.47 | - | - | - | |
80 | 1 | 40 | 83 | 54 | 0.28 | 0.00 | 0.02 | - | * | * |
2 | 63 | 64 | 38 | 0.00 | 0.00 | 0.08 | * | * | - | |
3 | 81 | 46 | 50 | 0.00 | 0.43 | 0.04 | * | - | * | |
40 | 1 | 167 | 84 | 67 | 0.00 | 0.00 | 0.00 | * | * | * |
2 | 105 | 62 | 95 | 0.00 | 0.00 | 0.00 | * | * | * | |
3 | 101 | 60 | 67 | 0.00 | 0.00 | 0.00 | * | * | * |
Comparison of cluster analysis methods
The two-level approach proposed in this work (SOM and linkage) is the best performing: the S (DB) scores are generally higher (lower) for all the number of clusters and the profiles are more stable. This suggests that the SOM preclustering is reducing the noise and generating more compact and well separated clusters. This is also confirmed by the per-cluster values of S (see Additional file 5). According to S, only the two-level approach shows reasonable evidence of correct cluster structure (S > 0.5) for most values of the number of clusters range.
While from the analysis of the two global S and DB indices both the linkage methods appear to be good choices for our two-level approach, a more detailed analysis of the performances on the ALL dataset revealed some differences. In particular, the complete linkage provides clearer and better separated clusters for conformations related to the most characteristic motions of the domain exploited by specific mutants (see the discussion of Figure 7 in the next section). Accordingly, the complete linkage algorithm was used in association with SOM for conformational and functional analysis of the SH3 domains.
The S and DB indices are often employed to identify the optimal number of clusters. In this case they show a rather even profile for the two-level approach, but it has to be noted that the significant values for the SOM and complete linkage approach range from 3 to 8. This is compatible with the optimal value of 5 that was selected by the Mojena's rule for the ALL dataset.
Additional File 6 shows a 2D RMSD matrix of the conformations in the ALL data set, where the lower triangular part reports the best 5 cluster solutions with GROMOS. While the RMSD maps are powerful exploratory tools for single trajectories, they become increasingly complex when several simulations are reported. Additionally the image shows a limitation of GROMOS nearest neighbour approach: a tendency to generate singletons. This is also a frequent event in single linkage clustering [5] that for this reason was not included in the test.
Conformational and functional analysis of the SH3 domains by SOM and complete linkage clustering
The proposed approach (SOM and complete linkage) with the SOM parameters and the sampling rate values set in the optimization stage was used to perform a conformational analysis of the MD trajectories of single SH3 domains independently as well as groups of simulations to detect similarities and differences in flexibilities among the systems. We will present the results for the following ensembles: R21G, WT and R21G, WT and six mutants.
dRMSD values between points A, B, C, D, in Figure 8.
dRMSD (Å) | ||||
---|---|---|---|---|
ABD | ACD | ABC | ABCD | |
Cluster 1 | 1.7 | 0.9 | 1.4 | 1.4 |
Cluster 2 | 1.0 | 1.5 | 1.6 | 1.4 |
Cluster 3 | 2.2 | 2.3 | 2.7 | 2.4 |
Cluster 4 | 2.0 | 2.1 | 1.6 | 2.2 |
Cluster 5 | 3.8 | 2.9 | 2.4 | 3.5 |
A closer look at the cluster composition shows the ability of the SOM to group conformations common to all domains, as well as to correctly separate the typical dynamics of each of the three domains with higher flexibility. The contributions of each mutant ensemble to the five clusters (Figure 7A) highlights that clusters 1 and 2 are populated by conformations from all the mutants. The larger contributions to cluster 1 are from the WT SH3 and the mutants with reduced conformational flexibility (R21A, N47A and A56G), while cluster 2 is more representative of R21G and N47G ensembles. Each of the remaining three clusters is dominated by one contribution: cluster 3 mainly by A56S, cluster 4 by N47G and cluster 5 is almost completely populated by conformations from the R21G ensemble.
An interesting feature arises from the topological nature of the SOM. Conformational transitions that occur in consecutive times along the MD trajectory involve conformations assigned to neighbour clusters on the map. An example for three trajectories is reported in Additional File 7 where the obtained clusters are annotated on the RMSD plots by colours. In the first part of the WT RMSD plot, frequent transitions occur between conformations in the green and blue clusters that are neighbours in the SOM. More clearly, in the A56S plot, transitions between conformations in the green and yellow clusters (that in the SOM are separated by the violet cluster) always occur in the trajectory through sampling of conformations in the violet cluster. In the R21G plot, while some transitions occur between conformations in neighbour clusters (blue and green), others (green to yellow or red to yellow) are separated by the violet cluster in the map and occur only through brief sampling of conformations of this type.
Previous studies suggested a hypothesis of the role of conformational flexibility in reducing the binding affinity [30]. To verify this we studied the effects of flexibility on the binding site geometry. As previously described, the binding pocket of the p41 peptide is flanked by the RT and n-src loops (Figure 3B), whose dynamics affects the binding specificity. Therefore, the inter-residue distances in the binding site of the SH3:p41 complex (all the heavy atoms in the residues interacting with p41, shown in Figure 3B) in the representative conformation of each cluster were compared to those in the experimental structure. The dRMSD values confirmed that the increase of conformational freedom induced by the N47G and R21G mutations (mainly described by clusters 4 and 5) produces a significant distortion of the binding site geometry (dRMSD = 1.7 and 1.8 Å), while other mutations did not produce comparable effects (dRMSD around 1 Å). These results correlate well with the binding affinity data reported for the N47G and R21G mutants [30]. In fact, even if the site of the former mutation is distant from the binding site, a significant decrease in binding affinity is produced. Moreover the R21G mutant that, on the basis of the favorable changes in the binding site stereo-electronic properties, was expected to give a significant increase in binding affinity to p41 (as observed for the R21A mutant), indeed produced a slight reduction in affinity.
Discussion
We have presented a strategy to compare conformational ensembles of protein domains with the goal of highlighting similarities and differences in functional motions. This strategy uses SOMs, that have recently been shown to be suitable for the analysis of individual MD trajectories [5], but takes also advantage of a two-level approach [26] where complete linkage clustering is performed on the SOM prototype vectors. To provide the users with optimal SOM parameters for this specific type of data, we developed a protocol by which we identified a small number of important parameters and calculated their optimal values (Table 5).
When dealing with large data sets of conformations, a major issue is the computational cost of the analysis. A possible strategy to overcome this problem is to compare only the average geometrical properties of subgroups of data, for example the ensemble RMSF on atom positions. However in this analysis small fluctuations are difficult to detect, the direction of motion associated to each peak is not considered, and comparisons are only pair-wise. These limitations appeared clearly in our analysis of the RMSF profiles of the Spc-SH3 domains. On the contrary, the proposed SOM approach retains high sensitivity, is able to differentiate motions with similar average fluctuation and is not restricted to pair-wise comparisons. A second strategy to reduce computational costs involves using geometrical clustering methods in a two-stage or sieved approach [5] by initially clustering only part of the data and then in a second step by adding the missing ones into existing clusters. However, this can lead to the loss or distortion of the topological relations among the original data and eventually to a biased grouping, if the selection at the first stage is not representative. Similarly during the SOM training, each data vector is compared only to the neuron vectors representing all the data already presented to the map, but the topological relations are intrinsically recorded and the representative geometries are dynamically updated avoiding a bias. This comes at increased computational cost in the training stage, but once a map is trained on a group of representative protein domains, it can be used for fast classification of conformational ensembles of similar systems.
In the proposed two-level approach, where complete linkage is applied on the SOM prototype vectors, the number of neurons identified by experimental design (Map size = 100) allows the SOM model to optimally capture the relevant features of a large input space in a smaller set of prototype vectors. These can be interpreted as an intermediate topological reference ("protoclusters" [26]), which allows efficient use of clustering algorithms to divide the prototypes into groups. The number of clusters has not to be specified in advance but it is identified only after hierarchical clustering on the prototype vectors by the Mojena's stopping rule. This generates a significant reduction in computational time and in noise compared to the direct application of linkage algorithms to the input data [26]. This also differs significantly from the approach adopted in a previous study by other authors [5] where only a SOM was used for clustering the original data and the number of clusters had to be selected in advance by the user, i.e. before the SOM learning takes place.
The results for the SH3 domains highlighted the specific advantages of a SOM in conformational and functional analysis. The major benefit is the possibility of providing a topological mapping of the conformational space embedded in a simple 2D visualisation that simplifies the identification of differences in the conformational dynamics (see Figures 5, 6, 7). Moreover, the map can adapt to record differences in both large and small fluctuations, as well as to group conformations associated with different directions of the same motion.
The proposed approach resulted in a very efficient comparison of multiple trajectories: low fluctuations, large concerted motions and intermediate dynamic perturbations were clearly and correctly detected (Figures 7 and 8, and Table 7). The comparison of inter-residue distances in the binding site [65] among the cluster representative conformations led to a functional interpretation of the observed differences. The increase of conformational freedom induced by the N47G and R21G mutations induces a distortion of the binding site geometry that explains the decreased ligand binding ability, while other mutations do not produce comparable effects.
The comparison of the two-level approach to other clustering methods showed that adding a SOM preclustering stage has the benefit of reducing the noise and generating more compact and well separated clusters (see Figure 4 and Additional file 5).
Our approach can be reliably used on other study-cases. Indeed, the number of statistically significant SOM parameters is small while their optimal values are expected to be suitable for conformational spaces with similar complexity, i.e. similar number of unique conformations sampled, to the SH3 study case. In particular, the optimal values of the three parameters associated with the SOM training process (Radius = 3, Training length = 5000 and Neighbour function = Gaussian) are expected to generalize quite well. The Map size parameter is more critical; the optimal value we found (100) allows the summarization of a great number of unique features, but it could be necessary to re-set it specifically for cases with a greater variance in the training data.
The approach presented here is independent from the method used to generate the structural ensemble and is reliable to describe both small and large differences. It is therefore suitable to also analyse combinations of ensembles from computational methods with a more extended sampling of the conformational space and from experiments (NMR ensembles or multiple X-ray depositions of the same structure).
The test-case presented here includes mutants of a single domain that can easily be aligned. A future development of this study is the identification of alternative representations of protein conformations that do not require the preliminary definition of structurally equivalent positions by structural alignment. This will allow an extension to the comparison of different domains, including distant homologous proteins.
Conclusions
The novel approach presented here to analyse conformational ensembles of protein domains resulted in a very efficient way of comparing multiple trajectories of the test case and analysing the role of conformational flexibility in modulating the domain function. The approach can easily be extended to other study cases. Further applications may include the treatment of ensemble of conformations derived from other sources of dynamical data, and the use of a trained map to classify conformational ensembles of similar systems. Future directions involve the extension of the approach to the comparison of protein flexibilities also for distant homologous proteins.
Declarations
Acknowledgements
We like to thank A Fornili, F Fraternali, J Kleinjung, M Sadowski and JW Saldanha for insightful discussions and critical revision of the manuscript. This research was supported by the National Institutes of Health [NIEHS Grant ES007685 to LB and DF], by a Marie Curie Intra European Fellowship within the 7th European Community Framework Programme [PIEF-GA-2008-220256 to AP] and partially by the Medical Research Council [Ref. U117581331 to AP].
Authors’ Affiliations
References
- Henzler-Wildman K, Kern D: Dynamic personalities of proteins. Nature 2007, 450: 964–72. 10.1038/nature06522View ArticlePubMedGoogle Scholar
- Tokuriki N, Tawfik DS: Protein dynamism and evolvability. Science 2009, 324: 203–7. 10.1126/science.1169375View ArticlePubMedGoogle Scholar
- Karplus M, Kuriyan J: Molecular dynamics and protein function. Proc Natl Acad Sci USA 2005, 102: 6679–85. 10.1073/pnas.0408930102PubMed CentralView ArticlePubMedGoogle Scholar
- van Gunsteren WF, Bakowies D, Baron R, Chandrasekhar I, Christen M, Daura X, Gee P, Geerke DP, Glättli A, Hünenberger PH, Kastenholz MA, Oostenbrink C, Schenk M, Trzesniak D, van der Vegt NFA, Yu HB: Biomolecular modeling: Goals, problems, perspectives. Angew Chem Int Ed Engl 2006, 45: 4064–92. 10.1002/anie.200502655View ArticlePubMedGoogle Scholar
- Shao J, Tanner SW, Thompson N, Cheatham TE: Clustering molecular dynamics trajectories: 1. Characterizing the Performance of different clustering algorithms. J Chem Theory Comput 2007, 3: 2312–2334. 10.1021/ct700119mView ArticleGoogle Scholar
- Chodera JD, Singhal N, Pande VS, Dill K, Swope W: Automatic discovery of metastable states for the construction of Markov models of macromolecular conformational dynamics. J Chem Phys 2007, 126: 155101. 10.1063/1.2714538View ArticlePubMedGoogle Scholar
- Noé F, Horenko I, Schütte C, Smith JC: Hierarchical Analysis of Conformational Dynamics in Biomolecules: Transition Networks of Metastable States. J Chem Phys 2007, 126: 155102. 10.1063/1.2714539View ArticlePubMedGoogle Scholar
- Muff S, Caflisch A: Kinetic analysis of molecular dynamics simulations reveals changes in the denatured state and switch of folding pathways upon single-point mutation of a b-sheet miniprotein. Proteins 2008, 70: 1185–1195.View ArticlePubMedGoogle Scholar
- Huisinga W, Best C, Roitzsch R, Schütte C, Cordes F: From Simulation Data to Conformational Ensembles: Structure and Dynamics based Methods. J Comput Chem 1999, 20: 1760–1774. 10.1002/(SICI)1096-987X(199912)20:16<1760::AID-JCC8>3.0.CO;2-2View ArticleGoogle Scholar
- Deuflhard P, Huisinga W, Fischer A, Schutte C: Identification of Almost Invariant Aggregates in Reversible Nearly Uncoupled Markov Chains. Linear Alg Appl 2000, 315: 39–59. 10.1016/S0024-3795(00)00095-1View ArticleGoogle Scholar
- Keller B, Daura X, van Gunsteren WF: Comparing geometric and kinetic cluster algorithms for molecular simulation data. J Chem Phys 2010, 132: 074110. 10.1063/1.3301140View ArticlePubMedGoogle Scholar
- Karpen ME, Tobias DJ, Brooks CL: Statistical clustering techniques for the analysis of long molecular dynamics trajectories: analysis of 2.2-ns trajectories of YPGDV. Biochemistry 1993, 32: 412–420. 10.1021/bi00053a005View ArticlePubMedGoogle Scholar
- Shenkin PS, McDonald DQ: Cluster analysis of molecular conformations. J Comput Chem 1994, 15: 899–916. 10.1002/jcc.540150811View ArticleGoogle Scholar
- Torda AE, van Gunsteren WF: Algorithms for clustering molecular dynamics configurations. J Comput Chem 1994, 15: 1331–1340. 10.1002/jcc.540151203View ArticleGoogle Scholar
- Pandini A, Bonati L: Conservation and specialization in PAS domain dynamics. Protein Eng Des Sel 2005, 18: 127–137. 10.1093/protein/gzi017View ArticlePubMedGoogle Scholar
- Pang A, Arinaminpathy Y, Sansom MS, Biggin PC: Comparative molecular dynamics--similar folds and similar motions? Proteins 2005, 61: 809–22. 10.1002/prot.20672View ArticlePubMedGoogle Scholar
- Cascella M, Micheletti C, Rothlisberger U, Carloni P: Evolutionarily conserved functional mechanics across pepsin-like and retroviral aspartic proteases. J Am Chem Soc 2005, 127: 3734–42. 10.1021/ja044608+View ArticlePubMedGoogle Scholar
- Pandini A, Mauri G, Bordogna A, Bonati L: Detecting similarities among distant homologous proteins by comparison of domain flexibilities. Protein Eng Des Sel 2007, 20: 285–299. 10.1093/protein/gzm021View ArticlePubMedGoogle Scholar
- Maguid S, Fernandez-Alberti S, Echave J: Evolutionary conservation of protein vibrational dynamics. Gene 2008, 422: 7–13. 10.1016/j.gene.2008.06.002View ArticlePubMedGoogle Scholar
- Velázquez-Muriel JA, Rueda M, Cuesta I, Pascual-Montano A, Orozco M, Carazo JM: Comparison of molecular dynamics and superfamily spaces of protein domain deformation. BMC Struct Biol 2009, 9: 6. 10.1186/1472-6807-9-6PubMed CentralView ArticlePubMedGoogle Scholar
- Raimondi F, Orozco M, Fanelli F: Deciphering the deformation modes associated with function retention and specialization in members of the Ras superfamily. Structure 2010, 18: 402–14. 10.1016/j.str.2009.12.015View ArticlePubMedGoogle Scholar
- Kohonen T: The self-organizing map. Proceedings of the Institute of Electrical and Electronics Engineers 1990, 78: 1464–1480.View ArticleGoogle Scholar
- Hyvönen MT, Hiltunen Y, El-Deredy W, Ojala T, Vaara J, Kovanen PT, Ala-Korpela M: Application of self-organizing maps in conformational analysis of lipids. J Am Chem Soc 2001, 123: 810–806. 10.1021/ja0025853View ArticlePubMedGoogle Scholar
- Murtola T, Kupiainen M, Falck E, Vattulainen I: Conformational analysis of lipid molecules by self-organizing maps. J Chem Phys 2007, 126: 054707. 10.1063/1.2429066View ArticlePubMedGoogle Scholar
- Bouvier G, Evrard-Todeschi N, Girault J-P, Bertho G: Automatic clustering of docking poses in virtual screening process using self-organizing map. Bioinformatics 2010, 26: 53–60. 10.1093/bioinformatics/btp623View ArticlePubMedGoogle Scholar
- Vesanto J, Alhoniemi E: Clustering of the Self-Organizing Map. IEEE T Neural Network 2000, 11: 586–600. 10.1109/72.846731View ArticleGoogle Scholar
- Xu R, Wunsch DC II: Clustering. Hoboken, New Jersey, John Wiley and Sons; 2009.Google Scholar
- Li SS: Specificity and versatility of SH3 and other proline-recognition domains: structural basis and implications for cellular signal transduction. Biochem J 2005, 390: 641–653.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang C, Pawley NH, Nicholson LK: The role of backbone motions in ligand binding to the c-Src SH3 domain. J Mol Biol 2001, 313: 873–887. 10.1006/jmbi.2001.5083View ArticlePubMedGoogle Scholar
- Casares S, López-Mayorga O, Vega MC, Cámara-Artigas A, Conejero-Lara F: Cooperative propagation of local stability changes from low-stability and high-stability regions in a SH3 domain. Proteins 2007, 67: 531–547. 10.1002/prot.21284View ArticlePubMedGoogle Scholar
- Berendsen HJC, van der Spoel D, van Drunen R: GROMACS: A message-passing parallel molecular dynamics implementation. Comp Phys Comm 1995, 91: 43–56. 10.1016/0010-4655(95)00042-EView ArticleGoogle Scholar
- Lindahl E, Hess B, van der Spoel D: Gromacs 3.0: A package for molecular simulation and trajectory analysis. J Mol Mod 2001, 7: 306–317.Google Scholar
- van der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJC: GROMACS: Fast, Flexible and Free. J Comput Chem 2005, 26: 1701–1718. 10.1002/jcc.20291View ArticlePubMedGoogle Scholar
- Berendsen HJC, Postma JPM, van Gunsteren WF, Hermans J: Intermolecular Forces. Edited by: Pullman B. Dordrecht, Reidel; 1981:331–342.View ArticleGoogle Scholar
- Darden T, York D, Pedersen L: Particle mesh Ewald: An N-log(N) method for Ewald sums in large Systems. J Chem Phys 1993, 98: 10089–10092. 10.1063/1.464397View ArticleGoogle Scholar
- Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG: A smooth particle mesh ewald potential. J Chem Phys 1995, 103: 8577–8592. 10.1063/1.470117View ArticleGoogle Scholar
- Miyamoto S, Kollman PA: Settle: An analytical version of the shake and Rattle algorithms for rigid water models. J Comput Chem 1992, 13: 952–962. 10.1002/jcc.540130805View ArticleGoogle Scholar
- Hess B, Bekker H, Berendsen HJC, Fraaije JGEM: LINCS: A linear constraint solver for molecular simulations. J Comput Chem 1997, 18: 1463–1472. 10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-HView ArticleGoogle Scholar
- Hess B: Similarities between principal components of protein dynamics and random diffusion. Phys Rev E 2000, 62: 8438–8448. 10.1103/PhysRevE.62.8438View ArticleGoogle Scholar
- Hess B: Convergence of sampling in protein simulations. Phys Rev E 2002, 65: 031910/1–031910/10.View ArticleGoogle Scholar
- Amadei A, Linssen ABM, Berendsen HJC: Essential dynamics of proteins. Proteins 1993, 17: 412–425. 10.1002/prot.340170408View ArticlePubMedGoogle Scholar
- Haykin S: Neural Networks: a comprehensive foundation. Upper Saddle River, New Jersey, Prentice-Hall Inc; 1999.Google Scholar
- Vesanto J, Himberg J, Alhonemi E, Parhankangas J: SOM Toolbox for Matlab 5. Espoo, Libella Oy; 2000.Google Scholar
- Wehrens R, Buydens LMC: Self- and Super-organising Maps in R: the kohonen package. J Stat Softw 2007., 21(5):Google Scholar
- Kohonen T: Self-Organizing maps. 3rd edition. Gottingen, Springer; 2001.View ArticleGoogle Scholar
- Chan CKK, Hsu AL, Tang SL, Halgamuge SK: Using Growing Self-Organizing Maps to Improve the Binning Process in Environmental Whole-Genome. J Biomed Biotechnol 2008, 2008: 513701.PubMed CentralView ArticlePubMedGoogle Scholar
- Newman AM, Cooper JB: AutoSOME: a clustering method for identifying gene expression modules without prior knowledge of cluster number. BMC Bioinformatics 2010, 11: 117. 10.1186/1471-2105-11-117PubMed CentralView ArticlePubMedGoogle Scholar
- Raskutti B, Leckie C: An Evaluation of Criteria for Measuring the Quality of Clusters. Proceedings of the Sixteenth International Joint Conference on Artificial Intelligence 1999, 905–910.Google Scholar
- Clarke B, Fokoué E, Zhang HH: Principles and theory for data mining and machine learning. Dordrecht, Springer; 2009.View ArticleGoogle Scholar
- Douglas C, Montgomery DC: Design and Analysis of Experiments, Student Solutions Manual. Wiley; 2005.Google Scholar
- Box EP, Hunter WG, Hunter JS: Statistics for Experimenters: An Introduction to Design, Data Analysis, and Model Building. John Wiley and Sons; 1978.Google Scholar
- Juran J: Juran's Quality Control Handbook. McGraw-Hill; 1988.Google Scholar
- Roy R: Design of Experiments Using The Taguchi Approach: 16 Steps to Product and Process Improvement. John Wiley and Sons; 2001.Google Scholar
- Myers RH, Montgomery DC: Response Surface Methodology: Process and Product Optimization Using Designed Experiments. Wiley Series in Probability and Statistics; 1995.Google Scholar
- SAS Institute Inc: JMP, Version 7 Cary, NC. 1989.Google Scholar
- Cressie N, Read TRC: Multinomial goodness-of-fit tests. J R Stat Soc Ser B 1984, 46: 440–464.Google Scholar
- John A Schinka, Wayne F Velicer (Eds): Handbook of Psychology, Research Methods in Psychology. Volume 2. John Wiley and Sons; 2003.Google Scholar
- Mojena R: Hierarchical grouping methods and stopping rules: An evaluation. Comput J 1977, 4: 359–363.View ArticleGoogle Scholar
- The Mathworks Inc: Statistics Toolbox 7 User's Guide Natick, MA. 1993.Google Scholar
- Daura X, Gademann K, Jaun B, Seebach D, van Gunsteren WF, Mark AE: Peptide Folding: When Simulation Meets Experiment. Angew Chem Int Ed 1999, 38: 236–240. 10.1002/(SICI)1521-3773(19990115)38:1/2<236::AID-ANIE236>3.0.CO;2-MView ArticleGoogle Scholar
- Rousseeuw PJ: Silhouettes: a Graphical Aid to the Interpretation and Validation of Cluster Analysis. Journal of Comp Appl Math 1987, 20: 53–65. 10.1016/0377-0427(87)90125-7View ArticleGoogle Scholar
- Davies D, Bouldin D: A cluster separation measure. IEEE Trans Pattern Anal Mach Intelligence 1979, 1: 224–227.View ArticleGoogle Scholar
- Kaufman L, Rousseeuw P: Finding Groups in Data. Wiley, New York; 1990.View ArticleGoogle Scholar
- Musacchio A, Noble M, Pauptit R, Wierenga R, Saraste M: Crystal structure of a Src-homology 3 (SH3) domain. Nature 1992, 359: 851–855. 10.1038/359851a0View ArticlePubMedGoogle Scholar
- Casares S, AB E, Eshuis H, López-Mayorga O, van Nuland NAJ, Conejero-Lara F: The high-resolution NMR structure of the R21A Spc-SH3:P41 complex: Understanding the determinants of binding affinity by comparison with Abl-SH3. BMC Struct Biol 2007, 7: 22. 10.1186/1472-6807-7-22PubMed CentralView ArticlePubMedGoogle Scholar
- Fernandez-Ballester G, Beltrao P, Gonzales JM, Song Y, Wilmanns M, Valencia A, Serrano L: Structure-based prediction of the saccharomyces cerevisiae SH3-ligand interactions. J Mol Biol 2009, 388: 902–916. 10.1016/j.jmb.2009.03.038View ArticlePubMedGoogle Scholar
- Vega MC, Martinez JC, Serrano L: Thermodynamic and structural characterization of Asn and Ala residues in the disallowed II' region of the Ramachandran plot. Protein Sci 2000, 9: 2322–2328. 10.1110/ps.9.12.2322PubMed CentralView ArticlePubMedGoogle Scholar
- Faraldo-Gomez JD, Forrest LR, Baaden M, Bond PJ, Domene C, Patargias G, Cuthbertson J, Sansom MSP: Conformational Sampling and Dynamics of Membrane Proteins From 10-Nanosecond Computer Simulations. Proteins 2004, 57: 783–791. 10.1002/prot.20257View ArticlePubMedGoogle Scholar
- Casares S, Sadqi M, López-Mayorga O, Martìnez JC, Conejero-Lara F: Structural cooperativity in the SH3 domain studied by site-directed mutagenesis and amide hydrogen exchange. FEBS Lett 2003, 539: 125–130. 10.1016/S0014-5793(03)00212-6View ArticlePubMedGoogle Scholar
Copyright
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.