Biomedical subjects

Celeste Sagui

Publications and source records attributed to Celeste Sagui.

12 recordsLinked to original sources

PubMed2006

The free energy landscape of small peptides as obtained from metadynamics with umbrella sampling corrections.

There is considerable interest in developing methodologies for the accurate evaluation of free energies, especially in the context of biomolecular simulations. Here, we report on a reexamination of the recently developed metadynamics method, which is explicitly designed to probe "rare events" and areas of phase space that are typically difficult to access with a molecular dynamics simulation. Specifically, we show that the accuracy of the free energy landscape calculated with the metadynamics method may be considerably improved when combined with umbrella sampling techniques. As test cases, we have studied the folding free energy landscape of two prototypical peptides: Ace-(Gly)(2)-Pro-(Gly)(3)-Nme in vacuo and trialanine solvated by both implicit and explicit water. The method has been implemented in the classical biomolecular code AMBER and is to be distributed in the next scheduled release of the code.

Algorithms↗

PubMed2006

Molecular dynamics simulations of DNA with polarizable force fields: convergence of an ideal B-DNA structure to the crystallographic structure.

We have investigated to what extent molecular dynamics (MD) simulations can reproduce DNA sequence-specific features, given different electrostatic descriptions and different cell environments. For this purpose, we have carried out multiple unrestrained MD simulations of the DNA duplex d(CCAACGTTGG)2. With respect to the electrostatic descriptions, two different force fields are studied: a traditional description based on atomic point charges and a polarizable force field. With respect to the cell environment, the difference between crystal and solution environments is emphasized, as well as the structural importance of divalent ions. By imposing the correct experimental unit cell environment, an initial configuration with two ideal B-DNA duplexes in the unit cell is shown to converge to the crystallographic structure. This convergence is measured by the appearance of sequence-dependent features that very closely resemble the crystallographic ones as well as by the decay of the all-atom root-mean-squared coordinates deviations (RMSD) with respect to the crystallographic structure. Given the appropriate crystallographic constraints, this is the first example of multiple nanosecond molecular dynamics trajectory that shows an ideal B-DNA model converging to an experimental structure, with a significant decay of RMSD.

Base Sequence↗

PubMed2006

Deprotonation of solvated formic acid: Car-Parrinello and metadynamics simulations.

The deprotonation of solvated formic acid was investigated theoretically with ab initio simulations. With the Car-Parrinello method, deprotonation and reprotonation by means of a proton wire were observed. The microscopics of these reactions were analyzed, and reveal the key role played by nearby water molecules in catalyzing the reactions. A constrained molecular dynamics calculation was carried out to estimate the dissociation free energy. Deprotonation of formic acid was further investigated with the recently developed metadynamics method using the formic acid oxygen coordination numbers as the collective variables. The determined free-energy landscape gives barriers similar to that obtained with the constrained free-energy calculation.

Computational Biology↗

PubMed2005

Quantum simulations of the structure and binding of glycopeptide antibiotic aglycons to cell wall analogues.

The recent rise of vancomycin-resistant enterococci (VRE) and vancomycin-resistant Staphylococcus aureus (VRSA) has given new impetus to the study of the binding between glycopeptide antibiotics and bacterial cell wall termini. Here, we report on an extensive first principles investigation of the binding of vancomycin, avoparcin, teicoplanin, and ristocetin aglycons with dipetides, Ac-d-Ala-X, where X = d-Lac and d-Ser (characteristic of VREs) and X = d-Ala, Gly (characteristic of non-VREs), and a model "methylated d-Ala" CH(2)CH(CH(3))COO(-), in liquid as well as gas phase. The gas-phase ordering of the binding, from strongest to weakest, is Gly, d-Ala, d-Ser, CH(2)CH(CH(3))COO(-), and d-Lac. Calculations show that the order of the Gly and d-Ala binding is reversed in solution. The results are in good agreement with recent experimental findings.

Anti-Bacterial Agents↗

PubMed2005

Exploring intramolecular reactions in complex systems with metadynamics: the case of the malonate anions.

We have determined the optimized structures, relative energies and intramolecular reactions for two anionic forms of malonic acid, anion malonate(-1) (HO(2)CCH(2)CO(2)(-)) and malonate(-2) ((-)O(2)CCH(2)CO(2)(-)). For this purpose we employed accurate quantum chemistry calculations using second-order Möller-Plesset perturbation theory and Density Functional Theory with an aug'-cc-p-VTZ basis set to determine the structures and energies, and a novel metadynamics method based on Car-Parrinello molecular dynamics for the thermal reactions. For both malonates, we found new isomers (keto and enol structures) characterized by CO(2) rotations and intramolecular proton transfers. These proton transfers characterize the keto-enol tautomerism that takes place both in the monoanion and dianion. In all cases, the keto tautomer is the more stable configuration. The metadynamics method allows the system to explore the potential energy surface in a few picoseconds, crossing activation barriers of 20-50 kcal/mol.

Journal Article↗

PubMed2005

Self-assembled patterns and strain-induced instabilities for modulated systems.

The self-assembled domain patterns of modulated systems are characteristic of a wide variety of chemical and physical systems, and are the result of competing interactions. From a technological point of view, there is considerable interest in these domain patterns, as they form suitable templates for the fabrication of nanostructures. We have analyzed the domains and instabilities that form in modulated systems, and show that a large variety of patterns--based on long-lived metastable or glassy states--may be formed as a compromise between the required equilibrium modulation period and the strain present in the system. The strain results from topologically constrained trajectories in phase space, that effectively preclude the equilibrium configuration.

Journal Article↗

PubMed2005

Surface solvation for an ion in a water cluster.

We have used molecular dynamics simulations to study the structural, dynamical, and thermodynamical properties of ions in water clusters. Careful evaluations of the free energy, internal energy, and entropy are used to address controversial or unresolved issues, related to the underlying physical cause of surface solvation, and the basic assumptions that go with it. Our main conclusions are the following. (i) The main cause of surface solvation of a single ion in a water cluster is both water and ion polarization, coupled to the charge and size of the ion. Interestingly, the total energy of the ion increases near the cluster surface, while the total energy of water decreases. Also, our analysis clearly shows that the cause of surface solvation is not the size of the total water dipole (unless this is too small). (ii) The entropic contribution is the same order of magnitude as the energetic contribution, and therefore cannot be neglected for quantitative results. (iii) A pure energetic analysis can give a qualitative description of the ion position at room temperature. (iv) We have observed surface solvation of a large positive iodinelike ion in a polarizable water cluster, but not in a nonpolarizable water cluster.

Journal Article↗

PubMed2005

New and exotic self-organized patterns for modulated nanoscale systems.

The self-assembled domain patterns of modulated systems are the result of competing short-range attractive and long-range repulsive interactions found in diverse physical and chemical systems. From an application point of view, there is considerable interest in these domain patterns, as they form templates suitable for the fabrication of nanostructures. In this work we have generated a variety of new and exotic patterns, which represent either metastable or glassy states. These patterns arise as a compromise between the required equilibrium modulation period and the strain resulting from topologically constrained trajectories in phase space that effectively preclude the equilibrium configuration.

Computer Simulation↗

PubMed2004

Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer in crystal environment: comparison of atomic point-charge, extra-point, and polarizable force fields.

Molecular dynamics simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure in a crystal environment. Three different force fields were used: a traditional description based on atomic point charges, a polarizable force field, and an "extra-point" force field (with additional charges on extranuclear sites). It is found that all the force fields reproduce fairly well the sequence-dependent features of the experimental structure. The polarizable force field, however, provides the most accurate representation of the crystal structure and the sequence-dependent effects observed in the experiment. These results point out to the need of the inclusion of polarization for accurate descriptions of DNA.

Computer Simulation↗

PubMed2004

First principles investigation of vancomycin and teicoplanin binding to bacterial cell wall termini.

The recent rise of vancomycin-resistant enterococci (VRE) has given new impetus to the study of the binding between glycopeptide antibiotics and bacterial cell wall termini. Here, we report on an extensive first principles investigation of the binding of vancomycin and teicoplanin with d-Ala-d-Lac (characteristic of VREs) and d-Ala-d-Ala (characteristic of non-VREs). Binding of both antibiotics to d-Ala-d-Ala was found to be stronger by about 3-5 kcal/mol and due primarily to the oxygen-oxygen lone-pair repulsion characteristic of the antibiotic/d-Ala-d-Lac complex. These results are in good agreement with recent experimental findings.

Alanine↗

PubMed2004

Ab initio calculation of electrostatic multipoles with Wannier functions for large-scale biomolecular simulations.

It has long been known that accurate electrostatics is a key issue for improving current force fields for large-scale biomolecular simulations. Typically, this calls for an improved and more accurate description of the molecular electrostatic potential, which eliminates the artifacts associated with current point charge-based descriptions. In turn, this involves the partitioning of the extended molecular charge distribution, so that charges and multipole moments can be assigned to different atoms. As an alternate to current approaches, we have investigated a charge partitioning scheme that is based on the maximally localized Wannier functions. This has the advantage of partitioning the charge, and placing it around the molecule in a chemically meaningful manner. Moreover, higher order multipoles may all be calculated without any undue numerical difficulties. Tests on isolated molecules and water dimers, show that the molecular electrostatic potentials generated by such a Wannier-function based approach are in excellent agreement with the density functional-based calculations.

Algorithms↗

PubMed2004

Towards an accurate representation of electrostatics in classical force fields: efficient implementation of multipolar interactions in biomolecular simulations.

The accurate simulation of biologically active macromolecules faces serious limitations that originate in the treatment of electrostatics in the empirical force fields. The current use of "partial charges" is a significant source of errors, since these vary widely with different conformations. By contrast, the molecular electrostatic potential (MEP) obtained through the use of a distributed multipole moment description, has been shown to converge to the quantum MEP outside the van der Waals surface, when higher order multipoles are used. However, in spite of the considerable improvement to the representation of the electronic cloud, higher order multipoles are not part of current classical biomolecular force fields due to the excessive computational cost. In this paper we present an efficient formalism for the treatment of higher order multipoles in Cartesian tensor formalism. The Ewald "direct sum" is evaluated through a McMurchie-Davidson formalism [L. McMurchie and E. Davidson, J. Comput. Phys. 26, 218 (1978)]. The "reciprocal sum" has been implemented in three different ways: using an Ewald scheme, a particle mesh Ewald (PME) method, and a multigrid-based approach. We find that even though the use of the McMurchie-Davidson formalism considerably reduces the cost of the calculation with respect to the standard matrix implementation of multipole interactions, the calculation in direct space remains expensive. When most of the calculation is moved to reciprocal space via the PME method, the cost of a calculation where all multipolar interactions (up to hexadecapole-hexadecapole) are included is only about 8.5 times more expensive than a regular AMBER 7 [D. A. Pearlman et al., Comput. Phys. Commun. 91, 1 (1995)] implementation with only charge-charge interactions. The multigrid implementation is slower but shows very promising results for parallelization. It provides a natural way to interface with continuous, Gaussian-based electrostatics in the future. It is hoped that this new formalism will facilitate the systematic implementation of higher order multipoles in classical biomolecular force fields.

Algorithms↗