Models of biomolecular structure and dynamics are often obtained by combining simulation or prediction approaches [e.g., comparative modeling, Molecular Dynamics (MD) simulations, Normal Mode Analysis (NMA), etc.] with experimental approaches [e.g., Nuclear Magnetic Resonance (NMR), X-ray crystallography, Small-Angle X-ray Scattering (SAXS), Electron Microscopy (EM), etc.] (Sali et al., ) (Figure 1). Such hybrid modeling extends the capabilities of experimental techniques, by enriching structural information and facilitating dynamics studies of biomolecules. This e-Book contains articles on methodological developments, applications, and challenges of hybrid biomolecular modeling that have been collected in the framework of the Frontiers Research Topic entitled “Hybrid Biomolecular Modeling.”
Figure 1
An example of hybrid modeling is fitting of structures of protein domains obtained by X-ray crystallography, NMR, or structure prediction into EM density maps of protein complexes (Kawabata, ; Birmanns et al., ; Tjioe et al., ; Yang et al., ). This allows obtaining high-resolution models of complexes when this cannot be achieved using a single experimental technique, as is often the case with large and flexible complexes (Cottevieille et al., ; Ciferri et al., ; Brown et al., ). This problem is addressed in the article by Habeck. The article is focused on a Bayesian inference approach to integrative biomolecular modeling by combining X-ray crystallography and cryo-EM data, but Habeck also discusses the computational challenges of this approach in a more general context of integrating other experimental data such as cross-linking/mass spectrometry and solid-state NMR data. The proposed approach is based on probabilistic models for cryo-EM maps and Markov chain Monte Carlo sampling of model structures from the posterior distribution.
Computational methods have been developed to predict the interactions between the protein subunits based on their shape complementary, electrostatic interactions, solvation energy, and statistical potential energy derived from the structural databases. This is known as molecular docking and one of its main challenges is the design of a reliable scoring function to assess the model quality. Inspired by the application of X-ray Free-Electron Lasers (XEFL) data in scoring models of conformational changes of complexes (Tokuhisa et al., ), Wang and Liu propose to use single-particle XFEL data for a more reliable scoring of models obtained by docking methods.
Computational approaches based on NMA or MD simulations have been developed to explore the conformational space of a model and identify the conformation (in this space) that best agrees with experimental data (Trabuco et al., ; Gorba and Tama, ; Jin et al., ). Devaurs et al. address this problem in the context of modeling based on experimental hydrogen/deuterium exchange (HDX) data. HDR data is often interpreted using an X-ray crystallography structure or a conformational ensemble obtained by MD simulations, though their correspondence with the HDR data is often not enough satisfactory. Devaurs et al. propose to select a single conformation that best fits the HDX data, from the conformational ensemble obtained with an extensive coarse-grained conformational sampling (of the given X-ray crystallography structure) that is biased with the information on the protein regions that produce the largest discrepancies with the HDX data.
Prischi and Pastore review an integrative structural modeling methodology that they have developed to determine the structure of weakly interacting molecular complexes. It combines NMR, SAXS, site directed mutagenesis, molecular docking, and MD simulations, and has been used by the authors and other groups to gain structural information on several iron-sulfur cluster (ISC) biogenesis complexes. The authors review these applications and discuss the advantages and limitations of this methodology as well as the future directions to improve it.
Woods et al. show a new application of an approach combining MD simulations, evolutionary sequence analysis, and Terahertz spectroscopy that they have developed to probe dynamics and allostery in rhodopsin. They show how the binding of the chlorophyll derivative, chlorin-e6 (Ce6) allosterically excites evolutionarily conserved communication pathways in rhodopsin that connect the ligand-binding site and the rest of the receptor.
Hsieh et al. present a NMA approach to analyze the dynamics of Dengue and Zika virus capsids based on their high-resolution cryo-EM models. They relate the differences identified in the dynamics of the E proteins in the two capsids to the differences observed in the two high-resolution models. They discuss the work that should be done in the future in order to fully characterize the dynamics of the two viruses.
Intrinsically disordered peptides and proteins present a challenge to experimental characterization of their functional conformations. Olson explores simulation techniques that could be used to build a computational framework for capturing conformational ensembles of such peptides and proteins. He explores temperature-based replica exchange methods for conformational ensemble sampling with implicit solvent models, as well as, explicit/implicit solvent hybrid replica exchange methods to capture the conformational ensemble of an intrinsically disordered peptide derived from the Ebola virus protein VP35. The author points out that intrinsically disordered peptides and proteins can be used as benchmarks to develop accurate methods for modeling conformational transitions.
The permeability of a cell membrane can be increased under the influence of an electric field of sufficient magnitude, which is known as membrane electroporation. Wriggers et al. address the problem of experimental and theoretical investigation of membrane electroporation. They extend, to the context of lipid bilayers and solvents, a statistical approach that they have originally developed for detecting allosteric signatures in MD simulations of well-structured proteins. This method is based on transforming time-domain information from MD trajectories into spatial heat maps that can be visualized on 3D molecular structures or in the form of interaction networks. The method is multiscale in the time domain and uses a mutual information approach for statistical bridging between the fast (local variables recorded by MD) and slow (global rate of change) time series. The mutual information method used with proteins was adapted to lipids and solvents by developing a new approach to probability density function estimation of random variables, which was described in a separate article (Kovacs et al.) in this e-Book.
We hope that this e-Book will be useful to experimentalists and method developers and that it will stimulate further use and development of hybrid biomolecular modeling methods. We thank all authors, co-authors, and reviewers for their contribution to this Research Topic and acknowledge the support from Frontiers Team members.
Statements
Author contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
BirmannsS.RusuM.WriggersW. (2011). Using Sculptor and Situs for simultaneous assembly of atomic components into low-resolution shapes. J. Struct. Biol.173, 428–435. 10.1016/j.jsb.2010.11.002
2
BrownA.AmuntsA.BaiX. C.SugimotoY.EdwardsP. C.MurshudovG.et al. (2014). Structure of the large ribosomal subunit from human mitochondria. Science346, 718–722. 10.1126/science.1258026
3
CiferriC.LanderG. C.MaiolicaA.HerzogF.AebersoldR.NogalesE. (2012). Molecular architecture of human polycomb repressive complex 2. Elife1:e00005. 10.7554/eLife.00005
4
CottevieilleM.LarquetE.JonicS.PetoukhovM. V.CapriniG.ParavisiS.et al. (2008). The subnanometer resolution structure of the glutamate synthase 1.2-MDa hexamer by cryoelectron microscopy and its oligomerization behavior in solution: functional implications. J. Biol. Chem.283, 8237–8249. 10.1074/jbc.M708529200
5
GorbaC.TamaF. (2010). Normal mode flexible fitting of high-resolution structures of biological molecules toward SAXS data. Bioinform. Biol. Insights4, 43–54. 10.4137/BBI.S4551
6
JinQ.SorzanoC. O.de la Rosa-TrevínJ. M.Bilbao-CastroJ. R.Núñez-RamírezR.LlorcaO.et al. (2014). Iterative elastic 3D-to-2D alignment method using normal modes for studying structural dynamics of large macromolecular complexes. Structure22, 496–506. 10.1016/j.str.2014.01.004
7
KawabataT. (2008). Multiple subunit fitting into a low-resolution density map of a macromolecular complex using a gaussian mixture model. Biophys. J.95, 4643–4658. 10.1529/biophysj.108.137125
8
SaliA.BermanH. M.SchwedeT.TrewhellaJ.KleywegtG.BurleyS. K.et al. (2015). Outcome of the first wwPDB hybrid/integrative methods task force workshop. Structure23, 1156–1167. 10.1016/j.str.2015.05.013
9
TjioeE.LaskerK.WebbB.WolfsonH. J.SaliA. (2011). MultiFit: a web server for fitting multiple protein structures into their electron microscopy density map. Nucleic Acids Res.39, W167–W170. 10.1093/nar/gkr490
10
TokuhisaA.JonicS.TamaF.MiyashitaO. (2016). Hybrid approach for structural modeling of biological systems from X-ray free electron laser diffraction patterns. J. Struct. Biol.194, 325–336. 10.1016/j.jsb.2016.03.009
11
TrabucoL. G.VillaE.MitraK.FrankJ.SchultenK. (2008). Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics. Structure16, 673–683. 10.1016/j.str.2008.03.005
12
YangZ.LaskerK.Schneidman-DuhovnyD.WebbB.HuangC. C.PettersenE. F.et al. (2012). UCSF Chimera, MODELLER, and IMP: an integrated modeling system. J. Struct. Biol.179, 269–278. 10.1016/j.jsb.2011.09.006
Summary
Keywords
structure, dynamics, Conformational variability, hybrid models, biomolecules
Citation
Jonic S, Miyashita O and Callebaut I (2018) Editorial: Hybrid Biomolecular Modeling. Front. Mol. Biosci. 5:98. doi: 10.3389/fmolb.2018.00098
Received
20 September 2018
Accepted
23 October 2018
Published
09 November 2018
Volume
5 - 2018
Edited by
Dieter W. Heermann, Universität Heidelberg, Germany
Reviewed by
Peter Virnau, Johannes Gutenberg-Universität Mainz, Germany
Updates
Copyright
© 2018 Jonic, Miyashita and Callebaut.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Slavica Jonic slavica.jonic@upmc.fr
This article was submitted to Biophysics, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.