Search

A computational technique is proposed which combines the string method with a sampling technique to determine minimum free energy paths. The technique only requires to compute the mean force and another conditional expectation locally along the string, and therefore can be applied even if the number of collective variables kept in the free energy calculation is large. This is in contrast with other free energy sampling techniques which aim at mapping the full free energy landscape and whose cost increases exponentially with the number of collective variables kept in the free energy. Provided that the number of collective variables is large enough, the new technique captures the mechanism of transition in that it allows to determine the committor function for the reaction and, in particular, the transition state region. The new technique is illustrated on the example of alanine dipeptide, in which we compute the minimum free energy path for the isomerization transition using either two or four dihedral angles as collective variables. It is shown that the mechanism of transition can be captured using the four dihedral angles, but it cannot be captured using only two of them.

In this work, we studied the poisoning of a nickel surface due to carbon. Performing ab initio simulations, within the framework of density functional theory, we computed the surface energy of the nickel (111) surface as a function of carbon coverage. On the basis of these results, we can assert that the stable state of the nickel/carbon surface is either a clean nickel surface or a fully carbon-covered nickel surface, which has a graphitic configuration. The relative stability of the two states depends on the temperature and partial pressure of the carbon gas. At fixed nominal carbon coverage, the most stable configurations are those forming carbon clusters. However, the nickel sites hosting these clusters change from hexagonal close packed/face centered cubic (hcp/fcc) sites to on-top sites when going toward larger clusters. This indicates that poisoning due to graphitic patches occurs on on-top sites.

Biological experiments were combined with molecular dynamics simulations to understand the importance of amino acidic residues present in the bromodomain of the yeast histone acetyltransferase Gcn5p. It was found that residue Pro371 plays an important role in the molecular recognition of the acetylated histone H4 tail by Gcn5p bromodomain. Crystallographic analysis of the complex showed that this residue does not directly interact with the histone substrate. It has been demonstrated that a double mutation Pro371Thr and Met372Ala in the Gcn5p bromodomain impairs chromatin remodeling activity. It is demonstrated here that, in this double mutant and in the fully deleted bromodomain strain, there is lower growth under amino acid deprivation conditions. By in vitro surface plasmon resonance (Biacore) experiments it is shown that the binding affinity of the double mutation to acetyl lysine 16 histone H4 peptide decreases. Molecular dynamics simulations were used to explain this loss in acetyl lysine-Gcn5p bromodomain affinity, in the double mutant. By comparing nanosecond molecular dynamics trajectories of the native as well as the single and doubly mutated bromodomain, it is concluded that the presence of Pro371 is important to the functionality of the Gcn5p bromodomain. In the simulation a point mutation involving this highly conserved residue induced an increase in the flexibility of the ZA loop, which in turn modulated the exposure of the binding pocket to the acetyl lysine. The combined double mutations (Pro371Thr-Met372Ala) not only markedly perturb the motion of the ZA loop but also destabilize the entire structure of the bromodomain., (Copyright 2005 John Wiley & Sons, Ltd.)

We give a new formula expressing the components of the mean force in terms of a conditional expectation which can be computed by Blue Moon sampling. This generalizes to the vectorial case a formula first derived by Ruiz-Montero et al. for a scalar reaction coordinate. We also discuss how to compute this conditional average by means of constrained stochastic dynamics which, unlike the usual constrained molecular dynamics, introduces no bias. Finally, we give a new perspective on bias removal by using constrained molecular dynamics.

We performed a room temperature molecular dynamics (MD) simulation on a system containing 1 carboxy-myoglobin (MbCO) molecule in a sucrose-water matrix of identical composition (89% [sucrose/(sucrose + water)] w/w) as for a previous trehalose-water-MbCO simulation (Cottone et al., Biophys J 2001;80:931-938). Results show that, as for trehalose, the amplitude of protein atomic mean-square fluctuations, on the nanosecond timescale, is reduced with respect to aqueous solutions also in sucrose. A detailed comparison as a function of residue number evidences mobility differences along the protein backbone, which can be related to a different efficacy in bioprotection. Different heme pocket structures are observed in the 2 systems. The joint distribution of the magnitude of the electric field at the CO oxygen atom and of the angle between the field and the CO unit vector shows a secondary maximum in sucrose, absent in trehalose. This can explain the CO stretching band profile (A substates distribution) differences evidenced by infrared spectroscopy in sucrose- and trehalose-coated MbCO (Giuffrida et al., J Phys Chem B 2004;108:15415-15421), and in particular the appearance of a further substate in sucrose. Analysis of hydrogen bonds at the protein-solvent interface shows that the fraction of water molecules shared between the protein and the sugar is lower in sucrose than in trehalose, in spite of a larger number of water molecules bound to the protein in the former system, thus indicating a lower protein-matrix coupling, as recently observed by Fourier transform infrared (FTIR) experiments (Giuffrida et al., J Phys Chem B 2004;108:15415-15421).

We study a two-dimensional gas of inelastic smooth hard dimers. Since the collisions between dimers are dissipative, being characterized by a coefficient of restitution alpha<1, and no external driving force is present, the energy of the system decreases in time and no stationary state is achieved. However, the resulting nonequilibrium state of the system displays several interesting properties in close analogy with systems of inelastic hard spheres, whose relaxational dynamics has been thoroughly explored. We generalize to inelastic systems a recently method introduced [G. Ciccotti and G. Kalibaeva, J. Stat. Phys. 115, 701 (2004)] to study the dynamics of rigid elastic bodies made up of different spheres held together by rigid bonds. Each dimer consists of two hard disks of diameter d, whose centers are separated by a fixed distance a. By describing the rigid bonds by means of holonomic constraints and deriving the appropriate collision rules between dimers, we reduce the dynamics to a set of equations which can be solved by means of event-driven simulation. After deriving the algorithm we study the decay of the total kinetic energy, and of the ratio between the rotational and the translational kinetic energy of inelastic dimers. We show numerically that the celebrated Haff's homogeneous cooling law t(-2), describing how the kinetic energy of an inelastic hard-sphere system with a constant coefficient of restitution decreases in time, holds even in the case of these nonspherical particles. We fully characterize this homogeneous decay process in terms of appropriate decay constants and confirm numerically the scaling behavior of the velocity distributions.

We generalize the linearized path integral approach to evaluate quantum time correlation functions for systems best described by a set of nuclear and electronic degrees of freedom, restricting ourselves to the adiabatic approximation. If the operators in the correlation function are nondiagonal in the electronic states, then this adiabatic linearized path integral approximation for the thermal averaged quantum dynamics presents interesting and distinctive features, which we derive and explore in this paper. The capability of these approximations to accurately reproduce the behavior of physical systems is demonstrated by calculating the diffusion constant for an excess electron in a metal-molten salt solution.

We discuss the general philosophy underlying the statistical behaviour, the dynamical evolution and the integration of the equations of motion for systems subject to constraints. We also show how all this is related to the treatment of general non-Hamiltonian systems. Then we introduce a family of algorithms derivable from approximations of the evolution operator obtained via the Trotter formula. Generalizing the treatment to time-dependent force fields we also show how one can adapt those algorithms to ordinary stochastic differential equations., (Copyright 2004 The Royal Society)

Information on protein internal motions is usually obtained through the analysis of atomic mean-square displacements, which are a measure of variability of the atomic positions distribution functions. We report a statistical approach to analyze molecular dynamics data on these displacements that is based on probability distribution functions. Using a technique inspired by the analysis of variance, we compute unbiased, reliable mean-square displacements of the atoms and analyze them statistically. We applied this procedure to characterize protein thermostability by comparing the results for a thermophilic enzyme and a mesophilic homolog. In agreement with previous experimental observations, our analysis suggests that the proteins surface regions can play a role in the different thermal behavior.

The Soret and the other transport coefficients, characterizing the heat and mass transport in binary mixtures, have been obtained by equilibrium and nonequilibrium molecular dynamics (EMD and NEMD, respectively). Two state points of the argon-krypton mixture are considered, for which experimental values of the Soret coefficient are available. To attempt a comparison between simulations and experiments the common enthalpy-diffusion-free expression for the heat flux has been chosen. The comparison of the simulations with the experiments shows a remarkable agreement, for all the several utilized EMD and NEMD techniques (dynamical and stationary). The techniques, used over 0.3 micros of total simulation time span, are slow convergent but have comparable performances.

We report on a molecular dynamics (MD) simulation of carboxy-myoglobin (MbCO) embedded in a water-trehalose system. The mean square fluctuations of protein atoms, calculated at different temperatures in the 100-300 K range, are compared with those from a previous MD simulation on an H2O-solvated MbCO and with experimental data from Mössbauer spectroscopy and incoherent elastic neutron scattering on trehalose-coated MbCO. The results show that, for almost all the atomic classes, the amplitude of the nonharmonic motions stemming from the interconversion among the protein's conformational substates is reduced with respect to the H2O-solvated system, and their onset is shifted toward higher temperature. Moreover, our simulation shows that, at 300 K, the heme performs confined diffusive motions as a whole, leaving the underlying harmonic vibrations unaltered.

We study by molecular dynamics simulations the hydration of beta-cyclodextrin. Our simulations show that within these barrel-shaped molecules hydrophobicity dominates, while at the top and bottom sides of the barrel interactions with water are mostly hydrophilic in nature. These results agree with crystallographic data at 120 K and, in particular, with the spontaneous hydration process of a cyclodextrin crystal in wet atmosphere. The predicted structure of the hydration shells is discussed and compared with previous molecular mechanics calculations which report an overall hydrophobic behavior. Moreover, the temperature dependence of the hydration process is discussed.

Molecular dynamics simulations, low temperature visible absorption spectroscopy, and resonance Raman spectroscopy have been performed on a mutant of the Scapharca inaequivalvis homodimeric hemoglobin, where residue threonine 72, at the subunit interface, has been substituted by isoleucine. Molecular dynamics simulation indicates that in the Thr-72-->Ile mutant several residues that have been shown to play a role in ligand binding fluctuate around orientations and distances similar to those observed in the x-ray structure of the CO derivative of the native hemoglobin, although the overall structure remains in the T state. Visible absorption spectroscopy data indicate that in the deoxy form the Soret band is less asymmetric in the mutant than in the native protein, suggesting a more planar heme structure; moreover, these data suggest a similar heme-solvent interaction in both the liganded and unliganded states of the mutant protein, at variance with that observed in the native protein. The "conformation sensitive" band III of the deoxy mutant protein is shifted to lower energy by >100 cm-1 with respect to the native one, about one-half of that observed in the low temperature photoproducts of both proteins, indicating a less polar or more hydrophobic heme environment. Resonance Raman spectroscopy data show a slight shift of the iron-proximal histidine stretching mode of the deoxy mutant toward lower frequency with respect to the native protein, which can be interpreted in terms of either a change in packing of the phenyl ring of Phe-97, as also observed from the simulation, or a loss of water in the heme pocket. In line with this latter interpretation, the number of water molecules that dynamically enters the intersubunit interface, as calculated by the molecular dynamics simulation, is lower in the mutant than in the native protein. The 10-ns photoproduct for the carbonmonoxy mutant derivative has a higher iron-proximal histidine stretching frequency than does the native protein. This suggests a subnanosecond relaxation that is slowed in the mutant, consistent with a stabilization of the R structure. Taken together, the molecular dynamics and the spectroscopic data indicate that the higher oxygen affinity displayed by the Thr-72-->Ile mutant is mainly due to a local perturbation in the dimer interface that propagates to the heme region, perturbing the polarity of the heme environment and propionate interactions. These changes are consistent with a destabilization of the T state and a stabilization of the R state in the mutant relative to the native protein.