A theoretical solvation model of peptides and proteins that mimics the heterogeneous membrane-water system was proposed. Our approach is based on the combined use of atomic parameters of solvation for water and hydrocarbons, which approximates the hydrated polar groups and acyl chains of lipids, respectively. This model was tested in simulations of several peptides: a nonpolar 20-mer polyleucine, a hydrophobic peptide with terminal polar groups, and a strongly amphiphilic peptide. The conformational space of the peptides in the presence of the membrane was studied by the Monte Carlo method. Unlike a polar solvent and vacuum, the membrane-like environment was shown to stabilize the alpha-helical conformation: low-energy structures have a helicity index of 100% in all cases. At the same time, the energetically most favorable orientations of the peptides relative to the membrane depend on their hydrophobic properties: nonpolar polyleucine is entirely immersed in the bilayer and the hydrophobic peptide with polar groups at the termini adopts a transbilayer orientation, whereas the amphiphilic peptide lies at the interface parallel to the membrane plane. The results of the simulations agree well with the available experimental data for these systems. In the following communications of this series, we plan to describe applications of the solvation model to membrane-bound proteins and peptides with biologically important functional activities.

Download full-text PDF

Source

Publication Analysis

Top Keywords

solvation model
12
polar groups
12
peptides proteins
8
hydrophobic peptide
8
amphiphilic peptide
8
peptides
5
[modeling peptides
4
membrane
4
proteins membrane
4
membrane environment
4

Similar Publications

Extending the MST Model to Large Biomolecular Systems: Parametrization of the ddCOSMO-MST Continuum Solvation Model.

J Comput Chem

January 2025

Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l'Alimentació, Universitat de Barcelona (UB), Barcelona, Spain.

Continuum solvation models such as the polarizable continuum model and the conductor-like screening model are widely used in quantum chemistry, but their application to large biosystems is hampered by their computational cost. Here, we report the parametrization of the Miertus-Scrocco-Tomasi (MST) model for the prediction of hydration free energies of neutral and ionic molecules based on the domain decomposition formulation of COSMO (ddCOSMO), which allows a drastic reduction of the computational cost by several orders of magnitude. We also introduce several novelties in MST, like a new definition of atom types based on hybridization and an automatic setup of the cavity for charged regions.

View Article and Find Full Text PDF

This study aims to shed light on the mechanism and kinetics of 1,4-dioxane degradation by hydroxyl radical (OH) across various solvation conditions to evaluate electronic and structural properties at the MP2/aug-cc-pVTZ level. Transition states (TS) structures determined in the gas phase and SMD solvation model reveal similar hydrogen abstraction patterns. In contrast, the explicit solvation model (ES) introduces significant changes, suggesting a kinetic preference for axial pathways.

View Article and Find Full Text PDF

The Search for the Optimal Methodology for Predicting Fluorinated Cathinone Drugs NMR Chemical Shifts.

Molecules

December 2024

Institute of Earth Sciences, Faculty of Natural Sciences, University of Silesia in Katowice, 60, Będzińska, 41-200 Sosnowiec, Poland.

Cathinone and its synthetic derivatives belong to organic compounds with narcotic properties. Their structural diversity and massive illegal use create the need to develop new analytical methods for their identification in different matrices. NMR spectroscopy is one of the most versatile methods for identifying the structure of organic substances.

View Article and Find Full Text PDF

The separation and structural identification of glycans are of great bioanalytical importance. To obtain a good understanding of the structural flexibility of glycans, replica exchange molecular dynamics (REMD) simulations were used based on AMBER force field calculations to create ensembles of glycan structures. Nonpolar surface area (NPSA) calculations based on continuum solvation (CS) models (Dhakal, R.

View Article and Find Full Text PDF

With the advancement of genetic code expansion, the field is progressing towards incorporating multiple non-canonical amino acids (ncAAs). The specificity of aminoacyl-tRNA synthetases (aaRSs) towards ncAAs is a critical factor, as engineered aaRSs frequently show polyspecificity, complicating the precise incorporation of multiple ncAAs. To address this challenge, predicting binding affinity can be beneficial.

View Article and Find Full Text PDF

Want AI Summaries of new PubMed Abstracts delivered to your In-box?

Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!