Loss of a water-mediated network results in reduced agonist affinity in a β2-adrenergic receptor clinical variant

State-dependent dynamics of extramembrane domains in the β2 -adrenergic receptor

G protein-coupled receptors (GPCRs) are membrane-bound signaling proteins that play an essential role in cellular signaling processes. Due to their intrinsic function of transmitting internal signals in response to external cues, these receptors are adapted to be highly dynamic in nature. The β2 -adrenergic receptor (β2 AR) is a representative member of the family that has been extensively analyzed in terms of its structure and activation. Although the structure of the transmembrane domain has been characterized in the different functional states of the receptor, the conformational dynamics of the extramembrane domains, especially the intrinsically disordered regions are still emerging. In this study, we analyze the state-dependent dynamics of extramembrane domains of β2 AR using atomistic molecular dynamics simulations. We introduce a parameter, the residue excess dynamics that allows us to better quantify receptor dynamics. Using this measure, we show that the dynamics of the extramembrane domains are sensitive to the receptor state. Interestingly, the ligand-bound intermediateR ' state shows the maximal dynamics compared to either the active R*G or inactive R states. Ligand binding appears to be correlated with high residue excess dynamics that are dampened upon G protein coupling. The intracellular loop-3 (ICL3) domain has a tendency to flip towards the membrane upon ligand binding, which could contribute to receptor "priming." We highlight an important ICL1-helix-8 interplay that is broken in the ligand-bound state but is retained in the active state. Overall, our study highlights the importance of characterizing the functional dynamics of the GPCR loop domains.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

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.

Assessing Protein-Protein Docking Protocols: Case Studies of G-Protein-Coupled Receptor Interactions

The interactions of G-protein-coupled receptors (GPCRs) with other proteins are critical in several cellular processes but resolving their structural dynamics remains challenging. An increasing number of GPCR complexes have been experimentally resolved but others including receptor variants are yet to be characterized, necessitating computational predictions of their interactions. Although integrative approaches with multi-scale simulations would provide rigorous estimates of their conformational dynamics, protein-protein docking remains a first tool of choice of many researchers due to the availability of open-source programs and easy to use web servers with reasonable predictive power. Protein-protein docking algorithms have limited ability to consider protein flexibility, environment effects, and entropy contributions and are usually a first step towards more integrative approaches. The two critical steps of docking: the sampling and scoring algorithms have improved considerably and their performance has been validated against experimental data. In this chapter, we provide an overview and generalized protocol of a few docking protocols using GPCRs as test cases. In particular, we demonstrate the interactions of GPCRs with extracellular protein ligands and an intracellular protein effectors (G-protein) predicted from docking approaches and test their limitations. The current chapter will help researchers critically assess docking protocols and predict experimentally testable structures of GPCR complexes.

Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment

The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

Molecular determinants of GPCR pharmacogenetics: Deconstructing the population variants in β2-adrenergic receptor

G protein-coupled receptors (GPCRs) are membrane proteins that play a central role in cell signaling and constitute one of the largest classes of drug targets. The molecular mechanisms underlying GPCR function have been characterized by several experimental and computational methods and provide an understanding of their role in physiology and disease. Population variants arising from nsSNPs affect the native function of GPCRs and have been implicated in differential drug response. In this chapter, we provide an overview on GPCR structure and activation, with a special focus on the β₂-adrenergic receptor (β₂-AR). First, we discuss the current understanding of the structural and dynamic features of the wildtype receptor. Subsequently, the population variants identified in this receptor from clinical and large-scale genomic studies are described. We show how computational approaches such as bioinformatics tools and molecular dynamics simulations can be used to characterize the variant receptors in comparison to the wildtype receptor. In particular, we discuss three examples of clinically important variants and discuss how the structure and function of these variants differ from the wildtype receptor at a molecular level. Overall, the chapter provides an overview of structure and function of GPCR variants and is a step towards the study of inter-individual differences and personalized medicine.

Hydration Dynamics in Biological Membranes: Emerging Applications of Terahertz Spectroscopy

Water drives the spontaneous self-assembly of lipids and proteins into quasi two-dimensional biological membranes that act as catalytic scaffolds for numerous processes central to life. However, the functional relevance of hydration in membrane biology is only beginning to be addressed, predominantly because of challenges associated with direct measurements of hydration microstructure and dynamics in a biological milieu. Our recent work on the novel interplay of membrane electrostatics and crowding in shaping membrane hydration dynamics utilizing terahertz (THz) spectroscopy represents an important step in this context. In this Perspective, we provide a glimpse into the ever-broadening functional landscape of hydration dynamics in biological membranes in the backdrop of the unique physical chemistry of water molecules. We further highlight the immense (and largely untapped) potential of the THz toolbox in addressing contemporary problems in membrane biology, while emphasizing the adaptability of the analytical framework reported recently by us to such studies.