Biomedical subjects

David van der Spoel

Publications and source records attributed to David van der Spoel.

11 recordsLinked to original sources

PubMed2006

Temperature and structural changes of water clusters in vacuum due to evaporation.

This paper presents a study on evaporation of pure water clusters. Molecular dynamics simulations between 20 ns and 3 micros of clusters ranging from 125 to 4096 molecules in vacuum were performed. Three different models (SPC, TIP4P, and TIP5P) were used to simulate water, starting at temperatures of 250, 275, and 300 K. We monitored the temperature, the number of hydrogen bonds, the tetrahedral order, the evaporation, the radial distribution functions, and the diffusion coefficients. The three models behave very similarly as far as temperature and evaporation are concerned. Clusters starting at a higher temperature show a higher initial evaporation rate and therefore reach the point where evaporation stop (around 240 K) sooner. The radius of the clusters is decreased by 0.16-0.22 nm after 0.5 micros (larger clusters tend to decrease their radius slightly more), which corresponds to around one evaporated molecule per nm(2). The cluster temperature seems to converge towards 215 K independent of cluster size, when starting at 275 K. We observe only small structural changes, but the clusters modeled by TIP5P show a larger percentage of molecules with low diffusion coefficient as t-->infinity, than those using the two other water models. TIP4P seems to be more structured and more hydrogen bonds are formed than in the other models as the temperature falls. The cooling rates are in good agreement with experimental results, and evaporation rates agree well with a phenomenological expression based on experimental observations.

Computer Simulation↗

PubMed2006

Lifting a wet glass from a table: a microscopic picture.

Why is it so hard to lift a wet glass from a table? Is it easier when there is whiskey between the glass and the table? Macroscopically, the picture is quite simple: two surfaces have to be disrupted that are connected indirectly through hydrogen bonds and/or van der Waals forces. In the beginning, a surface has to be created leading to surface tension, and after that a liquid bridge has to be broken. Here we study the phenomenon at the microscopic level using molecular dynamics simulations. The effective force between two quartz plates is measured at different distances and with different alcohol/water mixtures between them. This allows us to compute the total work necessary to "lift the glass from the table". Different aspects of the process, such as clustering and liquid ordering are discussed. We compare the structure of the liquid/glass interface to that of a liquid/vapor interface, for which we present simulation results, like surface tension, as well. On the basis of the simulations, we are able to provide a detailed description of the energetics during the separation process as a function of alcohol concentration. It is shown that there is a net entropy loss upon separating two plates with water or a 10% MeOH solution between them, whereas for higher alcohol concentrations, there is net entropy gain. These findings increase our understanding of the properties of colloid suspensions which is important for process technology.

Journal Article↗

PubMed2006

Protein folding kinetics and thermodynamics from atomistic simulations.

Determining protein folding kinetics and thermodynamics from all-atom molecular dynamics (MD) simulations without using experimental data represents a formidable scientific challenge because simulations can easily get trapped in local minima on rough free energy landscapes. This necessitates the computation of multiple simulation trajectories, which can be independent from each other or coupled in some manner, as, for example, in the replica exchange MD method. Here we present results obtained with a new analysis tool that allows the deduction of faithful kinetics data from a heterogeneous ensemble of simulation trajectories. The method is demonstrated on the decapeptide Chignolin for which we predict folding and unfolding time constants of 1.0 +/- 0.3 and 2.6 +/- 0.4 micros, respectively. We also derive the energetics of folding, and calculate a realistic melting curve for Chignolin.

Computer Simulation↗

PubMed2006

Thermodynamics of hydrogen bonding in hydrophilic and hydrophobic media.

The thermodynamics of hydrogen bond breaking and formation was studied in solutions of alcohol (methanol, ethanol, 1-propanol) molecules. An extensive series of over 400 molecular dynamics simulations with an aggregate length of over 900 ns was analyzed using an analysis technique in which hydrogen bond (HB) breaking is interpreted as an Eyring process, for which the Gibbs energy of activation DeltaG can be determined from the HB lifetime. By performing simulations at different temperatures, we were able to determine the enthalpy of activation DeltaH and the entropy of activation TDeltaS for this process from the Van't Hoff relation. The equilibrium thermodynamics was determined separately, based on the number of donor hydrogens that are involved in hydrogen bonds. Results (DeltaH) are compared to experimental data from Raman spectroscopy and found to be in good agreement for pure water and methanol. The DeltaG as well as the DeltaG are smooth functions of the composition of the mixtures. The main result of the calculations is that DeltaG is essentially independent of the environment (around 5 kJ/mol), suggesting that buried hydrogen bonds (e.g., in proteins) do not contribute significantly to protein stability. Enthalpically HB formation is a downhill process in all substances; however, for the alcohols there is an entropic barrier of 6-7 kJ/mol, at 298.15 K, which cannot be detected in pure water.

1-Propanol↗

PubMed2006

Blind docking of drug-sized compounds to proteins with up to a thousand residues.

Blind docking was introduced for the detection of possible binding sites and modes of peptide ligands by scanning the entire surface of protein targets. In the present study, the method is tested on a group of drug-sized compounds and proteins with up to a thousand amino acid residues. Both proteins from complex structures and ligand-free proteins were used as targets. Robustness, limitations and future perspectives of the method are discussed. It is concluded that blind docking can be used for unbiased mapping of the binding patterns of drug candidates.

Binding Sites↗

PubMed2005

Reproducible polypeptide folding and structure prediction using molecular dynamics simulations.

The folding of a polypeptide from an extended state to a well-defined conformation is studied using microsecond classical molecular dynamics (MD) simulations and replica exchange molecular dynamics (REMD) simulations in explicit solvent and in vacuo. It is shown that the solvated peptide folds many times in the REMD simulations but only a few times in the conventional simulations. From the folding events in the classical simulations we estimate an approximate folding time of 1-2 micros. The REMD simulations allow enough sampling to deduce a detailed Gibbs free energy landscape in three dimensions. The global minimum of the energy landscape corresponds to the native state of the peptide as determined previously by nuclear magnetic resonance (NMR) experiments. Starting from an extended state it takes about 50 ns before the native structure appears in the REMD simulations, about an order of magnitude faster than conventional MD. The calculated melting curve is in good qualitative agreement with experiment. In vacuo, the peptide collapses rapidly to a conformation that is substantially different from the native state in solvent.

Computer Simulation↗

PubMed2004

Model for the dynamics of a water cluster in an x-ray free electron laser beam.

A microscopic sample placed into a focused x-ray free electron laser beam will explode due to strong ionization on a femtosecond time scale. The dynamics of this Coulomb explosion has been modeled by Neutze et al. [Nature (London) 406, 752 (2000)] for a protein, using computer simulations. The results suggest that by using ultrashort exposures, structural information may be collected before the sample is destroyed due to radiation damage. In this paper a method is presented to include the effect of screening by free electrons in the sample in a molecular dynamics simulation. The electrons are approximated by a classical gas, and the electron distribution is calculated iteratively from the Poisson-Boltzmann equation. Test simulations of water clusters reveal the details of the explosion dynamics, as well as the evolution of the free electron gas during the beam exposure. We find that inclusion of the electron gas in the model slows down the Coulomb explosion. The hydrogen atoms leave the sample faster than the oxygen atoms, leading to a double layer of positive ions. A considerable electron density is located between these two layers. The fact that the hydrogens are found to explode much faster than the oxygens means that the diffracting part of the sample stays intact somewhat longer than the sample as a whole.

Computer Simulation↗

PubMed2004

Molecular dynamics simulations of peptides from the central domain of smooth muscle caldesmon.

The central domain of smooth muscle caldesmon contains a highly charged region consisting of ten 13-residue repeats. Experimental evidence obtained from the intact protein and fragments thereof suggests that this entire region forms a single stretch of stable alpha-helix. We have carried out molecular dynamics simulations on peptides consisting of one, two and three repeats to examine the mechanism of alpha-helical stability of the central domain at the atomic level. All three peptides show high helical stability on the timescale of the MD simulations. Deviations from alpha-helical structure in all the simulations arise mainly from the formation of long stretches of pi-helix. Interconversion between alpha-helical and pi-helical conformations occurs through insertion of water molecules into alpha-helical hydrogen bonds and subsequent formation of reverse turns. The alpha-helical structure is stabilized by electrostatic interactions (salt bridges) between oppositely charged sidechains with i,i+4 spacings, while the pi-helix is stabilized by i,i+5 salt bridge interactions. Possible i,i+3 salt bridges are of minor importance. There is a strong preference for salt bridges with a Glu residue N-terminal to a basic sidechain as compared to the opposite orientation. In the double and triple repeat peptides, strong i,i+4 salt bridges exist between the last Glu residue of one repeat and the first Lys residue of the next. This demonstrates a relationship between the repetitive nature of the central domain sequence and its ability to form very long stretches of alpha-helical structure.

Calmodulin-Binding Proteins↗

PubMed2003

Deformation of helix C in the low temperature L-intermediate of bacteriorhodopsin.

X-ray and electron diffraction studies of specific reaction intermediates, or reaction intermediate analogues, have produced a consistent picture of the structural mechanism of light-driven proton pumping by bacteriorhodopsin. Of central importance within this picture is the structure of the L-intermediate, which follows the retinal all-trans to 13-cis photoisomerization step of the K-intermediate and sets the stage for the primary proton transfer event from the positively charged Schiff base to the negatively charged Asp-85. Here we report the structural changes in bacteriorhodopsin following red light illumination at 150 K. Single crystal microspectrophotometry showed that only the L-intermediate is populated in three-dimensional crystals under these conditions. The experimental difference Fourier electron density map and refined crystallographic structure were consistent with those previously presented (Royant, A., Edman, K., Ursby, T., Pebay-Peyroula, E., Landau, E. M., and Neutze, R. (2000) Nature 406, 645-648; Royant, A., Edman, K., Ursby, T., Pebay-Peyroula, E., Landau, E. M., and Neutze, R. (2001) Photochem. Photobiol. 74, 794-804). Based on the refined crystallographic structures, molecular dynamic simulations were used to examine the influence of the conformational change of the protein that is associated with the K-to-L transition on retinal dynamics. Implications regarding the structural mechanism for proton pumping by bacteriorhodopsin are discussed.

Bacteriorhodopsins↗

PubMed2003

Modeling of enzyme-substrate complexes for the metalloproteases MMP-3, ADAM-9 and ADAM-10.

The matrix metalloproteases (MMPs) and the ADAMs (A Disintegrin And Metalloprotease domain) are proteolytic enzyme families containing a catalytic zinc ion, that are implicated in a variety of normal and pathological processes involving tissue remodeling and cancer. Synthetic MMP inhibitors have been designed for applications in pathological situations. However, a greater understanding of substrate binding and the catalytic mechanism is required so that more effective and selective inhibitors may be developed for both experimental and clinical purposes. By modeling a natural substrate spanning P4-P4' in complex with the catalytic domains, we aim to compare substrate-specificities between Stromelysin-1 (MMP-3), ADAM-9 and ADAM-10, with the aid of molecular dynamics simulations. Our results show that the substrate retains a favourable antiparallel beta-sheet conformation on the P-side in addition to the well-known orientation of the P'-region of the scissile bond, and that the primary substrate selectivity is dominated by the sidechains in the S1' pocket and the S2/S3 region. ADAM-9 has a hydrophobic residue as the central determinant in the S1' pocket, while ADAM-10 has an amphiphilic residue, which suggests a different primary specificity. The S2/S3 pocket is largely hydrophobic in all three enzymes. Inspired by our molecular dynamics calculations and supported by a large body of literature, we propose a novel, hypothetical, catalytic mechanism where the Zn-ion polarizes the oxygens from the catalytic glutamate to form a nucleophile, leading to a tetrahedral oxyanion anhydride transition state.

ADAM Proteins↗

PubMed2002

Efficient docking of peptides to proteins without prior knowledge of the binding site.

Reliability in docking of ligand molecules to proteins or other targets is an important challenge for molecular modeling. Applications of the docking technique include not only prediction of the binding mode of novel drugs, but also other problems like the study of protein-protein interactions. Here we present a study on the reliability of the results obtained with the popular AutoDock program. We have performed systematical studies to test the ability of AutoDock to reproduce eight different protein/ligand complexes for which the structure was known, without prior knowledge of the binding site. More specifically, we look at factors influencing the accuracy of the final structure, such as the number of torsional degrees of freedom in the ligand. We conclude that the Autodock program package is able to select the correct complexes based on the energy without prior knowledge of the binding site. We named this application blind docking, as the docking algorithm is not able to "see" the binding site but can still find it. The success of blind docking represents an important finding in the era of structural genomics.

Algorithms↗