Biomedical subjects

Tim N Heinz

Publications and source records attributed to Tim N Heinz.

4 recordsLinked to original sources

PubMed2005

Combining the lattice-sum and reaction-field approaches for evaluating long-range electrostatic interactions in molecular simulations.

A new scheme, the lattice-sum-emulated reaction-field (LSERF) method, is presented that combines the lattice-sum (LS) and reaction-field (RF) approaches for evaluating electrostatic interactions in molecular simulations. More precisely, the LSERF scheme emulates a RF calculation (based on an atomic cutoff) via the LS machinery. This is achieved by changing the form of the electrostatic interactions in a standard LS calculation (Coulombic) to the form corresponding to RF electrostatics (Coulombic plus quadratic reaction-field correction term, truncated at the cutoff distance). It is shown (both analytically and numerically) that in the limit of infinite reciprocal-space accuracy, (i) the LSERF scheme with a finite reaction-field cutoff and a given reaction-field permittivity is identical to the RF scheme with the same parameters (and an atomic cutoff), and (ii) the LSERF scheme is identical to the LS scheme in the limit of an infinite reaction-field cutoff, irrespective of the reaction-field permittivity. This new scheme offers two key advantages: (i) from a conceptual point of view, it shows that there is a continuity between the RF and LS schemes and unifies them into a common framework; (ii) from a practical point of view, it allows us to perform RF calculations with arbitrarily large reaction-field cutoff distances for the same computational costs as a corresponding LS calculation. The optimal choice for the cutoff will be the one that achieves the best compromise between artifacts arising from the dielectric heterogeneity of the system (short cutoff) and its artificial periodicity (long cutoff). The implementation of the LSERF method is extremely easy, requiring only very limited modifications of any standard LS code. For practical applications to biomolecular systems, the use of the LSERF scheme with large reaction-field cutoff distances is expected to represent a significant improvement over the current RF simulations involving comparatively much shorter cutoffs.

Journal Article↗

PubMed2005

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↗

PubMed2005

Measuring 1H-1H and 1H-13C RDCs in methyl groups: example of pulse sequences with numerically optimized coherence transfer schemes.

The optimization of coherence-transfer pulse-sequence elements (CTEs) is the most challenging step in the construction of heteronuclear correlation NMR experiments achieving sensitivity close to its theoretical maximum (in the absence of relaxation) in the shortest possible experimental time and featuring active suppression of undesired signals. As reported in the present article, this complex optimization problem in a space of high dimensionality turns out to be numerically tractable. Based on the application of molecular dynamics in the space of pulse-sequence variables, a general method is proposed for constructing optimized CTEs capable of transferring an arbitrary (generally non-Hermitian) spin operator encoding the chemical shift of heteronuclear spins to an arbitrary spin operator suitable for signal detection. The CTEs constructed in this way are evaluated against benchmarks provided by the theoretical unitary bound for coherence transfer and the minimal required transfer time (when available). This approach is used to design a set of NMR experiments enabling direct and selective observation of individual (1)H-transitions in (13)C-labeled methyl spin systems close to optimal sensitivity and using a minimal number of spectra. As an illustrative application of the method, optimized CTEs are used to quantitatively measure (1)H-(1)H and (1)H-(13)C residual dipolar couplings (RDCs) in a 17 kDa protein weakly aligned by means of Pf1 phages.

Bacterial Outer Membrane Proteins↗

PubMed2004

A fast pairlist-construction algorithm for molecular simulations under periodic boundary conditions.

A new grid-cell algorithm is presented that permits the fast construction of cutoff-based nonbonded pairlists in molecular simulations under periodic boundary conditions based on an arbitrary box shape. The key features of the method are (1) the use of a one-dimensional mask array (to determine which grid cells contain interacting atoms) that incorporates the effect of periodicity, and (2) the grouping of adjacent interacting cells of the mask array into stripes, which permits the handling of empty cells with a very low computational overhead. Testing of the algorithm on water systems of different sizes (containing about 2000 to 11,000 molecules) shows that the method (1) is about an order of magnitude more efficient compared to a standard (double-loop) algorithm, (2) achieves quasi-linear scaling in the number of atoms, (3) is weakly sensitive in terms of efficiency to the chosen number of grid cells, and (4) can be easily parallelized.

Journal Article↗