The Dynameomics rotamer library: amino acid side chain conformations and dynamics from comprehensive molecular dynamics simulations in water

We have recently completed systematic molecular dynamics simulations of 807 different proteins representing 95% of the known autonomous protein folds in an effort we refer to as Dynameomics. Here we focus on the analysis of side chain conformations and dynamics to create a dynamic rotamer library. Overall this library is derived from 31,000 occurrences of each of 86,217 different residues, or 2.7 × 10(9) rotamers. This dynamic library has 74% overlap of rotamer distributions with rotamer libraries derived from static high-resolution crystal structures. Seventy-five percent of the residues had an assignable primary conformation, and 68% of the residues had at least one significant alternate conformation. The average correlation time for switching between rotamers ranged from 22 ps for Met to over 8 ns for Cys; this time decreased 20-fold on the surface of the protein and modestly for dihedral angles further from the main chain. Side chain S(2) axis order parameters were calculated and they correlated well with those derived from NMR relaxation experiments (R = 0.9). Relationships relating the S(2) axis order parameters to rotamer occupancy were derived. Overall the Dynameomics rotamer library offers a comprehensive depiction of side chain rotamer preferences and dynamics in solution, and more realistic distributions for dynamic proteins in solution at ambient temperature than libraries derived from crystal structures, in particular charged surface residues are better represented. Details of the rotamer library are presented here and the library itself can be downloaded at http://www.dynameomics.org.     

Molecular mechanisms underlying deoxy‐ADP.Pi activation of pre‐powerstroke myosin

Myosin activation is a viable approach to treat systolic heart failure. We previously demonstrated that striated muscle myosin is a promiscuous ATPase that can use most nucleoside triphosphates as energy substrates for contraction. When 2‐deoxy ATP (dATP) is used, it acts as a myosin activator, enhancing cross‐bridge binding and cycling. In vivo, we have demonstrated that elevated dATP levels increase basal cardiac function and rescues function of infarcted rodent and pig hearts. Here we investigate the molecular mechanism underlying this physiological effect. We show with molecular dynamics simulations that the binding of dADP.Pi (dATP hydrolysis products) to myosin alters the structure and dynamics of the nucleotide binding pocket, myosin cleft conformation, and actin binding sites, which collectively yield a myosin conformation that we predict favors weak, electrostatic binding to actin. In vitro motility assays at high ionic strength were conducted to test this prediction and we found that dATP increased motility. These results highlight alterations to myosin that enhance cross‐bridge formation and reveal a potential mechanism that may underlie dATP‐induced improvements in cardiac function.     

Tumorigenic p53 mutants undergo common structural disruptions including conversion to α‐sheet structure

The p53 protein is a commonly studied cancer target because of its role in tumor suppression. Unfortunately, it is susceptible to mutation‐associated loss of function; approximately 50% of cancers are associated with mutations to p53, the majority of which are located in the central DNA‐binding domain. Here, we report molecular dynamics simulations of wild‐type (WT) p53 and 20 different mutants, including a stabilized pseudo‐WT mutant. Our findings indicate that p53 mutants tend to exacerbate latent structural‐disruption tendencies, or vulnerabilities, already present in the WT protein, suggesting that it may be possible to develop cancer therapies by targeting a relatively small set of structural‐disruption motifs rather than a multitude of effects specific to each mutant. In addition, α‐sheet secondary structure formed in almost all of the proteins. α‐Sheet has been hypothesized and recently demonstrated to play a role in amyloidogenesis, and its presence in the reported p53 simulations coincides with the recent re‐consideration of cancer as an amyloid disease.     

Structural and Dynamic Properties of the Human Prion Protein

Prion diseases involve the conformational conversion of the cellular prion protein (PrP^C) to its misfolded pathogenic form (PrP^Sc). To better understand the structural mechanism of this conversion, we performed extensive all-atom, explicit-solvent molecular-dynamics simulations for three structures of the wild-type human PrP (huPrP) at different pH values and temperatures. Residue 129 is polymorphic, being either Met or Val. Two of the three structures have Met in position 129 and the other has Val. Lowering the pH or raising the temperature induced large conformational changes of the C-terminal globular domain and increased exposure of its hydrophobic core. In some simulations, HA and its preceding S1-HA loop underwent large displacements. The C-terminus of HB was unstable and sometimes partially unfolded. Two hydrophobic residues, Phe-198 and Met-134, frequently became exposed to solvent. These conformational changes became more dramatic at lower pH or higher temperature. Furthermore, Tyr-169 and the S2-HB loop, or the X-loop, were different in the starting structures but converged to common conformations in the simulations for the Met-129, but not the Val-129, protein. α-Strands and β-strands formed in the initially unstructured N-terminus. α-Strand propensity in the N-terminus was different between the Met-129 and Val129 proteins, but β-strand propensity was similar. This study reveals detailed structural and dynamic properties of huPrP, providing insight into the mechanism of the conversion of PrP^C to PrP^Sc.     

Sensitivity of the folding/unfolding transition state ensemble of chymotrypsin inhibitor 2 to changes in temperature and solvent

To better characterize the transition state for folding/unfolding and its sensitivity to environmental changes, we have run multiple molecular dynamics simulations of chymotrypsin inhibitor 2 (CI2) under varying solvent conditions and temperature. The transition state structures agree well with experiment, and are similar under all of the conditions investigated here. Increasing the temperature leads to some movement in the position of the transition state along several reaction coordinates, as measured by changes in properties of the transition state structures. These structural changes are in the direction of a more native-like transition state as denaturation conditions become more severe, as expected for a Hammond effect. These structural changes are not, however, reflected in the global structure as measured by the total number of contacts or the average S-values. These results suggest that the small changes in average Phi-values with temperature seen by experiment may be due to an increase in the sensitivity of the transition state to mutation rather than a change in the average structure of the transition state. A simple analysis of the rates of unfolding indicates that the free energy barrier to unfolding decreases with increasing temperature, but even in our very high temperature simulations there is a small free energy barrier.     

Dystrophin is required for the formation of stable muscle attachments in the zebrafish embryo

A class of recessive lethal zebrafish mutations has been identified in which normal skeletal muscle differentiation is followed by a tissue-specific degeneration that is reminiscent of the human muscular dystrophies. Here, we show that one of these mutations, sapje, disrupts the zebrafish orthologue of the X-linked human Duchenne muscular dystrophy (DMD) gene. Mutations in this locus cause Duchenne or Becker muscular dystrophies in human patients and are thought to result in a dystrophic pathology through disconnecting the cytoskeleton from the extracellular matrix in skeletal muscle by reducing the level of dystrophin protein at the sarcolemma. This is thought to allow tearing of this membrane, which in turn leads to cell death. Surprisingly, we have found that the progressive muscle degeneration phenotype of sapje mutant zebrafish embryos is caused by the failure of embryonic muscle end attachments. Although a role for dystrophin in maintaining vertebrate myotendinous junctions (MTJs) has been postulated previously and MTJ structural abnormalities have been identified in the Dystrophin-deficient mdx mouse model, in vivo evidence of pathology based on muscle attachment failure has thus far been lacking. This zebrafish mutation may therefore provide a model for a novel pathological mechanism of Duchenne muscular dystrophy and other muscle diseases.     

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long‐range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.     

The hydration of amides in helices; a comprehensive picture from molecular dynamics,IR, and NMR

We examined the hydration of amides of alpha(3)D, a simple, designed three-helix bundle protein. Molecular dynamics calculations show that the amide carbonyls on the surface of the protein tilt away from the helical axis to interact with solvent water, resulting in a lengthening of the hydrogen bonds on this face of the helix. Water molecules are bonded to these carbonyl groups with partial occupancy ( approximately 50%-70%), and their interaction geometries show a large variation in their hydrogen bond lengths and angles on the nsec time scale. This heterogeneity is reflected in the carbonyl stretching vibration (amide I' band) of a group of surface Ala residues. The surface-exposed amides are broad, and shift to lower frequency (reflecting strengthening of the hydrogen bonds) as the temperature is decreased. By contrast, the amide I' bands of the buried (13)C-labeled Leu residues are significantly sharper and their frequencies are consistent with the formation of strong hydrogen bonds, independent of temperature. The rates of hydrogen-deuterium exchange and the proton NMR chemical shifts of the helical amide groups also depend on environment. The partial occupancy of the hydration sites on the surface of helices suggests that the interaction is relatively weak, on the order of thermal energy at room temperature. One unexpected feature that emerged from the dynamics calculations was that a Thr side chain subtly disrupted the helical geometry 4-7 residues N-terminal in sequence, which was reflected in the proton chemical shifts and the rates of amide proton exchange for several amides that engage in a mixed 3(10)/alpha/pi-helical conformation.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

A consensus view of fold space: combining SCOP, CATH, and the Dali Domain Dictionary

We have determined consensus protein-fold classifications on the basis of three classification methods, SCOP, CATH, and Dali. These classifications make use of different methods of defining and categorizing protein folds that lead to different views of protein-fold space. Pairwise comparisons of domains on the basis of their fold classifications show that much of the disagreement between the classification systems is due to differing domain definitions rather than assigning the same domain to different folds. However, there are significant differences in the fold assignments between the three systems. These remaining differences can be explained primarily in terms of the breadth of the fold classifications. Many structures may be defined as having one fold in one system, whereas far fewer are defined as having the analogous fold in another system. By comparing these folds for a nonredundant set of proteins, the consensus method breaks up broad fold classifications and combines restrictive fold classifications into metafolds, creating, in effect, an averaged view of fold space. This averaged view requires that the structural similarities between proteins having the same metafold be recognized by multiple classification systems. Thus, the consensus map is useful for researchers looking for fold similarities that are relatively independent of the method used to compare proteins. The 30 most populated metafolds, representing the folds of about half of a nonredundant subset of the PDB, are presented here. The full list of metafolds is presented on the Web.