Atomic-Level Characterization of the Activation Mechanism of SERCA by Calcium

We have performed molecular dynamics (MD) simulations to elucidate, in atomic detail, the mechanism by which the sarcoplasmic reticulum Ca2+-ATPase (SERCA) is activated by Ca2+. Crystal structures suggest that activation of SERCA occurs when the cytoplasmic head-piece, in an open (E1) conformation stabilized by Ca2+, undergoes a large-scale open-to-closed (E1 to E2) transition that is induced by ATP binding. However, spectroscopic measurements in solution suggest that these structural states (E1 and E2) are not tightly coupled to biochemical states (defined by bound ligands); the closed E2 state predominates even in the absence of ATP, in both the presence and absence of Ca2+. How is this loose coupling consistent with the high efficiency of energy transduction in the Ca2+-ATPase? To provide insight into this question, we performed long (500 ns) all-atom MD simulations starting from the open crystal structure, including a lipid bilayer and water. In both the presence and absence of Ca2+, we observed a large-scale open-to-closed conformational transition within 400 ns, supporting the weak coupling between structural and biochemical states. However, upon closer inspection, it is clear that Ca2+ is necessary and sufficient for SERCA to reach the precise geometrical arrangement necessary for activation of ATP hydrolysis. Contrary to suggestions from crystal structures, but in agreement with solution spectroscopy, the presence of ATP is not required for this activating transition. Principal component analysis showed that Ca2+ reshapes the free energy landscape of SERCA to create a path between the open conformation and the activated closed conformation. Thus the malleability of the free energy landscape is essential for SERCA efficiency, ensuring that ATP hydrolysis is tightly coupled to Ca2+ transport. These results demonstrate the importance of real-time dynamics in the formation of catalytically competent conformations of SERCA, with broad implications for understanding enzymatic catalysis in atomic detail

Molecular Basis of Ligand Dissociation in β-Adrenergic Receptors

The important and diverse biological functions of β-adrenergic receptors (βARs) have promoted the search for compounds to stimulate or inhibit their activity. In this regard, unraveling the molecular basis of ligand binding/unbinding events is essential to understand the pharmacological properties of these G protein-coupled receptors. In this study, we use the steered molecular dynamics simulation method to describe, in atomic detail, the unbinding process of two inverse agonists, which have been recently co-crystallized with β1 and β2ARs subtypes, along four different channels. Our results indicate that this type of compounds likely accesses the orthosteric binding site of βARs from the extracellular water environment. Importantly, reconstruction of forces and energies from the simulations of the dissociation process suggests, for the first time, the presence of secondary binding sites located in the extracellular loops 2 and 3 and transmembrane helix 7, where ligands are transiently retained by electrostatic and Van der Waals interactions. Comparison of the residues that form these new transient allosteric binding sites in both βARs subtypes reveals the importance of non-conserved electrostatic interactions as well as conserved aromatic contacts in the early steps of the binding process