Multiscale modeling of nucleosome dynamics

Energy landscape analyses of disordered histone tails reveal special organization of their conformational dynamics

Histone tails are highly flexible N- or C-terminal protrusions of histone proteins which facilitate the compaction of DNA into dense superstructures known as chromatin. On a molecular scale histone tails are polyelectrolytes with high degree of conformational disorder which allows them to function as biomolecular "switches", regulating various genetic processes. Unfortunately, their intrinsically disordered nature creates obstacles for comprehensive experimental investigation of both the structural and dynamical aspects of histone tails, because of which their conformational behaviors are still not well understood. In this work we have carried out ∼3 microsecond long all atom replica exchange molecular dynamics (REMD) simulations for each of four histone tails, H4, H3, H2B, and H2A, and probed their intrinsic conformational preferences. Our subsequent free energy landscape analysis demonstrated that most tails are not fully disordered, but show distinct conformational organization, containing specific flickering secondary structural elements. In particular, H4 forms β-hairpins, H3 and H2B adopt α-helical elements, while H2A is fully disordered. We rationalized observed patterns of conformational dynamics of various histone tails using ideas from physics of polyelectrolytes and disordered systems. We also discovered an intriguing re-entrant contraction-expansion of the tails upon heating, which is caused by subtle interplay between ionic screening and chain entropy.

Hierarchies in eukaryotic genome organization: Insights from polymer theory and simulations

Eukaryotic genomes possess an elaborate and dynamic higher-order structure within the limiting confines of the cell nucleus. Knowledge of the physical principles and the molecular machinery that govern the 3D organization of this structure and its regulation are key to understanding the relationship between genome structure and function. Elegant microscopy and chromosome conformation capture techniques supported by analysis based on polymer models are important steps in this direction. Here, we review results from these efforts and provide some additional insights that elucidate the relationship between structure and function at different hierarchical levels of genome organization.

The crystal structure of a self-activating G protein alpha subunit reveals its distinct mechanism of signal initiation

In animals, heterotrimeric guanine nucleotide-binding protein (G protein) signaling is initiated by G protein-coupled receptors (GPCRs), which activate G protein α subunits; however, the plant Arabidopsis thaliana lacks canonical GPCRs, and its G protein α subunit (AtGPA1) is self-activating. To investigate how AtGPA1 becomes activated, we determined its crystal structure. AtGPA1 is structurally similar to animal G protein α subunits, but our crystallographic and biophysical studies revealed that it had distinct properties. Notably, the helical domain of AtGPA1 displayed pronounced intrinsic disorder and a tendency to disengage from the Ras domain of the protein. Domain substitution experiments showed that the helical domain of AtGPA1 was necessary for self-activation and sufficient to confer self-activation to an animal G protein α subunit. These findings reveal the structural basis for a mechanism for G protein activation in Arabidopsis that is distinct from the well-established mechanism found in animals.

Force field parameters for the simulation of modified histone tails

We describe the development of force field parameters for methylated lysines and arginines, and acetylated lysine for the CHARMM all-atom force field. We also describe a CHARMM united-atom force field for modified sidechains suitable for use with fragment-based docking methods. The development of these parameters is based on results of ab initio quantum mechanics calculations of model compounds with subsequent refinement and validation by molecular mechanics and molecular dynamics simulations. The united-atom parameters are tested by fragment docking to target proteins using the MCSS procedure. The all-atom force field is validated by molecular dynamics simulations of multiple experimental structures. In both sets of calculations, the computational predictions using the force field were compared to the corresponding experimental structures. We show that the parameters yield an accurate reproduction of experimental structures. Together with the existing CHARMM force field, these parameters will enable the general modeling of post-translational modifications of histone tails.

Thermodynamic stability of histone H3 is a necessary but not sufficient driving force for its evolutionary conservation

Determining the forces that conserve amino acid positions in proteins across species is a fundamental pursuit of molecular evolution. Evolutionary conservation is driven by either a protein's function or its thermodynamic stability. Highly conserved histone proteins offer a platform to evaluate these driving forces. While the conservation of histone H3 and H4 “tail” domains and surface residues are driven by functional importance, the driving force behind the conservation of buried histone residues has not been examined. Using a computational approach, we determined the thermodynamically preferred amino acids at each buried position in H3 and H4. In agreement with what is normally observed in proteins, we find a significant correlation between thermodynamic stability and evolutionary conservation in the buried residues in H4. In striking contrast, we find that thermodynamic stability of buried H3 residues does not correlate with evolutionary conservation. Given that these H3 residues are not post-translationally modified and only regulate H3-H3 and H3-H4 stabilizing interactions, our data imply an unknown function responsible for driving conservation of these buried H3 residues.

Most proteins fold to a well-defined, three-dimensional structure, which can be delineated into the protein surface and its buried core. When comparing amino acid sequences of the same protein from different organisms, we would expect to find certain residue positions conserved due to the importance of that position in either maintaining the protein's function or its three-dimensional structure. In this study, we looked at residues in the buried core domains of histone proteins H3 and H4, which have no known function other than maintaining the three-dimensional structure of the protein. We find that perturbing protein stability (which is a measure of maintenance of the protein's structure) by mutating these residues compromises survival fitness in yeast. However, the stability conferred by buried amino acids of H3 alone cannot account for their evolutionary conservation, which is in striking contrast to other proteins where stability has been shown to be the driving force for sequence conservation. This conservation of H3 thus points to either new additional functions of H3 that have not been uncovered or a unique conservation mechanism that goes beyond survival pressure. These data therefore reveal a highly conserved domain that is distinct in its evolutionary conservation.

Stability tests on known and misfolded structures with discrete and all atom molecular dynamics simulations

The brevity of molecular dynamics simulations often limits their utility in developing and evaluating structural models of proteins. The duration of simulations can be increased greatly using discrete molecular dynamics (DMD). However, the trade off is that coarse graining, implicit solvent, and other time-saving procedures reduce the accuracy of DMD simulations. Here we address some of these issues by comparing results of DMD and conventional all atom MD simulations on proteins of known structure and misfolded proteins. DMD simulations were performed at a range of temperatures to identify a 'physiological' temperature for DMD that mimicked molecular motions of conventional MD simulations at 310K. We also compared results obtained with a new implicit solvent model developed here based on Miyazawa-Jernigan interaction pair potential to those obtained with a previously used model based on Kyte-Doolittle hydropathy scale. We compared DMD and all atom molecular dynamics with explicit water by simulating both correctly and incorrectly folded structures, and monomeric and dimeric α β-barrel structures to analyze the ability of these procedures to distinguish between good and bad models. Deviations from the correct structures were substantially greater with DMD, as would be expected from coarse-graining and longer simulation time. Deviations were smallest for β-strands and greatest for coiled loops. Structures of the incorrectly folded models were very poorly preserved during the DMD simulations; but both methods were able to distinguish between the correct and the incorrect structures based on differences in the magnitudes of the root mean squared deviation (RMSD) from the starting conformation.

Genome3D: a viewer-model framework for integrating and visualizing multi-scale epigenomic information within a three-dimensional genome

Background

New technologies are enabling the measurement of many types of genomic and epigenomic information at scales ranging from the atomic to nuclear. Much of this new data is increasingly structural in nature, and is often difficult to coordinate with other data sets. There is a legitimate need for integrating and visualizing these disparate data sets to reveal structural relationships not apparent when looking at these data in isolation.

Results

We have applied object-oriented technology to develop a downloadable visualization tool, Genome3D, for integrating and displaying epigenomic data within a prescribed three-dimensional physical model of the human genome. In order to integrate and visualize large volume of data, novel statistical and mathematical approaches have been developed to reduce the size of the data. To our knowledge, this is the first such tool developed that can visualize human genome in three-dimension. We describe here the major features of Genome3D and discuss our multi-scale data framework using a representative basic physical model. We then demonstrate many of the issues and benefits of multi-resolution data integration.

Conclusions

Genome3D is a software visualization tool that explores a wide range of structural genomic and epigenetic data. Data from various sources of differing scales can be integrated within a hierarchical framework that is easily adapted to new developments concerning the structure of the physical genome. In addition, our tool has a simple annotation mechanism to incorporate non-structural information. Genome3D is unique is its ability to manipulate large amounts of multi-resolution data from diverse sources to uncover complex and new structural relationships within the genome.

A didactic model of macromolecular crowding effects on protein folding

A didactic model is presented to illustrate how the effect of macromolecular crowding on protein folding and association is modeled using current analytical theory and discrete molecular dynamics. While analytical treatments of crowding may consider the effect as a potential of average force acting to compress a polypeptide chain into a compact state, the use of simulations enables the presence of crowding reagents to be treated explicitly. Using an analytically solvable toy model for protein folding, an approximate statistical thermodynamic method is directly compared to simulation in order to gauge the effectiveness of current analytical crowding descriptions. Both methodologies are in quantitative agreement under most conditions, indication that both current theory and simulation methods are capable of recapitulating aspects of protein folding even by utilizing a simplistic protein model.

Engineered allosteric activation of kinases in living cells

Studies of cellular and tissue dynamics benefit greatly from tools that can control protein activity with specificity and precise timing in living systems. We describe here a new approach to confer allosteric regulation specifically on the catalytic activity of kinases. A highly conserved portion of the kinase catalytic domain is modified with a small protein insert that inactivates catalytic activity, but does not affect other protein interactions. Catalytic activity is restored by addition of rapamycin or non-immunosuppresive analogs (Fig. 1A). We demonstrate the approach by specifically activating focal adhesion kinase (FAK) within minutes in living cells, thereby demonstrating a novel role for FAK in regulation of membrane dynamics. Molecular modeling and mutagenesis indicate that the protein insert reduces activity by increasing the flexibility of the catalytic domain. Drug binding restores activity by increasing rigidity. Successful regulation of Src and p38 suggest that modification of this highly conserved site will be applicable to other kinases.

Regulatory insertion removal restores maturation, stability and function of DeltaF508 CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) epithelial anion channel is a large multidomain membrane protein that matures inefficiently during biosynthesis. Its assembly is further perturbed by the deletion of F508 from the first nucleotide-binding domain (NBD1) responsible for most cystic fibrosis. The mutant polypeptide is recognized by cellular quality control systems and is proteolyzed. CFTR NBD1 contains a 32-residue segment termed the regulatory insertion (RI) not present in other ATP-binding cassette transporters. We report here that RI deletion enabled F508 CFTR to mature and traffic to the cell surface where it mediated regulated anion efflux and exhibited robust single chloride channel activity. Long-term pulse-chase experiments showed that the mature DeltaRI/DeltaF508 had a T(1/2) of approximately 14 h in cells, similar to the wild type. RI deletion restored ATP occlusion by NBD1 of DeltaF508 CFTR and had a strong thermostabilizing influence on the channel with gating up to at least 40 degrees C. None of these effects of RI removal were achieved by deletion of only portions of RI. Discrete molecular dynamics simulations of NBD1 indicated that RI might indirectly influence the interaction of NBD1 with the rest of the protein by attenuating the coupling of the F508-containing loop with the F1-like ATP-binding core subdomain so that RI removal overcame the perturbations caused by F508 deletion. Restriction of RI to a particular conformational state may ameliorate the impact of the disease-causing mutation.

PRIMO/PRIMONA: a coarse-grained model for proteins and nucleic acids that preserves near-atomistic accuracy

The new coarse graining model PRIMO/PRIMONA for proteins and nucleic acids is proposed. This model combines one to several heavy atoms into coarse-grained sites that are chosen to allow an analytical, high-resolution reconstruction of all-atom models based on molecular bonding geometry constraints. The accuracy of proposed reconstruction method in terms of structure and energetics is tested and compared with other popular reconstruction methods for a variety of protein and nucleic acid test sets.

Macromolecular crowding induces polypeptide compaction and decreases folding cooperativity

A cell's interior is comprised of macromolecules that can occupy up to 40% of its available volume. Such crowded environments can influence the stability of proteins and their rates of reaction. Using discrete molecular dynamics simulations, we investigate how both the size and number of neighboring crowding reagents affect the thermodynamic and folding properties of structurally diverse proteins. We find that crowding induces higher compaction of proteins. We also find that folding becomes less cooperative with the introduction of crowders into the system. The crowders may induce alternative non-native protein conformations, thus creating barriers for protein folding in highly crowded media.

Protein flexibility from discrete molecular dynamics simulations using quasi-physical potentials

We have applied all atoms discrete molecular dynamics (DMD) based on a quasi-physical potential to study the flexibility of an extended set of proteins for which atomistic MD simulations are available. The method uses pure physical potentials supplemented by information on secondary structure and despite its simplicity is able to reproduce with good accuracy the dynamics of proteins in solution. The method presents a clear improvement with respect to coarse-grained methods based on structure potentials and opens the possibility to explore dynamics of proteins out from the equilibrium and to trace conformational changes induced by interaction of proteins with both small and macromolecular ligands.

Allosteric Conformational Transition in Adenylate Kinase: Dynamic Correlations and Implication for Allostery

An essential aspect of protein science is to determine the deductive relationship between structure, dynamics, and various sets of functions. The role of dynamics is currently challenging our understanding of protein functions, both experimentally and theoretically. To verify the internal fluctuations and dynamics correlations in an enzyme protein undergoing conformational transitions, we have applied a coarse-grained dynamics algorithm using the elastic network model for adenylate kinase. Normal mode analysis reveals possible dynamical and allosteric pathways for the transition between the open and the closed states of adenylate kinase. As the ligands binding induces significant flexibility changes of the nucleotides monophosphate (NMP) domain and adenosine triphosphate (ATP) domain, the diagonalized correlation between different structural transition states shows that most correlated motions occur between the NMP domain and the helices surrounding the ATP domain. The simultaneous existence of positive and negative correlations indicates that the conformational changes of adenylate kinase take place in an allosteric manner. Analyses of the cumulated normal mode overlap coefficients and long-range correlated motion provide new insights of operating mechanisms and dynamics of adenylate kinase. They also suggest a quantitative dynamics criterion for determining the allosteric cooperativity, which may be applicable to other proteins.

Relaxation dynamics of nucleosomal DNA

Recent experimental and theoretical evidence demonstrates that proteins and water in the hydration layer can follow complex stretched exponential or power law relaxation dynamics. Here, we report on a 50 ns all atom molecular dynamics (MD) simulation of the yeast nucleosome, where the interactions between DNA, histones, surrounding water and ions are explicitly included. DNA interacts with the histone core in 14 locations, approximately every 10.4 base pairs. We demonstrate that all sites of interaction exhibit anomalously slow power law relaxation, extending up to 10 ns, while fast exponential relaxation dynamics of hundreds of picoseconds applies to DNA regions outside these locations. The appearance of 1/f(alpha) noise or pink noise in DNA dynamics is ubiquitous. For histone-bound nucleotide dynamics alpha --> 1 and is a signature of complexity of the protein-DNA interactions. For control purposes two additional DNA simulations free of protein are conducted. Both utilize the same sequence of DNA, as found the in the nucleosome. In one simulation the initial conformation of the double helix is a straight B-form. In the other, the initial conformation is super helical. Neither of these simulations exhibits the variation of alpha as a function of position, the measure of power law for dynamical behavior, which we observe in the nucleosome simulation. The unique correspondence (high alpha to DNA-histone interaction sites, low alpha to free DNA sites), suggests that alpha may be an important and new quantification of protein-DNA interactions for future experiments.

Multiscale approaches for studying energy transduction in dynein

Cytoplasmic dynein is an important motor that drives all minus-end directed movement along microtubules. Dynein is a complex motor whose processive motion is driven by ATP-hydrolysis. Dynein's run length has been measured to be several millimetres with typical velocities in the order of a few nanometres per second. Therefore, the average time between steps is a fraction of a second. When this time scale is compared with typical time scales for protein side chain and backbone movements (approximately 10(-9) s and approximately 10(-5) s, respectively), it becomes clear that a multi-timescale modelling approach is required to understand energy transduction in this protein. Here, we review recent efforts to use computational and mathematical modelling to understand various aspects of dynein's chemomechanical cycle. First, we describe a structural model of dynein's motor unit showing a heptameric organization of the motor subunits. Second, we describe our molecular dynamics simulations of the motor unit that are used to investigate the dynamics of the various motor domains. Third, we present a kinetic model of the coordination between the two dynein heads. Lastly, we investigate the various potential geometries of the dimer during its hydrolytic and stepping cycle.

Molecular renormalization group coarse-graining of polymer chains: application to double-stranded DNA

Coarse-graining of atomistic force fields allows us to investigate complex biological problems, occurring at long timescales and large length scales. In this work, we have developed an accurate coarse-grained model for double-stranded DNA chain, derived systematically from atomistic simulations. Our approach is based on matching correlators obtained from atomistic and coarse-grained simulations, for observables that explicitly enter the coarse-grained Hamiltonian. We show that this requirement leads to equivalency of the corresponding partition functions, resulting in a one-step renormalization. Compared to prior works exploiting similar ideas, the main novelty of this work is the introduction of a highly compact set of Hamiltonian basis functions, based on molecular interaction potentials. We demonstrate that such compactification allows us to reproduce many-body effects, generated by one-step renormalization, at low computational cost. In addition, compact Hamiltonians greatly increase the likelihood of finding unique solutions for the coarse-grained force-field parameter values. By successfully applying our molecular renormalization group coarse-graining technique to double-stranded DNA, we solved, for the first time, a long-standing problem in coarse-graining polymer systems, namely, how to accurately capture the correlations among various polymeric degrees of freedom. Excellent agreement is found among atomistic and coarse-grained distribution functions for various structural observables, including those not included in the Hamiltonian. We also suggest higher-order generalization of this method, which may allow capturing more subtle correlations in biopolymer dynamics.

Determinants of histone H4 N-terminal domain function during nucleosomal array oligomerization: roles of amino acid sequence, domain length, and charge density

Mg(2+)-dependent oligomerization of nucleosomal arrays is correlated with higher order folding transitions that stabilize chromosome structure beyond the 30-nm diameter fiber. In the present studies, we have employed a novel mutagenesis-based approach to identify the macromolecular determinants that control H4 N-terminal domain (NTD) function during oligomerization. Core histones were engineered in which 1) the H2A, H2B, and H3 NTDs were swapped onto the H4 histone fold; 2) the length of the H4 NTD and the H2A NTD on the H4 histone fold, were increased; 3) the charge density of the NTDs on the H4 histone fold was increased or decreased; and 4) the H4 NTD was placed on the H2B histone fold. Model nucleosomal arrays were assembled from wild type and mutant core histone octamers, and Mg(2+)-dependent oligomerization was characterized. The results demonstrated that the H2B and H3 NTDs could replace the H4 NTD, as could the H2A NTD if it was duplicated to the length of the native H4 NTD. Arrays oligomerized at lower salt concentrations as the length of the NTD on the H4 histone fold was increased. Mutations that decreased the NTD charge density required more Mg(2+) to oligomerize, whereas mutants that increased the charge density required less salt. Finally, the H4 NTD functioned differently when attached to the H2B histone fold than the H4 histone fold. These studies have revealed new insights into the biochemical basis for H4 NTD effects on genome architecture as well as the protein chemistry that underlies the function of the intrinsically disordered H4 NTD.

Coarse-grained dynamics of the receiver domain of NtrC: fluctuations, correlations and implications for allosteric cooperativity

Receiver domains are key molecular switches in bacterial signaling. Structural studies have shown that the receiver domain of the nitrogen regulatory protein C (NtrC) exists in a conformational equilibrium encompassing both inactive and active states, with phosphorylation of Asp54 allosterically shifting the equilibrium towards the active state. To analyze dynamical fluctuations and correlations in NtrC as it undergoes activation, we have applied a coarse-grained dynamics algorithm using elastic network models. Normal mode analysis reveals possible dynamical pathways for the transition of NtrC from the inactive state to the active state. The diagonalized correlation between the inactive and the active (phosphorylated) state shows that most correlated motions occur around the active site of Asp54 and in the region Thr82 to Tyr101. This indicates a coupled correlation of dynamics in the "Thr82-Tyr101" motion. With phosphorylation inducing significant flexibility changes around the active site and alpha3 and alpha4 helices, we find that this activation makes the active-site region and the loops of alpha3/beta4 and alpha4/beta5 more stable. This means that phosphorylation entropically favors the receiver domain in its active state, and the induced conformational changes occur in an allosteric manner. Analyses of the local flexibility and long-range correlated motion also suggest a dynamics criterion for determining the allosteric cooperativity of NtrC, and may be applicable to other proteins.

DNA sequence-dependent variation in nucleosome structure, stability, and dynamics detected by a FRET-based analysisThis paper is one of a selection of papers published in this Special Issue, entitled 29th Annual International Asilomar Chromatin and Chromosomes Conference, and has undergone the Journal’s usual peer review process

Förster resonance energy transfer (FRET) techniques provide powerful and sensitive methods for the study of conformational features in biomolecules. Here, we review FRET-based studies of nucleosomes, focusing particularly on our work comparing the widely used nucleosome standard, 5S rDNA, and 2 promoter-derived regulatory element-containing nucleosomes, mouse mammary tumor virus (MMTV)-B and GAL10. Using several FRET approaches, we detected significant DNA sequence-dependent structure, stability, and dynamics differences among the three. In particular, 5S nucleosomes and 5S H2A/H2B-depleted nucleosomal particles have enhanced stability and diminished DNA dynamics, compared with MMTV-B and GAL10 nucleosomes and particles. H2A/H2B-depleted nucleosomes are of interest because they are produced by the activities of many transcription-associated complexes. Significant location-dependent (intranucleosomal) stability and dynamics variations were also observed. These also vary among nucleosome types. Nucleosomes restrict regulatory factor access to DNA, thereby impeding genetic processes. Eukaryotic cells possess mechanisms to alter nucleosome structure, to generate DNA access, but alterations often must be targeted to specific nucleosomes on critical regulatory DNA elements. By endowing specific nucleosomes with intrinsically higher DNA accessibility and (or) enhanced facility for conformational transitions, DNA sequence-dependent nucleosome dynamics and stability variations have the potential to facilitate nucleosome recognition and, thus, aid in the crucial targeting process. This and other nucleosome structure and function conclusions from FRET analyses are discussed.

Epigenetics and the control of epithelial sodium channel expression in collecting duct

In eukaryotic nuclei, genomic DNA is compacted with histone and nonhistone proteins into a dynamic polymer termed chromatin. Reorganization of chromatin structure through histone modifications, the action of chromatin factors, or DNA methylation, can profoundly change gene expression. These epigenetic modifications allow heritable and potentially reversible changes in gene functioning to occur without altering the DNA sequence, thus extending the information potential of the genetic code. This review provides an introduction to epigenetic concepts for renal investigators and an overview of our work detailing an epigenetic pathway for aldosterone signaling and the control of epithelial Na(+) channel-alpha (ENaCalpha) subunit gene expression in the collecting duct. This new pathway involves a nuclear repressor complex, consisting of histone H3 Lys-79 methyltransferase disruptor of telomeric silencing-1a (Dot1a), ALL1 fused gene from chromosome 9 (Af9), a sequence-specific DNA-binding protein that binds the ENaCalpha promoter, and potentially other nuclear proteins. This complex regulates targeted histone H3 Lys-79 methylation of chromatin associated with the ENaCalpha promoter, thereby suppressing its transcriptional activity. Aldosterone disrupts the Dot1a-Af9 interaction by serum- and glucocorticoid-induced kinase-1 phosphorylation of Af9, and inhibits Dot1a and Af9 expression, resulting in histone H3 Lys-79 hypomethylation at specific subregions, and derepression of the ENaCalpha promoter. The Dot1a-Af9 pathway may also be involved in the control of genes implicated in renal fibrosis and hypertension.

Coarse-grained force field for the nucleosome from self-consistent multiscaling

A coarse-grained simulation model for the nucleosome is developed, using a methodology modified from previous work on the ribosome. Protein residues and DNA nucleotides are represented as beads, interacting through harmonic (for neighboring) or Morse (for nonbonded) potentials. Force-field parameters were estimated by Boltzmann inversion of the corresponding radial distribution functions obtained from a 5-ns all-atom molecular dynamics (MD) simulation, and were refined to produce agreement with the all-atom MD simulation. This self-consistent multiscale approach yields a coarse-grained model that is capable of reproducing equilibrium structural properties calculated from a 50-ns all-atom MD simulation. This coarse-grained model speeds up nucleosome simulations by a factor of 10(3) and is expected to be useful in examining biologically relevant dynamical nucleosome phenomena on the microsecond timescale and beyond.

Active nuclear receptors exhibit highly correlated AF-2 domain motions

Nuclear receptor ligand binding domains (LBDs) convert ligand binding events into changes in gene expression by recruiting transcriptional coregulators to a conserved activation function-2 (AF-2) surface. While most nuclear receptor LBDs form homo- or heterodimers, the human nuclear receptor pregnane X receptor (PXR) forms a unique and essential homodimer and is proposed to assemble into a functional heterotetramer with the retinoid X receptor (RXR). How the homodimer interface, which is located 30 A from the AF-2, would affect function at this critical surface has remained unclear. By using 20- to 30-ns molecular dynamics simulations on PXR in various oligomerization states, we observed a remarkably high degree of correlated motion in the PXR-RXR heterotetramer, most notably in the four helices that create the AF-2 domain. The function of such correlation may be to create "active-capable" receptor complexes that are ready to bind to transcriptional coactivators. Indeed, we found in additional simulations that active-capable receptor complexes involving other orphan or steroid nuclear receptors also exhibit highly correlated AF-2 domain motions. We further propose a mechanism for the transmission of long-range motions through the nuclear receptor LBD to the AF-2 surface. Taken together, our findings indicate that long-range motions within the LBD scaffold are critical to nuclear receptor function by promoting a mobile AF-2 state ready to bind coactivators.

Ab initio RNA folding by discrete molecular dynamics: from structure prediction to folding mechanisms

RNA molecules with novel functions have revived interest in the accurate prediction of RNA three-dimensional (3D) structure and folding dynamics. However, existing methods are inefficient in automated 3D structure prediction. Here, we report a robust computational approach for rapid folding of RNA molecules. We develop a simplified RNA model for discrete molecular dynamics (DMD) simulations, incorporating base-pairing and base-stacking interactions. We demonstrate correct folding of 150 structurally diverse RNA sequences. The majority of DMD-predicted 3D structures have <4 A deviations from experimental structures. The secondary structures corresponding to the predicted 3D structures consist of 94% native base-pair interactions. Folding thermodynamics and kinetics of tRNA(Phe), pseudoknots, and mRNA fragments in DMD simulations are in agreement with previous experimental findings. Folding of RNA molecules features transient, non-native conformations, suggesting non-hierarchical RNA folding. Our method allows rapid conformational sampling of RNA folding, with computational time increasing linearly with RNA length. We envision this approach as a promising tool for RNA structural and functional analyses.

Exploring the suitability of coarse-grained techniques for the representation of protein dynamics

A systematic study of two coarse-grained techniques for the description of protein dynamics is presented. The two techniques exploit either Brownian or discrete molecular dynamics algorithms applied in the context of simple C(alpha)-C(alpha) potentials, like those used in coarse-grained normal mode analysis. Coarse-grained simulations of the flexibility of protein metafolds are compared to those computed with fully atomistic molecular dynamics simulations using state-of-the-art physical potentials and explicit solvent. Both coarse-grained models efficiently capture critical features of the protein dynamics.

Probing protein aggregation using discrete molecular dynamics

Understanding the role of biomolecular dynamics in cellular processes leading to human diseases and the ability to rationally manipulate these processes is of fundamental importance in scientific research. The last decade has witnessed significant progress in probing biophysical behavior of proteins. However, we are still limited in understanding how changes in protein dynamics and inter-protein interactions occurring in short length- and time-scales lead to aberrations in their biological function. Bridging this gap in biology probed using computer simulations marks a challenging frontier in computational biology. Here we examine hypothesis-driven simplified protein models in conjunction with discrete molecular dynamics in the study of protein aggregation, implicated in series of neurodegenerative diseases, such as Alzheimer's and Huntington's diseases. Discrete molecular dynamics simulations of simplified protein models have emerged as a powerful methodology with its ability to bridge the gap in time and length scales from protein dynamics to aggregation, and provide an indispensable tool for probing protein aggregation.

Coarse-grained models reveal functional dynamics--II. Molecular dynamics simulation at the coarse-grained level--theories and biological applications

Molecular dynamics (MD) simulation has remained the most indispensable tool in studying equilibrium/non-equilibrium conformational dynamics since its advent 30 years ago. With advances in spectroscopy accompanying solved biocomplexes in growing sizes, sampling their dynamics that occur at biologically interesting spatial/temporal scales becomes computationally intractable; this motivated the use of coarse-grained (CG) approaches. CG-MD models are used to study folding and conformational transitions in reduced resolution and can employ enlarged time steps due to the absence of some of the fastest motions in the system. The Boltzmann-Inversion technique, heavily used in parameterizing these models, provides a smoothed-out effective potential on which molecular conformation evolves at a faster pace thus stretching simulations into tens of microseconds. As a result, a complete catalytic cycle of HIV-1 protease or the assembly of lipid-protein mixtures could be investigated by CG-MD to gain biological insights. In this review, we survey the theories developed in recent years, which are categorized into Folding-based and Molecular-Mechanics-based. In addition, physical bases in the selection of CG beads/time-step, the choice of effective potentials, representation of solvent, and restoration of molecular representations back to their atomic details are systematically discussed.

Recent advances in the study of nucleic acid flexibility by molecular dynamics

The recent use of molecular dynamics (MD) simulations to study flexibility of nucleic acids has been reviewed from an analysis of the publications appearing in the past two years (from 2005 till date). Despite the existence of some unsolved problems in the methodologies, these years have been witness to major advances in the field. Based on a critical review of the most recent contributions, excitement exists on the expected evolution of the field in the next years.

Diminished self-chaperoning activity of the DeltaF508 mutant of CFTR results in protein misfolding

The absence of a functional ATP Binding Cassette (ABC) protein called the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) from apical membranes of epithelial cells is responsible for cystic fibrosis (CF). Over 90% of CF patients carry at least one mutant allele with deletion of phenylalanine at position 508 located in the N-terminal nucleotide binding domain (NBD1). Biochemical and cell biological studies show that the ΔF508 mutant exhibits inefficient biosynthetic maturation and susceptibility to degradation probably due to misfolding of NBD1 and the resultant misassembly of other domains. However, little is known about the direct effect of the Phe508 deletion on the NBD1 folding, which is essential for rational design strategies of cystic fibrosis treatment. Here we show that the deletion of Phe508 alters the folding dynamics and kinetics of NBD1, thus possibly affecting the assembly of the complete CFTR. Using molecular dynamics simulations, we find that meta-stable intermediate states appearing on wild type and mutant folding pathways are populated differently and that their kinetic accessibilities are distinct. The structural basis of the increased misfolding propensity of the ΔF508 NBD1 mutant is the perturbation of interactions in residue pairs Q493/P574 and F575/F578 found in loop S7-H6. As a proof-of-principle that the S7-H6 loop conformation can modulate the folding kinetics of NBD1, we virtually design rescue mutations in the identified critical interactions to force the S7-H6 loop into the wild type conformation. Two redesigned NBD1-ΔF508 variants exhibited significantly higher folding probabilities than the original NBD1-ΔF508, thereby partially rescuing folding ability of the NBD1-ΔF508 mutant. We propose that these observed defects in folding kinetics of mutant NBD1 may also be modulated by structures separate from the 508 site. The identified structural determinants of increased misfolding propensity of NBD1-ΔF508 are essential information in correcting this pathogenic mutant.

Deletion of a single residue, phenylalanine at position 508, in the first nucleotide binding domain (NBD1) of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) is present in approximately 90% of cystic fibrosis (CF) patients. Experiments show that this mutant protein exhibits inefficient biosynthetic maturation and susceptibility to degradation probably due to misfolding of NBD1 and the resultant incorrect interactions of other domains. However, little is known about the direct effect of the Phe508 deletion on NBD1 folding. Here, using molecular dynamics simulations of NBD1-WT, NBD1-F508A, and NBD1-ΔF508, we show that the deletion of Phe508 indeed alters the kinetics of NBD1 folding. We also find that the intermediate states appearing on wild type and mutant folding pathways are populated differently and that their kinetic accessibilities are distinct. Moreover, we identified critical interactions not necessarily localized near position 508, such as Q493/P574 and F575/F587, to be significant structural elements influencing the kinetic difference between wild type and mutant NBD1. We propose that these observed alterations in folding kinetics of mutant NBD1 result in misassembly of the whole multi-domain protein, thereby causing its premature degradation.

Protein folding: then and now

Over the past three decades the protein folding field has undergone monumental changes. Originally a purely academic question, how a protein folds has now become vital in understanding diseases and our abilities to rationally manipulate cellular life by engineering protein folding pathways. We review and contrast past and recent developments in the protein folding field. Specifically, we discuss the progress in our understanding of protein folding thermodynamics and kinetics, the properties of evasive intermediates, and unfolded states. We also discuss how some abnormalities in protein folding lead to protein aggregation and human diseases.

Sequence-dependent variations associated with H2A/H2B depletion of nucleosomes

Mechanisms that can alter nucleosome structure to enhance DNA accessibility are of great interest because of their potential involvement in genomic processes. One such mechanism is H2A/H2B release from nucleosomes; it occurs in vivo and is involved in the in vitro activities of several transcription-associated complexes. Using fluorescence approaches based on Förster resonance energy transfer, we previously detected sequence-dependent structure/stability variations between 5S and two types of promoter nucleosomes (from yeast GAL10 or mouse mammary tumor virus promoters). Those variations included differing responses when nucleosomes were diluted to concentrations (sub-nM) known to produce H2A/H2B loss. Here, we show that treatment of these same three types of nucleosomes with the histone chaperone yNAP-1, which causes H2A/H2B release from nucleosomes in vitro, produces the same differential Förster resonance energy transfer responses, again demonstrating sequence-dependent variations associated with conditions that produce H2A/H2B loss. Single-molecule population data indicate that DNA dynamics on the particles produced by diluting nucleosomes to sub-nM concentrations follow two-state behavior. Rate information (determined by fluorescence correlation spectroscopy) suggests that these dynamics are enhanced in MMTV-B or GAL10 compared to 5S particles. Taken together, the results indicate that H2A/H2B loss has differing effects on 5S compared to these two promoter nucleosomes and the differences reflect sequence-dependent structure/stability variations in the depleted particles.

Global dynamics of newly constructed oligonucleosomes of conventional and variant H2A.Z histone

Background

Complexes of nucleosomes, which often occur in the gene promoter areas, are one of the fundamental levels of chromatin organization and thus are important for transcription regulation. Investigating the dynamic structure of a single nucleosome as well as nucleosome complexes is important for understanding transcription within chromatin. In a previous work, we highlighted the influence of histone variants on the functional dynamics of a single nucleosome using normal mode analysis developed by Bahar et al. The present work further analyzes the dynamics of nucleosome complexes (nucleosome oligomers or oligonucleosomes) such as dimer, trimer and tetramer (beads on a string model) with conventional core histones as well as with the H2A.Z histone variant using normal mode analysis.

Results

The global dynamics of oligonucleosomes reveal larger amplitude of motion within the nucleosomes that contain the H2A.Z variant with in-planar and out-of-planar fluctuations as the common mode of relaxation. The docking region of H2A.Z and the L1:L1 interactions between H2A.Z monomers of nucleosome (that are responsible for the highly stable nucleosome containing variant H2A.Z-histone) are highly dynamic throughout the first two dynamic modes.

Conclusion

Dissection of the dynamics of oligonucleosomes discloses in-plane as well as out-of-plane fluctuations as the common mode of relaxation throughout the global motions. The dynamics of individual nucleosomes and the combination of the relaxation mechanisms expressed by the individual nucleosome are quite interesting and highly dependent on the number of nucleosome fragments present in the complexes. Distortions generated by the non-planar dynamics influence the DNA conformation, and hence the histone-DNA interactions significantly alter the dynamics of the DNA. The variant H2A.Z histone is a major source of weaker intra- and inter-molecular correlations resulting in more disordered motions.

Atomic force microscopy imaging of SWI/SNF action: mapping the nucleosome remodeling and sliding

We propose a combined experimental (atomic force microscopy) and theoretical study of the structural and dynamical properties of nucleosomes. In contrast to biochemical approaches, this method allows us to determine simultaneously the DNA-complexed length distribution and nucleosome position in various contexts. First, we show that differences in the nucleoproteic structure observed between conventional H2A and H2A.Bbd variant nucleosomes induce quantitative changes in the length distribution of DNA-complexed with histones. Then, the sliding action of remodeling complex SWI/SNF is characterized through the evolution of the nucleosome position and wrapped DNA length mapping. Using a linear energetic model for the distribution of DNA-complexed length, we extract the net-wrapping energy of DNA onto the histone octamer and compare it to previous studies.

Structural flexibility of the nucleosome core particle at atomic resolution studied by molecular dynamics simulation

Comparative explicit solvent molecular dynamics (MD) simulations have been performed on a complete nucleosome core particle with and without N-terminal histone tails for more than 20 ns. Main purpose of the simulations was to study the dynamics of mobile elements such as histone N-terminal tails and how packing and DNA-bending influences the fine structure and dynamics of DNA. Except for the tails, histone and DNA molecules stayed on average close to the crystallographic start structure supporting the quality of the current force field approach. Despite the packing strain, no increase of transitions to noncanonical nucleic acid backbone conformations compared to regular B-DNA was observed. The pattern of kinks and bends along the DNA remained close to the experiment overall. In addition to the local dynamics, the simulations allowed the analysis of the superhelical mobility indicating a limited relative mobility of DNA segments separated by one superhelical turn (mean relative displacement of approximately +/-0.2 nm, mainly along the superhelical axis). An even higher rigidity was found for relative motions (distance fluctuations) of segments separated by half a superhelical turn (approximately +/-0.1 nm). The N-terminal tails underwent dramatic conformational rearrangements on the nanosecond time scale toward partially and transiently wrapped states around the DNA. Many of the histone tail changes corresponded to coupled association and folding events from fully solvent-exposed states toward complexes with the major and minor grooves of DNA. The simulations indicate that the rapid conformational changes of the tails can modulate the DNA accessibility within a few nanoseconds.