Search PubMed⌕ Search

Biomedical subjects

David J Wales

Publications and source records attributed to David J Wales.

At least 19 recordsLinked to original sources

Potential energy and free energy landscapes.

Familiar concepts for small molecules may require reinterpretation for larger systems. For example, rearrangements between geometrical isomers are usually considered in terms of transitions between the corresponding local minima on the underlying potential energy surface, V. However, transitions between bulk phases such as solid and liquid, or between the denatured and native states of a protein, are normally addressed in terms of free energy minima. To reestablish a connection with the potential energy surface we must think in terms of representative samples of local minima of V, from which a free energy surface is projected by averaging over most of the coordinates. The present contribution outlines how this connection can be developed into a tool for quantitative calculations. In particular, stepping between the local minima of V provides powerful methods for locating the global potential energy minimum, and for calculating global thermodynamic properties. When the transition states that link local minima are also sampled we can exploit statistical rate theory to obtain insight into global dynamics and rare events. Visualizing the potential energy landscape helps to explain how the network of local minima and transition states determines properties such as heat capacity features, which signify transitions between free energy minima. The organization of the landscape also reveals how certain systems can reliably locate particular structures on the experimental time scale from among an exponentially large number of local minima. Such directed searches not only enable proteins to overcome Levinthal's paradox but may also underlie the formation of "magic numbers" in molecular beams, the self-assembly of macromolecular structures, and crystallization.

Journal Article↗

Graph transformation method for calculating waiting times in Markov chains.

We describe an exact approach for calculating transition probabilities and waiting times in finite-state discrete-time Markov processes. All the states and the rules for transitions between them must be known in advance. We can then calculate averages over a given ensemble of paths for both additive and multiplicative properties in a nonstochastic and noniterative fashion. In particular, we can calculate the mean first-passage time between arbitrary groups of stationary points for discrete path sampling databases, and hence extract phenomenological rate constants. We present a number of examples to demonstrate the efficiency and robustness of this approach.

Journal Article↗

Theory of C2Hx species on Pt{110}(1x2): structure, stability, and thermal chemistry.

The adsorption of C2Hx (x=0-5) hydrocarbon fragments on Pt{110}(1x2) has been investigated using calculations based on density functional theory. For all the species, the most stable adsorption site identified completes the tetravalency of each carbon atom and involves the maximum possible number of Pt atoms subject to that constraint. The most stable adsorption sites for C2Hx fragments of stoichiometry x=2-5 involve ridge atoms, while trough sites stabilize C2H and C2 species. The relative stability of the fragments involved is compared via a free energy picture.

Journal Article↗

Super-Arrhenius diffusion in an undercooled binary Lennard-Jones liquid results from a quantifiable correlation effect.

On short time scales an underlying Arrhenius temperature dependence of the diffusion constant can be extracted from the fragile, super-Arrhenius diffusion of a binary Lennard-Jones mixture. This Arrhenius diffusion is related to the true super-Arrhenius behavior by a factor that depends on the average angle between steps in successive time windows. The correction factor accounts for the fact that on average, successive displacements are negatively correlated, and this effect can therefore be linked directly with the higher apparent activation energy for diffusion at low temperature.

Journal Article↗

Equilibrium thermodynamics from basin-sampling.

We present a "basin-sampling" approach for calculation of the potential energy density of states for classical statistical models. It combines a Wang-Landau-type uniform sampling of local minima and a novel approach for approximating the relative contributions from local minima in terms of the volumes of basins of attraction. We have employed basin-sampling to study phase changes in atomic clusters modeled by the Lennard-Jones potential and for ionic clusters. The approach proves to be efficient for systems involving broken ergodicity and has allowed us to calculate converged heat capacity curves for systems that could previously only be treated using the harmonic superposition approximation. Benchmarks are also provided by comparison with parallel tempering and Wang-Landau simulations, where these proved feasible.

Journal Article↗

Novel structural motifs in clusters of dipolar spheres: knots, links, and coils.

We present the structures of putative global potential energy minima for clusters bound by the Stockmayer (Lennard-Jones plus point dipole) potential. A rich variety of structures is revealed as the cluster size and dipole strength are varied. Most remarkable are groups of closed-loop structures with the topology of knots and links. Despite the large number of possibilities, energetically optimal structures exhibit only a few such topologies.

Journal Article↗

Global optimization and folding pathways of selected alpha-helical proteins.

The results of basin-hopping global optimization simulations are presented for four small, alpha-helical proteins described by a coarse-grained potential. A step-taking scheme that incorporates the local conformational preferences extracted from a large number of high-resolution protein structures is compared with an unbiased scheme. In addition, the discrete path sampling method is used to investigate the folding of one of the proteins, namely, the villin headpiece subdomain. Folding times from kinetic Monte Carlo simulations and iterative calculations based on a Markovian first-step analysis for the resulting stationary-point database are in good mutual agreement, but differ significantly from the experimental values, probably because the native state is not the global free energy minimum for the potential employed.

Algorithms↗

Energy landscapes and properties of biomolecules.

Thermodynamic and dynamic properties of biomolecules can be calculated using a coarse-grained approach based upon sampling stationary points of the underlying potential energy surface. The superposition approximation provides an overall partition function as a sum of contributions from the local minima, and hence functions such as internal energy, entropy, free energy and the heat capacity. To obtain rates we must also sample transition states that link the local minima, and the discrete path sampling method provides a systematic means to achieve this goal. A coarse-grained picture is also helpful in locating the global minimum using the basin-hopping approach. Here we can exploit a fictitious dynamics between the basins of attraction of local minima, since the objective is to find the lowest minimum, rather than to reproduce the thermodynamics or dynamics.

Algorithms↗

Diagnosing broken ergodicity using an energy fluctuation metric.

The Mountain and Thirumalai energy fluctuation metric, Omega(t), has been used to study the effective ergodicity of 60- and 256-atom binary Lennard-Jones mixtures in order to determine the reliability of the calculated diffusion constants at different energies. A plot of Omega(t) against 1time allows the identification of two distinct regimes: ergodic supercooled liquids, where Omega(t) approaches zero, and nonergodic glassy states, where Omega(t) asymptotically approaches a nonzero value on the molecular dynamics time scale. This approach seems to be more appropriate than attempting to define a threshold value for Omega(t)/Omega(0). The behavior of systems between these two limits, which are nonergodic on the time scale considered but may be approaching ergodicity, was examined for a range of simulation times. The calculated diffusion constants change as effective ergodicity is approached, moving closer to the Vogel-Tammann-Fulcher fit defined by higher-energy systems that are already considered to be effectively ergodic. Using the form of the decay of the metric as a measure of ergodicity, we have been able to reproduce the trend in fragility obtained by Sastry for a 256-atom system [Nature (London) 409, 164 (2001)], correcting some of our earlier results [J. Chem. Phys. 120, 8314 (2004)].

Journal Article↗

Influence of vibrational energy flow on isomerization of flexible molecules: incorporating non-Rice-Ramsperger-Kassel-Marcus kinetics in the simulation of dipeptide isomerization.

The conformational isomerization of a dipeptide, N-acetyl-tryptophan methyl amide (NATMA), is studied computationally by including important dynamical corrections to Rice-Ramsperger-Kassel-Marcus (RRKM) theory for the transition rate between pairs of isomers. The dynamical corrections arise from incomplete or sluggish vibrational energy flow in the dipeptide, a property suggested by the mode-selective chemistry that has been observed by Dian et al. [J. Chem. Phys. 120, 133 (2004)]. We compute the extent and rate of vibrational energy flow in NATMA quantum mechanically using local random matrix theory, which we then use to correct the RRKM theory rates. The latter rates are then introduced into a master equation to study the population dynamics of the dipeptide. Incomplete or slow vibrational energy flow is found to enhance the conformational selectivity of NATMA over RRKM estimates.

Chemistry, Physical↗

An ab initio study of tunneling splittings in the water trimer.

Tunneling splittings in the water trimer have been determined by the semiclassical WKB method, based on pathways characterized at the computational level of second-order Møller-Plesset theory with basis sets of aug-cc-pVTZ quality. This calculation takes into account the single-flip and bifurcation tunneling rearrangements of the water trimer. The predicted splittings are 37.93 cm(-1) for the flip and 6.50x10(-3) cm(-1) for bifurcation, which agree quite well with the corresponding experimental values of 43.52 cm(-1) and 9.63x10(-3) cm(-1).

Journal Article↗

Finding pathways between distant local minima.

We report a new algorithm for constructing pathways between local minima that involve a large number of intervening transition states on the potential energy surface. A significant improvement in efficiency has been achieved by changing the strategy for choosing successive pairs of local minima that serve as endpoints for the next search. We employ Dijkstra's algorithm [E. W. Dijkstra, Numer. Math. 1, 269 (1959)] to identify the "shortest" path corresponding to missing connections within an evolving database of local minima and the transition states that connect them. The metric employed to determine the shortest missing connection is a function of the minimized Euclidean distance. We present applications to the formation of buckminsterfullerene and to the folding of various biomolecules: the B1 domain of protein G, tryptophan zippers, and the villin headpiece subdomain. The corresponding pathways contain up to 163 transition states and will be used in future discrete path sampling calculations.

Journal Article↗

Protonated water clusters described by an empirical valence bond potential.

The properties of low-lying stationary points on the potential energy surfaces of singly protonated water clusters (H(2)O)(n)H(+), are investigated using an empirical valence bond potential. Candidate global minima are reported for n=2-4, 8, and 20-22. For n=8, the variation in the energies and structures of low-lying minima with the number of valence bond states included in the model is studied. For n=4 and 8, disconnectivity graphs are also reported and are compared to results for the equivalent neutral water clusters as described by the rigid TIP3P potential. For the larger clusters, n=20-22, the structural properties of the low energy minima are compared with recently published spectroscopic data on these systems. The observed differences between the n=20 and n=21 systems are qualitatively reproduced by the model potential, but the similarities between the n=21 and n=22 systems are not.

Journal Article↗

How many dimensions are required to approximate the potential energy landscape of a model protein?

A scheme to approximate the multidimensional potential energy landscape in terms of a minimal number of degrees of freedom is proposed using a linear transformation of the original atomic Cartesian coordinates. For one particular off-lattice model protein the inherent frustration can only be reproduced satisfactorily when a relatively large number of coordinates are employed. However, when this frustration is removed in a Go-type model, the number of coordinates required is significantly lower, especially around the global potential energy minimum. To aid our interpretation of the results we consider modified disconnectivity graphs where a measure of the structural diversity and a metric relation between the stationary points are incorporated.

Computer Simulation↗

The energy landscape as a unifying theme in molecular science.

The potential energy surface (PES) underlies most calculations of structure, dynamics and thermodynamics in molecular science. In this contribution connections between the topology of the PES and observable properties are developed for a coarse-grained model of virus capsid self-assembly. The model predicts that a thermodynamically stable, kinetically accessible icosahedral shell exists for pentameric building blocks of the right shape: not too flat and not too spiky. The structure of the corresponding PES is probably common to other systems where directed searches avoid Levinthal's paradox, such as 'magic number' clusters, protein folding and crystallization.

Binding Sites↗

Archetypal energy landscapes: dynamical diagnosis.

Recent studies have identified several motifs for potential energy surfaces corresponding to distinct dynamic and thermodynamic properties. The corresponding disconnectivity graphs were identified as "palm tree," "willow tree," and "banyan tree" patterns. In the present contribution we present a quantitative analysis of the relation between the topography and dynamics for each of these motifs. For the palm tree and willow tree forms we find that the arrangement of the stationary points in the monotonic sequences with respect to the global minimum is the most important factor in establishing the kinetic properties. However, the results are somewhat different for motifs involving a rough surface with several deep basins (banyan tree motif), with large barriers relative to the energy differences between minima. Here it is the size of the barrier for escape from the region relative to the barriers at the bottom that is most important. The present results may be helpful in distinguishing between the dynamics of "structure seeking" and "glass forming" systems.

Journal Article↗

Effect of salt bridges on the energy landscape of a model protein.

The effect of introducing salt bridges (gatekeepers) into an off-lattice three-color, 46-bead model protein is investigated in terms of the effect on global optimization statistics. The global minima for all the gatekeepers that exhibited faster folding in previous molecular dynamics studies are located more rapidly than for the original potential, although the global minimum itself may change. Visualization of the underlying potential energy surface using disconnectivity graphs reveals that the gatekeepers exhibit structure intermediate between the original potential and a Go model. Competition between low-lying minima and the global minimum is reduced in the gatekeepers compared to the original potential, and interconversion barriers are generally smaller.

Journal Article↗

Analysis of cooperativity and localization for atomic rearrangements.

We propose measures of localization and cooperativity for the analysis of atomic rearrangements. We show that for both clusters and bulk material cooperative rearrangements usually have significantly lower barriers than uncooperative ones, irrespective of the degree of localization. We also find that previous methods used to sample stationary points are biased towards rearrangements of particular types. Linear interpolation between local minima in double-ended transition state searches tends to produce cooperative rearrangements, while random perturbations of all the coordinates, as sometimes used in single-ended searches, have the opposite effect.

Journal Article↗