Determinants of Neutral Antagonism and Inverse Agonism in the ß2-Adrenergic Receptor.

Free-energy profiles for the activation/deactivation of the ß2-adrenergic receptor (ADRB2) with neutral antagonist and inverse agonist ligands have been determined with well-tempered multiple-walker (MW) metadynamics simulations. The inverse agonists carazolol and ICI118551 clearly favor single inactive conformational minima in both the binary and ternary ligand-receptor-G-protein complexes, in accord with the inverse-agonist activity of the ligands. The behavior of neutral antagonists is more complex, as they seem also to affect the recruitment of the G-protein. The results are analyzed in terms of the conformational states of the well-known microswitches that have been proposed as indicators of receptor activity.

Extended Metadynamics Protocol for Binding/Unbinding Free Energies of Peptide Ligands to Class A G-Protein-Coupled Receptors.

A metadynamics protocol is presented to characterize the binding and unbinding of peptide ligands to class A G-protein-coupled receptors (GPCRs). The protocol expands on the one previously presented for binding and unbinding small-molecule ligands to class A GPCRs and accounts for the more demanding nature of the peptide binding-unbinding process. It applies to almost all class A GPCRs. Exemplary simulations are described for subtypes Y1R, Y2R, and Y4R of the neuropeptide Y receptor family, vasopressin binding to the vasopressin V2 receptor (V2R), and oxytocin binding to the oxytocin receptor (OTR). Binding free energies and the positions of alternative binding sites are presented and, where possible, compared with the experiment.

Stereochemical Stability of Planar-Chiral Benzazepine Tricyclics: Inversion Energies of P- and S-Alkene Ligands.

Inversion barriers ΔG‡ for planar chiral phosphine-alkene and sulfonamide-alkene hybrid ligands based on phenyl-dibenz[b,f]azepine have been determined by density-functional theory calculations. Analysis of the structural and electronic characteristics of the minima and transition states explains the magnitudes of ΔG‡ and the geometrical changes during the inversion process. The steric repulsion caused by bulky substituents attached to the azepine nitrogen atom has a pronounced effect on the ΔG‡ value, explaining, inter alia, the stereochemical stability of the P- and S-alkene ligands when compared to the fluxional parent compound where X = H.

Activation/Deactivation Free-Energy Profiles for the ß2-Adrenergic Receptor: Ligand Modes of Action.

We use enhanced-sampling simulations with an effective collective variable to study the activation of the ß2-adrenergic receptor in the presence of ligands with different efficacy. The free-energy profiles are computed for the ligand-free (apo) receptor and binary (apo-receptor + G-protein α-subunit and receptor + ligand) and ternary complexes. The results are not only compatible with available experiments but also allow unprecedented structural insight into the nature of GPCR conformations along the activation pathway and their role in the activation mechanism. In particular, the simulations reveal an unexpected mode of action of partial agonists such as salmeterol and salbutamol that arises already in the binary complex without the G-protein. Specific differences in the polar interactions with residues in TM5, which are required to stabilize an optimal TM6 conformation that facilitates G-protein binding and receptor activation, play a major role in differentiating them from full agonists.

Stereochemistry-Driven Interactions of α,γ-Peptide Ligands with the Neuropeptide Y Y4-Receptor.

The G-protein-coupled Y4-receptor (Y4R) and its endogenous ligand, pancreatic polypeptide (PP), suppress appetite in response to food intake and, thus, are attractive drug targets for body-weight control. The C-terminus of human PP (hPP), T32-R33-P34-R35-Y36-NH2, penetrates deep into the binding pocket with its tyrosine-amide and di-arginine motif. Here, we present two C-terminally amidated α,γ-hexapeptides (1a/b) with sequence Ac-R31-γ-CBAA32-R33-L34-R35-Y36-NH2, where γ-CBAA is the (1R,2S,3R)-configured 2-(aminomethyl)-3-phenylcyclobutanecarboxyl moiety (1a) or its mirror image (1b). Both peptides bind the Y4R (Ki of 1a/b: 0.66/12 nM) and act as partial agonists (intrinsic activity of 1a/b: 50/39%). Their induced-fit binding poses in the Y4R pocket are unique and build ligand-receptor contacts distinct from those of the C-terminus of the endogenous ligand hPP. We conclude that energetically favorable interactions, although they do not match those of the native ligand hPP, still guarantee high binding affinity (with 1a rivaling hPP) but not the maximum receptor activation.

General Metadynamics Protocol To Simulate Activation/Deactivation of Class A GPCRs: Proof of Principle for the Serotonin Receptor.

We present a generally applicable metadynamics protocol for characterizing the activation free-energy profiles of class A G-protein coupled receptors and a proof-of-principle study for the 5HT1A-receptor. The almost universal A100 activation index, which depends on five inter-helix distances, is used as the single collective variable in well-tempered multiple-walker metadynamics simulations. Here, we show free-energy profiles for the serotonin receptor as binary (apo-receptor + G-protein-α-subunit and receptor + ligand) and ternary complexes with two prototypical orthosteric ligands: the full agonist serotonin and the partial agonist aripiprazole. Our results are not only compatible with previously reported experimental and computational data, but they also allow differences between active and inactive conformations to be determined in unprecedented atomic detail, and with respect to the so-called microswitches that have been suggested as determinants of activation, giving insight into their role in the activation mechanism.