Performance of the Cell processor for biomolecular simulations

The new Cell processor represents a turning point for computing intensive applications. Here, I show that for molecular dynamics it is possible to reach an impressive sustained performance in excess of 30 Gflops with a peak of 45 Gflops for the non-bonded force calculations, over one order of magnitude faster than a single core standard processor

Theoretical prediction of spectral and optical properties of bacteriochlorophylls in thermally disordered LH2 antenna complexes

A general approach for calculating spectral and optical properties of pigment-protein complexes of known atomic structure is presented. The method, that combines molecular dynamics simulations, quantum chemistry calculations and statistical mechanical modeling, is demonstrated by calculating the absorption and circular dichroism spectra of the B800-B850 BChls of the LH2 antenna complex from Rs. molischianum at room temperature. The calculated spectra are found to be in good agreement with the available experimental results. The calculations reveal that the broadening of the B800 band is mainly caused by the interactions with the polar protein environment, while the broadening of the B850 band is due to the excitonic interactions. Since it contains no fitting parameters, in principle, the proposed method can be used to predict optical spectra of arbitrary pigment-protein complexes of known structure.Comment: ReVTeX4, 11 pages, 9 figures, submitted to J. Chem. Phy

Cation‐π Interactions between Methylated Ammonium Groups and Tryptophan in the CHARMM36 Additive Force Field

Cation-π interactions between tryptophan and choline or trimethylated lysines are vital for many biological processes. The performance of the additive CHARMM36 force field against target quantum mechanical data is shown to reproduce QM equilibrium geometries but required modified Lennard-Jones potentials to accurately reproduce the QM interaction energies. The modified parameter set allows accurate modeling, including free energies, of cation-π indole-choline and indole-trimethylated lysines interactions relevant for protein–ligand, protein–membrane, and protein–protein interfaces.acceptedVersio

Cooperative binding of DNA and CBFβ to the Runt domain of the CBFα studied via MD simulations

The Runt domain (RD) is the DNA-binding region of the Runx genes. A related protein, known as core binding factor β (CBFβ) also binds to the RD to enhance RD–DNA interaction by 6- to 10-fold. Here, we report results from molecular dynamics (MD) simulations of RD alone, as a dimer in complexes with DNA and CBFβ and in a ternary complex with DNA and CBFβ. Consistent with the experimental findings, in the presence of CBFβ the estimated free energy of binding of RD to the DNA is more favorable, which is shown to be due to more favorable intermolecular interactions and desolvation contributions. Also contributing to the enhanced binding are favorable intramolecular interactions between the ‘wing’ residues (RD residues 139–145) and the ‘wing1’ residues (RD residues 104–116). The simulation studies also indicate that the RD–CBFβ binding is more favorable in the presence of DNA due to a more favorable RD–CBFβ interaction energy. In addition, it is predicted that long-range interactions involving ionic residues contribute to binding cooperativity. Results from the MD calculations are used to interpret a variety of experimental mutagenesis data. A novel role for RD Glu116 to the RD–CBFβ interaction is predicted

Intrinsic Conformational Energetics Associated with the Glycosyl Torsion in DNA: A Quantum Mechanical Study

AbstractThe glycosyl torsion (χ) in nucleic acids has long been recognized to be a major determinant of their conformational properties. χ torsional energetics were systematically mapped in deoxyribonucleosides using high-level quantum mechanical methods, for north and south sugar puckers and with γ in the g+ and trans conformations. In all cases, the syn conformation is found higher in energy than the anti. When γ is changed from g+ to trans, the anti orientation of the base is strongly destabilized, and the energy difference and barrier between anti and syn are significantly decreased. The barrier between anti and syn in deoxyribonucleosides is found to be less than 10kcal/mol and tends to be lower with purines than with pyrimidines. With γ=g+/χ=anti, a south sugar yields a significantly broader energy well than a north sugar with no energy barrier between χ values typical of A or B DNA. Contrary to the prevailing view, the syn orientation is not more stable with south puckers than with north puckers. The syn conformation is significantly more energetically accessible with guanine than with adenine in 5-nucleotides but not in nucleosides. Analysis of nucleic acid crystal structures shows that γ=trans/χ=anti is a minor but not negligible conformation. Overall, χ appears to be a very malleable structural parameter with the experimental χ distributions reflecting, to a large extent, the associated intrinsic torsional energetics

Dynamics of ions in the selectivity filter of the KcsA channel

The statistical and dynamical properties of ions in the selectivity filter of the KcsA ion channel are considered on the basis of molecular dynamics (MD) simulations of the KcsA protein embedded in a lipid membrane surrounded by an ionic solution. A new approach to the derivation of a Brownian dynamics (BD) model of ion permeation through the filter is discussed, based on unbiased MD simulations. It is shown that depending on additional assumptions, ion’s dynamics can be described either by under-damped Langevin equation with constant damping and white noise or by Langevin equation with a fractional memory kernel. A comparison of the potential of the mean force derived from unbiased MD simulations with the potential produced by the umbrella sampling method demonstrates significant differences in these potentials. The origin of these differences is an open question that requires further clarifications

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

Using Caenorhabditis elegans as a model organism for evaluating extracellular signal-regulated kinase docking domain inhibitors

We have recently identified several novel ATP-independent inhibitors that target the extracellular signal-regulated kinase-2 (ERK2) protein and inhibit substrate phosphorylation. To further characterize these compounds, we describe the use of C. elegans as a model organism. C. elegans is recognized as a versatile and cost effective model for use in drug development. These studies take advantage of the well characterized process of vulva development and egg laying, which requires MPK-1, the homolog to human ERK2. It is shown that treatment of C. elegans eggs or larvae prior to vulva formation with a previously identified lead compound (76) caused up to 50% reduction in the number of eggs produced from the adult worm. In contrast, compound 76 had no effect on egg laying in young adult or adult worms with fully formed vulva. The reduction in egg laying by the test compound was not due to effects on C. elegans life span, general toxicity, or non-specific stress. However, compound 76 did show selective inhibition of phosphorylation of LIN-1, a MPK-1 substrate essential for vulva precursor cell formation. Moreover, compound 76 inhibited cell fusion necessary for vulva maturation and reduced the multivulva phenotype in LET-60 (Ras) mutant worms that have constitutive activation of MPK-1. These findings support the use of C. elegans as a model organism to evaluate the selectivity and specificity of novel ERK targeted compounds