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. The path observed for benzene egress is a multistep ligand migration from the buried cavity to ultimate release through an opening between the F/G-, H-, and I-helices and requires a number of cooperative multiresidue and secondary-structure rearrangements within the C-terminal domain of L99A. These rearrangements are identical to those observed along the ground state to excited state transitions characterized by molecular dynamic simulations run on the Anton supercomputer. Analyses of the molecular properties of the residues lining the egress path suggest that protein surface electrostatic potential may play a role in the release mechanism. Simulations of wild-type T4 lysozyme also reveal that benzene-egress-associated dynamics in the L99A mutant are potentially exaggerations of the substrate-processivity-related dynamics of the wild type.

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.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

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.

Drug screening strategy for human membrane proteins: from NMR protein backbone structure to in silica- and NMR-screened hits.

About 8000 genes encode membrane proteins in the human genome. The information about their druggability will be very useful to facilitate drug discovery and development. The main problem, however, consists of limited structural and functional information about these proteins because they are difficult to produce biochemically and to study. In this paper we describe the strategy that combines Cell-free protein expression, NMR spectroscopy, and molecular DYnamics simulation (CNDY) techniques. Results of a pilot CNDY experiment provide us with a guiding light towards expedited identification of the hit compounds against a new uncharacterized membrane protein as a potentially druggable target. These hits can then be further characterized and optimized to develop the initial lead compound quicker. We illustrate such "omics" approach for drug discovery with the CNDY strategy applied to two example proteins: hypoxia-induced genes HIGD1A and HIGD1B.

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. The SH2 domain of Csk is an essential component for the down-regulation of all SFKs. A unique feature of the SH2 domain of Csk is the tight turn in place of the canonical CD loop in a surface site far removed from kinase domain interactions. In this study, we used a combination of experimental and computational methods to probe the importance of this difference by constructing a Csk variant with a longer SH2 CD loop to mimic the flexibility found in homologous kinase SH2 domains. Our results indicate that while the fold and function of the isolated domain and the full-length kinase are not affected by loop elongation, native protein dynamics that are essential for efficient catalysis are perturbed. We also identify key motifs and routes through which the distal SH2 site might influence catalysis at the active site. This study underscores the sensitivity of intramolecular signaling and catalysis to native protein dynamics that arise from modest changes in allosteric regions while providing a potential strategy to alter intrinsic activity and signaling modulation.

Influence of N-myristylation and ligand binding on the flexibility of the catalytic subunit of protein kinase A.

The catalytic (C) subunit of protein kinase A is regulated in part by cotranslational N-myristylation and ligand binding. Using a combination of time-resolved fluorescence anisotropy and molecular dynamics (MD) simulations, we characterized the effect of N-myristylation and ligand binding on C-subunit dynamics. Five single-site cysteine-substitution mutants of the C-subunit were engineered with and without N-terminal myristylation and labeled with fluorescein maleimide, and time-resolved fluorescence anisotropy decays were measured to assess the flexibility of the labeled regions in the presence and absence of ligands. A parallel set of in silico experiments were performed to complement the experimental findings. These experiments showed that myristylation produces both local and global effects on C-subunit dynamics. The local effects include stabilization of the N-terminus and myristate pocket, and the global effects include small increases in mobility along the C-tail at residue C343. Additionally, ligand binding was associated with an increase in mobility of the myristate binding pocket for both the myristylated and nonmyristylated enzyme on the basis of both the experimental and MD results. Also, MD simulations suggest that the myristylated protein exhibits increased dynamics when bound to ligands compared to the nonmyristylated protein.

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.

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. RESULTS: Using 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. CONCLUSION: Our 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.

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. Using Pol II ChIP-seq, GRO-seq, NET-seq, and mNET-seq data, we arrive at a robust set of criteria for Pol II pausing, independent of annotation, and find that a highly stable secondary structure is likely to form about 10-50 nucleotides upstream of a Pol II pausing site. Structure probing data confirm the existence of DNA secondary structures enriched at the promoter-proximal regions of paused genes in human cells. Using an in vitro transcription assay, we demonstrate that Pol II pausing at HSPA1B, a human heat shock gene, is affected by manipulating DNA secondary structure upstream of the pausing site. CONCLUSIONS: Our results indicate alternative DNA secondary structure formation as a mechanism for how GC-rich sequences regulate RNA Pol II promoter-proximal pausing genome-wide.

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. However, the mechanism for DNA breakage in patients without previous radiation exposure is unclear. Here, we investigate the role of non-cytotoxic levels of environmental factors, benzene, and diethylnitrosamine (DEN), and chemotherapeutic agents, etoposide, and doxorubicin, in generating DNA breakage at the patient breakpoint hotspots of the MLL and CBFB genes in human HSPCs. These conditions represent exposure to chemicals encountered daily or residual doses from chemotherapeutic drugs. Exposure of HSPCs to non-cytotoxic levels of environmental chemicals or chemotherapeutic agents causes DNA breakage at preferential sites in the human genome, including the leukemia-related genes MLL and CBFB. Though benzene, etoposide, and doxorubicin have previously been linked to leukemia formation, this is the first study to demonstrate a role for DEN in the generation of DNA breakage at leukemia-specific sites. These chemical-induced DNA breakpoints coincide with sites of predicted topoisomerase II cleavage. The distribution of breakpoints by exposure to non-cytotoxic levels of chemicals showed a similar pattern to fusion breakpoints in leukemia patients. Our findings demonstrate that HSPCs exposed to non-cytotoxic levels of environmental chemicals and chemotherapeutic agents are prone to topoisomerase II-mediated DNA damage at the leukemia-associated genes MLL and CBFB. These data suggest a role for long-term environmental chemical or residual chemotherapeutic drug exposure in generation of DNA breakage at sites with a propensity to form leukemia-causing gene rearrangements.

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. We use a combination of accelerated Molecular Dynamics (aMD) simulations and experimental methods to provide a new view of functional motions in the Csk scaffold. These computational studies suggest that high frequency vibrations in the SH2 domain are coupled through the N-terminal lobe of the kinase domain to motions in the SH3 domain. The effects of these reflexive movements on the kinase domain can be viewed using both Deuterium Exchange Mass Spectrometry (DXMS) and steady-state kinetic methods. Removal of several contacts, including a crystallographically unobserved N-terminal segment, between the SH3 and kinase domains short-circuit these coupled motions leading to reduced catalytic efficiency and stability of N-lobe motifs within the kinase domain. The data expands the model of Csk's activation whereby separate domains productively interact with two diametrically opposed surfaces of the kinase domain. Such reversible transitions may organize the active structure of the tyrosine kinase domain of Csk.

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.

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. Recently, our group presented an alternative aMD implementation where simulations are accelerated by lowering energy barriers of the potential energy surface. When coupled with thermodynamic integration simulations, this implementation showed very promising results. However, when applied to large systems, such as proteins, the simulation tends to be biased to high energy regions of the potential landscape. The reason for this behavior lies in the boost equation used since the highest energy barriers are dramatically more affected than the lower ones. To address this issue, in this work, we present a new boost equation that prevents oversampling of unfavorable high energy conformational states. The new boost potential provides not only better recovery of statistics throughout the simulation but also enhanced sampling of statistically relevant regions in explicit solvent MD simulations.

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. The T. brucei galactose synthesis pathway is one potential therapeutic target. Although galactose is essential for T. brucei survival, the parasite lacks the transporters required to intake galactose from the environment. UDP-galactose 4'-epimerase (TbGalE) is responsible for the epimerization of UDP-glucose to UDP-galactose and is therefore of great interest to medicinal chemists. Using molecular dynamics simulations, we investigate the atomistic motions of TbGalE in both the apo and holo states. The sampled conformations and protein dynamics depend not only on the presence of a UDP-sugar ligand, but also on the chirality of the UDP-sugar C4 atom. This dependence provides important insights into TbGalE function and may help guide future computer-aided drug discovery efforts targeting this protein.

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.

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. Using acetone in water as an example, we demonstrate that the enhanced configurational space sampling may be used to accurately and efficiently describe both the prominent features and line-shapes of absorption and emission spectra.

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 × 10(5), while allowing for the reconstruction of the correct canonical distribution of low-energy states; (ii) We then show, for a water box of 64 H(2)O 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.

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.