Por favor, use este identificador para citar o enlazar a este item:
http://hdl.handle.net/10261/41143
COMPARTIR / EXPORTAR:
SHARE BASE | |
Visualizar otros formatos: MARC | Dublin Core | RDF | ORE | MODS | METS | DIDL | DATACITE | |
Título: | Exact and efficient calculation of lagrange multipliers in biological polymers with constrained bond lengths and bond angles: Proteins and nucleic acids as example cases |
Autor: | García-Risueño, Pablo CSIC; Echenique, Pablo CSIC ORCID; Alonso, José Luis CSIC ORCID | Palabras clave: | Constraints Lagrange multipliers Bond angles proteins Nucleic acids |
Fecha de publicación: | 8-ago-2011 | Editor: | John Wiley & Sons | Citación: | Journal of Computational Chemistry 32(14): 3039-3046 (2011) | Resumen: | To accelerate molecular dynamics simulations, it is common to impose holonomic constraints on the hardest degrees of freedom. In this way, the time step used to integrate the equations of motion can be increased, thereby allowing longer total simulation times. The imposition of such constraints results in an aditional set of N equations (the equations of constraint) and unknowns (their associated Lagrange multipliers), whose solution is closely related to any algorithm implementing the constraints in Euclidean coordinates. In this work, it is shown that, due to the essentially linear structure of typical biological polymers the algebraic equations that need to be solved involve a matrix which is not only sparse, but also banded if the constraints are indexed in a skilful way. This allows the Lagrange multipliers to be obtained through a noniterative procedure, which can be considered exact up to machine precision, and which takes O(N) operations, instead of the usual O(N3) for generic molecular systems. We develop the formalism, and describe the appropriate indexing for a number of model molecules. Finally, we provide a numerical example of the technique in a series of polyalanine peptides of different lengths. Although a use of the Lagrange multipliers without any modification in the solution of the underlying ordinary differential equations yields unstable integration algorithms, the central role of these quantities makes their efficient calculation useful for the improvement of methods that correctly enforce the exact satisfaction of the constraints at each time step. We provide several examples of this. | Descripción: | This is a preprint version of a published work at http://onlinelibrary.wiley.com/doi/10.1002/jcc.21885/abstract | Versión del editor: | http://onlinelibrary.wiley.com/doi/10.1002/jcc.21885/abstract | URI: | http://hdl.handle.net/10261/41143 | ISSN: | 10.1002/jcc.21885 |
Aparece en las colecciones: | (IQF) Artículos |
Ficheros en este ítem:
Fichero | Descripción | Tamaño | Formato | |
---|---|---|---|---|
2011_LagrangeMultipliers.pdf | 240,15 kB | Adobe PDF | Visualizar/Abrir |
CORE Recommender
Page view(s)
348
checked on 22-abr-2024
Download(s)
253
checked on 22-abr-2024
Google ScholarTM
Check
NOTA: Los ítems de Digital.CSIC están protegidos por copyright, con todos los derechos reservados, a menos que se indique lo contrario.