scholarly journals Refining amino acid hydrophobicity for dynamics simulation of membrane proteins

PeerJ ◽  
2018 ◽  
Vol 6 ◽  
pp. e4230 ◽  
Author(s):  
Ronald D. Hills, Jr

Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems.

2017 ◽  
Author(s):  
Ronald D Hills, Jr

Coarse-grained simulations enable the study of membrane proteins in the context of their native environment but require reliable parameters. The CgProt force field is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. The reassignment of polar sidechain sites was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlate with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These data points can be used to further discriminate between alternate force field parameters. Available partitioning data was also used to reparameterize the representation of the polar peptide backbone for non-alanine residues. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in molecular dynamics simulations of a variety of membrane protein systems.


2017 ◽  
Author(s):  
Ronald D Hills, Jr

Coarse-grained simulations enable the study of membrane proteins in the context of their native environment but require reliable parameters. The CgProt force field is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. The reassignment of polar sidechain sites was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlate with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These data points can be used to further discriminate between alternate force field parameters. Available partitioning data was also used to reparameterize the representation of the polar peptide backbone for non-alanine residues. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in molecular dynamics simulations of a variety of membrane protein systems.


2019 ◽  
Vol 123 (37) ◽  
pp. 7829-7839
Author(s):  
Karolina Ziȩba ◽  
Magdalena Ślusarz ◽  
Rafał Ślusarz ◽  
Adam Liwo ◽  
Cezary Czaplewski ◽  
...  

2021 ◽  
Author(s):  
Christos Lamprakis ◽  
Ioannis Andreadelis ◽  
John Manchester ◽  
Camilo Velez-Vega ◽  
José S. Duca ◽  
...  

<p>Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3 we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.</p>


2020 ◽  
Vol 22 (41) ◽  
pp. 24028-24040
Author(s):  
Maryam S. Sadeghi ◽  
Mohammad Reza Moghbeli ◽  
William A. Goddard

Application of the novel quantum based coarse grained force field (CGq FF) for formation of a micelle from 250 chains of 2000 Dalton CG-MePEG23-b-PCL9 block copolymer in water at 310.15 K.


Sign in / Sign up

Export Citation Format

Share Document