Mechanisms for Benzene Dissociation through the Excited State of T4 Lysozyme L99A Mutant.

The atomic-level mechanisms that coordinate ligand release from protein pockets are only known for a handful of proteins. Here, we report results from accelerated molecular dynamics simulations for benzene dissociation from the buried cavity of the T4 lysozyme Leu99Ala mutant (L99A). In these simulations, benzene is released through a previously characterized, sparsely populated room-temperature excited state of the mutant, explaining the coincidence for experimentally measured benzene off rate and apo protein slow-timescale NMR relaxation rates between ground and excited states.

Alternative DNA secondary structure formation affects RNA polymerase II promoter-proximal pausing in human.

Background: Alternative DNA secondary structures can arise from single-stranded DNA when duplex DNA is unwound during DNA processes such as transcription, resulting in the regulation or perturbation of these processes. We identify sites of high propensity to form stable DNA secondary structure across the human genome using Mfold and ViennaRNA programs with parameters for analyzing DNA.

Results: The promoter-proximal regions of genes with paused transcription are significantly and energetically more favorable to form DNA secondary structure than non-paused genes or genes without RNA polymerase II (Pol II) binding.

Deep sequencing of 10,000 human genomes.

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.

Environmental and chemotherapeutic agents induce breakage at genes involved in leukemia-causing gene rearrangements in human hematopoietic stem/progenitor cells.

Hematopoietic stem and progenitor cells (HSPCs) give rise to all of the cells that make up the hematopoietic system in the human body, making their stability and resilience especially important. Damage to these cells can severely impact cell development and has the potential to cause diseases, such as leukemia. Leukemia-causing chromosomal rearrangements have largely been studied in the context of radiation exposure and are formed by a multi-step process, including an initial DNA breakage and fusion of the free DNA ends.

Theoretical Insights Reveal Novel Motions in Csk's SH3 Domain That Control Kinase Activation.

The Src family of tyrosine kinases (SFKs) regulate numerous aspects of cell growth and differentiation and are under the principal control of the C-terminal Src Kinase (Csk). Although Csk and SFKs share conserved kinase, SH2 and SH3 domains, they differ considerably in three-dimensional structure, regulatory mechanism, and the intrinsic kinase activities. Although the SH2 and SH3 domains are known to up- or down-regulate tyrosine kinase function, little is known about the global motions in the full-length kinase that govern these catalytic variations.

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.

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.

Distal loop flexibility of a regulatory domain modulates dynamics and activity of C-terminal SRC kinase (csk).

The Src family of tyrosine kinases (SFKs) regulate numerous aspects of cell growth and differentiation and are under the principal control of the C-terminal Src Kinase (Csk). Csk and SFKs share a modular design with the kinase domain downstream of the N-terminal SH2 and SH3 domains that regulate catalytic function and membrane localization. While the function of interfacial segments in these multidomain kinases are well-investigated, little is known about how surface sites and long-range, allosteric coupling control protein dynamics and catalytic function.

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.

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.

Exploring the Photophysical Properties of Molecular Systems Using Excited State Accelerated ab Initio Molecular Dynamics.

In the present work, we employ excited state accelerated ab initio molecular dynamics (A-AIMD) to efficiently study the excited state energy landscape and photophysical topology of a variety of molecular systems. In particular, we focus on two important challenges for the modeling of excited electronic states: (i) the identification and characterization of conical intersections and crossing seams, in order to predict different and often competing radiationless decay mechanisms, and (ii) the description of the solvent effect on the absorption and emission spectra of chemical species in solution. In particular, using as examples the Schiff bases formaldimine and salicylidenaniline, we show that A-AIMD can be readily employed to explore the conformational space around crossing seams in molecular systems with very different photochemistry.

The molecular dynamics of Trypanosoma brucei UDP-galactose 4'-epimerase: a drug target for African sleeping sickness.

During the past century, several epidemics of human African trypanosomiasis, a deadly disease caused by the protist Trypanosoma brucei, have afflicted sub-Saharan Africa. Over 10 000 new victims are reported each year, with hundreds of thousands more at risk. As current drug treatments are either highly toxic or ineffective, novel trypanocides are urgently needed.

Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations.

Molecular dynamics (MD) is one of the most common tools in computational chemistry. Recently, our group has employed accelerated molecular dynamics (aMD) to improve the conformational sampling over conventional molecular dynamics techniques. In the original aMD implementation, sampling is greatly improved by raising energy wells below a predefined energy level.

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.

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.

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.

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

Background: Gene 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.