Dataset Information

Unfavorable regions in the ramachandran plot: Is it really steric hindrance? The interacting quantum atoms perspective.

ABSTRACT: Accurate description of the intrinsic preferences of amino acids is important to consider when developing a biomolecular force field. In this study, we use a modern energy partitioning approach called Interacting Quantum Atoms to inspect the cause of the ? and ? torsional preferences of three dipeptides (Gly, Val, and Ile). Repeating energy trends at each of the molecular, functional group, and atomic levels are observed across both (1) the three amino acids and (2) the ?/? scans in Ramachandran plots. At the molecular level, it is surprisingly electrostatic destabilization that causes the high-energy regions in the Ramachandran plot, not molecular steric hindrance (related to the intra-atomic energy). At the functional group and atomic levels, the importance of key peptide atoms (O_i_-1 , C_i , N_i , N_i₊₁ ) and some sidechain hydrogen atoms (H_? ) are identified as responsible for the destabilization seen in the energetically disfavored Ramachandran regions. Consistently, the O_i_-1 atoms are particularly important for the explanation of dipeptide intrinsic behavior, where electrostatic and steric destabilization unusually complement one another. The findings suggest that, at least for these dipeptides, it is the peptide group atoms that dominate the intrinsic behavior, more so than the sidechain atoms. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

SUBMITTER: Maxwell PI

PROVIDER: S-EPMC5659141 | biostudies-literature | 2017 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Unfavorable regions in the ramachandran plot: Is it really steric hindrance? The interacting quantum atoms perspective.

Maxwell Peter I PI Popelier Paul L A PLA

Journal of computational chemistry 20170825 29

Accurate description of the intrinsic preferences of amino acids is important to consider when developing a biomolecular force field. In this study, we use a modern energy partitioning approach called Interacting Quantum Atoms to inspect the cause of the φ and ψ torsional preferences of three dipeptides (Gly, Val, and Ile). Repeating energy trends at each of the molecular, functional group, and atomic levels are observed across both (1) the three amino acids and (2) the φ/ψ scans in Ramachandran ...[more]

PMID: 28841241

Similar Datasets

Project description:Energy decomposition analysis (EDA) is a well-established approach to dissect the interaction energy into chemically sound components. Despite the inherent requirement of reference states has been a long-standing object of debate, the direct relation with the molecular orbital analysis helps in building up predictive models. The alternative molecular energy decomposition schemes that decompose the total energy into atomic and diatomic contributions, such as the interacting quantum atoms (IQA), has no external reference requirements and also the intra- and intermolecular interactions are treated on equal footing. However, a connection with heuristic chemical models are limited, bringing about a somewhat narrower predictive power. While efforts to reconcile the bonding picture obtained by both methodologies have been discussed in the past, a synergic combination of them has not been tackled yet. Herein, we present the use of IQA decomposition of the individual terms arising from the EDA in the context of intermolecular interactions, henceforth EDA-IQA. The method is applied to a molecular set covering a wide range of interaction types, including hydrogen bonding, charge-dipole, π-π and halogen interactions. We find that the electrostatic energy from EDA, entirely seen as intermolecular, leads to meaningful and non-negligible intra-fragment contributions upon IQA decomposition, originated from charge penetration. EDA-IQA also affords the decomposition of the Pauli repulsion term into intra- and inter-fragment contributions. The intra-fragment term is destabilizing, particularly for the moieties that are net acceptors of charge, while the inter-fragment Pauli term is actually stabilizing. In the case of the orbital interaction term, the sign and magnitude of the intra-fragment contribution at equilibrium geometries is largely driven by the amount of charge transfer, while the inter-fragment contribution is clearly stabilizing. EDA-IQA terms show a smooth behavior along the intermolecular dissociation path of selected systems. The new EDA-IQA methodology provides a richer energy decomposition scheme that aims at bridging the gap between the two main distinct real-space and Hilbert-space methodologies. Via this approach, the partitioning can be used directionally on all the EDA terms aiding in identifying the causal effects on geometries and/or reactivity.

Dataset Information

Unfavorable regions in the ramachandran plot: Is it really steric hindrance? The interacting quantum atoms perspective.

Publications

Unfavorable regions in the ramachandran plot: Is it really steric hindrance? The interacting quantum atoms perspective.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets