title: Monte Carlo Sampling with Linear Inverse Kinematics for Simulation of Protein Flexible Regions.
abstract: A Monte Carlo linear inverse-kinematics method for the simulation of protein chains with fixed ends is introduced. It includes backbone bond-angle bending and simultaneous loop and ring closure to allow full proline ring flexibility. An obstacle to linear null-space methods is the eventual drift of the end group. Maintenance of the end group at its initial position by occasional reset is performed in a way that is consistent with the overall methodology and minimally disruptive to the current conformation. The implementation permitted multiple rigid regions within the chain, enabling the simulation of domain movements where domains are rigid bodies connected by flexible interdomain regions. The method was tested on polyalanine, polyglycine, loop 6 of triosephosphate isomerase, and glutamine binding protein. Simulations of glutamine binding protein, where only 11 of the 226 residues at the interdomain bending regions were flexible, accurately reproduced the experimentally determined domain movement.
authors.affiliation: D'Arcy Thompson Centre for Computational Biology, School of Computing Sciences, University of East Anglia , Norwich NR4 7TJ, U.K.
authors.forename: Steven
authors.initials: S
authors.lastname: Hayward
authors.affiliation: Institute of Molecular and Cellular Biosciences, The University of Tokyo , 1-1-1 Yayoi, Bunkyo, Tokyo 113-0032, Japan.
authors.forename: Akio
authors.initials: A
authors.lastname: Kitao
mesh_terms: D001696:Biomechanical Phenomena; D009010:Monte Carlo Method; D010455:Peptides / Q000737:chemistry / Q000378:metabolism; D017434:Protein Structure, Tertiary; D011506:Proteins / Q000737:chemistry* / Q000378:metabolism; D014305:Triose-Phosphate Isomerase / Q000737:chemistry / Q000378:metabolism
chemical_list: D010455:Peptides; D011506:Proteins; C019529:polyalanine; C011080:polyglycine; D014305:Triose-Phosphate Isomerase