Biomedical subjects

Giovanni Ciccotti

Publications and source records attributed to Giovanni Ciccotti.

8 recordsLinked to original sources

PubMed2006

An adiabatic linearized path integral approach for quantum time-correlation functions II: a cumulant expansion method for improving convergence.

Linearized mixed quantum-classical simulations are a promising approach for calculating time-correlation functions. At the moment, however, they suffer from some numerical problems that may compromise their efficiency and reliability in applications to realistic condensed-phase systems. In this paper, we present a method that improves upon the convergence properties of the standard algorithm for linearized calculations by implementing a cumulant expansion of the relevant averages. The effectiveness of the new approach is tested by applying it to the challenging computation of the diffusion of an excess electron in a metal-molten salt solution.

Journal Article↗

PubMed2006

String method in collective variables: minimum free energy paths and isocommittor surfaces.

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.

Journal Article↗

PubMed2005

Blue moon sampling, vectorial reaction coordinates, and unbiased constrained dynamics.

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.

Journal Article↗

PubMed2005

The inelastic hard dimer gas: a nonspherical model for granular matter.

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.

Journal Article↗

PubMed2005

An adiabatic linearized path integral approach for quantum time correlation functions: electronic transport in metal-molten salt solutions.

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.

Journal Article↗

PubMed2004

Deterministic and stochastic algorithms for mechanical systems under constraints.

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.

Journal Article↗

PubMed2004

Atomic mean-square displacements in proteins by molecular dynamics: a case for analysis of variance.

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.

Algorithms↗

PubMedDate not supplied

The role of loop ZA and Pro371 in the function of yeast Gcn5p bromodomain revealed through molecular dynamics and experiment.

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.

Acetylation↗