Search

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

The chemokine receptor CXCR4 and its cognate ligand CXCL12 mediate pathways that lead to cell migration and chemotaxis. Although the structural details of related receptor-ligand complexes have been resolved, the roles of the N-terminal domain of the receptor and post-translational sulfation that are determinants of ligand selectivity and affinity remain unclear. Here, we analyze the structural dynamics induced by receptor sulfation by combining molecular dynamics, docking and network analysis. The sulfotyrosine residues, 7Ys N-term , 12Ys N-term and 21Ys N-term allow the N-terminal domain of the apo-sulfated receptor to adopt an "open" conformation that appears to facilitate ligand binding. The overall topology of the CXCR4-CXCL12 complex is independent of the sulfation state, but an extensive network of protein-protein interactions characterizes the sulfated receptor, in line with its increased ligand affinity. The altered interactions of sulfotyrosine residues, such as 21Ys N-term -47R CXCL12 replacing the 21Y N-term -13F CXCL12 interaction, propagate via allosteric pathways towards the receptor lumen. In particular, our results suggest that the experimentally-reported receptor-ligand interactions 262D 6.58 -8R CXCL12 and 277E 7.28 -12R CXCL12 could be dependent on the sulfation state of the receptor and need to be carefully analyzed. Our work is an important step in understanding chemokine-receptor interactions and how post-translational modifications could modulate receptor-ligand complexes., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Dynamic interactions between membrane proteins span a range of spatio-temporal scales and determine several cellular outcomes. Experimental methods in structure determination are able to resolve static protein-protein complexes at the membrane, but lack the resolution required for disordered, flexible domains and dynamic interactions. Computational approaches could bridge the resolution gap and help to unravel molecular details underlying these crucial interactions. Here, we review current approaches to predict dynamic membrane-protein complexes, with a focus on G protein-coupled receptors (GPCRs). Ensemble coarse-grain simulations have captured the conformational heterogeneity of several membrane receptor complexes. In conjunction, the conformational plasticity of protein interfaces especially encompassing unstructured domains is well represented by atomistic simulations. A combined integrative approach will pave the way forward to understand the molecular details of these dynamic complexes., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

The CXC chemokine receptor 1 (CXCR1) is an important member of the G protein-coupled receptor (GPCR) family in which the extracellular N-terminal domain has been implicated in ligand binding and selectivity. The structure of this domain has not yet been elucidated due to its inherent dynamics, but experimental evidence points toward membrane-dependent organization and dynamics. To gain molecular insight into the interaction of the N-terminal domain with the membrane bilayer, we performed a series of microsecond time scale atomistic simulations of the N-terminal domain of CXCR1 in the presence and absence of POPC bilayers. Our results show that the peptide displays a high propensity to adopt a β-sheet conformation in the presence of the membrane bilayer. The interaction of the peptide with the membrane bilayer was found to be transient in our simulations. Interestingly, a scrambled peptide, containing the same residues in a randomly varying sequence, did not exhibit membrane-modulated structural dynamics. These results suggest that sequence-dependent electrostatics, modulated by the membrane, could play an important role in folding of the N-terminal domain. We believe that our results reinforce the emerging paradigm that cellular membranes could be important modulators of function of G protein-coupled receptors such as CXCR1., (Copyright © 2017 Elsevier B.V. All rights reserved.)