51
|
Schuyler AD, Jernigan RL, Qasba PK, Ramakrishnan B, Chirikjian GS. Iterative cluster-NMA: A tool for generating conformational transitions in proteins. Proteins 2009; 74:760-76. [PMID: 18712827 DOI: 10.1002/prot.22200] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Computational models provide insight into the structure-function relationship in proteins. These approaches, especially those based on normal mode analysis, can identify the accessible motion space around a given equilibrium structure. The large magnitude, collective motions identified by these methods are often well aligned with the general direction of the expected conformational transitions. However, these motions cannot realistically be extrapolated beyond the local neighborhood of the starting conformation. In this article, the iterative cluster-NMA (icNMA) method is presented for traversing the energy landscape from a starting conformation to a desired goal conformation. This is accomplished by allowing the evolving geometry of the intermediate structures to define the local accessible motion space, and thus produce an appropriate displacement. Following the derivation of the icNMA method, a set of sample simulations are performed to probe the robustness of the model. A detailed analysis of beta1,4-galactosyltransferase-T1 is also given, to highlight many of the capabilities of icNMA. Remarkably, during the transition, a helix is seen to be extended by an additional turn, emphasizing a new unknown role for secondary structures to absorb slack during transitions. The transition pathway for adenylate kinase, which has been frequently studied in the literature, is also discussed.
Collapse
Affiliation(s)
- Adam D Schuyler
- Department of Neurology, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | | | | | | |
Collapse
|
52
|
Flexible fitting of high-resolution x-ray structures into cryoelectron microscopy maps using biased molecular dynamics simulations. Biophys J 2008; 95:5692-705. [PMID: 18849406 DOI: 10.1529/biophysj.108.139451] [Citation(s) in RCA: 80] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
A methodology for flexible fitting of all-atom high-resolution structures into low-resolution cryoelectron microscopy (cryo-EM) maps is presented. Flexibility of the modeled structure is simulated by classical molecular dynamics and an additional effective potential is introduced to enhance the fitting process. The additional potential is proportional to the correlation coefficient between the experimental cryo-EM map and a synthetic map generated for an all-atom structure being fitted to the map. The additional forces are calculated as a gradient of the correlation coefficient. During the molecular dynamics simulations under the additional forces, the molecule undergoes a conformational transition that maximizes the correlation coefficient, which results in a high-accuracy fit of all-atom structure into a cryo-EM map. Using five test proteins that exhibit structural rearrangement during their biological activity, we demonstrate performance of our method. We also test our method on the experimental cryo-EM of elongation factor G and show that the model obtained is comparable to previous studies. In addition, we show that overfitting can be avoided by assessing the quality of the fitted model in terms of correlation coefficient and secondary structure preservation.
Collapse
|
53
|
Matsumoto A, Kamata T, Takagi J, Iwasaki K, Yura K. Key interactions in integrin ectodomain responsible for global conformational change detected by elastic network normal-mode analysis. Biophys J 2008; 95:2895-908. [PMID: 18515366 PMCID: PMC2527288 DOI: 10.1529/biophysj.108.131045] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2008] [Accepted: 05/06/2008] [Indexed: 01/03/2023] Open
Abstract
Integrin, a membrane protein with a huge extracellular domain, participates in cell-cell and cell-extracellular-matrix interactions for metazoan. A group of integrins is known to perform a large-scale structural change when the protein is activated, but the activation mechanism and generality of the conformational change remain to be elucidated. We performed normal-mode analysis of the elastic network model on integrin alpha(V)beta(3) ectodomain in the bent form and identified key residues that influenced molecular motions. Iterative normal-mode calculations demonstrated that the specific nonbonded interactions involving the key residues work as a snap to keep integrin in the bent form. The importance of the key residues for the conformational change was further verified by mutation experiments, in which integrin alpha(IIb)beta(3) was used. The conservation pattern of amino acid residues among the integrin family showed that the characteristic pattern of residues seen around these key residues is found in the limited groups of integrin beta-chains. This conservation pattern suggests that the molecular mechanism of the conformational change relying on the interactions found in integrin alpha(V)beta(3) is unique to the limited types of integrins.
Collapse
Affiliation(s)
- Atsushi Matsumoto
- Quantum Bioinformatics Team, Center for Computational Science and Engineering, Japan Atomic Energy Agency, 8-1 Umemidai, Kizugawa, Kyoto 619-0215, Japan.
| | | | | | | | | |
Collapse
|
54
|
Multiple subunit fitting into a low-resolution density map of a macromolecular complex using a gaussian mixture model. Biophys J 2008; 95:4643-58. [PMID: 18708469 PMCID: PMC2576401 DOI: 10.1529/biophysj.108.137125] [Citation(s) in RCA: 82] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
Recently, electron microscopy measurement of single particles has enabled us to reconstruct a low-resolution 3D density map of large biomolecular complexes. If structures of the complex subunits can be solved by x-ray crystallography at atomic resolution, fitting these models into the 3D density map can generate an atomic resolution model of the entire large complex. The fitting of multiple subunits, however, generally requires large computational costs; therefore, development of an efficient algorithm is required. We developed a fast fitting program, “gmfit”, which employs a Gaussian mixture model (GMM) to represent approximated shapes of the 3D density map and the atomic models. A GMM is a distribution function composed by adding together several 3D Gaussian density functions. Because our model analytically provides an integral of a product of two distribution functions, it enables us to quickly calculate the fitness of the density map and the atomic models. Using the integral, two types of potential energy function are introduced: the attraction potential energy between a 3D density map and each subunit, and the repulsion potential energy between subunits. The restraint energy for symmetry is also employed to build symmetrical origomeric complexes. To find the optimal configuration of subunits, we randomly generated initial configurations of subunit models, and performed a steepest-descent method using forces and torques of the three potential energies. Comparison between an original density map and its GMM showed that the required number of Gaussian distribution functions for a given accuracy depended on both resolution and molecular size. We then performed test fitting calculations for simulated low-resolution density maps of atomic models of homodimer, trimer, and hexamer, using different search parameters. The results indicated that our method was able to rebuild atomic models of a complex even for maps of 30 Å resolution if sufficient numbers (eight or more) of Gaussian distribution functions were employed for each subunit, and the symmetric restraints were assigned for complexes with more than three subunits. As a more realistic test, we tried to build an atomic model of the GroEL/ES complex by fitting 21-subunit atomic models into the 3D density map obtained by cryoelectron microscopy using the C7 symmetric restraints. A model with low root mean-square deviations (14.7 Å) was obtained as the lowest-energy model, showing that our fitting method was reasonably accurate. Inclusion of other restraints from biological and biochemical experiments could further enhance the accuracy.
Collapse
|
55
|
Abstract
In fitting atomic structures into EM maps, it often happens that the map corresponds to a different conformation of the structure. We have developed a new methodology to handle these situations that preserves the covalent geometry of the structure and allows the modeling of large deformations. The first goal is achieved by working in generalized coordinates (positional and internal coordinates), and the second by avoiding harmonic potentials. Instead, we use dampers (shock absorbers) between every pair of atoms, combined with a force field that attracts the atomic structure toward incompletely occupied regions of the EM map. The trajectory obtained by integrating the resulting equations of motion converges to a conformation that, in our validation cases, was very close to the target atomic structure. Compared to current methods, our approach is more efficient and robust against wrong solutions and to overfitting, and does not require user intervention or subjective decisions. Applications to the computation of transition pathways between known conformers, homology and loop modeling, as well as protein docking, are also discussed.
Collapse
|
56
|
X-ray solution scattering (SAXS) combined with crystallography and computation: defining accurate macromolecular structures, conformations and assemblies in solution. Q Rev Biophys 2008; 40:191-285. [PMID: 18078545 DOI: 10.1017/s0033583507004635] [Citation(s) in RCA: 845] [Impact Index Per Article: 52.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Crystallography supplies unparalleled detail on structural information critical for mechanistic analyses; however, it is restricted to describing low energy conformations of macromolecules within crystal lattices. Small angle X-ray scattering (SAXS) offers complementary information about macromolecular folding, unfolding, aggregation, extended conformations, flexibly linked domains, shape, conformation, and assembly state in solution, albeit at the lower resolution range of about 50 A to 10 A resolution, but without the size limitations inherent in NMR and electron microscopy studies. Together these techniques can allow multi-scale modeling to create complete and accurate images of macromolecules for modeling allosteric mechanisms, supramolecular complexes, and dynamic molecular machines acting in diverse processes ranging from eukaryotic DNA replication, recombination and repair to microbial membrane secretion and assembly systems. This review addresses both theoretical and practical concepts, concerns and considerations for using these techniques in conjunction with computational methods to productively combine solution scattering data with high-resolution structures. Detailed aspects of SAXS experimental results are considered with a focus on data interpretation tools suitable to model protein and nucleic acid macromolecular structures, including membrane protein, RNA, DNA, and protein-nucleic acid complexes. The methods discussed provide the basis to examine molecular interactions in solution and to study macromolecular flexibility and conformational changes that have become increasingly relevant for accurate understanding, simulation, and prediction of mechanisms in structural cell biology and nanotechnology.
Collapse
|
57
|
Yang LW, Chng CP. Coarse-grained models reveal functional dynamics--I. Elastic network models--theories, comparisons and perspectives. Bioinform Biol Insights 2008; 2:25-45. [PMID: 19812764 PMCID: PMC2735964 DOI: 10.4137/bbi.s460] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
In this review, we summarize the progress on coarse-grained elastic network models (CG-ENMs) in the past decade. Theories were formulated to allow study of conformational dynamics in time/space frames of biological interest. Several highlighted models and their underlined hypotheses are introduced in physical depth. Important ENM offshoots, motivated to reproduce experimental data as well as to address the slow-mode-encoded configurational transitions, are also introduced. With the theoretical developments, computational cost is significantly reduced due to simplified potentials and coarse-grained schemes. Accumulating wealth of data suggest that ENMs agree equally well with experiment in describing equilibrium dynamics despite their distinct potentials and levels of coarse-graining. They however do differ in the slowest motional components that are essential to address large conformational changes of functional significance. The difference stems from the dissimilar curvatures of the harmonic energy wells described for each model. We also provide our views on the predictability of 'open to close' (open-->close) transitions of biomolecules on the basis of conformational selection theory. Lastly, we address the limitations of the ENM formalism which are partially alleviated by the complementary CG-MD approach, to be introduced in the second paper of this two-part series.
Collapse
Affiliation(s)
- Lee-Wei Yang
- Institute of Molecular and Cellular Biosciences, University of Tokyo, Tokyo 113-0032, Japan.
| | | |
Collapse
|
58
|
Schröder GF, Brunger AT, Levitt M. Combining efficient conformational sampling with a deformable elastic network model facilitates structure refinement at low resolution. Structure 2008; 15:1630-41. [PMID: 18073112 DOI: 10.1016/j.str.2007.09.021] [Citation(s) in RCA: 198] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2007] [Revised: 09/11/2007] [Accepted: 09/26/2007] [Indexed: 11/19/2022]
Abstract
Structural studies of large proteins and protein assemblies are a difficult and pressing challenge in molecular biology. Experiments often yield only low-resolution or sparse data that are not sufficient to fully determine atomistic structures. We have developed a general geometry-based algorithm that efficiently samples conformational space under constraints imposed by low-resolution density maps obtained from electron microscopy or X-ray crystallography experiments. A deformable elastic network (DEN) is used to restrain the sampling to prior knowledge of an approximate structure. The DEN restraints dramatically reduce over-fitting, especially at low resolution. Cross-validation is used to optimally weight the structural information and experimental data. Our algorithm is robust even for noise-added density maps and has a large radius of convergence for our test case. The DEN restraints can also be used to enhance reciprocal space simulated annealing refinement.
Collapse
Affiliation(s)
- Gunnar F Schröder
- Department of Structural Biology, Stanford University Stanford, CA 94305, USA.
| | | | | |
Collapse
|
59
|
Delarue M. Dealing with structural variability in molecular replacement and crystallographic refinement through normal-mode analysis. ACTA CRYSTALLOGRAPHICA SECTION D: BIOLOGICAL CRYSTALLOGRAPHY 2007; 64:40-8. [PMID: 18094466 PMCID: PMC2394787 DOI: 10.1107/s0907444907053516] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/03/2007] [Accepted: 10/26/2007] [Indexed: 01/08/2023]
Abstract
Normal-mode analysis (NMA) can be used to generate multiple structural variants of a given template model, thereby increasing the chance of finding the molecular-replacement solution. Here, it is shown that it is also possible to directly refine the amplitudes of the normal modes against experimental data (X-ray or cryo-EM), generalizing rigid-body refinement methods by adding just a few additional degrees of freedom that sample collective and large-amplitude movements. It is also argued that the situation where several (conformations of) models are present simultaneously in the crystal can be studied with adjustable occupancies using techniques derived from statistical thermodynamics and already used in molecular modelling.
Collapse
Affiliation(s)
- Marc Delarue
- Unité de Dynamique Structurale des Macromolecules, Departement de Biologie Structurale et Chimie, URA 2185 du CNRS, Institut Pasteur, 25 Rue du Dr Roux, 75015 Paris, France.
| |
Collapse
|
60
|
Fitting low-resolution cryo-EM maps of proteins using constrained geometric simulations. Biophys J 2007; 94:1613-21. [PMID: 17993504 DOI: 10.1529/biophysj.107.115949] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Recent experimental advances in producing density maps from cryo-electron microscopy (cryo-EM) have challenged theorists to develop improved techniques to provide structural models that are consistent with the data and that preserve all the local stereochemistry associated with the biomolecule. We develop a new technique that maintains the local geometry and chemistry at each stage of the fitting procedure. A geometric simulation is used to drive the structure from some appropriate starting point (a nearby experimental structure or a modeled structure) toward the experimental density, via a set of small incremental motions. Structural motifs such as alpha-helices can be held rigid during the fitting procedure as the starting structure is brought into alignment with the experimental density. After validating this procedure on simulated data for adenylate kinase and lactoferrin, we show how cryo-EM data for two different GroEL structures can be fit using a starting x-ray crystal structure. We show that by incorporating the correct local stereochemistry in the modeling, structures can be obtained with effective resolution that is significantly higher than might be expected from the nominal cryo-EM resolution.
Collapse
|
61
|
Normal-mode flexible fitting of high-resolution structure of biological molecules toward one-dimensional low-resolution data. Biophys J 2007; 94:1589-99. [PMID: 17993489 DOI: 10.1529/biophysj.107.122218] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We present a method for reconstructing a 3D structure from a pair distribution function by flexibly fitting known x-ray structures toward a conformation that agrees with the low-resolution data. This method uses a linear combination of low-frequency normal modes from elastic-network description of the molecule in an iterative manner to deform the structure optimally to conform to the target pair distribution function. A simple function, pair distance distribution function between atoms, is chosen as a test model to establish computational algorithms-optimization algorithm and scoring function-that can utilize low-resolution 1D data. To select a correct structural model based on less information, we developed a scoring function that takes into account a characteristic of pair distribution functions. In addition, we employ a new optimization algorithm, the trusted region method, that relies on both first and second derivatives of the scoring function. Illustrative results of our studies on simulated 1D data from five different proteins, for which large conformational changes are known to occur, are presented.
Collapse
|
62
|
Lacapère JJ, Pebay-Peyroula E, Neumann JM, Etchebest C. Determining membrane protein structures: still a challenge! Trends Biochem Sci 2007; 32:259-70. [PMID: 17481903 DOI: 10.1016/j.tibs.2007.04.001] [Citation(s) in RCA: 133] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2006] [Revised: 03/07/2007] [Accepted: 04/13/2007] [Indexed: 11/20/2022]
Abstract
Determination of structures and dynamics events of transmembrane proteins is important for the understanding of their function. Analysis of such events requires high-resolution 3D structures of the different conformations coupled with molecular dynamics analyses describing the conformational pathways. However, the solution of 3D structures of transmembrane proteins at atomic level remains a particular challenge for structural biochemists--the need for purified and functional transmembrane proteins causes a 'bottleneck'. There are various ways to obtain 3D structures: X-ray diffraction, electron microscopy, NMR and modelling; these methods are not used exclusively of each other, and the chosen combination depends on several criteria. Progress in this field will improve knowledge of ligand-induced activation and inhibition of membrane proteins in addition to aiding the design of membrane-protein-targeted drugs.
Collapse
Affiliation(s)
- Jean-Jacques Lacapère
- INSERM, U773, Centre de Recherche Biomédicale Bichat Beaujon CRB3, Faculté de Médecine X. Bichat, Université Paris 7, BP 416, F-75018, Paris, France.
| | | | | | | |
Collapse
|
63
|
Jernigan RL, Kloczkowski A. Packing regularities in biological structures relate to their dynamics. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2006; 350:251-76. [PMID: 16957327 PMCID: PMC2039702 DOI: 10.1385/1-59745-189-4:251] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/28/2023]
Abstract
The high packing density inside proteins leads to certain geometric regularities and also is one of the most important contributors to the high extent of cooperativity manifested by proteins in their cohesive domain motions. The orientations between neighboring nonbonded residues in proteins substantially follow the similar geometric regularities, regardless of whether the residues are on the surface or buried, a direct result of hydrophobicity forces. These orientations are relatively fixed and correspond closely to small deformations from those of the face-centered cubic lattice, which is the way in which identical spheres pack at the highest density. Packing density also is related to the extent of conservation of residues, and we show this relationship for residue packing densities by averaging over a large sample or residue packings. There are three regimes: (1) over a broad range of packing densities the relationship between sequence entropy and inverse packing density is nearly linear, (2) over a limited range of low packing densities the sequence entropy is nearly constant, and (3) at extremely low packing densities the sequence entropy is highly variable. These packing results provide important justification for the simple elastic network models that have been shown for a large number of proteins to represent protein dynamics so successfully, even when the models are extremely coarse grained. Elastic network models for polymeric chains are simple and could be combined with these protein elastic networks to represent partially denatured parts of proteins. Finally, we show results of applications of the elastic network model to study the functional motions of the ribosome, based on its known structure. These results indicate expected correlations among its components for the step-wise processing steps in protein synthesis, and suggest ways to use these elastic network models to develop more detailed mechanisms, an important possibility because most experiments yield only static structures.
Collapse
Affiliation(s)
- Robert L Jernigan
- Department of Biochemistry, Biophysics, and Molecular Biology, Laurence H. Baker Center for Bioinformatics and Biological Statistics, Iowa State University, Ames, IA, USA
| | | |
Collapse
|
64
|
Almeida WI, Martins OB, Carvalho-Alves PC. Self-association of isolated large cytoplasmic domain of plasma membrane H+ -ATPase from Saccharomyces cerevisiae: role of the phosphorylation domain in a general dimeric model for P-ATPases. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2006; 1758:1768-76. [PMID: 17026955 DOI: 10.1016/j.bbamem.2006.08.018] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/10/2006] [Revised: 08/21/2006] [Accepted: 08/22/2006] [Indexed: 10/24/2022]
Abstract
Large cytoplasmic domain (LCD) plasma membrane H+ -ATPase from S. cerevisiae was expressed as two fusion polypeptides in E. coli: a DNA sequence coding for Leu353-Ileu674 (LCDh), comprising both nucleotide (N) and phosphorylation (P) domains, and a DNA sequence coding for Leu353-Thr543 (LCDDeltah, lacking the C-terminus of P domain), were inserted in expression vectors pDEST-17, yielding the respective recombinant plasmids. Overexpressed fusion polypeptides were solubilized with 6 M urea and purified on affinity columns, and urea was removed by dialysis. Their predicted secondary structure contents were confirmed by CD spectra. In addition, both recombinant polypeptides exhibited high-affinity 2',3'-O-(2,4,6-trinitrophenyl)adenosine-5'-triphosphate (TNP-ATP) binding (Kd = 1.9 microM and 2.9 microM for LCDh and LCDDeltah, respectively), suggesting that they have native-like folding. The gel filtration profile (HPLC) of purified LCDh showed two main peaks, with molecular weights of 95 kDa and 39 kDa, compatible with dimeric and monomeric forms, respectively. However, a single elution peak was observed for purified LCDDeltah, with an estimated molecular weight of 29 kDa, as expected for a monomer. Together, these data suggest that LCDh exist in monomer-dimer equilibrium, and that the C-terminus of P domain is necessary for self-association. We propose that such association is due to interaction between vicinal P domains, which may be of functional relevance for H+ -ATPase in native membranes. We discuss a general dimeric model for P-ATPases with interacting P domains, based on published crystallography and cryo-electron microscopy evidence.
Collapse
Affiliation(s)
- W I Almeida
- Instituto de Bioquímica Médica, Centro de Ciências da Saúde, Universidade Federal do Rio de Janeiro, Rio de Janeiro, 21941-590, Brazil.
| | | | | |
Collapse
|
65
|
Lindahl E, Azuara C, Koehl P, Delarue M. NOMAD-Ref: visualization, deformation and refinement of macromolecular structures based on all-atom normal mode analysis. Nucleic Acids Res 2006; 34:W52-6. [PMID: 16845062 PMCID: PMC1538881 DOI: 10.1093/nar/gkl082] [Citation(s) in RCA: 248] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
Normal mode analysis (NMA) is an efficient way to study collective motions in biomolecules that bypasses the computational costs and many limitations associated with full dynamics simulations. The NOMAD-Ref web server presented here provides tools for online calculation of the normal modes of large molecules (up to 100 000 atoms) maintaining a full all-atom representation of their structures, as well as access to a number of programs that utilize these collective motions for deformation and refinement of biomolecular structures. Applications include the generation of sets of decoys with correct stereochemistry but arbitrary large amplitude movements, the quantification of the overlap between alternative conformations of a molecule, refinement of structures against experimental data, such as X-ray diffraction structure factors or Cryo-EM maps and optimization of docked complexes by modeling receptor/ligand flexibility through normal mode motions. The server can be accessed at the URL .
Collapse
Affiliation(s)
- Erik Lindahl
- Unite de Dynamique Structurale des Macromolécules, URA 2185 du C.N.R.S., Institut Pasteur75015 Paris, France
- Stockholm Bioinformatics Center, Stockholm University106 91 Stockholm, Sweden
| | - Cyril Azuara
- Unite de Dynamique Structurale des Macromolécules, URA 2185 du C.N.R.S., Institut Pasteur75015 Paris, France
| | - Patrice Koehl
- Computer Science Department and Genome Center, University of CaliforniaDavis, CA 95616, USA
| | - Marc Delarue
- Unite de Dynamique Structurale des Macromolécules, URA 2185 du C.N.R.S., Institut Pasteur75015 Paris, France
- To whom correspondence should be addressed. Tel: +33 1 45 68 86 05; Fax: +33 1 45 68 86 04;
| |
Collapse
|
66
|
Abstract
Computational studies of large macromolecular assemblages have come a long way during the past 10 years. With the explosion of computer power and parallel computing, timescales of molecular dynamics simulations have been extended far beyond the hundreds of picoseconds timescale. However, limitations remain for studies of large-scale conformational changes occurring on timescales beyond nanoseconds, especially for large macromolecules. In this review, we describe recent methods based on normal mode analysis that have enabled us to study dynamics on the microsecond timescale for large macromolecules using different levels of coarse graining, from atomically detailed models to those employing only low-resolution structural information. Emerging from such studies is a control principle for robustness in Nature's machines. We discuss this idea in the context of large-scale functional reorganization of the ribosome, virus particles, and the muscle protein myosin.
Collapse
Affiliation(s)
- Florence Tama
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, California 92037, USA
| | | |
Collapse
|
67
|
Jeong JI, Lattman EE, Chirikjian GS. A method for finding candidate conformations for molecular replacement using relative rotation between domains of a known structure. ACTA CRYSTALLOGRAPHICA. SECTION D, BIOLOGICAL CRYSTALLOGRAPHY 2006; 62:398-409. [PMID: 16552141 PMCID: PMC2836325 DOI: 10.1107/s0907444906002204] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/02/2005] [Accepted: 01/18/2006] [Indexed: 11/11/2022]
Abstract
This paper presents a methodology to obtain candidate conformations of multidomain proteins for use in molecular replacement. For each separate domain, the orientational relationship between the template and the target structure is obtained using standard molecular replacement. The orientational relationships of the domains are then used to calculate the relative rotation between the domains in the target conformation by using pose-estimation techniques from the field of robotics and computer vision. With the angle of relative rotation between the domains as a cost function, iterative normal-mode analysis is used to drive the template structure to a candidate conformation that matches the X-ray crystallographic data obtained for the target conformation. The selection of the correct intra-protein domain orientations from among the many spurious maxima in the rotation function (including orientations obtained from domains in symmetry mates rather than within the same copy of the protein) presents a challenge. This problem is resolved by checking R factors of each domain, measuring the absolute value of relative rotation between domains, and evaluating the cost value after each candidate conformation is driven to convergence with iterative NMA. As a validation, the proposed method is applied to three test proteins: ribose-binding protein, lactoferrin and calcium ATPase. In each test case, the orientation and translation of the final candidate conformation in the unit cell are generated correctly from the suggested procedure. The results show that the proposed method can yield viable candidate conformations for use in molecular replacement and can reveal the structural details and pose of the target conformation in the crystallographic unit cell.
Collapse
Affiliation(s)
- Jay I. Jeong
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218 USA
| | - Eaton E. Lattman
- Department of Biophysics, Johns Hopkins University, Baltimore, Maryland 21218
| | - Gregory S. Chirikjian
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218 USA
| |
Collapse
|
68
|
Bonvin AMJJ. Flexible protein–protein docking. Curr Opin Struct Biol 2006; 16:194-200. [PMID: 16488145 DOI: 10.1016/j.sbi.2006.02.002] [Citation(s) in RCA: 216] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2005] [Revised: 01/11/2006] [Accepted: 02/06/2006] [Indexed: 10/25/2022]
Abstract
Predicting the structure of protein-protein complexes using docking approaches is a difficult problem whose major challenges include identifying correct solutions, and properly dealing with molecular flexibility and conformational changes. Flexibility can be addressed at several levels: implicitly, by smoothing the protein surfaces or allowing some degree of interpenetration (soft docking) or by performing multiple docking runs from various conformations (cross or ensemble docking); or explicitly, by allowing sidechain and/or backbone flexibility. Although significant improvements have been achieved in the modeling of sidechains, methods for the explicit inclusion of backbone flexibility in docking are still being developed. A few novel approaches have emerged involving collective degrees of motion, multicopy representations and multibody docking, which should allow larger conformational changes to be modeled.
Collapse
Affiliation(s)
- Alexandre M J J Bonvin
- Bijvoet Center for Biomolecular Research, Science Faculty, Utrecht University, NL-3584 CH, Utrecht, The Netherlands.
| |
Collapse
|
69
|
Chu JW, Voth GA. Coarse-grained modeling of the actin filament derived from atomistic-scale simulations. Biophys J 2006; 90:1572-82. [PMID: 16361345 PMCID: PMC1367308 DOI: 10.1529/biophysj.105.073924] [Citation(s) in RCA: 143] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2005] [Accepted: 11/22/2005] [Indexed: 11/18/2022] Open
Abstract
A coarse-grained (CG) procedure that incorporates the information obtained from all-atom molecular dynamics (MD) simulations is presented and applied to actin filaments (F-actin). This procedure matches the averaged values and fluctuations of the effective internal coordinates that are used to define a CG model to the values extracted from atomistic MD simulations. The fluctuations of effective internal coordinates in a CG model are computed via normal-mode analysis (NMA), and the computed fluctuations are matched with the atomistic MD results in a self-consistent manner. Each actin monomer (G-actin) is coarse-grained into four sites, and each site corresponds to one of the subdomains of G-actin. The potential energy of a CG G-actin contains three bonds, two angles, and one dihedral angle; effective harmonic bonds are used to describe the intermonomer interactions in a CG F-actin. The persistence length of a CG F-actin was found to be sensitive to the cut-off distance of assigning intermonomer bonds. Effective harmonic bonds for a monomer with its third nearest neighboring monomers are found to be necessary to reproduce the values of persistence length obtained from all-atom MD simulations. Compared to the elastic network model, incorporating the information of internal coordinate fluctuations enhances the accuracy and robustness for a CG model to describe the shapes of low-frequency vibrational modes. Combining the fluctuation-matching CG procedure and NMA, the achievable time- and length scales of modeling actin filaments can be greatly enhanced. In particular, a method is described to compute the force-extension curve using the CG model developed in this work and NMA. It was found that F-actin is easily buckled under compressive deformation, and a writhing mode is developed as a result. In addition to the bending and twisting modes, this novel writhing mode of F-actin could also play important roles in the interactions of F-actin with actin-binding proteins and in the force-generation process via polymerization.
Collapse
Affiliation(s)
- Jhih-Wei Chu
- Center for Biophysical Modeling and Simulation and Department of Chemistry, University of Utah, Salt Lake City, Utah 84112-0850, USA
| | | |
Collapse
|
70
|
Nicolay S, Sanejouand YH. Functional modes of proteins are among the most robust. PHYSICAL REVIEW LETTERS 2006; 96:078104. [PMID: 16606146 DOI: 10.1103/physrevlett.96.078104] [Citation(s) in RCA: 68] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/14/2005] [Indexed: 05/08/2023]
Abstract
It is shown that a small subset of modes which are likely to be involved in protein functional motions of large amplitude can be determined by retaining the most robust normal modes obtained using different protein models. This result should prove helpful in the context of several applications proposed recently, like for solving difficult molecular replacement problems or for fitting atomic structures into low-resolution electron density maps. It may also pave the way for the development of methods allowing us to predict such motions accurately.
Collapse
Affiliation(s)
- S Nicolay
- Laboratoire de Physique, Ecole Normale Supérieure, 46 allées d Italie, 69364 Lyon Cedex 07, France
| | | |
Collapse
|
71
|
Bahar I, Rader AJ. Coarse-grained normal mode analysis in structural biology. Curr Opin Struct Biol 2006; 15:586-92. [PMID: 16143512 PMCID: PMC1482533 DOI: 10.1016/j.sbi.2005.08.007] [Citation(s) in RCA: 531] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2005] [Revised: 07/09/2005] [Accepted: 08/24/2005] [Indexed: 10/25/2022]
Abstract
The realization that experimentally observed functional motions of proteins can be predicted by coarse-grained normal mode analysis has renewed interest in applications to structural biology. Notable applications include the prediction of biologically relevant motions of proteins and supramolecular structures driven by their structure-encoded collective dynamics; the refinement of low-resolution structures, including those determined by cryo-electron microscopy; and the identification of conserved dynamic patterns and mechanically key regions within protein families. Additionally, hybrid methods that couple atomic simulations with deformations derived from coarse-grained normal mode analysis are able to sample collective motions beyond the range of conventional molecular dynamics simulations. Such applications have provided great insight into the underlying principles linking protein structures to their dynamics and their dynamics to their functions.
Collapse
Affiliation(s)
- Ivet Bahar
- Department of Computational Biology, University of Pittsburgh, W1043 Biomedical Science Tower, 200 Lothrop Street, Pittsburgh, PA 15261, USA.
| | | |
Collapse
|
72
|
Chennubhotla C, Rader AJ, Yang LW, Bahar I. Elastic network models for understanding biomolecular machinery: from enzymes to supramolecular assemblies. Phys Biol 2005; 2:S173-80. [PMID: 16280623 DOI: 10.1088/1478-3975/2/4/s12] [Citation(s) in RCA: 127] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
With advances in structure genomics, it is now recognized that knowledge of structure alone is insufficient to understand and control the mechanisms of biomolecular function. Additional information in the form of dynamics is needed. As demonstrated in a large number of studies, the machinery of proteins and their complexes can be understood to a good approximation by adopting Gaussian (or elastic) network models (GNM) for simplified normal mode analyses. While this approximation lacks chemical details, it provides us with a means for assessing the collective motions of large structures/assemblies and perform a comparative analysis of a series of proteins, thus providing insights into the mechanical aspects of biomolecular dynamics. In this paper, we discuss recent applications of GNM to a series of enzymes as well as large structures such as the HK97 bacteriophage viral capsids. Understanding the dynamics of large protein structures can be computationally challenging. To this end, we introduce a new approach for building a hierarchical, reduced rank representation of the protein topology and consequently the fluctuation dynamics.
Collapse
Affiliation(s)
- Chakra Chennubhotla
- Department of Computational Biology, School of Medicine, University of Pittsburgh, W1040 BST 200 Lothrop Street, Pittsburgh, PA 15261, USA
| | | | | | | |
Collapse
|
73
|
Llorca O. Introduction to 3D reconstruction of macromolecules using single particle electron microscopy. Acta Pharmacol Sin 2005; 26:1153-64. [PMID: 16174429 DOI: 10.1111/j.1745-7254.2005.00203.x] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
Abstract
Single-particle electron microscopy has now reached maturity, becoming a commonly used method in the examination of macromolecular structure. Using a small amount of purified protein, isolated molecules are observed under the electron microscope and the data collected can be averaged into a 3D reconstruction. Single-particle electron microscopy is an appropriate tool for the analysis of proteins that can only be obtained in modest quantities, like many of the large complexes currently of interest in biomedicine. Whilst the use of electron microscopy expands, new methods are being developed and improved to deal with further challenges, such as reaching higher resolutions and the combination of information at different levels of structural detail. More importantly, present methodology is still not robust enough when studying certain tricky proteins like those displaying extensive conformational flexibility and a great deal of user expertise is required, posing a threat to the consistency of the final structure. This mini review describes a brief outline of the methods currently used in the 3D analysis of macromolecules using single-particle electron microscopy, intended for those first approaching this field. A summary of methods, techniques, software, and some recent work is presented. The spectacular improvements to the technique in recent years, its advantages and limitations compared to other structural methods, and its future developments are discussed.
Collapse
Affiliation(s)
- Oscar Llorca
- Centro de Investigaciones Biologicas, Consejo Superior de Investigaciones Cientificas, Ramiro de Maeztu, 9 Campus Universidad Complutense, 28040 Madrid, Spain.
| |
Collapse
|
74
|
Lindahl E, Delarue M. Refinement of docked protein-ligand and protein-DNA structures using low frequency normal mode amplitude optimization. Nucleic Acids Res 2005; 33:4496-506. [PMID: 16087736 PMCID: PMC1183489 DOI: 10.1093/nar/gki730] [Citation(s) in RCA: 53] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Prediction of structural changes resulting from complex formation, both in ligands and receptors, is an important and unsolved problem in structural biology. In this work, we use all-atom normal modes calculated with the Elastic Network Model as a basis set to model structural flexibility during formation of macromolecular complexes and refine the non-bonded intermolecular energy between the two partners (protein–ligand or protein–DNA) along 5–10 of the lowest frequency normal mode directions. The method handles motions unrelated to the docking transparently by first applying the modes that improve non-bonded energy most and optionally restraining amplitudes; in addition, the method can correct small errors in the ligand position when the first six rigid-body modes are switched on. For a test set of six protein receptors that show an open-to-close transition when binding small ligands, our refinement scheme reduces the protein coordinate cRMS by 0.3–3.2 Å. For two test cases of DNA structures interacting with proteins, the program correctly refines the docked B-DNA starting form into the expected bent DNA, reducing the DNA cRMS from 8.4 to 4.8 Å and from 8.7 to 5.4 Å, respectively. A public web server implementation of the refinement method is available at .
Collapse
Affiliation(s)
- Erik Lindahl
- Unité de Biochimie Structurale, URA 2185 du CNRS, Institut Pasteur25 Rue du Dr Roux, F-75015 Paris, France
- Stockholm Bioinformatics Center, Stockholm UniversitySE-17156 Stockholm, Sweden
| | - Marc Delarue
- Unité de Biochimie Structurale, URA 2185 du CNRS, Institut Pasteur25 Rue du Dr Roux, F-75015 Paris, France
- To whom correspondence should be addressed. Tel: +33 1 45 68 8605; Fax: +33 1 45 68 8604;
| |
Collapse
|