Desensitization dynamics of the AMPA receptor.

To perform their physiological functions, amino methyl propionic acid receptors (AMPARs) cycle through active, resting, and desensitized states, and dysfunction in AMPAR activity is associated with various neurological disorders. Transitions among AMPAR functional states, however, are largely uncharacterized at atomic resolution and are difficult to examine experimentally. Here, we report long-timescale molecular dynamics simulations of dimerized AMPAR ligand-binding domains (LBDs), whose conformational changes are tightly coupled to changes in AMPAR functional states, in which we observed LBD dimer activation and deactivation upon ligand binding and unbinding at atomic resolution.

Molecular Basis of Small-Molecule Binding to α-Synuclein.

Intrinsically disordered proteins (IDPs) are implicated in many human diseases. They have generally not been amenable to conventional structure-based drug design, however, because their intrinsic conformational variability has precluded an atomic-level understanding of their binding to small molecules. Here we present long-time-scale, atomic-level molecular dynamics (MD) simulations of monomeric α-synuclein (an IDP whose aggregation is associated with Parkinson's disease) binding the small-molecule drug fasudil in which the observed protein-ligand interactions were found to be in good agreement with previously reported NMR chemical shift data.

Structural mechanism for Bruton's tyrosine kinase activation at the cell membrane.

Bruton's tyrosine kinase (Btk) is critical for B cell proliferation and activation, and the development of Btk inhibitors is a vigorously pursued strategy for the treatment of various B cell malignancies. A detailed mechanistic understanding of Btk activation has, however, been lacking. Here, inspired by a previous suggestion that Btk activation might depend on dimerization of its lipid-binding PH-TH module on the cell membrane, we performed long-timescale molecular dynamics simulations of membrane-bound PH-TH modules and observed that they dimerized into a single predominant conformation.

Atomic-level characterization of protein-protein association.

Despite the biological importance of protein-protein complexes, determining their structures and association mechanisms remains an outstanding challenge. Here, we report the results of atomic-level simulations in which we observed five protein-protein pairs repeatedly associate to, and dissociate from, their experimentally determined native complexes using a molecular dynamics (MD)-based sampling approach that does not make use of any prior structural information about the complexes. To study association mechanisms, we performed additional, conventional MD simulations, in which we observed numerous spontaneous association events.

Entry from the Lipid Bilayer: A Possible Pathway for Inhibition of a Peptide G Protein-Coupled Receptor by a Lipophilic Small Molecule.

The pathways that G protein-coupled receptor (GPCR) ligands follow as they bind to or dissociate from their receptors are largely unknown. Protease-activated receptor-1 (PAR1) is a GPCR activated by intramolecular binding of a tethered agonist peptide that is exposed by thrombin cleavage. By contrast, the PAR1 antagonist vorapaxar is a lipophilic drug that binds in a pocket almost entirely occluded from the extracellular solvent.

Atomic structure of a toxic, oligomeric segment of SOD1 linked to amyotrophic lateral sclerosis (ALS).

Fibrils and oligomers are the aggregated protein agents of neuronal dysfunction in ALS diseases. Whereas we now know much about fibril architecture, atomic structures of disease-related oligomers have eluded determination. Here, we determine the corkscrew-like structure of a cytotoxic segment of superoxide dismutase 1 (SOD1) in its oligomeric state.

Quantitative Characterization of the Binding and Unbinding of Millimolar Drug Fragments with Molecular Dynamics Simulations.

A quantitative characterization of the binding properties of drug fragments to a target protein is an important component of a fragment-based drug discovery program. Fragments typically have a weak binding affinity, however, making it challenging to experimentally characterize key binding properties, including binding sites, poses, and affinities. Direct simulation of the binding equilibrium by molecular dynamics (MD) simulations can provide a computational route to characterize fragment binding, but this approach is so computationally intensive that it has thus far remained relatively unexplored.

Molecular Basis of Ligand Dissociation from the Adenosine A2A Receptor.

How drugs dissociate from their targets is largely unknown. We investigated the molecular basis of this process in the adenosine A2Areceptor (A2AR), a prototypical G protein-coupled receptor (GPCR). Through kinetic radioligand binding experiments, we characterized mutant receptors selected based on molecular dynamic simulations of the antagonist ZM241385 dissociating from the A2AR.

Demonstrating an Order-of-Magnitude Sampling Enhancement in Molecular Dynamics Simulations of Complex Protein Systems.

Molecular dynamics (MD) simulations can describe protein motions in atomic detail, but transitions between protein conformational states sometimes take place on time scales that are infeasible or very expensive to reach by direct simulation. Enhanced sampling methods, the aim of which is to increase the sampling efficiency of MD simulations, have thus been extensively employed. The effectiveness of such methods when applied to complex biological systems like proteins, however, has been difficult to establish because even enhanced sampling simulations of such systems do not typically reach time scales at which convergence is extensive enough to reliably quantify sampling efficiency.

Assessing the Accuracy of Two Enhanced Sampling Methods Using EGFR Kinase Transition Pathways: The Influence of Collective Variable Choice.

Structurally elucidating transition pathways between protein conformations gives deep mechanistic insight into protein behavior but is typically difficult. Unbiased molecular dynamics (MD) simulations provide one solution, but their computational expense is often prohibitive, motivating the development of enhanced sampling methods that accelerate conformational changes in a given direction, embodied in a collective variable. The accuracy of such methods is unclear for complex protein transitions, because obtaining unbiased MD data for comparison is difficult.

Structural basis for modulation of a G-protein-coupled receptor by allosteric drugs.

The design of G-protein-coupled receptor (GPCR) allosteric modulators, an active area of modern pharmaceutical research, has proved challenging because neither the binding modes nor the molecular mechanisms of such drugs are known. Here we determine binding sites, bound conformations and specific drug-receptor interactions for several allosteric modulators of the M2 muscarinic acetylcholine receptor (M2 receptor), a prototypical family A GPCR, using atomic-level simulations in which the modulators spontaneously associate with the receptor. Despite substantial structural diversity, all modulators form cation-π interactions with clusters of aromatic residues in the receptor extracellular vestibule, approximately 15 Å from the classical, 'orthosteric' ligand-binding site.

Recovery from slow inactivation in K+ channels is controlled by water molecules.

Application of a specific stimulus opens the intracellular gate of a K(+) channel (activation), yielding a transient period of ion conduction until the selectivity filter spontaneously undergoes a conformational change towards a non-conductive state (inactivation). Removal of the stimulus closes the gate and allows the selectivity filter to interconvert back to its conductive conformation (recovery). Given that the structural differences between the conductive and inactivated filter are very small, it is unclear why the recovery process can take up to several seconds.

Transitions to catalytically inactive conformations in EGFR kinase.

The epidermal growth factor receptor (EGFR) is a key protein in cellular signaling, and its kinase domain (EGFR kinase) is an intensely pursued target of small-molecule drugs. Although both catalytically active and inactive conformations of EGFR kinase have been resolved crystallographically, experimental characterization of the transitions between these conformations remains difficult. Using unbiased, all-atom molecular dynamics simulations, we observed EGFR kinase spontaneously transition from the active to the so-called "Src-like inactive" conformation by way of two sets of intermediate conformations: One corresponds to a previously identified locally disordered state and the other to previously undescribed "extended" conformations, marked by the opening of the ATP-binding site between the two lobes of the kinase domain.

Molecular determinants of drug-receptor binding kinetics.

It is increasingly appreciated that the rates at which drugs associate with and dissociate from receptors--the binding kinetics--directly impact drug efficacy and safety. The molecular determinants of drug-receptor binding kinetics remain poorly understood, however, especially when compared with the well-known factors that affect binding affinity. The rational modulation of kinetics during lead optimization thus remains challenging.

The dynamic process of β(2)-adrenergic receptor activation.

G-protein-coupled receptors (GPCRs) can modulate diverse signaling pathways, often in a ligand-specific manner. The full range of functionally relevant GPCR conformations is poorly understood. Here, we use NMR spectroscopy to characterize the conformational dynamics of the transmembrane core of the β(2)-adrenergic receptor (β(2)AR), a prototypical GPCR.

Mechanism of Cd2+ coordination during slow inactivation in potassium channels.

In K+ channels, rearrangements of the pore outer vestibule have been associated with C-type inactivation gating. Paradoxically, the crystal structure of Open/C-type inactivated KcsA suggests these movements to be modest in magnitude. In this study, we show that under physiological conditions, the KcsA outer vestibule undergoes relatively large dynamic rearrangements upon inactivation.

Structure and dynamics of the M3 muscarinic acetylcholine receptor.

Acetylcholine, the first neurotransmitter to be identified, exerts many of its physiological actions via activation of a family of G-protein-coupled receptors (GPCRs) known as muscarinic acetylcholine receptors (mAChRs). Although the five mAChR subtypes (M1-M5) share a high degree of sequence homology, they show pronounced differences in G-protein coupling preference and the physiological responses they mediate. Unfortunately, despite decades of effort, no therapeutic agents endowed with clear mAChR subtype selectivity have been developed to exploit these differences.

Thermodynamic coupling between activation and inactivation gating in potassium channels revealed by free energy molecular dynamics simulations.

The amount of ionic current flowing through K(+) channels is determined by the interplay between two separate time-dependent processes: activation and inactivation gating. Activation is concerned with the stimulus-dependent opening of the main intracellular gate, whereas inactivation is a spontaneous conformational transition of the selectivity filter toward a nonconductive state occurring on a variety of timescales. A recent analysis of multiple x-ray structures of open and partially open KcsA channels revealed the mechanism by which movements of the inner activation gate, formed by the inner helices from the four subunits of the pore domain, bias the conformational changes at the selectivity filter toward a nonconductive inactivated state.

Activation mechanism of the β2-adrenergic receptor.

A third of marketed drugs act by binding to a G-protein-coupled receptor (GPCR) and either triggering or preventing receptor activation. Although recent crystal structures have provided snapshots of both active and inactive functional states of GPCRs, these structures do not reveal the mechanism by which GPCRs transition between these states. Here we propose an activation mechanism for the β(2)-adrenergic receptor, a prototypical GPCR, based on atomic-level simulations in which an agonist-bound receptor transitions spontaneously from the active to the inactive crystallographically observed conformation.

Pathway and mechanism of drug binding to G-protein-coupled receptors.

How drugs bind to their receptors--from initial association, through drug entry into the binding pocket, to adoption of the final bound conformation, or "pose"--has remained unknown, even for G-protein-coupled receptor modulators, which constitute one-third of all marketed drugs. We captured this pharmaceutically critical process in atomic detail using the first unbiased molecular dynamics simulations in which drug molecules spontaneously associate with G-protein-coupled receptors to achieve final poses matching those determined crystallographically. We found that several beta blockers and a beta agonist all traverse the same well-defined, dominant pathway as they bind to the β(1)- and β(2)-adrenergic receptors, initially making contact with a vestibule on each receptor's extracellular surface.

On the structural basis of modal gating behavior in K(+) channels.

Modal-gating shifts represent an effective regulatory mechanism by which ion channels control the extent and time course of ionic fluxes. Under steady-state conditions, the K(+) channel KcsA shows three distinct gating modes, high-P(o), low-P(o) and a high-frequency flicker mode, each with about an order of magnitude difference in their mean open times. Here we show that in the absence of C-type inactivation, mutations at the pore-helix position Glu71 unmask a series of kinetically distinct modes of gating in a side chain-specific way.

Structural basis for the coupling between activation and inactivation gates in K(+) channels.

The coupled interplay between activation and inactivation gating is a functional hallmark of K(+) channels. This coupling has been experimentally demonstrated through ion interaction effects and cysteine accessibility, and is associated with a well defined boundary of energetically coupled residues. The structure of the K(+) channel KcsA in its fully open conformation, in addition to four other partial channel openings, richly illustrates the structural basis of activation-inactivation gating.

Building Markov state models along pathways to determine free energies and rates of transitions.

An efficient method is proposed for building Markov models with discrete states able to accurately describe the slow relaxation of a complex system with two stable conformations. First, the reaction pathway described by a set of collective variables between the two stable states is determined using the string method with swarms of trajectories. Then, short trajectories are initiated at different points along this pathway to build the state-to-state transition probability matrix.

Finding transition pathways using the string method with swarms of trajectories.

An approach to find transition pathways in complex systems is presented. The method, which is related to the string method in collective variables of Maragliano et al. (J.