Molecular dynamics simulation of the heart type fatty acid binding protein in a crystal environment
- 2 July 2021
- journal article
- research article
- Published by Informa UK Limited in Journal of Biomolecular Structure and Dynamics
- Vol. 39 (10), 3459-3468
- https://doi.org/10.1080/07391102.2020.1773315
Abstract
Crystallographic data comes from a space-time average over all the unit cells within the crystal, so dynamic phenomena do not contribute significantly to the diffraction data. Many efforts have been made to reconstitute the movement of the macromolecules and explore the microstates that the confined proteins can adopt in the crystalline network. We explored different strategies to simulate a heart fatty acid binding protein (H-FABP) crystal by means of Molecular Dynamics (MD) simulations. We evaluate the effect of introducing restraints according to experimental isotropic B-factors and we analyzed the H-FABP motions in the crystal using Principal Component Analysis (PCA), isotropic and anisotropic B-factors. We compared the behavior of the protein simulated in the crystal confinement versus in solution, and we observed the effect of that confinement in the mobility of the protein residues. Restraining one-third of C alpha atoms based on experimental B-factors produce lower B-factors than simulations without restraints, showing that the position restraint of the atoms with the lowest experimental B-factor is a good strategy to maintain the geometry of the crystal with an obvious decrease in the degrees of motion of the protein. PCA shows that, as position restraint reduces the conformational space explored by the system, the motion of the crystal is better recovered, for an essential subspace of the same size, in the simulations without restraints. Restraining only one C alpha seems to be a good balance between giving flexibility to the system and preserving its structure.Keywords
Funding Information
- Consejo Nacional de Investigaciones Científicas y Técnicas
- Universidad Nacional de La Plata
This publication has 49 references indexed in Scilit:
- GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputersSoftwareX, 2015
- The future of the protein data bankPeptide Science, 2012
- Optimization of parameters for molecular dynamics simulation using smooth particle‐mesh Ewald in GROMACS 4.5Journal of Computational Chemistry, 2011
- High temperatures enhance cooperative motions between CBM and catalytic domains of a thermostable cellulase: mechanism insights from essential dynamicsPhysical Chemistry Chemical Physics, 2011
- A large data set comparison of protein structures determined by crystallography and NMR: Statistical test for structural differences and the effect of crystal packingProteins-Structure Function and Bioinformatics, 2007
- Canonical sampling through velocity rescalingThe Journal of Chemical Physics, 2007
- The Protein Data BankNucleic Acids Research, 2000
- On the convergence of the conformational coordinates basis set obtained by the essential dynamics analysis of proteins' molecular dynamics simulationsProteins-Structure Function and Bioinformatics, 1999
- Essential dynamics of proteinsProteins-Structure Function and Bioinformatics, 1993
- The missing term in effective pair potentialsThe Journal of Physical Chemistry, 1987