A spectrin repeat unit was subject to extension using cyclic expansion nonequilibrium molecular dynamics. Periodic boundary conditions were used to examine the effects of the contiguous alpha-helical linker on the force response. The measured force-extension curve shows a linear increase in the forc

A spectrin repeat unit was subject to extension using cyclic expansion nonequilibrium molecular dynamics. Periodic boundary conditions were used to examine the effects of the contiguous alpha-helical linker on the force response. The measured force-extension curve shows a linear increase in the force response when the spectrin repeat unit is extended by approximately 0.4 nm. After that point, the force response peaks and subsequently declines. The peak in the force response marks the point where the spectrin repeat unit undergoes a change in its material properties from a strongly elastic material to a mostly viscous one, on the timescales of the simulations. The force peak is also correlated with rupture of the alpha-helical linker, and is likely the event responsible for the peaks in the sawtooth-pattern force-extension curves measured by atomic force microscopy experiments. Rupture of the linker involves simultaneously breaking approximately four hydrogen bonds that maintain the alpha-helical linker. After this initial rupture, the linker undergoes simple helix-to-coil transitions as the spectrin repeat unit continues to be extended. The implications of linker rupture in the interpretation of unfolding and atomic force microscopy experiments are also discussed.

Nonequilibrium molecular dynamics simulations were used to calculate the elastic properties of a spectrin repeat unit. A contiguous alpha-helical linker was constructed by employing periodic boundary conditions, allowing a novel scheme for evaluating the thermodynamic force as a function of extensio

Nonequilibrium molecular dynamics simulations were used to calculate the elastic properties of a spectrin repeat unit. A contiguous alpha-helical linker was constructed by employing periodic boundary conditions, allowing a novel scheme for evaluating the thermodynamic force as a function of extension. By measuring the force-extension response under small extensions, spectrin was observed to behave primarily as an elastic material with a spring constant of 1700 /- 100 pN/nm. The implications of this spring constant, in terms of the properties of the spectrin tetramer, are also discussed.

The unfolding pathways of multiple-repeat spectrin molecules were examined using steered molecular dynamics (SMD) simulations to forcibly unfold double- and triple-repeat spectrin molecules. Although SMD has previously been used to study other repeating-domain proteins, spectrin offers a unique chal

The unfolding pathways of multiple-repeat spectrin molecules were examined using steered molecular dynamics (SMD) simulations to forcibly unfold double- and triple-repeat spectrin molecules. Although SMD has previously been used to study other repeating-domain proteins, spectrin offers a unique challenge in that the linker connecting repeat units has a definite secondary structure, that of an alpha-helix. Therefore, the boundary conditions imposed on a double- or triple-repeat spectrin must be carefully considered if any relationship to the real system is to be deduced. This was accomplished by imposing additional forces on the system which ensure that the terminal alpha-helices behave as if there were no free noncontiguous helical ends. The results of the SMD simulations highlight the importance of the rupture of the alpha-helical linker on the subsequent unfolding events. Rupture of the linker propagates unfolding in the adjacent repeat units by destabilizing the tertiary structure, ultimately resulting in complete unfolding of the affected repeat unit. Two dominant classes of unfolding pathways are observed after the initial rupture of a linker which involve either rupture of another linker (possibly adjacent) or rupture of the basic tertiary structure of a repeat unit. The relationship between the force response observed on simulation timescales and those of experiment or physiological conditions is also discussed.

The fluctuation theorem describes the distribution of work done on small systems which have been pushed out of equilibrium in response to an external field. The theorem has recently been a subject of much interest for describing single-molecule experiments and simulations. In this communication, it

The fluctuation theorem describes the distribution of work done on small systems which have been pushed out of equilibrium in response to an external field. The theorem has recently been a subject of much interest for describing single-molecule experiments and simulations. In this communication, it is shown how the fluctuation theorem can be extended to describe fluctuations not only in the work done on a system, but also in a reaction coordinate. The extension explored in this work allows for a generalized derivation of Hummer and Szabo's expression (G. Hummer and A. Szabo, Proc. Natl. Acad. Sci. 98, 3658 (2001)) for reconstructing the potential of mean force from nonequilibrium trajectories. The derivation demonstrates how implementation of this expression can be more easily facilitated. Atomistic simulations of a biomolecular system are presented which support these results.

The synthetic atomic force microscopy (AFM) method is developed to simulate a periodically replicated atomistic system subject to force and length fluctuations characteristic of an AFM experiment. This new method is used to examine the forced-extension and subsequent rupture of the alpha-helical lin

The synthetic atomic force microscopy (AFM) method is developed to simulate a periodically replicated atomistic system subject to force and length fluctuations characteristic of an AFM experiment. This new method is used to examine the forced-extension and subsequent rupture of the alpha-helical linker connecting periodic images of a spectrin protein repeat unit. A two-dimensional potential of mean force (PMF) along the length and a reaction coordinate describing the state of the linker was calculated. This PMF reveals that the basic material properties of the spectrin repeat unit are sensitive to the state of linker, an important feature that cannot be accounted for in a one-dimensional PMF. Furthermore, nonequilibrium simulations were generated to examine the rupture event in the context of the fluctuation theorem. These atomistic simulations demonstrate that trajectories which are in apparent violation of the second law can overcome unfolding barriers at significantly reduced rupture forces.

The ff94 force field that is commonly associated with the Amber simulation package is one of the most widely used parameter sets for biomolecular simulation. After a decade of extensive use and testing, limitations in this force field, such as over-stabilization of alpha-helices, were reported by us

The ff94 force field that is commonly associated with the Amber simulation package is one of the most widely used parameter sets for biomolecular simulation. After a decade of extensive use and testing, limitations in this force field, such as over-stabilization of alpha-helices, were reported by us and other researchers. This led to a number of attempts to improve these parameters, resulting in a variety of "Amber" force fields and significant difficulty in determining which should be used for a particular application. We show that several of these continue to suffer from inadequate balance between different secondary structure elements. In addition, the approach used in most of these studies neglected to account for the existence in Amber of two sets of backbone phi/psi dihedral terms. This led to parameter sets that provide unreasonable conformational preferences for glycine. We report here an effort to improve the phi/psi dihedral terms in the ff99 energy function. Dihedral term parameters are based on fitting the energies of multiple conformations of glycine and alanine tetrapeptides from high level ab initio quantum mechanical calculations. The new parameters for backbone dihedrals replace those in the existing ff99 force field. This parameter set, which we denote ff99SB, achieves a better balance of secondary structure elements as judged by improved distribution of backbone dihedrals for glycine and alanine with respect to PDB survey data. It also accomplishes improved agreement with published experimental data for conformational preferences of short alanine peptides and better accord with experimental NMR relaxation data of test protein systems.

The thermodynamic properties and phase behavior of water in confined regions can vary significantly from that observed in the bulk. This is particularly true for systems in which the confinement is on the molecular-length scale. In this study, we use molecular dynamics simulations and a powerful sol

The thermodynamic properties and phase behavior of water in confined regions can vary significantly from that observed in the bulk. This is particularly true for systems in which the confinement is on the molecular-length scale. In this study, we use molecular dynamics simulations and a powerful solvent analysis technique based on inhomogenous solvation theory to investigate the properties of water molecules that solvate the confined regions of protein active sites. Our simulations and analysis indicate that the solvation of protein active sites that are characterized by hydrophobic enclosure and correlated hydrogen bonds induce atypical entropic and enthalpic penalties of hydration. These penalties apparently stabilize the protein-ligand complex with respect to the independently solvated ligand and protein, which leads to enhanced binding affinities. Our analysis elucidates several challenging cases, including the super affinity of the streptavidin-biotin system.

Understanding the underlying physics of the binding of small-molecule ligands to protein active sites is a key objective of computational chemistry and biology. It is widely believed that displacement of water molecules from the active site by the ligand is a principal (if not the dominant) source o

Understanding the underlying physics of the binding of small-molecule ligands to protein active sites is a key objective of computational chemistry and biology. It is widely believed that displacement of water molecules from the active site by the ligand is a principal (if not the dominant) source of binding free energy. Although continuum theories of hydration are routinely used to describe the contributions of the solvent to the binding affinity of the complex, it is still an unsettled question as to whether or not these continuum solvation theories describe the underlying molecular physics with sufficient accuracy to reliably rank the binding affinities of a set of ligands for a given protein. Here we develop a novel, computationally efficient descriptor of the contribution of the solvent to the binding free energy of a small molecule and its associated receptor that captures the effects of the ligand displacing the solvent from the protein active site with atomic detail. This descriptor quantitatively predicts (R(2) = 0.81) the binding free energy differences between congeneric ligand pairs for the test system factor Xa, elucidates physical properties of the active-site solvent that appear to be missing in most continuum theories of hydration, and identifies several features of the hydration of the factor Xa active site relevant to the structure-activity relationship of its inhibitors.

SU-8 is an epoxy-novolac resin and a well-established negative photoresist for microfabrication and microengineering. The photopolymerized resist is an extremely highly crosslinked polymer showing outstanding chemical and physical robustness with residual surface epoxy groups amenable for chemical f

SU-8 is an epoxy-novolac resin and a well-established negative photoresist for microfabrication and microengineering. The photopolymerized resist is an extremely highly crosslinked polymer showing outstanding chemical and physical robustness with residual surface epoxy groups amenable for chemical functionalization. In this paper we describe, for the first time, the preparation and surface modification of SU-8 particles shaped as microbars, the attachment of appropriate linkers, and the successful application of these particles to multistep solid-phase synthesis leading to oligonucleotides and peptides attached in an unambiguous manner to the support surface.

The current system of dissemination of scientific data and knowledge is far less efficient than it needs to be to facilitate improved collaborative science, especially considering current publication vehicles and infrastructure. There is a growing movement promoting more Open Science with the belief

The current system of dissemination of scientific data and knowledge is far less efficient than it needs to be to facilitate improved collaborative science, especially considering current publication vehicles and infrastructure. There is a growing movement promoting more Open Science with the belief that a more transparent scientific process can perform far more effectively. The logical extension of this concept is full transparency – exposing a researcher’s complete record of progress to the public in near real time. Not only will such a process enable ongoing data sharing it also provides an opportunity to develop collaborative communities of scientists and, at the conclusion of data acquisition, can enable communal extraction of conclusions when necessary. We have named this approach Open Notebook Science and have demonstrated its implementation and feasibility with the UsefulChem project, started in the summer of 2005, with the aim of synthesizing novel anti-malarial compounds. Our system currently uses free hosted services using general blog and wiki functions to facilitate replication across any scientific domains. These services are not chemically intelligent and are limited to text and graphic based data sharing only. For Open Notebook Chemistry the ability to intelligently manipulate, manage and search chemical structures and associated data is necessary and we have demonstrated proof of concept capabilities by integrating with the ChemSpider service, a free access online database managing chemical structures and focused on developing a structure centric community for chemists. This work will require the development of a chemically intelligent software platform to extend the capabilities of both the blog and the wiki environment for managing Open Notebook Science. The exposure of raw experimental procedures and data in a semantically rich format will enable the participation of both human and autonomous agents in the process of scientific discovery. This phenomenon of spontaneous group intelligence, referred to as “Crowdsourcing”, has proven valuable in several contexts. Already, productive collaborations have been forged within the UsefulChem project with groups from Indiana University, Nanyang Technological University, the National Cancer Institute and UC San Francisco.