Fractional Brownian dynamics in proteins.

Correlation functions describing relaxation processes in proteins and other complex molecular systems are known to exhibit a nonexponential decay. The simulation study presented here shows that fractional Brownian dynamics is a good model for the internal dynamics of a lysozyme molecule in solution. We show that both the dynamic structure factor and the associated memory function fit well the corresponding analytical functions calculated from the model. The numerical analysis is based on autoregressive modeling of time series.

nMoldyn: a program package for a neutron scattering oriented analysis of molecular dynamics simulations.

We present a new implementation of the program nMoldyn, which has been developed for the computation and decomposition of neutron scattering intensities from Molecular Dynamics trajectories (Comp. Phys. Commun 1995, 91, 191-214). The new implementation extends the functionality of the original version, provides a much more convenient user interface (both graphical/interactive and batch), and can be used as a tool set for implementing new analysis modules. This was made possible by the use of a high-level language, Python, and of modern object-oriented programming techniques. The quantities that can be calculated by nMoldyn are the mean-square displacement, the velocity autocorrelation function as well as its Fourier transform (the density of states) and its memory function, the angular velocity autocorrelation function and its Fourier transform, the reorientational correlation function, and several functions specific to neutron scattering: the coherent and incoherent intermediate scattering functions with their Fourier transforms, the memory function of the coherent scattering function, and the elastic incoherent structure factor. The possibility to compute memory function is a new and powerful feature that allows to relate simulation results to theoretical studies.

Influence of molecular flexibility on DNA radiosensitivity: a simulation study.

Radiation damage in DNA is caused mainly by hydroxyl radicals which are generated by ionizing radiation in water and removing hydrogen atoms from the DNA chain. This damage affects certain nucleotide sequences more than others due to differences in the local structure of the DNA chains. This sequence dependence has been analyzed experimentally and calculated theoretically for a rigid DNA model. In this paper we take into account the flexibility of the DNA chain and show how it modifies the strand breakage probabilities. We use a simple harmonic model for DNA flexibility which permits the study of a long (68 base pair) fragment with modest computational effort. The essential influence of flexibility is an increased breakage probability towards the ends of the fragment, which can also be identified in the experimental data.

Theoretical study of the conformation of the H-protein lipoamide arm as a function of its terminal group.

The glycine decarboxylase complex consists of four different proteins (the L-, P-, H-, and T-proteins). The H-protein plays a central role in communication among the other enzymes, as its lipoamide arm interacts successively with each of the components of the complex. The crystal structures of two states of the H-protein have been resolved: the oxidized form, Hox at 2 A and the methylamine-loaded form, Hmet at 2.2 A. However, the position of the arm for the reduced form, Hred, is still unknown. We have performed numerical free-energy calculations in order to better understand the differences in the structures and to elucidate the conformation of the arm in Hred. The results of the simulations are in agreement with the crystallographic results, as the minima of the free energy surface for Hox and Hmet correspond to the crystal structures. For Hred, we observe a single minimum in which the arm is on the surface of the H-protein, close to its position in the Hox structure. In all of our simulations, the lower, lysine portion of the arm remains bound to the protein, which substantially reduces the number of accessible arm configurations. An analysis of the stability of Hmet in the cavity shows that electrostatic interactions are crucial for locking the arm in the bottom of the cavity, especially near Glu14. In addition, the analysis shows that there is a water molecule, also observed in the crystallographic structure, that binds to the arm's terminal NH3+ group and helps to fix it in the cavity. In conclusion, because of the close agreement of the results of our calculations with the available experimental evidence, we are able to suggest a structural basis for the observed behavior. Proteins 1999;36: 228-237.

Tertiary and quaternary conformational changes in aspartate transcarbamylase: a normal mode study.

Aspartate transcarbamylase (ATCase) initiates the pyrimidine biosynthetic pathway in Escherichia coli. Binding of aspartate to this allosteric enzyme induces a cooperative transition between the tensed (T) and relaxed (R) states of the enzyme which involves large quaternary and tertiary rearrangements. The mechanisms of the transmission of the regulatory signal to the active site (60 A away) and that of the cooperative transition are not known in detail, although a large number of single, double, and triple site-specific mutants and chimeric forms of ATCase have been obtained and kinetically characterized. A previous analysis of the very low-frequency normal modes of both the T and R state structures of ATCase identified some of the large-amplitude motions mediating the intertrimer elongation and rotation that occur during the cooperative transition (Thomas et al., J. Mol. Biol. 257:1070-1087, 1996; Thomas et al., J. Mol. Biol. 261:490-506, 1996). As a complement to that study, the deformation of the quaternary and tertiary structure of ATCase by normal modes below 5 cm(-1) is investigated in this article. The ability of the modes to reproduce the domain motions occurring during the transition is analyzed, with special attention to the interdomain closure in the catalytic chain, which has been shown to be critical for homotropic cooperativity. The calculations show a coupling between the quaternary motions and more localized motions involving specific residues. The particular dynamic behavior of these residues is examined in the light of biochemical results to obtain insights into their role in the transmission of the allosteric signal.

Analysis of domain motions in large proteins.

We present a new approach for determining dynamical domains in large proteins, either based on a comparison of different experimental structures, or on a simplified normal mode calculation for a single conformation. In a first step, a deformation measure is evaluated for all residues in the protein; a high deformation indicates highly flexible interdomain regions. The sufficiently rigid parts of the protein are then classified into rigid domains and low-deformation interdomain regions on the basis of their global motion. We demonstrate the techniques on three proteins: citrate synthase, which has been the subject of earlier domain analyses, HIV-1 reverse transcriptase, which has a rather complex domain structure, and aspartate transcarbamylase as an example of a very large protein. These examples show that the comparison of conformations and the normal mode analysis lead to essentially the same domain identification, except for cases where the experimental conformations differ by the presence of a large ligand, such as a DNA strand. Normal mode analysis has the advantage of requiring only one experimental structure and of providing a more detailed picture of domain movements, e.g. the splitting of domains into subdomains at higher frequencies.

Analysis of domain motions by approximate normal mode calculations.

The identification of dynamical domains in proteins and the description of the low-frequency domain motions are one of the important applications of numerical simulation techniques. The application of these techniques to large proteins requires a substantial computational effort and therefore cannot be performed routinely, if at all. This article shows how physically motivated approximations permit the calculation of low-frequency normal modes in a few minutes on standard desktop computers. The technique is based on the observation that the low-frequency modes, which describe domain motions, are independent of force field details and can be obtained with simplified mechanical models. These models also provide a useful measure for rigidity in proteins, allowing the identification of quasi-rigid domains. The methods are validated by application to three well-studied proteins, crambin, lysozyme, and ATCase. In addition to being useful techniques for studying domain motions, the success of the approximations provides new insight into the relevance of normal mode calculations and the nature of the potential energy surface of proteins.