No Arabic abstract
Guanylate binding proteins (GBPs) are soluble dynamin-like proteins with structured domains that undergo a conformational transition for GTP-controlled oligomerization to exert their function as part of the innate immune system of mammalian cells - attacking intra-cellular parasites by disrupting their membranes. The structural basis and mechanism of this process is unknown. Therefore, we apply neutron spin echo, X-ray scattering, fluorescence, and EPR spectroscopy as techniques for integrative dynamic structural biology to human GBP1 (hGBP1). We mapped hGBP1s essential dynamics from nanoseconds to milliseconds by motional spectra of sub-domains. We find a GTP-independent flexibility of the C-terminal effector domain in the $mu$s-regime and structurally characterize conformers being essential that hGBP1 can open like a pocketknife for oligomerization. This unveils the intrinsic flexibility, a GTP-triggered association of the GTPase-domains and assembly-dependent GTP-hydrolysis as functional design principles of hGBP1 that control its reversible oligomerization in polar assemblies and the subsequent formation of condensates.
We outline recent developments in artificial intelligence (AI) and machine learning (ML) techniques for integrative structural biology of intrinsically disordered proteins (IDP) ensembles. IDPs challenge the traditional protein structure-function paradigm by adapting their conformations in response to specific binding partners leading them to mediate diverse, and often complex cellular functions such as biological signaling, self organization and compartmentalization. Obtaining mechanistic insights into their function can therefore be challenging for traditional structural determination techniques. Often, scientists have to rely on piecemeal evidence drawn from diverse experimental techniques to characterize their functional mechanisms. Multiscale simulations can help bridge critical knowledge gaps about IDP structure function relationships - however, these techniques also face challenges in resolving emergent phenomena within IDP conformational ensembles. We posit that scalable statistical inference techniques can effectively integrate information gleaned from multiple experimental techniques as well as from simulations, thus providing access to atomistic details of these emergent phenomena.
The elastic network (EN) is a prime model that describes the long-time dynamics of biomolecules. However, the use of harmonic potentials renders this model insufficient for studying large conformational changes of proteins (e.g. stretching of proteins, folding and thermal unfolding). Here, we extend the capabilities of the EN model by using a harmonic approximation described by Lennard-Jones (LJ) interactions for far contacts and native contacts obtained from the standard overlap criterion as in the case of Go-like models. While our model is validated against the EN model by reproducing the equilibrium properties for a number of proteins, we also show that the model is suitable for the study of large conformation changes by providing various examples. In particular, this is illustrated on the basis of pulling simulations that predict with high accuracy the experimental data on the rupture force of the studied proteins. Furthermore, in the case of DDFLN4 protein, our pulling simulations highlight the advantages of our model with respect to Go-like approaches, where the latter fail to reproduce previous results obtained by all-atom simulations that predict an additional characteristic peak for this protein. In addition, folding simulations of small peptides yield different folding times for alpha-helix and beta-hairpin, in agreement with experiment, in this way providing further opportunities for the application of our model in studying large conformational changes of proteins. In contrast to the EN model, our model is suitable for both normal mode analysis and molecular dynamics simulation. We anticipate that the proposed model will find applications in a broad range of problems in biology, including, among others, protein folding and thermal unfolding.
Previous studies of the flexibilities of ancestral proteins suggests that proteins evolve their function by altering their native state ensemble. Here we propose a more direct method of visualizing this by measuring the changes in the vibrational density of states (VDOS) of proteins as they evolve. Through analysis of VDOS profiles of ancestral and extant proteins we observe that $beta$-lactamase and thioredoxins evolve by altering their density of states in the terahertz region. Particularly, the shift in VDOS profiles between ancestral and extant proteins suggests that nature utilize dynamic allostery for functional evolution. Moreover, we also show that VDOS profile of individual position can be used to describe the flexibility changes, particularly those without any amino acid substitution.
The Virtual Institute for Integrative Biology (VIIB) is a Latin American initiative for achieving global collaborative e-Science in the areas of bioinformatics, genome biology, systems biology, metagenomics, medical applications and nanobiotechnolgy. The scientific agenda of VIIB includes: construction of databases for comparative genomics, the AlterORF database for alternate open reading frames discovery in genomes, bioinformatics services and protein simulations for biotechnological and medical applications. Human resource development has been promoted through co-sponsored students and shared teaching and seminars via video conferencing. E-Science challenges include: interoperability and connectivity concerns, high performance computing limitations, and the development of customized computational frameworks and flexible workflows to efficiently exploit shared resources without causing impediments to the user. Outreach programs include training workshops and classes for high school teachers and students and the new Adopt-a-Gene initiative. The VIIB has proved an effective way for small teams to transcend the critical mass problem, to overcome geographic limitations, to harness the power of large scale, collaborative science and improve the visibility of Latin American science It may provide a useful paradigm for developing further e-Science initiatives in Latin America and other emerging regions.
We propose a stochastic model for gene transcription coupled to DNA supercoiling, where we incorporate the experimental observation that polymerases create supercoiling as they unwind the DNA helix, and that these enzymes bind more favourably to regions where the genome is unwound. Within this model, we show that when the transcriptionally induced flux of supercoiling increases, there is a sharp crossover from a regime where torsional stresses relax quickly and gene transcription is random, to one where gene expression is highly correlated and tightly regulated by supercoiling. In the latter regime, the model displays transcriptional bursts, waves of supercoiling, and up-regulation of divergent or bidirectional genes. It also predicts that topological enzymes which relax twist and writhe should provide a pathway to down-regulate transcription. This article has been published in Physical Review Letters, May 2016.