Accurate modeling of DNA conformational flexibility by a multivariate Ising model.
Ontology highlight
ABSTRACT:
SUBMITTER: Liebl K
PROVIDER: S-EPMC8053917 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.
| S-EPMC3574579 | biostudies-literature
Project description:Kinase structures in the inactive "DFG-out" state provide a wealth of druggable binding site variants. The conformational plasticity of this state can be mainly described by different conformations of binding site-forming elements such as DFG motif, A-loop, P-loop, and ?C-helix. Compared to DFG-in structures, DFG-out structures are largely underrepresented in the Protein Data Bank (PDB). Thus, structure-based drug design efforts for DFG-out inhibitors may benefit from an efficient approach to generate an ensemble of DFG-out structures. Accordingly, the presented modeling pipeline systematically generates homology models of kinases in several DFG-out conformations based on a sophisticated creation of template structures that represent the major states of the flexible structural elements. Eighteen template classes were initially selected from all available kinase structures in the PDB and subsequently employed for modeling the entire kinome in different DFG-out variants by fusing individual structural elements to multiple chimeric template structures. Molecular dynamics simulations revealed that conformational transitions between the different DFG-out states generally do not occur within trajectories of a few hundred nanoseconds length. This underlines the benefits of the presented homology modeling pipeline to generate relevant conformations of "DFG-out" kinase structures for subsequent in silico screening or binding site analysis studies.
| S-EPMC6852311 | biostudies-literature
Project description:In this study we provide the first comprehensive map of DNA conformational flexibility in Saccharomyces cerevisiae complete genome. Flexibility plays a key role in DNA supercoiling and DNA/protein binding, regulating DNA transcription, replication or repair. Specific interest in flexibility analysis concerns its relationship with human genome instability. Enrichment in flexible sequences has been detected in unstable regions of human genome defined fragile sites, where genes map and carry frequent deletions and rearrangements in cancer. Flexible sequences have been suggested to be the determinants of fragile gene proneness to breakage; however, their actual role and properties remain elusive. Our in silico analysis carried out genome-wide via the StabFlex algorithm, shows the conserved presence of highly flexible regions in budding yeast genome as well as in genomes of other Saccharomyces sensu stricto species. Flexibile peaks in S. cerevisiae identify 175 ORFs mapping on their 3'UTR, a region affecting mRNA translation, localization and stability. (TA)n repeats of different extension shape the central structure of peaks and co-localize with polyadenylation efficiency element (EE) signals. ORFs with flexible peaks share common features. Transcripts are characterized by decreased half-life: this is considered peculiar of genes involved in regulatory systems with high turnover; consistently, their function affects biological processes such as cell cycle regulation or stress response. Our findings support the functional importance of flexibility peaks, suggesting that the flexible sequence may be derived by an expansion of canonical TAYRTA polyadenylation efficiency element. The flexible (TA)n repeat amplification could be the outcome of an evolutionary neofunctionalization leading to a differential 3'-end processing and expression regulation in genes with peculiar function. Our study provides a new support to the functional role of flexibility in genomes and a strategy for its characterization inside human fragile sites.
| S-EPMC4393290 | biostudies-literature
Project description:Single-stranded DNA (ssDNA) is an essential intermediate in various DNA metabolic processes and interacts with a large number of proteins. Due to its flexibility, the conformations of ssDNA in solution can only be described using statistical approaches, such as flexibly jointed or worm-like chain models. However, there is limited data available to assess such models quantitatively, especially for describing the flexibility of short ssDNA and RNA. To address this issue, we performed FRET studies of a series of oligodeoxythymidylates, (dT)N, over a wide range of salt concentrations and chain lengths (10 < or = N < or = 70 nucleotides), which provide systematic constraints for testing theoretical models. Unlike in mechanical studies where available ssDNA conformations are averaged out during the time it takes to perform measurements, fluorescence lifetimes may act here as an internal clock that influences fluorescence signals depending on how fast the ssDNA conformations fluctuate. A reasonably good agreement could be obtained between our data and the worm-like chain model provided that limited relaxations of the ssDNA conformations occur within the fluorescence lifetime of the donor. The persistence length thus estimated ranges from 1.5 nm in 2 M NaCl to 3 nm in 25 mM NaCl.
| S-EPMC1304100 | biostudies-literature
Project description:Molecular self-assembly is a promising approach to the preparation of nanostructures. DNA, in particular, shows great potential to be a superb molecular system. Synthetic DNA molecules have been programmed to assemble into a wide range of nanostructures. It is generally believed that rigidities of DNA nanomotifs (tiles) are essential for programmable self-assembly of well defined nanostructures. Recently, we have shown that adequate conformational flexibility could be exploited for assembling 3D objects, including tetrahedra, dodecahedra, and buckyballs, out of DNA three-point star motifs. In the current study, we have integrated tensegrity principle into this concept to assemble well defined, complex nanostructures in both 2D and 3D. A symmetric five-point-star motif (tile) has been designed to assemble into icosahedra or large nanocages depending on the concentration and flexibility of the DNA tiles. In both cases, the DNA tiles exhibit significant flexibilities and undergo substantial conformational changes, either symmetrically bending out of the plane or asymmetrically bending in the plane. In contrast to the complicated natures of the assembled structures, the approach presented here is simple and only requires three different component DNA strands. These results demonstrate that conformational flexibility could be explored to generate complex DNA nanostructures. The basic concept might be further extended to other biomacromolecular systems, such as RNA and proteins.
| S-EPMC2504817 | biostudies-literature
Project description:The translesion synthesis (TLS) DNA polymerases Rev1 and Polζ function together in DNA lesion bypass during DNA replication, acting as nucleotide inserter and extender polymerases, respectively. While the structural characterization of the Saccharomyces cerevisiae Polζ in its DNA-bound state has illuminated how this enzyme synthesizes DNA, a mechanistic understanding of TLS also requires probing conformational changes associated with DNA- and Rev1 binding. Here, we used single-particle cryo-electron microscopy to determine the structure of the apo Polζ holoenzyme. We show that compared with its DNA-bound state, apo Polζ displays enhanced flexibility that correlates with concerted motions associated with expansion of the Polζ DNA-binding channel upon DNA binding. We also identified a lysine residue that obstructs the DNA-binding channel in apo Polζ, suggesting a gating mechanism. The Polζ subunit Rev7 is a hub protein that directly binds Rev1 and is a component of several other protein complexes such as the shieldin DNA double-strand break repair complex. We analyzed the molecular interactions of budding yeast Rev7 in the context of Polζ and those of human Rev7 in the context of shieldin using a crystal structure of Rev7 bound to a fragment of the shieldin-3 protein. Overall, our study provides new insights into Polζ mechanism of action and the manner in which Rev7 recognizes partner proteins.
| S-EPMC8319531 | biostudies-literature
Project description:A central question in biophysics is whether DNA sequence affects its mechanical properties, which are thought to influence nucleosome positioning and gene expression. Previous attempts to answer this question have been hindered by an inability to resolve DNA structure and dynamics at the base-pair level. Here we use a model to measure the effects of sequence on the stability of DNA under bending. Sequence is shown to influence DNA's flexibility and its ability to form kinks, which arise when certain motifs slide past others to form non-native contacts. A mechanism for nucleosome positioning is proposed in which sequence influences DNA-histone binding by altering the local base-pair-level structure when subject to the curvature necessary for binding.
| S-EPMC3410732 | biostudies-literature
Project description:In the course of evolution, proteins show a remarkable conservation of their three-dimensional structure and their biological function, leading to strong evolutionary constraints on the sequence variability between homologous proteins. Our method aims at extracting such constraints from rapidly accumulating sequence data, and thereby at inferring protein structure and function from sequence information alone. Recently, global statistical inference methods (e.g. direct-coupling analysis, sparse inverse covariance estimation) have achieved a breakthrough towards this aim, and their predictions have been successfully implemented into tertiary and quaternary protein structure prediction methods. However, due to the discrete nature of the underlying variable (amino-acids), exact inference requires exponential time in the protein length, and efficient approximations are needed for practical applicability. Here we propose a very efficient multivariate Gaussian modeling approach as a variant of direct-coupling analysis: the discrete amino-acid variables are replaced by continuous Gaussian random variables. The resulting statistical inference problem is efficiently and exactly solvable. We show that the quality of inference is comparable or superior to the one achieved by mean-field approximations to inference with discrete variables, as done by direct-coupling analysis. This is true for (i) the prediction of residue-residue contacts in proteins, and (ii) the identification of protein-protein interaction partner in bacterial signal transduction. An implementation of our multivariate Gaussian approach is available at the website http://areeweb.polito.it/ricerca/cmp/code.
| S-EPMC3963956 | biostudies-literature
Project description:There has been a lot of work fitting Ising models to multivariate binary data in order to understand the conditional dependency relationships between the variables. However, additional covariates are frequently recorded together with the binary data, and may influence the dependence relationships. Motivated by such a dataset on genomic instability collected from tumor samples of several types, we propose a sparse covariate dependent Ising model to study both the conditional dependency within the binary data and its relationship with the additional covariates. This results in subject-specific Ising models, where the subject's covariates influence the strength of association between the genes. As in all exploratory data analysis, interpretability of results is important, and we use ?1 penalties to induce sparsity in the fitted graphs and in the number of selected covariates. Two algorithms to fit the model are proposed and compared on a set of simulated data, and asymptotic results are established. The results on the tumor dataset and their biological significance are discussed in detail.
| S-EPMC4425428 | biostudies-literature
Project description:By using the Ising model formulation for combinatorial optimization with 0-1 binary variables, we investigated the extent to which partisan gerrymandering is possible from a random but even distribution of supporters. Assuming that an electoral district consists of square subareas and that each subarea shares at least one edge with other subareas in the district, it was possible to find the most tilted assignment of seats in most cases. However, in cases where supporters' distribution included many enclaves, the maximum tilted assignment was usually found to fail. We also discussed the proposed algorithm is applicable to other fields such as the redistribution of delivery destinations.
| S-EPMC8655093 | biostudies-literature