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.

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.

Methylene blue binding to DNA with alternating GC base sequence: continuum treatment of salt effects.

Methylene blue (MB), an efficient singlet oxygen generating photoactive dye, binds to DNA and allows photosensitized reactions to be used for sequence-specific cleavage of the DNA backbone. Intercalation and groove binding are possible binding modes of the dye, depending on base sequences and environmental conditions. In a recent modeling study of methylene blue binding to a double stranded DNA decamer with an alternating GC sequence, six structural models for intercalation structures and for minor and major groove binding have been obtained.

Conformational deformability of RNA: a harmonic mode analysis.

The harmonic mode analysis method was used to characterize the conformational deformability of regular Watson-Crick paired, mismatch- and bulge-containing RNA. Good agreement between atomic Debye-Waller factors derived from x-ray crystallography of a regular RNA oligonucleotide and calculated atomic fluctuations was obtained. Calculated helical coordinate fluctuations showed a small sequence dependence of up to approximately 30-50%.

Force field based conformational analysis of RNA structural motifs: GNRA tetraloops and their pyrimidine relatives.

The protocol of conformational analysis applied here to ribonucleotide oligomers combines conformational search in the space of torsion angles and energy minimization using the AMBER4.1 force field with a continuum treatment of electrostatic solute-solvent interactions. RNA fragments with 5'-GGGCGNNAGCCU-3' sequences commonly fold into hairpins with four-membered loops.

Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model.

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A -form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure.

Analysis of the stability of looped-out and stacked-in conformations of an adenine bulge in DNA using a continuum model for solvent and ions.

A combination of conformational search, energy minimization, and energetic evaluation using a continuum solvent treatment has been employed to study the stability of various conformations of the DNA fragment d(CGCAGAA)/d(TTCGCG) containing a single adenine bulge. The extra-helical (looped-out) bulge conformation derived from a published x-ray structure and intra-helical (stacked bulge base) model structures partially based on nuclear magnetic resonance (NMR) data were used as start structures for the conformational search. Solvent-dependent contributions to the stability of the conformations were calculated from the solvent exposed molecular surface area and by using the finite difference Poisson-Boltzmann approach.

Combining structural analysis of DNA with search routines for the detection of transcription regulatory elements.

Motivation: Analysis of unannotated genomic sequences for regulatory regions depends on a reliable recognition of individual cis-acting elements. Since some of them have a very low conserved sequence pattern, additional criteria are required.

Results: Using molecular modelling techniques, we have created a complete database for the conversion of base sequences into profiles of structural parameters of DNA.

NMR studies and restrained-molecular-dynamics calculations of a long A+T-rich stretch in DNA. Effects of phosphate charge and solvent approximations.

The nonamer duplex d(GCAAAAACG).d(CGTTTTTGC) was studied by 1H-NMR at 500 MHz. With the exception of the H5' and H5" sugar protons, all protons were assigned by two-dimensional NMR experiments [NOE spectroscopy (NOESY), double-quantum-filtered J-correlated spectroscopy (DQF-COSY) and total correlation spectroscopy (TOCSY)].

Conformation of d(GGGATCCC)2 in crystals and in solution studied by X-ray diffraction, Raman spectroscopy and molecular modelling.

In the crystal, d(GGGATCCC)2 forms an A-DNA double helix as known from a single crystal X-ray diffraction study. Accordingly, in the Raman spectra of crystals the A-family marker bands at 664, 705, 807 and 1101 cm-1 and the spectral characteristics in the region 1200 to 1500 cm-1 clearly demonstrate the A-form as the dominant conformation. Bands at 691, 850, and 1080 cm-1, however, indicate that a minor fraction of the octamer molecules in the crystal is in an unusual, still not unequivocally identified conformation possibly belonging to the B-family.

Conformational and helicoidal analysis of the molecular dynamics of proteins: "curves," dials and windows for a 50 psec dynamic trajectory of BPTI.

A new procedure for the graphic analysis of molecular dynamics (MD) simulations on proteins is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the conformational and helicoidal parameters for each structure entering the analysis via the method "Curves," developed for proteins by Sklenar, Etchebest, and Lavery (Proteins: Structure, Function Genet. 6:46-60, 1989) followed by a novel computer graphic display of the results.

Conformational and helicoidal analysis of 30 PS of molecular dynamics on the d(CGCGAATTCGCG) double helix: "curves", dials and windows.

A new procedure for the analysis of the structure and molecular dynamics of duplex DNA is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the values of the conformational and helicoidal parameters for each structure entering the analysis using the method "Curves" developed by Lavery and Sklenar, J. Biomol.

Defining the structure of irregular nucleic acids: conventions and principles.

The algorithm "Curves", that we have recently presented in this journal (J. Biolmol. Str.

Describing protein structure: a general algorithm yielding complete helicoidal parameters and a unique overall axis.

We present a general and mathematically rigorous algorithm which allows the helicoidal structure of a protein to be calculated starting from the atomic coordinates of its peptide backbone. This algorithm yields a unique curved axis which quantifies the folding of the backbone and a full set of helicoidal parameters describing the location of each peptide unit. The parameters obtained form a complete and independent set and can therefore be used for analyzing, comparing, or reconstructing protein backbone geometry.

The definition of generalized helicoidal parameters and of axis curvature for irregular nucleic acids.

An algorithm is presented which solves the problem of obtaining a rigorous helicoidal description of an irregular nucleic acid segment. Central to this approach is the definition of a function describing simultaneously the curvature of the nucleic acid segment in question and the corresponding stepwise variation of helicoidal parameters along the segment. Minimisation of this function leads to an optimal distribution of the conformational irregularity of the segment between these two components.

The flexibility of the nucleic acids: (II). The calculation of internal energy and applications to mononucleotide repeat DNA.

Results concerning the flexibility of mononucleotide repeat DNA are presented using a novel methodology, denoted "SIR", to describe continuous changes in the structure of the nucleic acid. This methodology, combined with internal energy calculations and analytical energy gradients allows us to determine optimal conformations of poly(dG).poly(dC) and poly (dA).

The flexibility of the nucleic acids: (I). "SIR", a novel approach to the variation of polymer geometry in constrained systems.

A novel and powerful methodology is developed which allows the alteration of molecular structures subjected to constraints and its application to polynucleotides with mononucleotide repeat symmetry, including the treatment of the flexible sugar rings is described. In contrast to procedures proposed by other authors, the constraints are formulated as differential equations which are linear with respect to the differentials of the geometrical variables. These equations can be solved easily by stepwise numerical integration involving sucessive infinitesimal rotations (SIR).