PubMed Health⌕ Search

Biomedical subjects

Vincent Kräutler

Publications and source records attributed to Vincent Kräutler.

5 recordsLinked to original sources

A multiple time step algorithm compatible with a large number of distance classes and an arbitrary distance dependence of the time step size for the fast evaluation of nonbonded interactions in molecular simulations.

A new algorithm is introduced to perform the multiple time step integration of the equations of motion for a molecular system, based on the splitting of the nonbonded interactions into a series of distance classes. The interactions between particle pairs in successive classes are updated at a progressively decreasing frequency. Unlike previous multiple time-stepping schemes relying on distance classes, the present algorithm sorts interacting particle pairs by their next update times rather than by their update frequencies. For this reason, the proposed scheme is extremely flexible with respect to the number of classes that can be employed (up to hundred or more) and the distance dependence of the relative time step size (arbitrary integer function of the distance). It can also easily be adapted to classes defined based on a criterion other than the interparticle distance (e.g., interaction magnitude). Different variants of the algorithm are tested in terms of accuracy and efficiency for simulations of a pure water system (6167 molecules) under truncated-octahedral periodic boundary conditions, and compared to the twin-range method standardly used with GROMOS96 (short- and long-range cutoff distances of 0.8 and 1.4 nm, pair list and intermediate-range interactions updated every five steps). In particular, multiple time-stepping schemes with an accuracy comparable to that of the twin-range method can be designed, that permit to increase the effective (long-range) cutoff distance from 1.4 to 3.0 nm with a performance loss of only about a factor 2. This result is quite encouraging, considering the benefits of doubling the cutoff radius in the context of (bio-)molecular simulations.

Journal Article↗

Use of molecular dynamics in the design and structure determination of a photoinducible beta-hairpin.

The study presented here consists of three parts. In the first, the ability of a set of differently substituted diazobenzene-based linkers to act as photoswitchable beta-turn building blocks was assessed. A 12-residue peptide known to form beta-hairpins was taken as the basis for the modeling process. The central (beta-turn) residue pair was successively replaced by six symmetrically ((o,o), (m,m), or (p,p)) substituted (aminomethyl/carboxymethyl or aminoethyl/carboxyethyl) diazobenzene derivatives leading to a set of peptides with a photoswitchable backbone conformation. The folding behavior of each peptide was then investigated by performing molecular dynamics simulations in water (4 ns) and in methanol (10 ns) at room temperature. The simulations suggest that (o,o)- and (m,m)-substituted linkers with a single methylene spacer are significantly better suited to act as photoswitchable beta-turn building blocks than the other linkers examined in this study. The peptide containing the (m,m)-substituted linker was synthesized and characterized by NMR in its cis configuration. In the second part of this study, the structure of this peptide was refined using explicit-solvent simulations and NOE distance restraints, employing a variety of refinement protocols (instantaneous and time-averaged restraining as well as unrestrained simulations). We show that for this type of systems, even short simulations provide a significant improvement in our understanding of their structure if physically meaningful force fields are employed. In the third part, unrestrained explicit-solvent simulations starting from either the NMR model structure (75 ns) or a fully extended structure (25 ns) are shown to converge to a stable beta-hairpin. The resulting ensemble is in good agreement with experimental data, indicating successful structure prediction of the investigated hairpin by classical explicit-solvent molecular dynamics simulations.

Journal Article↗

A photoinducible beta-hairpin.

A photochromic azobenzene linker was incorporated as a turn element into an amino acid sequence known to fold into a beta-hairpin structure in aqueous solution. Oligomer formation when the linker was in its thermodynamically favored trans form prohibited structure determination. Light-induced conformational change of the linker to the cis form led to the formation of monomers which exhibited a well-defined beta-hairpin structure as determined by (1)H NMR. The rate of the light-induced cis-to-trans isomerization of the azobenzene-containing peptide was 30% slower compared to the unsubstituted chromophore. These results suggest that suitably substituted azobenzenes can be used as photoinducible turn elements to investigate and control the folding and stability of beta-sheets.

Azo Compounds↗

An improved nucleic acid parameter set for the GROMOS force field.

Over the past decades, the GROMOS force field for biomolecular simulation has primarily been developed for performing molecular dynamics (MD) simulations of polypeptides and, to a lesser extent, sugars. When applied to DNA, the 43A1 and 45A3 parameter sets of the years 1996 and 2001 produced rather flexible double-helical structures, in which the Watson-Crick hydrogen-bonding content was more limited than expected. To improve on the currently available parameter sets, the nucleotide backbone torsional-angle parameters and the charge distribution of the nucleotide bases are reconsidered based on quantum-chemical data. The new 45A4 parameter set resulting from this refinement appears to perform well in terms of reproducing solution NMR data and canonical hydrogen bonding. The deviation between simulated and experimental observables is now of the same order of magnitude as the uncertainty in the experimental values themselves.

Algorithms↗

The GROMOS software for biomolecular simulation: GROMOS05.

We present the latest version of the Groningen Molecular Simulation program package, GROMOS05. It has been developed for the dynamical modelling of (bio)molecules using the methods of molecular dynamics, stochastic dynamics, and energy minimization. An overview of GROMOS05 is given, highlighting features not present in the last major release, GROMOS96. The organization of the program package is outlined and the included analysis package GROMOS++ is described. Finally, some applications illustrating the various available functionalities are presented.

Algorithms↗