Make Illumina part of your DNA Sign up for PNAS Online eTocs
Link: Info for AuthorsLink: Editorial BoardLink: AboutLink: SubscribeLink: AdvertiseLink: ContactLink: Sitemap Link: PNAS Home
Proceedings of the National Academy of Sciences
Link: Current Issue "" Link: Archives "" Link: Online Submission ""  Link: Advanced Search

Published online on January 10, 2007, 10.1073/pnas.0605534104
PNAS | January 16, 2007 | vol. 104 | no. 3 | 796-801


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Supporting Information
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a colleague
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My File Cabinet
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via ISI Web of Science (10)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Rueda, M.
Right arrow Articles by Orozco, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Rueda, M.
Right arrow Articles by Orozco, M.
Social Bookmarking
Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg  
What's this?

 Previous Article  | Table of Contents |  Next Article 

PHYSICAL SCIENCES / BIOLOGICAL SCIENCES / COMPUTER SCIENCES / BIOPHYSICS
A consensus view of protein dynamics

Manuel Rueda*,{dagger}, Carles Ferrer-Costa*,{dagger}, Tim Meyer*,{dagger},{ddagger}, Alberto Pérez*,{dagger}, Jordi Camps{dagger},§, Adam Hospital*,{dagger},§, Josep Lluis Gelpí*,{dagger},{ddagger}, and Modesto Orozco*,{dagger},{ddagger},§,

*Molecular Modelling and Bioinformatics Unit and §Structural Biology Node, Institut de Recerca Biomèdica, Parc Científic de Barcelona, Josep Samitier 1-5, 08028 Barcelona, Spain; {dagger}Computational Biology Program, Barcelona Supercomputing Center, Jordi Girona 31, Edifici Nexus II, 08028 Barcelona, Spain; and {ddagger}Departament de Bioquímica i Biologia Molecular, Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 645, 08028 Barcelona, Spain

Edited by Harold A. Scheraga, Cornell University, Ithaca, NY, and approved November 10, 2006 (received for review July 6, 2006)


   Abstract
Top
 Abstract
Results and Discussion
 Materials and Methods
 Acknowledgements
 References
 
The dynamics of proteins in aqueous solution has been investigated through a massive approach based on "state of the art" molecular dynamics simulations performed for all protein metafolds using the four most popular force fields (OPLS, CHARMM, AMBER, and GROMOS). A detailed analysis of the massive database of trajectories (>1.5 terabytes of data obtained using {approx}50 years of CPU) allowed us to obtain a robust-consensus picture of protein dynamics in aqueous solution.

force field | molecular dynamics | molecular modeling | protein structure

Flexibility is a key determinant of the biological functionality of proteins. A significant percentage of proteins are unfolded in the absence of ligand, and many others change their conformation as a result of the presence of other molecules or changes in the environment (1–4). Experimental representation of flexibility, even a priori possible, is in general difficult, which makes atomistic simulation the only viable alternative to study this important phenomenon for many proteins. Among the different theoretical methods available for description of protein flexibility, molecular dynamics is probably the most powerful (5–13). Since it was first applied to proteins in the late 1970s, molecular dynamics (MD) has been largely used to study the dynamics of proteins (14, 15). Unfortunately, due to the cost of simulations and the diversity of force fields, a consensus view of protein dynamics has not been yet obtained using this technique. In this paper, we present a systematic study of the most populated protein metafolds [see supporting information (SI) Data Set] using state of the art atomistic MD simulation conditions and the four most widely used force fields [OPLS (O), CHARMM (C), AMBER (A), and GROMOS (G)] (16–25). The result of this massive supercomputer effort (>1.5 terabytes of data and a computational equivalent to 50 CPU years) is a consensus picture of protein dynamics under conditions close to the physiological ones.


   Results and Discussion
Top
 Abstract
 Results and Discussion
Materials and Methods
 Acknowledgements
 References
 
Supporting Information. For further details, see SI Data Set and SI Figs. 5–17

Force Field Convergence. Previous to any dynamic study, we need to determine whether force fields are providing a similar picture of protein structure and whether such a picture is similar to that derived experimentally. Analysis of collected samplings indicates an average divergence ({alpha} in Eq. 1) of {approx}2 Å between force fields (slightly larger deviations are found in G simulations), which, considering the thermal noise of the simulations ({Omega} in Eq. 2), suggests a similarity of {approx}70% between the four samplings and an average "effective distance ({Omega}–1) between them of only 1.4 Å (slightly larger values are always obtained for G-simulations; see Table 2). This finding indicates that all simulations are, in fact, sampling a similar region of the conformational space. This suggestion is supported by the analysis of the radii of gyration and solvent accessible surface (differences of {approx}1% in radii of gyration and 2% in solvent-accessible surface between the four force fields; see Fig. 1).


Figure 1
View larger version (47K):
[in this window]
[in a new window]

  		Fig. 1. Average values for rmsd (Upper Left), rRMSD (based on TM-score; Upper Right) with respect to experimental structures, radii of gyration (Lower Left), and solvent accessible surface (Lower Right) for A, C, O, and G force fields.

Not only are the structures sampled in the four simulations similar, but they are also close to the experimental conformations (see Fig. 1). Thus, differences in radii of gyration between MD samplings and experimental structures are <1%, and average differences in solvent accessible surface are {approx}5%. The average backbone rmsds between simulated and experimental structures are {approx}2.0 Å (A, 1.9; C, 2.0; O, 1.9; G, 2.5 Å), close to the thermal noise of MD simulations. Quite surprisingly, the presence or absence of disulfide bridges does not modify the deviation of our samplings from experimental structures, which is, however, dependent on the origin of the experimental structure. Thus, proteins showing the largest displacements with experiment correspond always to structures solved by NMR spectroscopy (see Fig. 1 and Table 1). The deviations in these cases are due to displacements in regions of large flexibility, where NMR structures are in general poorly defined. This is clear when calculations are repeated using the TM-score, a more robust measure that reduces the weight of flexible entities in the rms-fit. In this case, no system is found for which deviation from experimental data are >2.5 Å (see Fig. 1). Moderate deviations between MD trajectories and experimental model might then reflect either lack of quality in some parts of the experimental model, or just the intrinsic differences between simulation and experimental conditions. Clearly, the commonly made assumption that any deviation from experimental structure of MD simulations is always due to simulation errors is not correct and can ignore important physical characteristics of proteins.


View this table:
[in this window]
[in a new window]

  		Table 1. Structures representative of protein metafolds

 
Native protein–protein contacts are well preserved during simulations (in average, 90% correct predictions). The total number of intramolecular hydrogen bonds found in experimental structures is very well preserved (see Fig. 2), but obviously, hydrogen bond donor and acceptors are more promiscuous in MD than inferred from single experimental structures (see below). All MD simulations preserve very well (90%; see Fig. 2) secondary structure ({alpha}-helices and beta-sheets). In summary: (i) different force fields provide similar picture of protein structure, and (ii) MD simulations sample regions close to experimental structures. Divergences between simulations and experimental structure signal proteins with very flexible moieties whose exact conformation is often not well defined. The two prerequisites noted above are then fulfilled, and we can then safely analyze the dynamics properties of proteins.


Figure 2
View larger version (54K):
[in this window]
[in a new window]

  		Fig. 2. Coverage of contacts (Upper Left), conservation of secondary structure for {alpha}-helix (Upper Right), beta-sheet (Lower Left), and relation between experimental and MD hydrogen bonds (intra; Lower Right) for A, C, O, and G force fields.

B-factors are the simplest method to analyze local deformability and have the advantage to be accessible from x-ray data, which provides us with an additional test on the quality of MD-trajectories. As noted in Fig. 3, there is a very good agreement between experimental and MD-derived B-factors (see Fig. 3). Detailed analysis (see SI Fig. 5) show that such an agreement is especially good for A, C, and O simulations, whereas B-factor distribution is slightly displaced to higher values in G trajectories. In all of the cases, beta-sheet segments are the most rigid, whereas turns are in general the most flexible ones, in good agreement to what is deduced from experimental data. In fact, the only significant difference between MD and x-ray B-factors is that, although experimental B-factors above 60 Å2 are rare, they are not so uncommon in MD simulations. Inspection of B-factor distribution in the different proteins (see example in Fig. 3 and SI Fig. 10) show a good correlation (average Spearman's rank correlation of 0.7 between experimental and MD distributions). This means that proteins residues that appear as the most flexible ones in MD-simulations are those with the largest mobility in the x-ray structures. It seems that MD-simulation is just allowing the movement of flexible residues that were "frozen" by the crystal lattice. This hypothesis can be tested by computation of crystal effective temperatures (see Eq. 3). Considering harmonic residue oscillations as those with B-factors below 60 Å2 (see Fig. 3) we obtain an "effective temperature" of 290 K for the crystal. This demonstrates that the effect of crystal lattice is not, as often assumed, a general "effective cooling" of the general protein. The reduction in kinetic energy (i.e., the reduction of "effective temperature") induced by the lattice is localized in a few very flexible residues, whose mobility is severely reduced with respect to the situation found in diluted aqueous solution.


Figure 3
View larger version (25K):
[in this window]
[in a new window]

  		Fig. 3. B-factor distribution for x-ray structures and from MD simulations (Upper), and per residue B-factors for two proteins: 1CZT and 1LKI (Lower).

Lindemann's disorder index (see Eq. 4) provides a global picture of the flexibility of proteins compared with that of macroscopic solid or liquids (26). Quite independent of the force field and on the presence or absence of disulfide bridges MD simulations suggest that proteins as a whole behave like dense liquids ({Delta}L=0.28 ± 0.06). However, this is the combination of the "solid-like" properties of the buried interior {Delta}L=0.18 ± 0.04 (0.16 ± 0.04 if only buried backbones are considered) and of pure liquid properties of the exposed side-chains {Delta}L=0.38 ± 0.07. In summary, as suggested by experimental studies in particular systems (26, 27), proteins are "melted-solids" with strong differences between the near-solid interior and the full-liquid exterior.

Essential dynamics provides and additional test of the global flexibility of proteins. Interestingly, all of the force fields provide a similar picture of the type of movements that are more important in describing protein similarity ({Gamma} {approx} 85% see Table 2 and SI Figs. 6–15). Space dimensionality (i.e., the number of relevant deformation nodes, ref. 28) depends linearly with the number of residues of the protein (r2 > 0.94; see Fig. 4), with no clear changes induced by the presence of disulfide bridges. Given A, C, and O simulations, we can derive a consensus regression equation (dim = 32.5 + 0.58 x Nres; r2 = 0.93), which indicates that each residue adds only 0.6 dimensions to the conformational space, indicating that the global structure of the protein strongly limits, in a noncooperative manner, the number of accessible conformations of the constituting peptides.


View this table:
[in this window]
[in a new window]

  		Table 2. Similarity indexes

 


Figure 4
View larger version (7K):
[in this window]
[in a new window]

  		Fig. 4. Dimensionality versus the number of residues (Left) and Schlitter's entropy in kcal/mol–1·K–1 (Right) for A, C, O, and G force fields.

The first essential deformation movements have associated very low frequencies (<70 cm–1 for the first 25 essential modes in the C{alpha} space) indicating that significant deformations are expected at room temperature along these modes. Entropies derived from these frequencies (see Materials and Methods) depend in all cases linearly with the number of residues (r2 > 0.98; see Fig. 4), confirming the lack of cooperative effects in determining the accessible conformational space of mono-domain proteins. Perfect match is found between A, O, and C simulations, whereas G trajectories lead to more disordered structures. Interestingly, the presence of disulfide bridges does not introduce any sizeable effect in determining the entropy of the folded protein. Thus, we can conclude that, in contrast with the general believe, disulfide bridges do not especially restrict the conformational space of native proteins.

A single experimental structure presents a clearly defined pattern of stabilizing interactions (salt bridges, hydrogen bonds, and hydrophobic interactions). Obviously, when flexibility is allowed, such a pattern becomes more diffuse and interacting groups are more promiscuous. Thus, despite the total number of stabilizing interactions preserved in MD, there is, in all cases, a reduction in the number of "permanent" interactions (defined as those occurring for >80% of the trajectory) with respect to those in the experimental structure. Thus, the number of permanent hydrogen bonds and hydrophobic interactions are {approx}68 and 81% of those existing in experimental structures. The number of "permanent" salt bridges is {approx}77% the experimental ones for A, C, and O simulations, whereas a bigger loss of salt bridges is found in G simulations (see Table 3). Once "permanent" interactions are located, we can determine their contribution to the global stiffness of the protein by deriving the associated force-constant (see Eq. 6). Results in Table 2 show that, in agreement with chemical intuition and irrespective of the force field used, hydrogen bonds are individually stiffer (force constants in the range 22–29 kcal/mol Å2) than hydrophobic contacts (force constants in the range 3–6 kcal/mol Å2). More surprising is the soft nature of salt bridges (29), especially of those involving lysine, which in fact behave similarly (in terms of stiffness) to hydrophobic contacts. MD strongly suggests that individually, hydrogen bonds are the main responsible of protein stiffness, whereas salt bridges (especially with Lys as donor) are quite labile and promiscuous. It is then clear that none of the theoretically more stabilizing interactions disulfide and salt bridges are controlling protein dynamics.


View this table:
[in this window]
[in a new window]

  		Table 3. Force constants (in kcal/mol A2) associated with different types of protein interactions

 
Protein dynamics is complex; even in the hydrophobic core, side chains may undergo rotameric changes on a time scale longer than that studied in the present work, which is limited to study short time scale dynamics. Furthermore, classic force fields are rough approximations to the real potential energy function. However, our massive dynamic analysis provides strong evidence that most of the behavior detected during MD simulations is force field independent and that deviations from experimental structures, when they exist, are not spurious, but signal a physically meaningful process. Proteins behave like biphasic systems, with a solid core surrounded by external liquid layer, whose flexibility is mostly restricted by hydrogen bonds, not by disulfide or salt bridges. At the monodomain level, the complexity of the accessible conformational space increases slowly, but linearly with the size of the protein, with no obvious cooperativity effects. Finally, and quite surprisingly, crystal lattice does not reduce the kinetic energy of all atoms of the protein, but just of a small set of otherwise too-flexible atoms, those that, when the lattice is released, lead to the largest conformational changes.


   Materials and Methods
Top
 Abstract
 Results and Discussion
 Materials and Methods
Acknowledgements
 References
 
Model Selection. We built a list of proteins to simulate the 30 most populated folds according to SCOP (30, 31), CATH (32), Dali (33, 34), and Dagget's databases (35) (see SI Data Set). When several good three-dimensional models exist in a class, we favor structures without bound ligands.

System Setup. As in our MODEL database (http://mmb.pcb.ub.es/MODEL), we used a common, automatic setup procedure designed to guarantee reasonable ionization states, no electrostatic unbalances, and a good hydration before simulation starts. Thus, in all cases, experimental structures (for NMR, the first one deposited in the Protein Data Bank) were titrated to define the major ionic state a neutral pH, neutralized by ions, minimized, heated, and hydrated (with special care in introducing structural waters). These systems were then preequilibrated for 0.5 ns with parm99-AMBER force field, and then equilibrated (0.5 ns) with each force field (see below).

Simulation Details. Equilibrated structures were used as starting points for 10-ns production trajectories, performed at constant pressure (1 atm) and temperature (300 K) using standard coupling schemes (the same in all cases). Trajectories for three (Protein Data Bank ID codes 1CQY, 1OPC, and 1KTE) ultrarepresentative structures of {alpha}, beta, and {alpha}/beta folds (28) were extended to 100 (A, C, and O) or 50 (G) ns to evaluate the reliability of 10 ns trajectories. Results in SI Fig. 16 strongly suggest that 10 ns is long enough for simulations for many analysis purposes. Particle Mesh Ewald approaches were used to deal with long-range effects (36). Integration of motion equations was performed every 1 fs, the vibrations of bonds involving hydrogen being removed by SHAKE/RATTLE algorithm (37, 38). Four force fields were used in production runs: AMBER (17) (A), CHARMM (18, 19) (C), OPLS/AA (20–23) (O), and GROMOS (24, 25) (G). TIP3P (39, 40) was used as water model for A, C, and O simulations, whereas G calculations were done using the SPC model (41) for coherence with the force field. Simulations were done using parallel versions of AMBER8 (42) (A-simulations), NAMD (43, 44) (C and O simulations), and single processor version of GROMACS (45, 46) (G simulations).

Control simulations were performed to determine whether force fields were able to recognize unfolding conditions and lead to the destabilization of the folded form. For this purpose, 30-ns simulations of 1CQY, 1OPC, and 1KTE were performed using 8 M urea and a temperature of 368 K. OPLS (47) parameters for urea were used for O and A simulations, whereas specific CHARMM parameters (48) were used in C simulations. Due to the lack of specific parameters, these control simulations were not performed using GROMOS force field. Results in SI Fig. 16 confirm that force fields are able to recognize strong denaturing conditions starting the unfolding of the proteins (which might take micro- to millisecond to complete) in the 0- to 30-ns simulation window. The final set with control simulations (10 ns; A, O, C, and G trajectories) was performed for three proteins (Protein Data Bank ID codes 2GB1, 1LYS, and 1UBQ) for which residual dipolar couplings and S2 distributions have been reported (see references in the legend of SI Fig. 17) from NMR experiments. Comparison of MD- and NMR-derived measures of flexibility is quite satisfactory (see SI Fig. 17) giving confidence in the quality of our theoretical estimates of flexibility (it is not possible to perform this comparison for representative proteins in SI Table 4 due to the lack of experimental data).

Analysis of Trajectories. Trajectories were analyzed to obtain structural and dynamic properties. Structural descriptors include backbone rmsd, TM score (49), radii of gyration, solvent accessible surface for all heavy atoms (50), secondary structure (51), intramolecular hydrogen bonds, contact maps, Ramachandran plots (52, 53), solvent contact profiles, and others. The global similarity between the structures collected with the different force fields was obtained by computing the pair-cross rmsd (i.e., the rmsd between all of the snapshots collected in the two trajectories; see Eq. 1), and a related similarity index (see Eq. 2).

Formula 1

where N is the number of atoms and M is the number of frames

Formula 2

Local dynamic properties are represented by the B factors, which are compared with x-ray values to obtain a measure of "crystal effective" temperatures (see Eq. 3). Global flexibility was represented by Lindemann's indexes (26) (see Eq. 4) and by essential dynamics protocols considering either Cartesian or mass-weighted covariance matrices (54), which provided us with measures as protein entropy (55–57) or the dimensionality of the conformational space (28). The similarity between the nature of the essential deformation patterns of two trajectories (of the same protein) was determined by {Gamma} index (see Eq. 5; ref. 58), considering a small set of 25 eigenvectors (which represent {approx}80–85% of protein variance).

Formula 3


Formula 4

where a' is the most probable nonbonded near-neighbor distance, N is the number of atoms, and {Delta}ri stands for the fluctuation of atom i from its average position.

Formula 5

where n is the number of size of the important space (n = 25 here) and {nu} stands for one unitary eigenvector. Note that {Gamma} = 1 indicates that both trajectories are sampling the same type of essential movements, whereas {Gamma} = 0 means that they are orthogonal.

The flexibility related to key interactions (hydrogen bonds, hydrophobic contacts, and salt bridges) was computed by assuming the harmonic oscillator model (see Eq. 6). Because the use of harmonic model implies a Gaussian distribution of distances, we computed force constant (K in Eq. 6) for "stable" interactions considered as those found in >80% of the trajectory. A salt bridge was defined when the distance N{zeta}(K)/C{zeta} (R) – C{gamma}(D)/C{delta}(E) was <6.5 Å, hydrogen bonds were annotated by using standard PTRAJ criteria (42), and hydrophobic contacts were counted when two hydrophobic groups (A, V, P, F, M, I, L, W), which are not neighbors in the sequence, have their Cbeta atoms separated by <10 Å.

Formula 6

where kb is the Boltzmann constant, T is the absolute temperature, and {Delta}X is the oscillation in the interaction distance from average values.


   Acknowledgements
Top
 Abstract
 Results and Discussion
 Materials and Methods
 Acknowledgements
References
 
We are indebted to Prof. F. J. Luque and Dr. P. Bernadó for helpful discussions and A. Caballero-Herrera for CHARMM urea parameters. T.M. and A.P. are fellows of the Spanish and Catalan Ministries of Science. This work was supported by the Spanish Ministry of Education and Science (BIO2006-01602 and GEN2003-20642-CO9-07). All calculations were carried out on MareNostrum, a supercomputer with {approx}4,524 64-bit Myrinet-connected processors.


   Footnotes
 

Abbreviations: MD, molecular dynamics; O, OPLS; C, CHARMM; A, AMBER; G, GROMOS.

To whom correspondence should be addressed. E-mail: modesto{at}mmb.pcb.ub.es

Author contributions: M.R., C.F.-C., and T.M. contributed equally to this work; M.R., C.F.-C., A.H., J.L.G., and M.O. designed research; M.R., C.F.-C., T.M., A.P., and J.C. performed research; M.R., C.F.-C., T.M., and A.P. analyzed data; and M.O. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS direct submission.

This article contains supporting information online at www.pnas.org/cgi/content/full/0605534104/DC1.

© 2007 by The National Academy of Sciences of the USA


   References
Top
 Abstract
 Results and Discussion
 Materials and Methods
 Acknowledgements
 References
 

  1. Waldron, TT & Murphy, KP. (2003) Biochemistry 42, 5058–5064.[CrossRef][Medline]
  2. Lee-Wei, Y & Bahar, I. (2005) Structure (London) 13, 893–904.
  3. Sacquin-Mora, S & Lavery, R. (2006) Biophys J 90, 2706–2717.
  4. Remy, I, Wilson, IA & Michnick, SW. (1999) Science 283, 990–993.[Abstract/Free Full Text]
  5. McCammon, JA, Gelin, BR & Karplus, M. (1977) Nature 267, 585–590.[CrossRef][Medline]
  6. Karplus, M & McCammon, JA. (1986) Sci Am 254, 42–51.[ISI][Medline]
  7. Allen, MP & Tildesley, DJ. (1989) Computer Simulation of Liquids (Clarendon, Oxford,).
  8. Brooks, CL III, Karplus, M & Pettitt, BM. (1987) Proteins: A theoretical Perspective of Dynamics, Structure and Thermodynamics (Cambridge Univ Press, Cambridge, UK,).
  9. Warshel, A. (1976) Nature 260, 679–683.[CrossRef][Medline]
  10. Ponder, JW & Case, DA. (2003) Adv Protein Chem 66, 27–85.[ISI][Medline]
  11. van Gunsteren, WF, Berendsen, HJ, Hermans, J, Hol, WG & Postma, JP. (1983) Proc Natl Acad Sci USA 80, 4315–4319.[Abstract/Free Full Text]
  12. Tirado-Rives, J & Jorgensen, WL. (1991) Biochemistry 30, 3864–3871.[CrossRef][Medline]
  13. Tirado-Rives, J & Jorgensen, WL. (1990) J Am Chem Soc 112, 2773–2781.[CrossRef]
  14. Karplus, M & McCammon, JA. (2002) Nat Struct Biol 9, 646–652.[CrossRef][ISI][Medline]
  15. Karplus, M & Kuriyan, J. (2005) Proc Natl Acad Sci USA 102, 6679–6685.[Abstract/Free Full Text]
  16. Weiner, SJ, Kollman, PA, Case, DA, Singh, UC, Ghio, C, Alagona, G, Profeta, S Jr & Weiner, P. (1984) J Am Chem Soc 106, 765–784.[CrossRef]
  17. Cornell, WD, Cieplak, P, Bayly, CI, Gould, IR, Merz, KM Jr, Ferguson, DM, Spellmeyer, DC, Fox, T, Caldwell, JW & Kollman, PA. (1995) J Am Chem Soc 117, 5179–5197.[CrossRef]
  18. MacKerell, AD Jr, Bashford, D, Bellott, M, Dunbrack, RL, Evanseck, JD, Field, MJ, Fischer, S, Gao, J, Guo, H & Ha, S, et al. (1998) J Phys Chem B 102, 3586–3616.[CrossRef]
  19. MacKerell, AD Jr, Wiorkiewicz-Kuczera, J & Karplus, M. (1995) J Am Chem Soc 117, 11946–11975.[CrossRef]
  20. Damm, W, Frontera, A, Tirado-Rives, J & Jorgensen, WL. (1997) J Comput Chem 18, 1955–1970.[CrossRef]
  21. Jorgensen, WL, Maxwell, DS & Tirado-Rives, J. (1996) J Am Chem Soc 118, 11225–11236.[CrossRef]
  22. Kaminski, G, Duffy, EM, Matsui, T & Jorgensen, WL. (1994) J Phys Chem 98, 13077–13082.[CrossRef]
  23. Kaminski, GA, Friesner, RA, Tirado-Rives, J & Jorgensen, WL. (2001) J Phys Chem B 105, 6474–6487.[CrossRef]
  24. Ott, K-H & Meyer, B. (1996) J Comput Chem 17, 1068–1084.[CrossRef]
  25. Hermans, J, Berendsen, HJC, Van Gunsteren, WF & Postma, JPM. (1984) Biopolymers 23, 1513–1518.[CrossRef][ISI]
  26. Zhou, Y, Vitkup, D & Karplus, M. (1999) J Mol Biol 285, 1371–1375.[CrossRef][ISI][Medline]
  27. Lindorff-Larsen, K, Best, RB, Depristo, MA, Dobson, CM & Vendruscolo, M. (2005) Nature 433, 128–132.[CrossRef][Medline]
  28. Perez, A, Blas, JR, Rueda, M, López-Bes, JM, de La Cruz, X, Luque, FJ & Orozco, M. (2005) J Chem Theory Comput 1, 790–800.[CrossRef]
  29. Strop, P & Mayo, SL. (2000) Biochemistry 39, 1251–1255.[CrossRef][Medline]
  30. Andreeva, A, Howorth, D, Brenner, SE, Hubbard, TJ, Chothia, C & Murzin, AG. (2004) Nucleic Acids Res 32, D226–D229.[Abstract/Free Full Text]
  31. Murzin, AG, Brenner, SE, Hubbard, T & Chothia, C. (1995) J Mol Biol 247, 536–540.[CrossRef][ISI][Medline]
  32. Pearl, F, Todd, A, Sillitoe, I, Dibley, M, Redfern, O, Lewis, T, Bennett, C, Marsden, R, Grant, A & Lee, D, et al. (2005) Nucleic Acids Res 33, D247–D251.[Abstract/Free Full Text]
  33. Dietmann, S, Park, J, Notredame, C, Heger, A, Lappe, M & Holm, L. (2001) Nucleic Acids Res 29, 55–57.[Abstract/Free Full Text]
  34. Holm, L & Sander, C. (1995) Trends Biochem Sci 20, 478–480.[CrossRef][ISI][Medline]
  35. Day, R, Beck, DA, Armen, RS & Daggett, V. (2003) Protein Sci 12, 2150–2160.[Abstract/Free Full Text]
  36. Darden, TL, York, D & Pedersen, L. (1993) J Chem Phys 98, 10089–10092.[CrossRef]
  37. Ryckaert, JP, Ciccotti, G & Berendsen, HJC. (1977) J Comput Phys 23, 327–341.[CrossRef]
  38. Andersen, HC. (1983) J Comput Phys 52, 24–34.[CrossRef]
  39. Jorgensen, WL, Chandrasekhar, J, Madura, JD, Impey, RW & Klein, ML. (1983) J Chem Phys 79, 926–935.[CrossRef]
  40. Mahoney, MW & Jorgensen, WL. (2000) J Chem Phys 112, 8910–8922.[CrossRef]
  41. Berendsen, HJC, Postma, JPM, van Gunsteren, WF & Hermans, J. (1981) in Intermolecular Forces ed. Pullman, B. (Reidel, Dordrecht, The Netherlands,) pp. 331.
  42. Case, DA, Pearlman, DA, Caldwell, JW, Cheatham, TE III, Ross, WS, Simmerling, CL, Darden, TL, Marz, KM, Stanton, RV & Cheng, AL, et al. (2004) AMBER 8 (University of California, San Francisco,).
  43. Kale, L, Skeel, R, Bhandarkar, M, Brunner, R, Gursoy, A, Krawetz, N, Phillips, J, Shinozaki, A, Varadarajan, K & Schulten, K. (1999) J Comput Phys 151, 283–312.[CrossRef]
  44. Phillips, JC, Braun, R, Wang, W, Gumbart, J, Tajkhorshid, E, Villa, E, Chipot, C, Skeel, RD, Kale, L & Schulten, K. (2005) J Comput Chem 26, 1781–1802.[CrossRef][ISI][Medline]
  45. Berendsen, HJC, van der Spoel, D & van Drunen, R. (1995) Comput Phys Commun 91, 43–56.[CrossRef]
  46. Van Der Spoel, D, Lindahl, E, Hess, B, Groenhof, G, Mark, AE & Berendsen, HJ. (2005) J Comput Chem 26, 1701–1718.[CrossRef][ISI][Medline]
  47. Tirado-Rives, J, Orozco, M & Jorgensen, WL. (1997) Biochemistry 36, 7313–7329.[CrossRef][Medline]
  48. Caballero-Herrera, A & Nilsson, A. (2006) J Mol Struct 758, 139–148.
  49. Zhang, Y & Skolnick, J. (2004) Proteins 57, 702–710.[CrossRef][ISI][Medline]
  50. Hubbard, SJ & Thornton, JM. (1993) NACCESS (Dept of Biochemistry and Molecular Biology, University College, London,).
  51. Kabsch, W & Sander, C. (1983) Biopolymers 22, 2577–2637.[CrossRef][ISI][Medline]
  52. Laskowski, RA, MacArthur, MW, Moss, DS & Thornton, JM. (1993) J Appl Crystallogr 26, 283–291.[CrossRef]
  53. Morris, AL, MacArthur, MW, Hutchinson, EG & Thornton, JM. (1992) Proteins 12, 345–364.[CrossRef][ISI][Medline]
  54. Amadei, A, Linssen, AB & Berendsen, HJ. (1993) Proteins 17, 412–425.[CrossRef][ISI][Medline]
  55. Harris, SA, Gavathiotis, E, Searle, MS, Orozco, M & Laughton, CA. (2001) J Am Chem Soc 123, 12658–12663.[CrossRef][ISI][Medline]
  56. Schlitter, J. (1993) Chem Phys Lett 215, 617–621.[CrossRef]
  57. Andricioaei, I & Karplus, M. (2001) J Chem Phys 115, 6289–6292.[CrossRef]
  58. Hess, B. (2000) Phys Rev E Stat Nonlin Soft Matter Phys 62, 8438–8448.

Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg    What's this?


This article has been cited by other articles in HighWire Press-hosted journals:


Home page
 Nucleic Acids ResHome page
A. Perez, F. Lankas, F. J. Luque, and M. Orozco
Towards a molecular dynamics consensus view of B-DNA flexibility
Nucleic Acids Res., April 1, 2008; 36(7): 2379 - 2394.
[Abstract] [Full Text] [PDF]

Home page
 Brief BioinformHome page
X. Liu and H. A. Karimi
High-throughput modeling and analysis of protein structural dynamics
Brief Bioinform, November 1, 2007; 8(6): 432 - 445.
[Abstract] [Full Text] [PDF]

Home page
 J. Biol. Chem.Home page
J. Fort, L. R. de la Ballina, H. E. Burghardt, C. Ferrer-Costa, J. Turnay, C. Ferrer-Orta, I. Uson, A. Zorzano, J. Fernandez-Recio, M. Orozco, et al.
The Structure of Human 4F2hc Ectodomain Provides a Model for Homodimerization and Electrostatic Interaction with Plasma Membrane
J. Biol. Chem., October 26, 2007; 282(43): 31444 - 31452.
[Abstract] [Full Text] [PDF]

Home page
 BioinformaticsHome page
M. Seeber, M. Cecchini, F. Rao, G. Settanni, and A. Caflisch
Wordom: a program for efficient analysis of molecular dynamics simulations
Bioinformatics, October 1, 2007; 23(19): 2625 - 2627.
[Abstract] [Full Text] [PDF]

Home page
 J. Biol. Chem.Home page
B. Holst, J. Mokrosinski, M. Lang, E. Brandt, R. Nygaard, T. M. Frimurer, A. G. Beck-Sickinger, and T. W. Schwartz
Identification of an Efficacy Switch Region in the Ghrelin Receptor Responsible for Interchange between Agonism and Inverse Agonism
J. Biol. Chem., May 25, 2007; 282(21): 15799 - 15811.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Supporting Information
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a colleague
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My File Cabinet
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via ISI Web of Science (10)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Rueda, M.
Right arrow Articles by Orozco, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Rueda, M.
Right arrow Articles by Orozco, M.
Social Bookmarking
Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg  
What's this?