Conformational Sampling of a Biomolecular Rugged Energy Landscape

Author(s):  
Jakub Rydzewski ◽  
Rafal Jakubowski ◽  
Giuseppe Nicosia ◽  
Wieslaw Nowak
2016 ◽  
Vol 473 (12) ◽  
pp. 1651-1662 ◽  
Author(s):  
Shinji Iida ◽  
Haruki Nakamura ◽  
Junichi Higo

We introduce various, recently developed, generalized ensemble methods, which are useful to sample various molecular configurations emerging in the process of protein–protein or protein–ligand binding. The methods introduced here are those that have been or will be applied to biomolecular binding, where the biomolecules are treated as flexible molecules expressed by an all-atom model in an explicit solvent. Sampling produces an ensemble of conformations (snapshots) that are thermodynamically probable at room temperature. Then, projection of those conformations to an abstract low-dimensional space generates a free-energy landscape. As an example, we show a landscape of homo-dimer formation of an endothelin-1-like molecule computed using a generalized ensemble method. The lowest free-energy cluster at room temperature coincided precisely with the experimentally determined complex structure. Two minor clusters were also found in the landscape, which were largely different from the native complex form. Although those clusters were isolated at room temperature, with rising temperature a pathway emerged linking the lowest and second-lowest free-energy clusters, and a further temperature increment connected all the clusters. This exemplifies that the generalized ensemble method is a powerful tool for computing the free-energy landscape, by which one can discuss the thermodynamic stability of clusters and the temperature dependence of the cluster networks.


2016 ◽  
Vol 37 (31) ◽  
pp. 2687-2700 ◽  
Author(s):  
Shinji Iida ◽  
Tadaaki Mashimo ◽  
Takashi Kurosawa ◽  
Hironobu Hojo ◽  
Hiroya Muta ◽  
...  

2020 ◽  
Vol 48 (4) ◽  
pp. 1707-1724
Author(s):  
Jane R. Allison

Proteins are dynamic molecules that can transition between a potentially wide range of structures comprising their conformational ensemble. The nature of these conformations and their relative probabilities are described by a high-dimensional free energy landscape. While computer simulation techniques such as molecular dynamics simulations allow characterisation of the metastable conformational states and the transitions between them, and thus free energy landscapes, to be characterised, the barriers between states can be high, precluding efficient sampling without substantial computational resources. Over the past decades, a dizzying array of methods have emerged for enhancing conformational sampling, and for projecting the free energy landscape onto a reduced set of dimensions that allow conformational states to be distinguished, known as collective variables (CVs), along which sampling may be directed. Here, a brief description of what biomolecular simulation entails is followed by a more detailed exposition of the nature of CVs and methods for determining these, and, lastly, an overview of the myriad different approaches for enhancing conformational sampling, most of which rely upon CVs, including new advances in both CV determination and conformational sampling due to machine learning.


2021 ◽  
Author(s):  
Jose A Villegas ◽  
Tasneem M Vaid ◽  
Michael E Johnson ◽  
Terry W Moore

One of the principal difficulties in computational modeling of macromolecules is the vast conformational space that arises out of large numbers of atomic degrees of freedom. This problem is a familiar issue in the area of protein-protein docking, where models of protein complexes are generated from the monomeric subunits. Although restriction of molecular flexibility is a commonly used approximation that decreases the dimensionality of the problem, the seemingly endless number of possible ways two binding partners can interact generally necessitates the use of further approximations to explore the search space. Recently, growing interest in using computational tools to build predictive models of PROTAC-mediated complexes has led to the application of state-of-the-art protein-protein docking techniques to tackle this problem. Additionally, the atomic degrees of freedom introduced by flexibility of linkers used in the construction of PROTACs further expands the configurational search space, a problem that can be tackled with conformational sampling tools. However, repurposing existing tools to carry out protein-protein docking and linker conformer generation independently results in extensive sampling of structures incompatible with PROTAC-mediated complex formation. Here we show that it is possible to restrict the search to the space of protein-protein conformations that can be bridged by a PROTAC molecule with a given linker composition by using a cyclic coordinate descent algorithm to position PROTACs into complex-bound configurations. We use this methodology to construct a picture of the energy landscape of PROTAC-mediated interactions in a model test case, and show that the global minimum lies in the space of native-like conformations.


2004 ◽  
Vol 15 (07) ◽  
pp. 933-937 ◽  
Author(s):  
HANDAN ARKIN

A combination of replica exchange Monte Carlo sampling techniques and energy landscape paving approach is presented. This hybrid algorithm combines the features of the energy landscape paving (ELP) and replica exchange methods (REM). I have tested its performance in studying the low-energy conformations of the benchmark peptide Met-enkephalin.


2016 ◽  
Vol 18 (42) ◽  
pp. 29170-29182 ◽  
Author(s):  
Qiang Shao

A novel in silico approach (NMA–ITS) is introduced to rapidly and effectively sample the configuration space and give quantitative data for exploring the conformational changes of proteins.


2004 ◽  
Vol 71 ◽  
pp. 1-14
Author(s):  
David Leys ◽  
Jaswir Basran ◽  
François Talfournier ◽  
Kamaldeep K. Chohan ◽  
Andrew W. Munro ◽  
...  

TMADH (trimethylamine dehydrogenase) is a complex iron-sulphur flavoprotein that forms a soluble electron-transfer complex with ETF (electron-transferring flavoprotein). The mechanism of electron transfer between TMADH and ETF has been studied using stopped-flow kinetic and mutagenesis methods, and more recently by X-ray crystallography. Potentiometric methods have also been used to identify key residues involved in the stabilization of the flavin radical semiquinone species in ETF. These studies have demonstrated a key role for 'conformational sampling' in the electron-transfer complex, facilitated by two-site contact of ETF with TMADH. Exploration of three-dimensional space in the complex allows the FAD of ETF to find conformations compatible with enhanced electronic coupling with the 4Fe-4S centre of TMADH. This mechanism of electron transfer provides for a more robust and accessible design principle for interprotein electron transfer compared with simpler models that invoke the collision of redox partners followed by electron transfer. The structure of the TMADH-ETF complex confirms the role of key residues in electron transfer and molecular assembly, originally suggested from detailed kinetic studies in wild-type and mutant complexes, and from molecular modelling.


1997 ◽  
Vol 7 (3) ◽  
pp. 395-421 ◽  
Author(s):  
Jin Wang ◽  
Steven S. Plotkin ◽  
Peter G. Wolynes
Keyword(s):  

2020 ◽  
Author(s):  
Pia Vervoorts ◽  
Stefan Burger ◽  
Karina Hemmer ◽  
Gregor Kieslich

The zeolitic imidazolate frameworks ZIF-8 and ZIF-67 harbour a series of fascinating stimuli responsive properties. Looking at their responsitivity to hydrostatic pressure as stimulus, open questions exist regarding the isotropic compression with non-penetrating pressure transmitting media. By applying a state-of-the-art high-pressure powder X-ray diffraction setup, we revisit the high-pressure behaviour of ZIF-8 and ZIF-67 up to <i>p</i> = 0.4 GPa in small pressure increments. We observe a drastic, reversible change of high-pressure powder X-ray diffraction data at <i>p</i> = 0.3 GPa, discovering large volume structural flexibility in ZIF-8 and ZIF-67. Our results imply a shallow underlying energy landscape in ZIF-8 and ZIF-67, an observation that might point at rich polymorphism of ZIF-8 and ZIF-67, similar to ZIF-4(Zn).<br>


Sign in / Sign up

Export Citation Format

Share Document