Theoretical determination of conformational paths in citrate synthase
- 1 May 1992
- journal article
- research article
- Published by Wiley in Biopolymers
- Vol. 32 (5) , 561-574
- https://doi.org/10.1002/bip.360320512
Abstract
Two methods are developed for the theoretical determination of a conformational path between two well-documented forms, a closed form and the open form [Remington et al. (1982) J. Mol. Biol. 158, 111–152] of pig heart citrate synthase, a dimeric enzyme of 2 × 437 residues. The first method uses the minimization of the sum of the potential energies at a set of equidistant points, according to Elber and Karplus [(1987) Chem. Phys. Lett. 139, 375–380]. The initialization of the algorithm is modified to account for large-angle rotations of many groups by performing the interpolations in the space of internal polar coordinates of a set of generalized Jacobi vectors earlier introduced by Durup [(1991) J. Phys. Chem. 95, 1817–1829] and by carefully testing all choices of directions of rotation for determining the initialized midpoint between the known forms. The path includes intermediate points, created by successive splittings of each interval into two equal parts, with a partial energy minimization performed after each splitting. The minimization encounters the well-known local-minima problem, which here is handled by low-temperature molecular dynamics annealing. It is shown that the best ratio of potential energy decrease to rms deviation is achieved by running the dynamics at 50 K, as compared to 100 K and above. The main character of the path obtained is the occurrence of strong to-and-fro variations of some dihedral angles at specific stages along the path. The second method, which we name directed dynamics, uses only low-temperature molecular dynamics simulations by starting trajectories from each of the two known forms with initial velocities directed toward the other one. The procedure is iterated by restarting trajectory pairs after the points of closest approach of the preceding pair. The two half-paths thus built eventually meet after 70 iterations. This method provides a second path with strong similarities, as well as some differences, with respect to the path obtained by the first method.Keywords
This publication has 91 references indexed in Scilit:
- Molecular dynamics simulations with experimental restraintsAccounts of Chemical Research, 1991
- Refinement of the influenza virus hemagglutinin by simulated annealingJournal of Molecular Biology, 1990
- Motions of an α‐helical polypeptide: Comparison of molecular and harmonic dynamicsBiopolymers, 1990
- Revision of the methodology in enzyme kinetics: A fractal approachJournal of Theoretical Biology, 1989
- Calculations of electrostatic properties in proteinsJournal of Molecular Biology, 1987
- Crystal structure analysis and molecular model of a complex of citrate synthase with oxaloacetate and S-acetonyl-coenzyme AJournal of Molecular Biology, 1984
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983
- Molecular dynamics of ferrocytochrome cJournal of Molecular Biology, 1981
- The protein data bank: A computer-based archival file for macromolecular structuresJournal of Molecular Biology, 1977
- A simplified representation of protein conformations for rapid simulation of protein foldingJournal of Molecular Biology, 1976