4 research outputs found
Explicit factorization of external coordinates in constrained Statistical Mechanics models
If a macromolecule is described by curvilinear coordinates or rigid
constraints are imposed, the equilibrium probability density that must be
sampled in Monte Carlo simulations includes the determinants of different
mass-metric tensors. In this work, we explicitly write the determinant of the
mass-metric tensor G and of the reduced mass-metric tensor g, for any molecule,
general internal coordinates and arbitrary constraints, as a product of two
functions; one depending only on the external coordinates that describe the
overall translation and rotation of the system, and the other only on the
internal coordinates. This work extends previous results in the literature,
proving with full generality that one may integrate out the external
coordinates and perform Monte Carlo simulations in the internal conformational
space of macromolecules. In addition, we give a general mathematical argument
showing that the factorization is a consequence of the symmetries of the metric
tensors involved. Finally, the determinant of the mass-metric tensor G is
computed explicitly in a set of curvilinear coordinates specially well-suited
for general branched molecules.Comment: 22 pages, 2 figures, LaTeX, AMSTeX. v2: Introduccion slightly
extended. Version in arXiv is slightly larger than the published on
An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints
In this work, we introduce an algorithm to compute the derivatives of
physical observables along the constrained subspace when flexible constraints
are imposed on the system (i.e., constraints in which the hard coordinates are
fixed to configuration-dependent values). The presented scheme is exact, it
does not contain any tunable parameter, and it only requires the calculation
and inversion of a sub-block of the Hessian matrix of second derivatives of the
function through which the constraints are defined. We also present a practical
application to the case in which the sought observables are the Euclidean
coordinates of complex molecular systems, and the function whose minimization
defines the constraints is the potential energy. Finally, and in order to
validate the method, which, as far as we are aware, is the first of its kind in
the literature, we compare it to the natural and straightforward
finite-differences approach in three molecules of biological relevance:
methanol, N-methyl-acetamide and a tri-glycine peptideComment: 13 pages, 8 figures, published versio
Quantum mechanical calculation of the effects of stiff and rigid constraints in the conformational equilibrium of the Alanine dipeptide
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.Comment: 37 pages, 4 figures, LaTeX, BibTeX, AMSTe