Biomedical subjects

Mark E Tuckerman

Publications and source records attributed to Mark E Tuckerman.

24 records · Page 2Linked to original sources

PubMed2004

Reaction pathway of the [4 + 2] Diels-Alder adduct formation on Si(100)-2 x 1.

Despite a long history of experimental and theoretical investigation, the mechanism of the Diels-Alder (DA) reaction has been controversial since its discovery 80 years ago. From these investigations, two schools of thought have emerged, namely that the reaction can proceed via a concerted, symmetric or asymmetric mechanism or via a nonconcerted mechanism involving a zwitterion or diradical as an intermediate. Here, we employ finite temperature ab initio molecular dynamics simulations, employing forces computed "on the fly" from electronic structure calculations, to investigate the microscopic mechanism of DA adduct formation between 1,3-butadiene and the Si(100)-2x1 surface. Free energy profiles and nonequilibrium trajectories strongly suggest a nonconcerted mechanism that forms a zwitterionic intermediate state. This mechanism, which begins with a nucleophilic attack of the C=C double bond on the positive member of a charge-asymmetric buckled Si-Si dimer, was previously shown to be common to the formation of a wide range of adducts that can form on the surface.

Journal Article↗

PubMed2004

On-the-fly localization of electronic orbitals in Car-Parrinello molecular dynamics.

The ab initio molecular-dynamics formalism of Car and Parrinello is extended to preserve the locality of the orbitals. The supplementary term in the Lagrangian does not affect the nuclear dynamics, but ensures "on the fly" localization of the electronic orbitals within a periodic supercell in the Gamma-point approximation. The relationship between the resulting equations of motion and the formation of a gauge-invariant Lagrangian combined with a gauge-fixing procedure is briefly discussed. The equations of motion can be used to generate a very stable and easy to implement numerical integration algorithm. It is demonstrated that this algorithm can be used to compute the trajectory of the maximally localized orbitals, known as Wannier orbitals, in ab initio molecular dynamics with only a modest increase in the overall computer time. In the present paper, the new method is implemented within the generalized gradient approximation to Kohn-Sham density-functional theory employing plane wave basis sets and atomic pseudopotentials. In the course of the presentation, we briefly discuss how the present approach can be combined with localized basis sets to design fast linear scaling ab initio molecular-dynamics methods.

Journal Article↗

PubMed2004

Scalable fine-grained parallelization of plane-wave-based ab initio molecular dynamics for large supercomputers.

Many systems of great importance in material science, chemistry, solid-state physics, and biophysics require forces generated from an electronic structure calculation, as opposed to an empirically derived force law to describe their properties adequately. The use of such forces as input to Newton's equations of motion forms the basis of the ab initio molecular dynamics method, which is able to treat the dynamics of chemical bond-breaking and -forming events. However, a very large number of electronic structure calculations must be performed to compute an ab initio molecular dynamics trajectory, making the efficiency as well as the accuracy of the electronic structure representation critical issues. One efficient and accurate electronic structure method is the generalized gradient approximation to the Kohn-Sham density functional theory implemented using a plane-wave basis set and atomic pseudopotentials. The marriage of the gradient-corrected density functional approach with molecular dynamics, as pioneered by Car and Parrinello (R. Car and M. Parrinello, Phys Rev Lett 1985, 55, 2471), has been demonstrated to be capable of elucidating the atomic scale structure and dynamics underlying many complex systems at finite temperature. However, despite the relative efficiency of this approach, it has not been possible to obtain parallel scaling of the technique beyond several hundred processors on moderately sized systems using standard approaches. Consequently, the time scales that can be accessed and the degree of phase space sampling are severely limited. To take advantage of next generation computer platforms with thousands of processors such as IBM's BlueGene, a novel scalable parallelization strategy for Car-Parrinello molecular dynamics is developed using the concept of processor virtualization as embodied by the Charm++ parallel programming system. Charm++ allows the diverse elements of a Car-Parrinello molecular dynamics calculation to be interleaved with low latency such that unprecedented scaling is achieved. As a benchmark, a system of 32 water molecules, a common system size employed in the study of the aqueous solvation and chemistry of small molecules, is shown to scale on more than 1500 processors, which is impossible to achieve using standard approaches. This degree of parallel scaling is expected to open new opportunities for scientific inquiry.

Journal Article↗

PubMed2003

Molecular dynamics study of the connection between flap closing and binding of fullerene-based inhibitors of the HIV-1 protease.

The complementary spatial relationship between fullerene C(60) and the hydrophobic cavity region of the human immunodeficiency virus (HIV) protease, which houses the active site of the enzyme, has led to the suggestion that fullerene-based derivatives could have potential use as effective HIV protease inhibitors. The ability of such compounds to desolvate the cavity region leads to a strong hydrophobic interaction between the C(60) moiety and residues in the cavity region. In this study, the connection between the motion of the so-called flexible flaps of the cavity and favorable binding of a fullerene-based protease inhibitor is explored using multiple-time scale molecular dynamics simulations and free energy techniques. In addition, the effect of the interaction between the C(60) moiety and the residues in the cavity region on the water content of the cavity is also investigated. Conformational free energy profiles along a suitably chosen flap opening coordinate show a considerable barrier to flap opening in the presence of the inhibitor, while no such barrier exists for the protease alone. This result is interpreted in terms of a strong hydrophobic interaction between the C(60) moiety and the flexible flaps, which cause the latter to close tightly around the inhibitor, thereby expelling water from the cavity and leading to a favorable binding interaction. This interpretation is rationalized by direct analysis of the water content in the cavity in the presence and absence of the inhibitor.

Binding, Competitive↗

PubMed2002

The nature and transport mechanism of hydrated hydroxide ions in aqueous solution.

Compared to other ions, protons (H(+)) and hydroxide ions (OH(-)) exhibit anomalously high mobilities in aqueous solutions. On a qualitative level, this behaviour has long been explained by 'structural diffusion' the continuous interconversion between hydration complexes driven by fluctuations in the solvation shell of the hydrated ions. Detailed investigations have led to a clear understanding of the proton transport mechanism at the molecular level. In contrast, hydroxide ion mobility in basic solutions has received far less attention, even though bases and base catalysis play important roles in many organic and biochemical reactions and in the chemical industry. The reason for this may be attributed to the century-old notion that a hydrated OH(-) can be regarded as a water molecule missing a proton, and that the transport mechanism of such a 'proton hole' can be inferred from that of an excess proton by simply reversing hydrogen bond polarities. However, recent studies have identified OH(-) hydration complexes that bear little structural similarity to proton hydration complexes. Here we report the solution structures and transport mechanisms of hydrated hydroxide, which we obtained from first-principles computer simulations that explicitly treat quantum and thermal fluctuations of all nuclei. We find that the transport mechanism, which differs significantly from the proton hole picture, involves an interplay between the previously identified hydration complexes and is strongly influenced by nuclear quantum effects.

Journal Article↗

PubMed2002

Using novel variable transformations to enhance conformational sampling in molecular dynamics.

One of the computational "grand challenges" is to develop methodology capable of sampling conformational equilibria in systems with rough energy landscapes. Here, a significant step forward is made by combining molecular dynamics with a novel variable transformation designed to enhance sampling by reducing barriers without introducing bias and, thus, to preserve, perfectly, equilibrium properties.

Models, Chemical↗