Quantum mechanical calculation of the effects of stiff and rigid constraints in the conformational equilibrium of the Alanine dipeptide

Abstract

If constraints are imposed on a macromolecule, two inequivalent classical models may be used: the stiff and the rigid one. This work studies the effects of such constraints on the Conformational Equilibrium Distribution (CED) of the model dipeptide HCO-L-Ala-NH2 without any simplifying assumption. We use ab initio Quantum Mechanics calculations including electron correlation at the MP2 level to describe the system, and we measure the conformational dependence of all the correcting terms to the naive CED based in the Potential Energy Surface (PES) that appear when the constraints are considered. These terms are related to mass-metric tensors determinants and also occur in the Fixman's compensating potential. We show that some of the corrections are non-negligible if one is interested in the whole Ramachandran space. On the other hand, if only the energetically lower region, containing the principal secondary structure elements, is assumed to be relevant, then, all correcting terms may be neglected up to peptides of considerable length. This is the first time, as far as we know, that the analysis of the conformational dependence of these correcting terms is performed in a relevant biomolecule with a realistic potential energy function.

arXiv comment37 pages, 4 figures, LaTeX, BibTeX, AMSTeX

References in corpus (4)

• Definition of Systematic, Approximately Separable and Modular Internal Coordinates (SASMIC) for macromolecular simulation
• Robust ab initio calculation of condensed matter: transparent convergence through semicardinal multiresolution analysis
• A physically meaningful method for the comparison of potential energy functions
• Do theoretical physicists care about the protein-folding problem?

Cited by in corpus (6)

• The canonical equilibrium of constrained molecular models
• An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints
• A Quantum-mechanical Approach for Constrained Macromolecular Chains
• Efficient model chemistries for peptides. II. Basis set convergence in the B3LYP method
• From constrained stochastic processes to the nonlinear sigma model. Two old problems revisited
• Dynamical aspects of inextensible chains