Investigation of the conformational space of hydrophobic-polar heteropolymers by gyration tensor based parameters

2022 ◽  
Vol 552 ◽  
pp. 111372
Author(s):  
Buket Taşdizen Özyurt ◽  
Handan Arkın
Keyword(s):  
2020 ◽  
Author(s):  
Abhishek Singh ◽  
Reman K. Singh ◽  
G Naresh Patwari

The rational design of conformationally controlled foldable modules can lead to a deeper insight into the conformational space of complex biological molecules where non-covalent interactions such as hydrogen bonding and π-stacking are known to play a pivotal role. Squaramides are known to have excellent hydrogen bonding capabilities and hence, are ideal molecules for designing foldable modules that can mimic the secondary structures of bio-molecules. The π-stacking induced folding of bis-squaraines tethered using aliphatic primary and secondary-diamine linkers of varying length is explored with a simple strategy of invoking small perturbations involving the length linkers and degree of substitution. Solution phase NMR investigations in combination with molecular dynamics simulations suggest that bis-squaraines predominantly exist as extended conformations. Structures elucidated by X-ray crystallography confirmed a variety of folded and extended secondary conformations including hairpin turns and 𝛽-sheets which are determined by the hierarchy of π-stacking relative to N–H···O hydrogen bonds.


2020 ◽  
Author(s):  
Abhishek Singh ◽  
Reman K. Singh ◽  
G Naresh Patwari

The rational design of conformationally controlled foldable modules can lead to a deeper insight into the conformational space of complex biological molecules where non-covalent interactions such as hydrogen bonding and π-stacking are known to play a pivotal role. Squaramides are known to have excellent hydrogen bonding capabilities and hence, are ideal molecules for designing foldable modules that can mimic the secondary structures of bio-molecules. The π-stacking induced folding of bis-squaraines tethered using aliphatic primary and secondary-diamine linkers of varying length is explored with a simple strategy of invoking small perturbations involving the length linkers and degree of substitution. Solution phase NMR investigations in combination with molecular dynamics simulations suggest that bis-squaraines predominantly exist as extended conformations. Structures elucidated by X-ray crystallography confirmed a variety of folded and extended secondary conformations including hairpin turns and 𝛽-sheets which are determined by the hierarchy of π-stacking relative to N–H···O hydrogen bonds.


2020 ◽  
Author(s):  
Lim Heo ◽  
Collin Arbour ◽  
Michael Feig

Protein structures provide valuable information for understanding biological processes. Protein structures can be determined by experimental methods such as X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, or cryogenic electron microscopy. As an alternative, in silico methods can be used to predict protein structures. Those methods utilize protein structure databases for structure prediction via template-based modeling or for training machine-learning models to generate predictions. Structure prediction for proteins distant from proteins with known structures often results in lower accuracy with respect to the true physiological structures. Physics-based protein model refinement methods can be applied to improve model accuracy in the predicted models. Refinement methods rely on conformational sampling around the predicted structures, and if structures closer to the native states are sampled, improvements in the model quality become possible. Molecular dynamics simulations have been especially successful for improving model qualities but although consistent refinement can be achieved, the improvements in model qualities are still moderate. To extend the refinement performance of a simulation-based protocol, we explored new schemes that focus on an optimized use of biasing functions and the application of increased simulation temperatures. In addition, we tested the use of alternative initial models so that the simulations can explore conformational space more broadly. Based on the insight of this analysis we are proposing a new refinement protocol that significantly outperformed previous state-of-the-art molecular dynamics simulation-based protocols in the benchmark tests described here. <br>


2019 ◽  
Author(s):  
Chem Int

The full conformational space of N-formyl-L-alanine-amide was explored by the semi-empirical method AM1 coupled to the Multi Niche Crowding (MNC) genetic algorithm implemented in a package of programs developed in our laboratory. The structural and energy analysis of the resulting conformational space E(,ψ) exhibits 5 regions or minima ɣL, ɣD, ɛL, D and αD. The technique provides better detection of local and global minima within a reasonable time.


2020 ◽  
Vol 27 (4) ◽  
pp. 321-328 ◽  
Author(s):  
Yanru Li ◽  
Ying Zhang ◽  
Jun Lv

Background: Protein folding rate is mainly determined by the size of the conformational space to search, which in turn is dictated by factors such as size, structure and amino-acid sequence in a protein. It is important to integrate these factors effectively to form a more precisely description of conformation space. But there is no general paradigm to answer this question except some intuitions and empirical rules. Therefore, at the present stage, predictions of the folding rate can be improved through finding new factors, and some insights are given to the above question. Objective: Its purpose is to propose a new parameter that can describe the size of the conformational space to improve the prediction accuracy of protein folding rate. Method: Based on the optimal set of amino acids in a protein, an effective cumulative backbone torsion angles (CBTAeff) was proposed to describe the size of the conformational space. Linear regression model was used to predict protein folding rate with CBTAeff as a parameter. The degree of correlation was described by the coefficient of determination and the mean absolute error MAE between the predicted folding rates and experimental observations. Results: It achieved a high correlation (with the coefficient of determination of 0.70 and MAE of 1.88) between the logarithm of folding rates and the (CBTAeff)0.5 with experimental over 112 twoand multi-state folding proteins. Conclusion: The remarkable performance of our simplistic model demonstrates that CBTA based on optimal set was the major determinants of the conformation space of natural proteins.


2020 ◽  
Vol 22 (9) ◽  
pp. 635-648 ◽  
Author(s):  
Korosh Mashayekh ◽  
Shahrzad Sharifi ◽  
Tahereh Damghani ◽  
Maryam Elyasi ◽  
Mohammad S. Avestan ◽  
...  

Background: c-Met kinase plays a critical role in a myriad of human cancers, and a massive scientific work was devoted to design more potent inhibitors. Objective: In this study, 16 molecular dynamics simulations of different complexes of potent c-Met inhibitors with U-shaped binding mode were carried out regarding the dynamic ensembles to design novel potent inhibitors. Methods: A cluster analysis was performed, and the most representative frame of each complex was subjected to the structure-based pharmacophore screening. The GOLD docking program investigated the interaction energy and pattern of output hits from the virtual screening. The most promising hits with the highest scoring values that showed critical interactions with c-Met were presented for ADME/Tox analysis. Results: The screening yielded 45,324 hits that all of them were subjected to the docking studies and 10 of them with the highest-scoring values having diverse structures were presented for ADME/Tox analyses. Conclusion: The results indicated that all the hits shared critical Pi-Pi stacked and hydrogen bond interactions with Tyr1230 and Met1160 respectively.


Polymers ◽  
2020 ◽  
Vol 13 (1) ◽  
pp. 99
Author(s):  
Cristian Privat ◽  
Sergio Madurga ◽  
Francesc Mas ◽  
Jaime Rubio-Martínez

Solvent pH is an important property that defines the protonation state of the amino acids and, therefore, modulates the interactions and the conformational space of the biochemical systems. Generally, this thermodynamic variable is poorly considered in Molecular Dynamics (MD) simulations. Fortunately, this lack has been overcome by means of the Constant pH Molecular Dynamics (CPHMD) methods in the recent decades. Several studies have reported promising results from these approaches that include pH in simulations but focus on the prediction of the effective pKa of the amino acids. In this work, we want to shed some light on the CPHMD method and its implementation in the AMBER suitcase from a conformational point of view. To achieve this goal, we performed CPHMD and conventional MD (CMD) simulations of six protonatable amino acids in a blocked tripeptide structure to compare the conformational sampling and energy distributions of both methods. The results reveal strengths and weaknesses of the CPHMD method in the implementation of AMBER18 version. The change of the protonation state according to the chemical environment is presumably an improvement in the accuracy of the simulations. However, the simulations of the deprotonated forms are not consistent, which is related to an inaccurate assignment of the partial charges of the backbone atoms in the CPHMD residues. Therefore, we recommend the CPHMD methods of AMBER program but pointing out the need to compare structural properties with experimental data to bring reliability to the conformational sampling of the simulations.


2021 ◽  
Vol 22 (3) ◽  
pp. 1364
Author(s):  
V. V. Krishnan ◽  
Timothy Bentley ◽  
Alina Xiong ◽  
Kalyani Maitra

Both nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations are routinely used in understanding the conformational space sampled by peptides in the solution state. To investigate the role of single-residue change in the ensemble of conformations sampled by a set of heptapeptides, AEVXEVG with X = L, F, A, or G, comprehensive NMR, and MD simulations were performed. The rationale for selecting the particular model peptides is based on the high variability in the occurrence of tri-peptide E*L between the transmembrane β-barrel (TMB) than in globular proteins. The ensemble of conformations sampled by E*L was compared between the three sets of ensembles derived from NMR spectroscopy, MD simulations with explicit solvent, and the random coil conformations. In addition to the estimation of global determinants such as the radius of gyration of a large sample of structures, the ensembles were analyzed using principal component analysis (PCA). In general, the results suggest that the -EVL- peptide indeed adopts a conformational preference that is distinctly different not only from a random distribution but also from other peptides studied here. The relatively straightforward approach presented herein could help understand the conformational preferences of small peptides in the solution state.


2021 ◽  
Author(s):  
Théo Jaffrelot Inizan ◽  
Frédéric Célerse ◽  
Olivier Adjoua ◽  
Dina El Ahdab ◽  
Luc-Henri Jolly ◽  
...  

We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs).


Sign in / Sign up

Export Citation Format

Share Document