Levi C. T. Pierce

Routine Access to Millisecond Time Scale Events with Accelerated Molecular Dynamics

In this work, we critically assess the ability of the all-atom enhanced sampling method accelerated molecular dynamics (aMD) to investigate conformational changes in proteins that typically occur on the millisecond time scale. We combine aMD with the inherent power of graphics processor units (GPUs) and apply the implementation to the bovine pancreatic trypsin inhibitor (BPTI). A 500 ns aMD simulation is compared to a previous millisecond unbiased brute force MD simulation carried out on BPTI, showing that the same conformational space is sampled by both approaches. To our knowledge, this represents the first implementation of aMD on GPUs and also the longest aMD simulation of a biomolecule run to date. Our implementation is available to the community in the latest release of the Amber software suite (v12), providing routine access to millisecond events sampled from dynamics simulations using off the shelf hardware.

Over half of breakpoints in gene pairs involved in cancer-specific recurrent translocations are mapped to human chromosomal fragile sites.

BackgroundGene rearrangements such as chromosomal translocations have been shown to contribute to cancer development. Human chromosomal fragile sites are regions of the genome especially prone to breakage, and have been implicated in various chromosome abnormalities found in cancer. However, there has been no comprehensive and quantitative examination of the location of fragile sites in relation to all chromosomal aberrations.ResultsUsing up-to-date databases containing all cancer-specific recurrent translocations, we have examined 444 unique pairs of genes involved in these translocations to determine the correlation of translocation breakpoints and fragile sites in the gene pairs. We found that over half (52%) of translocation breakpoints in at least one gene of these gene pairs are mapped to fragile sites. Among these, we examined the DNA sequences within and flanking three randomly selected pairs of translocation-prone genes, and found that they exhibit characteristic features of fragile DNA, with frequent AT-rich flexibility islands and the potential of forming highly stable secondary structures.ConclusionOur study is the first to examine gene pairs involved in all recurrent chromosomal translocations observed in tumor cells, and to correlate the location of more than half of breakpoints to positions of known fragile sites. These results provide strong evidence to support a causative role for fragile sites in the generation of cancer-specific chromosomal rearrangements.

Adaptive Accelerated Molecular Dynamics (Ad-AMD) Revealing the Molecular Plasticity of P450cam

An extended accelerated molecular dynamics (AMD) methodology called adaptive AMD is presented. Adaptive AMD (Ad-AMD) is an efficient and robust conformational space sampling algorithm that is particularly-well suited to proteins with highly structured potential energy surfaces exhibiting complex, large-scale collective conformational transitions. Ad-AMD simulations of substrate-free P450cam reveal that this system exists in equilibrium between a fully and partially open conformational state. The mechanism for substrate binding depends on the size of the ligand. Larger ligands enter the P450cam binding pocket, and the resulting substrate-bound system is trapped in an open conformation via a population shift mechanism. Small ligands, which fully enter the binding pocket, cause an induced-fit mechanism, resulting in the formation of an energetically stable closed conformational state. These results are corroborated by recent experimental studies and potentially provide detailed insight into the functional dynamics and conformational behavior of the entire cytochrome-P450 superfamily.

On the Use of Accelerated Molecular Dynamics to Enhance Configurational Sampling in Ab Initio Simulations

We have implemented the accelerated molecular dynamics approach (Hamelberg, D.; Mongan, J.; McCammon, J. A. J. Chem. Phys. 2004, 120 (24), 11919) in the framework of ab initio MD (AIMD). Using three simple examples, we demonstrate that accelerated AIMD (A-AIMD) can be used to accelerate solvent relaxation in AIMD simulations and facilitate the detection of reaction coordinates: (i) We show, for one cyclohexane molecule in the gas phase, that the method can be used to accelerate the rate of the chair-to-chair interconversion by a factor of ∼1 × 105, while allowing for the reconstruction of the correct canonical distribution of low-energy states; (ii) We then show, for a water box of 64 H2O molecules, that A-AIMD can also be used in the condensed phase to accelerate the sampling of water conformations, without affecting the structural properties of the solvent; and (iii) The method is then used to compute the potential of mean force (PMF) for the dissociation of Na−Cl in water, accelerating the convergence by a factor of ∼3−4 compared to conventional AIMD simulations.(2) These results suggest that A-AIMD is a useful addition to existing methods for enhanced conformational and phase-space sampling in solution. While the method does not make the use of collective variables superfluous, it also does not require the user to define a set of collective variables that can capture all the low-energy minima on the potential energy surface. This property may prove very useful when dealing with highly complex multidimensional systems that require a quantum mechanical treatment.

Dynamics and calcium association to the N-terminal regulatory domain of human cardiac troponin C: a multiscale computational study.

Troponin C (TnC) is an important regulatory molecule in cardiomyocytes. Calcium binding to site II in TnC initiates a series of molecular events that result in muscle contraction. The most direct change upon Ca2+ binding is an opening motion of the molecule that exposes a hydrophobic patch on the surface allowing for Troponin I to bind. Molecular dynamics simulations were used to elucidate the dynamics of this crucial protein in three different states: apo, Ca2+-bound, and Ca2+-TnI-bound. Dynamics between the states are compared, and the Ca2+-bound system is investigated for opening motions. On the basis of the simulations, NMR chemical shifts and order parameters are calculated and compared with experimental observables. Agreement indicates that the simulations sample the relevant dynamics of the system. Brownian dynamics simulations are used to investigate the calcium association of TnC. We find that calcium binding gives rise to correlative motions involving the EF hand and collective motions conducive of formation of the TnI-binding interface. We furthermore indicate the essential role of electrostatic steering in facilitating diffusion-limited binding of Ca2+.

Role of DNA secondary structures in fragile site breakage along human chromosome 10

The formation of alternative DNA secondary structures can result in DNA breakage leading to cancer and other diseases. Chromosomal fragile sites, which are regions of the genome that exhibit chromosomal breakage under conditions of mild replication stress, are predicted to form stable DNA secondary structures. DNA breakage at fragile sites is associated with regions that are deleted, amplified or rearranged in cancer. Despite the correlation, unbiased examination of the ability to form secondary structures has not been evaluated in fragile sites. Here, using the Mfold program, we predict potential DNA secondary structure formation on the human chromosome 10 sequence, and utilize this analysis to compare fragile and non-fragile DNA. We found that aphidicolin (APH)-induced common fragile sites contain more sequence segments with potential high secondary structure-forming ability, and these segments clustered more densely than those in non-fragile DNA. Additionally, using a threshold of secondary structure-forming ability, we refined legitimate fragile sites within the cytogenetically defined boundaries, and identified potential fragile regions within non-fragile DNA. In vitro detection of alternative DNA structure formation and a DNA breakage cell assay were used to validate the computational predictions. Many of the regions identified by our analysis coincide with genes mutated in various diseases and regions of copy number alteration in cancer. This study supports the role of DNA secondary structures in common fragile site instability, provides a systematic method for their identification and suggests a mechanism by which DNA secondary structures can lead to human disease.

DNA secondary structure at chromosomal fragile sites in human disease.

DNA has the ability to form a variety of secondary structures that can interfere with normal cellular processes, and many of these structures have been associated with neurological diseases and cancer. Secondary structure-forming sequences are often found at chromosomal fragile sites, which are hotspots for sister chromatid exchange, chromosomal translocations, and deletions. Structures formed at fragile sites can lead to instability by disrupting normal cellular processes such as DNA replication and transcription. The instability caused by disruption of replication and transcription can lead to DNA breakage, resulting in gene rearrangements and deletions that cause disease. In this review, we discuss the role of DNA secondary structure at fragile sites in human disease.

Accelerating chemical reactions: Exploring reactive free-energy surfaces using accelerated ab initio molecular dynamics

A biased potential molecular dynamics simulation approach, accelerated molecular dynamics (AMD), has been implemented in the framework of ab initio molecular dynamics for the study of chemical reactions. Using two examples, the double proton transfer reaction in formic acid dimer and the hypothetical adiabatic ring opening and subsequent rearrangement reactions in methylenecyclopropane, it is demonstrated that ab initio AMD can be readily employed to efficiently explore the reactive potential energy surface, allowing the prediction of chemical reactions and the identification of metastable states. An adaptive variant of the AMD method is developed, which additionally affords an accurate representation of both the free-energy surface and the mechanism associated with the chemical reaction of interest and can also provide an estimate of the reaction rate.

Structural insight into the separate roles of inositol tetraphosphate and deacetylase-activating domain in activation of histone deacetylase 3

Histone deacetylases (HDACs) repress transcription by deacetylating acetyllysines on specific histone tails. HDAC3 is implicated in neurodegenerative diseases, certain leukemias, and even in disrupting HIV‐1 latency. A recent crystal structure of HDAC3 in complex with the deacetylase‐activating domain (DAD) of its corepressor complex revealed an inositol tetraphosphate (IP4) molecule at the protein–protein interface. IP4 was shown to play an important, yet mechanistically ambiguous, role in the activity of HDAC3. The goal of this article is to explore the conformational ensemble of HDAC3 in its inactive apo state and in the presence of each or both of DAD and IP4. Using triplicate, 100 ns molecular dynamic simulations, we study the apo, ternary, and intermediate DAD‐bound or IP4‐bound HDAC3 states. We find that a population‐shift effect is induced by the presence of each corepressor, and is most notable in the presence of both. Our results offer new insights into the change in dynamics necessary for the activation of HDAC3 and highlight the roles of IP4 and DAD in this process.

DNA Topoisomerases Participate in Fragility of the Oncogene RET

Fragile site breakage was previously shown to result in rearrangement of the RET oncogene, resembling the rearrangements found in thyroid cancer. Common fragile sites are specific regions of the genome with a high susceptibility to DNA breakage under conditions that partially inhibit DNA replication, and often coincide with genes deleted, amplified, or rearranged in cancer. While a substantial amount of work has been performed investigating DNA repair and cell cycle checkpoint proteins vital for maintaining stability at fragile sites, little is known about the initial events leading to DNA breakage at these sites. The purpose of this study was to investigate these initial events through the detection of aphidicolin (APH)-induced DNA breakage within the RET oncogene, in which 144 APH-induced DNA breakpoints were mapped on the nucleotide level in human thyroid cells within intron 11 of RET, the breakpoint cluster region found in patients. These breakpoints were located at or near DNA topoisomerase I and/or II predicted cleavage sites, as well as at DNA secondary structural features recognized and preferentially cleaved by DNA topoisomerases I and II. Co-treatment of thyroid cells with APH and the topoisomerase catalytic inhibitors, betulinic acid and merbarone, significantly decreased APH-induced fragile site breakage within RET intron 11 and within the common fragile site FRA3B. These data demonstrate that DNA topoisomerases I and II are involved in initiating APH-induced common fragile site breakage at RET, and may engage the recognition of DNA secondary structures formed during perturbed DNA replication.