Atomistic Monte Carlo simulation of lipid membranes.

Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules.

Solvation oscillations and excited-state dynamics of 2-amino- and 2-hydroxy-7-nitrofluorene and its 2'-deoxyriboside.

Push-pull substituted fluorenes are considered for use as dynamic solvation probes in polynucleotides. Their fluorescence band is predicted (by simulations) to show weak spectral oscillations on the subpicosecond time scale depending on the nucleotide sequence. The oscillations reflect the local far-infrared spectrum of the environment around the probe molecule.

Using internal and collective variables in Monte Carlo simulations of nucleic acid structures: chain breakage/closure algorithm and associated Jacobians.

This article describes a method for solving the geometric closure problem for simplified models of nucleic acid structures by using the constant bond lengths approximation. The resulting chain breakage/closure equations, formulated in the space of variable torsion and bond angles, are easy to solve, and have only two solutions. The analytical simplicity is in contrast with the high complexity of the closure problem in the torsion angle space with at most 16 solutions, which has been dealt with by several authors and was solved analytically by Wu and Deem (J.

Molecular flexibility in ab initio drug docking to DNA: binding-site and binding-mode transitions in all-atom Monte Carlo simulations.

The dynamics of biological processes depend on the structure and flexibility of the interacting molecules. In particular, the conformational diversity of DNA allows for large deformations upon binding. Drug-DNA interactions are of high pharmaceutical interest since the mode of action of anticancer, antiviral, antibacterial and other drugs is directly associated with their binding to DNA.

Structural and energetic origins of sequence-specific DNA bending: Monte Carlo simulations of papillomavirus E2-DNA binding sites.

DNA bending is an important structural feature for indirect readout in protein-DNA recognition. The binding of papillomavirus E2 transcription factors to their DNA binding sites is associated with DNA bending, providing an attractive model system to study the origins of sequence-specific DNA bending. The consensus E2 target is of the general form ACCGN(4)CGGT with a variable four base pair region.

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps.

Molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide basepair steps are reported. The objective is to obtain the calculated dynamical structure for at least two copies of each case, use the results to examine issues with regard to convergence and dynamical stability of MD on DNA, and determine the significance of sequence context effects on all unique dinucleotide steps. This information is essential to understand sequence effects on DNA structure and has implications on diverse problems in the structural biology of DNA.

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. I. Research design and results on d(CpG) steps.

We describe herein a computationally intensive project aimed at carrying out molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide base sequences. This initiative was undertaken by an international collaborative effort involving nine research groups, the "Ascona B-DNA Consortium" (ABC). Calculations were carried out on the 136 cases imbedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs.

Methylene blue binding to DNA with alternating AT base sequence: minor groove binding is favored over intercalation.

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes.