Summary The dynamic behavior of proteins is important for an understanding of their function and folding. We have performed molecular dynamics simulations of the native state and unfolding pathways of over 1000 proteins, representing the majority of folds in globular proteins. These data are stored and organized using an innovative database approach, which can be mined to obtain both general and specific information about the dynamics and folding/unfolding of proteins, relevant subsets thereof, and individual proteins. Here we describe the project in general terms and the type of information contained in the database. Then we provide examples of mining the database for information relevant to protein folding, structure building, the effect of single-nucleotide polymorphisms, and drug design. The native state simulation data and corresponding analyses for the 100 most populated metafolds, together with related resources, are publicly accessible through www.dynameomics.org.

The need for data-centric scientific tools is growing; domains like biology, chemistry, and physics are increasingly adopting computational approaches. As a result, scientists must now deal with the challenges of big data. To address these challenges, we built a visual analytics platform named DIVE: Data Intensive Visualization Engine. DIVE is a data-agnostic, ontologically-expressive software framework capable of streaming large datasets at interactive speeds. Here we present the technical details of the DIVE platform, multiple usage examples, and a case study from the Dynameomics molecular dynamics project. We specifically highlight our novel contributions to structured data model manipulation and high-throughput streaming of large, structured datasets.

The α-tocopherol transfer protein (α-TTP) is a liver protein that transfers α-tocopherol (vitamin E) to very-low-density lipoproteins (VLDLs). These VLDLs are then circulated throughout the body to maintain blood α-tocopherol levels. Mutations to the α-TTP gene are associated with ataxia with vitamin E deficiency (AVED), a disease characterized by peripheral nerve degeneration. In this study, molecular dynamics simulations of the E141K and R59W disease-associated mutants were performed. The mutants displayed disruptions in and around the ligand-binding pocket. Structural analysis and ligand docking to the mutant structures predicted a decreased affinity for α-tocopherol. To determine the detailed mechanism of the mutation-related changes, we developed a new tool called ContactWalker that analyzes contact differences between mutant and wild-type proteins and highlights pathways of altered contacts within the mutant proteins. Taken together, our findings are in agreement with experiment and suggest structural explanations for the reduced ability of the mutants to bind and carry α-tocopherol.

The p53 tumor suppressor protein performs a critical role in stimulating apoptosis and cell cycle arrest in response to oncogenic stress. The function of p53 can be compromised by mutation, leading to increased risk of cancer; approximately 50% of cancers are associated with mutations in the p53 gene, the majority of which are in the core DNA-binding domain. The Y220C mutation of p53, for example, destabilizes the core domain by 4 kcal/mol, leading to rapid denaturation and aggregation. The associated loss of tumor suppressor functionality is associated with approximately 75 000 new cancer cases every year. Destabilized p53 mutants can be 'rescued' and their function restored; binding of a small molecule into a pocket on the surface of mutant p53 can stabilize its wild-type structure and restore its function. Here, we describe an in silico algorithm for identifying potential rescue pockets, including the algorithm's integration with the Dynameomics molecular dynamics data warehouse and the DIVE visual analytics engine. We discuss the results of the application of the method to the Y220C p53 mutant, entailing finding a putative rescue pocket through MD simulations followed by an in silico search for stabilizing ligands that dock into the putative rescue pocket. The top three compounds from this search were tested experimentally and one of them bound in the pocket, as shown by nuclear magnetic resonance, and weakly stabilized the mutant.

Binaries and documentation are available at http://www.dynameomics.org/DIVE/DIVESetup.exe.