2,405 research outputs found
A Central Partition of Molecular Conformational Space.III. Combinatorial Determination of the Volume Spanned by a Molecular System
In the first work of this series [physics/0204035] it was shown that the
conformational space of a molecule could be described to a fair degree of
accuracy by means of a central hyperplane arrangement. The hyperplanes divide
the espace into a hierarchical set of cells that can be encoded by the face
lattice poset of the arrangement. The model however, lacked explicit rotational
symmetry which made impossible to distinguish rotated structures in
conformational space. This problem was solved in a second work
[physics/0404052] by sorting the elementary 3D components of the molecular
system into a set of morphological classes that can be properly oriented in a
standard 3D reference frame. This also made possible to find a solution to the
problem that is being adressed in the present work: for a molecular system
immersed in a heat bath we want to enumerate the subset of cells in
conformational space that are visited by the molecule in its thermal wandering.
If each visited cell is a vertex on a graph with edges to the adjacent cells,
here it is explained how such graph can be built
A Central Partition of Molecular Conformational Space. IV. Extracting information from the graph of cells
In previous works [physics/0204035, physics/0404052, physics/0509126] a
procedure was described for dividing the -dimensional
conformational space of a molecular system into a number of discrete cells,
this partition allowed the building of a combinatorial structure from data
sampled in molecular dynamics trajectories: the graph of cells, that encodes
the set of cells in conformational space that are visited by the system in its
thermal wandering. Here we outline a set of procedures for extracting useful
information from this structure: 1st) interesting regions in the volume
occupied by the system in conformational space can be bounded by a polyhedral
cone whose faces are determined empirically from a set of relations between the
coordinates of the molecule, 2nd) it is also shown that this cone can be
decomposed into a hierarchical set of smaller cones, 3rd) the set of cells in a
cone can be encoded by a simple combinatorial sequence.Comment: added an intrduction and reference
Ab initio study of alanine polypeptide chains twisting
We have investigated the potential energy surfaces for alanine chains
consisting of three and six amino acids. For these molecules we have calculated
potential energy surfaces as a function of the Ramachandran angles Phi and Psi,
which are widely used for the characterization of the polypeptide chains. These
particular degrees of freedom are essential for the characterization of
proteins folding process. Calculations have been carried out within ab initio
theoretical framework based on the density functional theory and accounting for
all the electrons in the system. We have determined stable conformations and
calculated the energy barriers for transitions between them. Using a
thermodynamic approach, we have estimated the times of characteristic
transitions between these conformations. The results of our calculations have
been compared with those obtained by other theoretical methods and with the
available experimental data extracted from the Protein Data Base. This
comparison demonstrates a reasonable correspondence of the most prominent
minima on the calculated potential energy surfaces to the experimentally
measured angles Phi and Psi for alanine chains appearing in native proteins. We
have also investigated the influence of the secondary structure of polypeptide
chains on the formation of the potential energy landscape. This analysis has
been performed for the sheet and the helix conformations of chains of six amino
acids.Comment: 24 pages, 10 figure
CHARMM36m: An improved force field for folded and intrinsically disordered proteins.
The all-atom additive CHARMM36 protein force field is widely used in molecular modeling and simulations. We present its refinement, CHARMM36m (http://mackerell.umaryland.edu/charmm_ff.shtml), with improved accuracy in generating polypeptide backbone conformational ensembles for intrinsically disordered peptides and proteins
Maximum Flux Transition Paths of Conformational Change
Given two metastable states A and B of a biomolecular system, the problem is
to calculate the likely paths of the transition from A to B. Such a calculation
is more informative and more manageable if done for a reduced set of collective
variables chosen so that paths cluster in collective variable space. The
computational task becomes that of computing the "center" of such a cluster. A
good way to define the center employs the concept of a committor, whose value
at a point in collective variable space is the probability that a trajectory at
that point will reach B before A. The committor "foliates" the transition
region into a set of isocommittors. The maximum flux transition path is defined
as a path that crosses each isocommittor at a point which (locally) has the
highest crossing rate of distinct reactive trajectories. (This path is
different from that of the MaxFlux method of Huo and Straub.) It is argued that
such a path is nearer to an ideal path than others that have been proposed with
the possible exception of the finite-temperature string method path. To make
the calculation tractable, three approximations are introduced, yielding a path
that is the solution of a nonsingular two-point boundary-value problem. For
such a problem, one can construct a simple and robust algorithm. One such
algorithm and its performance is discussed.Comment: 7 figure
Predicting Transcription Factor Specificity with All-Atom Models
The binding of a transcription factor (TF) to a DNA operator site can
initiate or repress the expression of a gene. Computational prediction of sites
recognized by a TF has traditionally relied upon knowledge of several cognate
sites, rather than an ab initio approach. Here, we examine the possibility of
using structure-based energy calculations that require no knowledge of bound
sites but rather start with the structure of a protein-DNA complex. We study
the PurR E. coli TF, and explore to which extent atomistic models of
protein-DNA complexes can be used to distinguish between cognate and
non-cognate DNA sites. Particular emphasis is placed on systematic evaluation
of this approach by comparing its performance with bioinformatic methods, by
testing it against random decoys and sites of homologous TFs. We also examine a
set of experimental mutations in both DNA and the protein. Using our explicit
estimates of energy, we show that the specificity for PurR is dominated by
direct protein-DNA interactions, and weakly influenced by bending of DNA.Comment: 26 pages, 3 figure
Secondary-Structure Design of Proteins by a Backbone Torsion Energy
We propose a new backbone-torsion-energy term in the force field for protein
systems. This torsion-energy term is represented by a double Fourier series in
two variables, the backbone dihedral angles phi and psi. It gives a natural
representation of the torsion energy in the Ramachandran space in the sense
that any two-dimensional energy surface periodic in both phi and psi can be
expanded by the double Fourier series. We can then easily control
secondary-structure-forming tendencies by modifying the torsion-energy surface.
For instance, we can increase/decrease the alpha-helix-forming-tendencies by
lowering/raising the torsion-energy surface in the alpha-helix region and
likewise increase/decrease the beta-sheet-forming tendencies by
lowering/raising the surface in the beta-sheet region in the Ramachandran
space. We applied our approach to AMBER parm94 and AMBER parm96 force fields
and demonstrated that our modifications of the torsion-energy terms resulted in
the expected changes of secondary-structure-forming-tendencies by performing
folding simulations of alpha-helical and beta-hairpin peptides.Comment: 13 pages, (Revtex4), 5 figure
Excitons in a Photosynthetic Light-Harvesting System: A Combined Molecular Dynamics/Quantum Chemistry and Polaron Model Study
The dynamics of pigment-pigment and pigment-protein interactions in
light-harvesting complexes is studied with a novel approach which combines
molecular dynamics (MD) simulations with quantum chemistry (QC) calculations.
The MD simulations of an LH-II complex, solvated and embedded in a lipid
bilayer at physiological conditions (with total system size of 87,055 atoms)
revealed a pathway of a water molecule into the B800 binding site, as well as
increased dimerization within the B850 BChl ring, as compared to the
dimerization found for the crystal structure. The fluctuations of pigment (B850
BChl) excitation energies, as a function of time, were determined via ab initio
QC calculations based on the geometries that emerged from the MD simulations.
From the results of these calculations we constructed a time-dependent
Hamiltonian of the B850 exciton system from which we determined the linear
absorption spectrum. Finally, a polaron model is introduced to describe quantum
mechanically both the excitonic and vibrational (phonon) degrees of freedom.
The exciton-phonon coupling that enters into the polaron model, and the
corresponding phonon spectral function are derived from the MD/QC simulations.
It is demonstrated that, in the framework of the polaron model, the absorption
spectrum of the B850 excitons can be calculated from the autocorrelation
function of the excitation energies of individual BChls, which is readily
available from the combined MD/QC simulations. The obtained result is in good
agreement with the experimentally measured absorption spectrum.Comment: REVTeX3.1, 23 pages, 13 (EPS) figures included. A high quality PDF
file of the paper is available at
http://www.ks.uiuc.edu/Publications/Papers/PDF/DAMJ2001/DAMJ2001.pd
Molecular dynamics study of accelerated ion-induced shock waves in biological media
We present a molecular dynamics study of the effects of carbon- and iron-ion induced shock waves in DNA duplexes in liquid water. We use the CHARMM force field implemented within the MBN Explorer simulation package to optimize and equilibrate DNA duplexes in liquid water boxes of different sizes and shapes. The translational and vibrational degrees of freedom of water molecules are excited according to the energy deposited by the ions and the subsequent shock waves in liquid water are simulated. The pressure waves generated are studied and compared with an analytical hydrodynamics model which serves as a benchmark for evaluating the suitability of the simulation boxes. The energy deposition in the DNA backbone bonds is also monitored as an estimation of biological damage, something which is not possible with the analytical model
- …