Molecular dynamics simulations of N-terminal peptides from a nucleotide binding protein

Abstract

Molecular dynamics (MD) simulations of N-terminal peptides from lactate dehydrogenase (LDH) with increasing length and individual secondary structure elements were used to study their stability in relation to folding. Ten simulations of 1-2 ns of different peptides in water starting from the coordinates of the crystal structure were performed. The stability of the peptides was compared qualitatively by analyzing the root mean square deviation (RMSD) from the crystal structure, radius of gyration, secondary and tertiary structure, and solvent accessible surface area. In agreement with earlier MD studies, relatively short (< 15 amino acids) peptides containing individual secondary structure elements were generally found to be unstable; the hydrophobic a,-helix of the nucleotide binding fold displayed a significantly higher stability, however. Our simulations further showed that the first pap supersecondary unit of the characteristic dinucleotide binding fold (Rossmann fold) of LDH is somewhat more stable than other units of similar length and that the a,-helix, which unfolds by itself, is stabilized by binding to this unit. This finding suggests that the first pap unit could function as a n N-terminal folding nucleus, upon which the remainder of the polypeptide chain can be assembled. Indeed, simulations with longer units ( p a f k~ a n d papapp) showed that all structural elements of these units are rather stable. The outcome of our studies is in line with suggestions that folding of the N-terminal portion of LDH in vivo can be a cotranslational process that takes place during the ribosomal peptide synthesis.

GROMACS: Fast, flexible, and free

“…The topic of protein folding has long been one of the challenges pursued by the GROMACS authors 130–133. More recently, longer simulations have become possible, in the range of 50 ns and longer.…”

Molecular dynamics simulation of the antimicrobial salivary peptide histatin‐5 in water and in trifluoroethanol: a microscopic description of the water destructuring effect

“…Histatin‐5 had a total final charge +5 e and the electrostatic interaction was treated considering a cut‐off distance of 12.0 Å. A similar cut‐off distance was used in recent MD studies dealing with peptides of similar length, or with protein having a comparable or higher charge (22–25), and has been shown to give reliable results. A cut‐off distance of 12.0 Å was chosen for the van der Waal's interactions.…”

Activation mechanism of Yarrowia lipolytica lipase immobilized on carbon nanotubes

“…42 ED analysis can identify functionally relevant displacements of groups of residues by concentrating on a subset of the principal eigenvalues and eigenvectors of the residue pair covariance matrix calculated from MD simulations. 43 The trajectories of the lipase from the MD simulations were projected onto the first eigenvector 42,44 and the residue-based root mean fluctuations along this main direction were calculated (Fig. 2).…”