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.

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.

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.

Architecture of the MKK6-p38α complex defines the basis of MAPK specificity and activation.

The mitogen-activated protein kinase (MAPK) p38α is a central component of signaling in inflammation and the immune response and is, therefore, an important drug target. Little is known about the molecular mechanism of its activation by double phosphorylation from MAPK kinases (MAP2Ks), because of the challenge of trapping a transient and dynamic heterokinase complex. We applied a multidisciplinary approach to generate a structural model of p38α in complex with its MAP2K, MKK6, and to understand the activation mechanism. Integrating cryo-electron microscopy with molecular dynamics simulations, hydrogen-deuterium exchange mass spectrometry, and experiments in cells, we demonstrate a dynamic, multistep phosphorylation mechanism, identify catalytically relevant interactions, and show that MAP2K-disordered amino termini determine pathway specificity. Our work captures a fundamental step of cell signaling: a kinase phosphorylating its downstream target kinase.

Machine Learning and Enhanced Sampling Simulations for Computing the Potential of Mean Force and Standard Binding Free Energy.

Computational capabilities are rapidly increasing, primarily because of the availability of GPU-based architectures. This creates unprecedented simulative possibilities for the systematic and robust computation of thermodynamic observables, including the free energy of a drug binding to a target. In contrast to calculations of relative binding free energy, which are nowadays widely exploited for drug discovery, we here push the boundary of computing the binding free energy and the potential of mean force. We introduce a novel protocol that leverages enhanced sampling, machine learning, and ad hoc algorithms to limit human intervention, computing time, and free parameters in free energy calculations. We first validate the method on a host-guest system, and then we apply the protocol to glycogen synthase kinase 3 beta, a protein kinase of pharmacological interest. Overall, we obtain a good correlation with experimental values in relative and absolute terms. While we focus on protein-ligand binding, the strategy is of broad applicability to any complex event that can be described with a path collective variable. We systematically discuss key details that influence the final result. The parameters and simulation settings are available at PLUMED-NEST to allow full reproducibility.

Investigating Drug-Target Residence Time in Kinases through Enhanced Sampling Simulations.

It is widely accepted that drug-target association and dissociation rates directly affect drug efficacy and safety. To rationally optimize drug binding kinetics, one must know the atomic arrangement of the protein-ligand complex during the binding/unbinding process in order to detect stable and metastable states. Whereas experimental approaches can determine kinetic constants with fairly good accuracy, computational approaches based on molecular dynamics (MD) simulations can deliver the atomistic details of the unbinding process. Furthermore, they can also be utilized prospectively to predict residence time (i.e., the inverse of unbinding kinetics constant, koff) with an acceptable level of accuracy. Here, we report a novel method based on adiabatic bias MD with an electrostatics-like collective variable (dubbed elABMD) for sampling protein-ligand dissociation events in two kinases. elABMD correctly ranked a ligand series on glucokinase, in agreement with experimental data and previous calculations. Subsequently, we applied the new method prospectively to a congeneric series of GSK-3ß inhibitors. For this series, new crystal structures were generated and the residence time was experimentally measured with surface plasmon resonance (SPR). There was good agreement between computational predictions and experimental measures, suggesting that elABMD is an innovative and efficient tool for calculating residence times.