scholarly journals Ensemble completeness in conformer sampling: the case of small macrocycles

2021 ◽  
Vol 13 (1) ◽  
Author(s):  
Lea Seep ◽  
Anne Bonin ◽  
Katharina Meier ◽  
Holger Diedam ◽  
Andreas H. Göller

AbstractIn this study we compare the three algorithms for the generation of conformer ensembles Biovia BEST, Schrödinger Prime macrocycle sampling (PMM) and Conformator (CONF) form the University of Hamburg, with ensembles derived for exhaustive molecular dynamics simulations applied to a dataset of 7 small macrocycles in two charge states and three solvents. Ensemble completeness is a prerequisite to allow for the selection of relevant diverse conformers for many applications in computational chemistry. We apply conformation maps using principal component analysis based on ring torsions. Our major finding critical for all applications of conformer ensembles in any computational study is that maps derived from MD with explicit solvent are significantly distinct between macrocycles, charge states and solvents, whereas the maps for post-optimized conformers using implicit solvent models from all generator algorithms are very similar independent of the solvent. We apply three metrics for the quantification of the relative covered ensemble space, namely cluster overlap, variance statistics, and a novel metric, Mahalanobis distance, showing that post-optimized MD ensembles cover a significantly larger conformational space than the generator ensembles, with the ranking PMM > BEST >> CONF. Furthermore, we find that the distributions of 3D polar surface areas are very similar for all macrocycles independent of charge state and solvent, except for the smaller and more strained compound 7, and that there is also no obvious correlation between 3D PSA and intramolecular hydrogen bond count distributions.

2019 ◽  
Author(s):  
Peng He ◽  
Sheila Sarkar ◽  
Emilio Gallicchio ◽  
Tom Kurtzman ◽  
Lauren Wickstrom

<p>This study investigates the role of hydration and its relationship to the conformational equilibrium of the host molecule β-cyclodextrin. Molecular dynamics simulations indicate that the unbound β-cyclodextrin exhibits two state behavior in explicit solvent due to the opening and closing of its cavity. In implicit solvent, these transitions are not observed and there is one dominant conformation of β-cyclodextrin with an open cavity. Based on these observations, we investigate the hypothesis that the expulsion of thermodynamically unfavorable water molecules into the bulk plays an important role in controlling the accessibility of the closed macrostate at room temperature. We compare the results of the molecular mechanics analytical generalized Born plus non-polar solvation approach to those obtained through Grid Inhomogeneous Solvation Theory analysis with explicit solvation to elucidate the thermodynamic forces at play. The calculations help to illustrate the deficiencies of continuum solvent models and demonstrate the key role of the thermodynamics of enclosed hydration in driving the conformational equilibrium of molecules in solution. </p>


2019 ◽  
Vol 92 (2) ◽  
pp. 259-267
Author(s):  
Karlo Sović ◽  
Tea Ostojić ◽  
Ines Primožič ◽  
Tomica Hrenar ◽  
Mirjana Skočibušić ◽  
...  

Full conformational space of cinchonine and cinchonidine has been investigated by means of statistical analysis of quantum chemical molecular dynamics simulations. Recently developed procedure comprising principal component analysis of molecular dynamics trajectories was applied on cinchonine and cinchonidine as well as on their protonated and methylated quaternary derivatives. The method for full conformational analysis includes Cartesian coordinates sampling through quantum chemical molecular dynamics simulations, reduction of dimensionality by principal component analysis, determination of probability distributions in a reduced space of Cartesian coordinates and search for all the strict extrema points in probability distribution functions. In order to gain crucial insight in the understanding of chirality induction of these alkaloids, comparison of the determined conformational spaces of pseudo-enantiomers has been made. It was shown that protonation of the quinuclidine nitrogen atom stabilizes the conformers with the intramolecular 1N–H∙∙∙9O hydrogen bond whereas methylation on the same position results in the reduction of the domain of internal coordinates responsible for the conformational space.


2021 ◽  
Author(s):  
Eric Lang ◽  
Emily Baker ◽  
Derek Woolfson ◽  
Adrian Mulholland

We test a range of standard implicit solvent models and protein forcefields for a set of 5 experimentally characterized, designed α-helical peptides. 65 combinations of forcefield and implicit solvent models are evaluated in >800 µs of molecular dynamics simulations. The data show that implicit solvent models generally fail to reproduce the experimentally observed secondary structure content, and none performs well for all 5 peptides. The results show that these models are not usefully predictive.


2018 ◽  
Author(s):  
Maud Jusot ◽  
Dirk Stratmann ◽  
Marc Vaisset ◽  
Jacques Chomilier ◽  
Juan Cortés

Small cyclic peptides represent a promising class of therapeutic molecules with unique chemical properties. However, the poor knowledge of their structural characteristics makes their computational design and structure prediction a real challenge. In order to better describe their conformational space, we developed a method, named EGSCyP, for the exhaustive exploration of the energy landscape of small head-to-tail cyclic peptides. The method can be summarized by (i) a global exploration of the conformational space based on a mechanistic representation of the peptide and the use of robotics-based algorithms to deal with the closure constraint, (ii) an all-atom refinement of the obtained conformations. EGSCyP can handle D-form residues and N-methylations. Two strategies for the side-chains placement were implemented and compared. To validate our approach, we applied it to a set of three variants of cyclic RGDFV pentapeptides, including the drug candidate Cilengitide. A comparative analysis was made with respect to replica exchange molecular dynamics simulations in implicit solvent. It results that the EGSCyP method provides a very complete characterization of the conformational space of small cyclic pentapeptides.


2019 ◽  
Author(s):  
Peng He ◽  
Sheila Sarkar ◽  
Emilio Gallicchio ◽  
Tom Kurtzman ◽  
Lauren Wickstrom

<p>This study investigates the role of hydration and its relationship to the conformational equilibrium of the host molecule β-cyclodextrin. Molecular dynamics simulations indicate that the unbound β-cyclodextrin exhibits two state behavior in explicit solvent due to the opening and closing of its cavity. In implicit solvent, these transitions are not observed and there is one dominant conformation of β-cyclodextrin with an open cavity. Based on these observations, we investigate the hypothesis that the expulsion of thermodynamically unfavorable water molecules into the bulk plays an important role in controlling the accessibility of the closed macrostate at room temperature. We compare the results of the molecular mechanics analytical generalized Born plus non-polar solvation approach to those obtained through Grid Inhomogeneous Solvation Theory analysis with explicit solvation to elucidate the thermodynamic forces at play. The calculations help to illustrate the deficiencies of continuum solvent models and demonstrate the key role of the thermodynamics of enclosed hydration in driving the conformational equilibrium of molecules in solution. </p>


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.


Sign in / Sign up

Export Citation Format

Share Document