Biomedical subjects

Philippe H Hünenberger

Publications and source records attributed to Philippe H Hünenberger.

At least 19 recordsLinked to original sources

PubMed2006

The transition between the B and Z conformations of DNA investigated by targeted molecular dynamics simulations with explicit solvation.

The transition between the B and Z conformations of double-helical deoxyribonucleic acid (DNA) belongs to the most complex and elusive conformational changes occurring in biomolecules. Since the accidental discovery of the left-handed Z-DNA form in the late 1970s, research on this DNA morphology has been engaged in resolving questions relative to its stability, occurrence, and function in biological processes. While the occurrence of Z-DNA in vivo is now widely recognized and the major factors influencing its thermodynamical stability are largely understood, the intricate conformational changes that take place during the B-to-Z transition are still unknown at the atomic level. In this article, we report simulations of this transition for the 3'-(CGCGCG)-5' hexamer duplex using targeted molecular dynamics with the GROMOS96 force field in explicit water under different ionic-strength conditions. The results suggest that for this oligomer length and sequence, the transition mechanism involves: 1), a stretched intermediate conformation, which provides a simple solution to the important sterical constraints involved in this transition; 2), the transient disruption of Watson-Crick hydrogen-bond pairing, partly compensated energetically by an increase in the number of solute-solvent hydrogen bonds; and 3), an asynchronous flipping of the bases compatible with a zipperlike progression mechanism.

Computer Simulation↗

PubMed2006

Interaction of the sugars trehalose, maltose and glucose with a phospholipid bilayer: a comparative molecular dynamics study.

Molecular dynamics simulations are used to investigate the interaction of the sugars trehalose, maltose, and glucose with a phospholipid bilayer at atomic resolution. Simulations of the bilayer in the absence or in the presence of sugar (2 molal concentration for the disaccharides, 4 molal for the monosaccharide) are carried out at 325 and 475 K. At 325 K, the three sugars are found to interact directly with the lipid headgroups through hydrogen bonds, replacing water at about one-fifth to one-quarter of the hydrogen-bonding sites provided by the membrane. Because of its small size and of the reduced topological constraints imposed on the hydroxyl group locations and orientations, glucose interacts more tightly (at identical effective hydroxyl group concentration) with the lipid headgroups when compared to the disaccharides. At high temperature, the three sugars are able to prevent the thermal disruption of the bilayer. This protective effect is correlated with a significant increase in the number of sugar-headgroups hydrogen bonds. For the disaccharides, this change is predominantly due to an increase in the number of sugar molecules bridging three or more lipid molecules. For glucose, it is primarily due to an increase in the number of sugar molecules bound to one or bridging two lipid molecules.

Computer Simulation↗

PubMed2006

Configurational entropies of lipids in pure and mixed bilayers from atomic-level and coarse-grained molecular dynamics simulations.

Single-chain and single-fragment configurational entropies of lipid tails in hydrated lipid bilayers are evaluated from molecular dynamics simulations using the quasi-harmonic approximation. The entropy distribution along individual acyl tails is obtained and compared to that of corresponding hydrocarbon chains in the liquid phase. We consider pure dipalmitoylphosphatidylcholine and mixed dioleoylphosphatidylcholine/dioleoylphosphatidylethanolamine bilayers. The systems are modeled at different levels of spatial resolution: In an atomic-level (AL) model all (heavy) atoms are explicitly simulated; in a coarse-grained (CG) model particles (beads) representing groups of covalently bound atoms are used, which map approximately four non-hydrogen atoms to one interaction site. Single-chain and single-fragment entropies and correlations between the motions of (single) acyl chains are compared. A good correspondence is found between the flexibility of the AL and CG models. The loss in configurational entropy due to the reduction in the number of degrees of freedom upon coarse-graining of the model is estimated. The CG model shows about 4 times faster convergence of the chain entropies than the more detailed AL model. Corrections to the quasi-harmonic entropy estimates were found to be small for the CG model. For the AL model, the correction due to mode anharmonicities is small, but the correction due to pairwise (supralinear) mode correlations is sizable.

Computer Simulation↗

PubMed2006

Biomolecular modeling: Goals, problems, perspectives.

Computation based on molecular models is playing an increasingly important role in biology, biological chemistry, and biophysics. Since only a very limited number of properties of biomolecular systems is actually accessible to measurement by experimental means, computer simulation can complement experiment by providing not only averages, but also distributions and time series of any definable quantity, for example, conformational distributions or interactions between parts of systems. Present day biomolecular modeling is limited in its application by four main problems: 1) the force-field problem, 2) the search (sampling) problem, 3) the ensemble (sampling) problem, and 4) the experimental problem. These four problems are discussed and illustrated by practical examples. Perspectives are also outlined for pushing forward the limitations of biomolecular modeling.

Computer Simulation↗

PubMed2006

Computation of methodology-independent ionic solvation free energies from molecular simulations. II. The hydration free energy of the sodium cation.

The raw ionic solvation free energies computed from atomistic (explicit-solvent) simulations are extremely sensitive to the boundary conditions (finite or periodic system, system shape, and size) and treatment of electrostatic interactions (Coulombic, lattice sum, or cutoff based) used during these simulations. In the present article, it is shown that correction terms can be derived for the effect of (A) an incorrect solvent polarization around the ion due to the use of an approximate (not strictly Coulombic) electrostatic scheme; (B) the finite size or artificial periodicity of the simulated system; (C) an improper summation scheme to evaluate the potential at the ion site and the possible presence of a liquid-vacuum interface in the simulated system. Taking the hydration free energy of the sodium cation as a test case, it is shown that the raw solvation free energies obtained using seven different types of boundary conditions and electrostatic schemes commonly used in explicit-solvent simulations (for a total of 72 simulations differing in the corresponding simulation parameters) can be corrected so as to obtain a consistent value for this quantity.

Cations↗

PubMed2006

Comparison of atomic-level and coarse-grained models for liquid hydrocarbons from molecular dynamics configurational entropy estimates.

Molecular liquids can be modeled at different levels of spatial resolution. In atomic-level (AL) models, all (heavy) atoms can be explicitly simulated. In coarse-grained (CG) models, particles (beads) that represent groups of covalently bound atoms are used as elementary units. Ideally, a CG model should reproduce the thermodynamic and structural properties of the corresponding AL model after mapping to the lower-resolution scale. In the present work, two such models are investigated: (i) the classical GROMOS atomic-level model; (ii) a CG model recently proposed by Marrink et al., which maps approximately four non-hydrogen atoms to one bead [J. Phys. Chem. B 2004, 108, 750]. The study is restricted to n-alkanes whose aliphatic fragments are abundantly found in lipids of biological interest. Additionally, cis-9-octadecene is included, as a template chain of the lipid dioleoylphosphatidylcholine (DOPC). The two representations of molecules in the liquid phase are compared in terms of average molecular structures, extent of configurational space sampled, and single-molecule entropies. An approximate method is used to estimate the rotational contributions to the absolute configurational entropy. Good correspondence between the AL and CG representations is found. The loss in configurational entropy due to the reduction in degrees of freedom upon coarse-graining of the model is estimated.

Algorithms↗

PubMed2006

Conformational and dynamical properties of disaccharides in water: a molecular dynamics study.

Explicit-solvent molecular dynamics simulations (50 ns, 300 K) of the eight reducing glucose disaccharides (kojibiose, sophorose, nigerose, laminarabiose, maltose, cellobiose, isomaltose, and gentiobiose) have been carried out using the GROMOS 45A4 force field (including a recently reoptimized carbohydrate parameter set), to investigate and compare their conformational preferences, intramolecular hydrogen-bonding patterns, torsional dynamics, and configurational entropies. The calculated average values of the glycosidic torsional angles agree well with available experimental data, providing validation for the force field and simulation methodology employed in this study. These simulations show in particular that: 1) (1-->6)-linked disaccharides are characterized by an increased flexibility, the absence of any persistent intramolecular hydrogen bond and a significantly higher configurational entropy (compared to the other disaccharides); 2) cellobiose presents a highly persistent interresidue hydrogen bond and a significantly lower configurational entropy (compared to the other disaccharides); 3) persistent hydrogen bonds are observed for all disaccharides (except (1-->6)-linked) and typically involve a hydrogen donor in the reducing residue and an acceptor in the nonreducing one; 4) the probability distributions associated with the glycosidic dihedral angles and psi are essentially unimodal for all disaccharides, and full rotation around these angles occurs at most once or twice for (never for psi) on the 50-ns timescale; and 5) the timescales associated with torsional transitions (except around and psi) range from approximately 30 ps (rotation of hydroxyl groups) to the nanosecond range (rotation of the lactol and hydroxymethyl groups, and around the omega-glycosidic dihedral angle in (1-->6)-linked disaccharides).

Carbohydrate Sequence↗

PubMed2006

Computation of methodology-independent ionic solvation free energies from molecular simulations. I. The electrostatic potential in molecular liquids.

The computation of ionic solvation free energies from atomistic simulations is a surprisingly difficult problem that has found no satisfactory solution for more than 15 years. The reason is that the charging free energies evaluated from such simulations are affected by very large errors. One of these is related to the choice of a specific convention for summing up the contributions of solvent charges to the electrostatic potential in the ionic cavity, namely, on the basis of point charges within entire solvent molecules (M scheme) or on the basis of individual point charges (P scheme). The use of an inappropriate convention may lead to a charge-independent offset in the calculated potential, which depends on the details of the summation scheme, on the quadrupole-moment trace of the solvent molecule, and on the approximate form used to represent electrostatic interactions in the system. However, whether the M or P scheme (if any) represents the appropriate convention is still a matter of on-going debate. The goal of the present article is to settle this long-standing controversy by carefully analyzing (both analytically and numerically) the properties of the electrostatic potential in molecular liquids (and inside cavities within them). Restricting the discussion to real liquids of "spherical" solvent molecules (represented by a classical solvent model with a single van der Waals interaction site), it is concluded that (i) for Coulombic (or straight-cutoff truncated) electrostatic interactions, the M scheme is the appropriate way of calculating the electrostatic potential; (ii) for non-Coulombic interactions deriving from a continuously differentiable function, both M and P schemes generally deliver an incorrect result (for which an analytical correction must be applied); and (iii) finite-temperature effects, including intermolecular orientation correlations and a preferential orientational structure in the neighborhood of a liquid-vacuum interface, must be taken into account. Applications of these results to the computation methodology-independent ionic solvation free energies from molecular simulations will be the scope of a forthcoming article.

Journal Article↗

PubMed2006

Development of a lattice-sum method emulating nonperiodic boundary conditions for the treatment of electrostatic interactions in molecular simulations: a continuum-electrostatics study.

Artifacts induced by the application of periodic boundary conditions and lattice-sum methods in explicit-solvent simulations of (bio-)molecular systems are nowadays a major concern in the computer-simulation community. The present article reports a first step toward the design of a modified lattice-sum algorithm emulating nonperiodic boundary conditions, and therefore exempt of such periodicity-induced artifacts. This result is achieved here in the (more simple) context of continuum electrostatics. It is shown that an appropriate modification of the periodic Poisson equation and of its boundary conditions leads to a continuum-electrostatics scheme, which, although applied under periodic boundary conditions, exactly mimics the nonperiodic situation. The possible extension of this scheme to explicit-solvent simulations is outlined and its practical implementation will be described in more details in a forthcoming article.

Algorithms↗

PubMed2006

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↗

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

Molecular dynamics simulations of phospholipid bilayers: Influence of artificial periodicity, system size, and simulation time.

This article investigates the convergence of structural and dynamical properties with system size and with time in molecular dynamics simulations of solvated phospholipid bilayers performed at constant volume under periodic boundary conditions using lattice-sum electrostatics. The electron density profile across the bilayer, the carbon-deuterium order parameters, and the surface tension are shown to be converged for a bilayer containing 36 lipids per leaflet and simulated over a period of 3-4 ns. Reasonable estimates for these properties can already be obtained from a system containing 16 lipids per leaflet. The convergence limit of 36 lipids per leaflet and the investigation of the correlation between lipid headgroup dipoles suggest a correlation length of about 3-5 nm in the lateral directions for a hydrated DPPC bilayer in the liquid-crystalline phase. Although these (relatively small) system sizes and (relatively short) time scales appear sufficient to obtain converged collective structural properties at constant volume, two restrictions should be kept in mind: (i) the relaxation times associated with the motion of individual lipids may be much longer and (ii) simulated properties converge significantly faster under constant volume conditions as compared to constant pressure conditions. Therefore, an accurate assessment of the dynamical properties of the system or of the relaxation of the bilayer under constant pressure conditions may require longer simulation time scales.

Computer Simulation↗

PubMed2005

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↗

PubMed2005

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↗

PubMed2005

A new GROMOS force field for hexopyranose-based carbohydrates.

A new parameter set (referred to as 45A4) is developed for the explicit-solvent simulation of hexopyranose-based carbohydrates. This set is compatible with the most recent version of the GROMOS force field for proteins, nucleic acids, and lipids, and the SPC water model. The parametrization procedure relies on: (1) reassigning the atomic partial charges based on a fit to the quantum-mechanical electrostatic potential around a trisaccharide; (2) refining the torsional potential parameters associated with the rotations of the hydroxymethyl, hydroxyl, and anomeric alkoxy groups by fitting to corresponding quantum-mechanical profiles for hexopyranosides; (3) adapting the torsional potential parameters determining the ring conformation so as to stabilize the (experimentally predominant) (4)C(1) chair conformation. The other (van der Waals and nontorsional covalent) parameters and the rules for third and excluded neighbors are taken directly from the most recent version of the GROMOS force field (except for one additional exclusion). The new set is general enough to define parameters for any (unbranched) hexopyranose-based mono-, di-, oligo- or polysaccharide. In the present article, this force field is validated for a limited set of monosaccharides (alpha- and beta-D-glucose, alpha- and beta-D-galactose) and disaccharides (trehalose, maltose, and cellobiose) in solution, by comparing the results of simulations to available experimental data. More extensive validation will be the scope of a forthcoming article. (c) 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1400-1412, 2005.

Carbohydrate Conformation↗

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

Trehalose-protein interaction in aqueous solution.

A variety of sugars are known to enhance the stability of biomaterials. Trehalose, a nonreducing disaccharide composed of two alpha, alpha(1 --> 1)-linked D-glucopyranose units, appears to be one of the most effective protectants. Both in vivo and in vitro, trehalose protects biostructures such as proteins and membranes from damage due to dehydration, heat, or cold. However, despite the significant amount of experimental data on this disaccharide, no clear picture of the molecular mechanism responsible for its stabilizing properties has emerged yet. Three major hypotheses (water-trehalose hydrogen-bond replacement, coating by a trapped water layer, and mechanical inhibition of the conformational fluctuations) have been proposed to explain the stabilizing effect of trehalose on proteins. To investigate the nature of protein-trehalose-water interactions in solution at the molecular level, two molecular dynamics simulations of the protein lysozyme in solution at room temperature have been carried out, one in the presence (about 0.5 M) and one in the absence of trehalose. The results show that the trehalose molecules cluster and move toward the protein, but neither completely expel water from the protein surface nor form hydrogen bonds with the protein. Furthermore, the coating by trehalose does not significantly reduce the conformational fluctuations of the protein compared to the trehalose-free system. Based on these observations, a model is proposed for the interaction of trehalose molecules with a protein in moderately concentrated solutions, at room temperature and on the nanosecond timescale.

Computer Simulation↗