Quantifying conformational changes in GPCRs : glimpse of a common functional mechanism

Background: G-protein-coupled receptors (GPCRs) are important drug targets and a better understanding of their molecular mechanisms would be desirable. The crystallization rate of GPCRs has accelerated in recent years as techniques have become more sophisticated, particularly with respect to Class A GPCRs interacting with G-proteins. These developments have made it possible for a quantitative analysis of GPCR geometrical features and binding-site conformations, including a statistical comparison between Class A GPCRs in active (agonist-bound) and inactive (antagonist-bound) states. Results: Here we implement algorithms for the analysis of interhelical angles, distances, interactions and binding-site volumes in the transmembrane domains of 25 Class A GPCRs (7 active and 18 inactive). Two interhelical angles change in a statistically significant way between average inactive and active states: TM3-TM6 (by -9°) and TM6-TM7 (by +12°). A third interhelical angle: TM5-TM6 shows a trend, changing by -9°. In the transition from inactive to active states, average van der Waals interactions between TM3 and TM7 significantly increase as the average distance between them decreases by >2 Å. Average H-bonding between TM3 and TM6 decreases but is seemingly compensated by an increase in H-bonding between TM5 and TM6. In five Class A GPCRs, crystallized in both active and inactive states, increased H-bonding of agonists to TM6 and TM7, relative to antagonists, is observed. These protein-agonist interactions likely favour a change in the TM6-TM7 angle, which creates a narrowing in the binding pocket of activated receptors and an average ~200 Å3 reduction in volume. Conclusions: In terms of similar conformational changes and agonist binding pattern, Class A GPCRs appear to share a common mechanism of activation, which can be exploited in future drug development

Statistics for the analysis of molecular dynamics simulations : providing P values for agonist-dependent GPCR activation

Molecular dynamics (MD) is the common computational technique for assessing efficacy of GPCR-bound ligands. Agonist efficacy measures the capability of the ligand-bound receptor of reaching the active state in comparison with the free receptor. In this respect, agonists, neutral antagonists and inverse agonists can be considered. A collection of MD simulations of both the ligand-bound and the free receptor are needed to provide reliable conclusions. Variability in the trajectories needs quantification and proper statistical tools for meaningful and non-subjective conclusions. Multiple-factor (time, ligand, lipid) ANOVA with repeated measurements on the time factor is proposed as a suitable statistical method for the analysis of agonist-dependent GPCR activation MD simulations. Inclusion of time factor in the ANOVA model is consistent with the time-dependent nature of MD. Ligand and lipid factors measure agonist and lipid influence on receptor activation. Previously reported MD simulations of adenosine A2a receptor (A2aR) are reanalyzed with this statistical method. TM6-TM3 and TM7-TM3 distances are selected as dependent variables in the ANOVA model. The ligand factor includes the presence or absence of adenosine whereas the lipid factor considers DOPC or DOPG lipids. Statistical analysis of MD simulations shows the efficacy of adenosine and the effect of the membrane lipid composition. Subsequent application of the statistical methodology to NECA A2aR agonist, with resulting P values in consistency with its pharmacological profile, suggests that the method is useful for ligand comparison and potentially for dynamic structure-based virtual screening