Structure of a GRK5-Calmodulin Complex Reveals Molecular Mechanism of GRK Activation and Substrate Targeting.

The phosphorylation of G protein-coupled receptors (GPCRs) by GPCR kinases (GRKs) facilitates arrestin binding and receptor desensitization. Although this process can be regulated by Ca2+-binding proteins such as calmodulin (CaM) and recoverin, the molecular mechanisms are poorly understood. Here, we report structural, computational, and biochemical analysis of a CaM complex with GRK5, revealing how CaM shapes GRK5 response to calcium. The CaM N and C domains bind independently to two helical regions at the GRK5 N and C termini to inhibit GPCR phosphorylation, though only the C domain interaction disrupts GRK5 membrane association, thereby facilitating cytoplasmic translocation. The CaM N domain strongly activates GRK5 via ordering of the amphipathic αN-helix of GRK5 and allosteric disruption of kinase-RH domain interaction for phosphorylation of cytoplasmic GRK5 substrates. These results provide a framework for understanding how two functional effects, GRK5 activation and localization, can cooperate under control of CaM for selective substrate targeting by GRK5.

A positively tuned voltage indicator for extended electrical recordings in the brain.

Genetically encoded voltage indicators (GEVIs) enable optical recording of electrical signals in the brain, providing subthreshold sensitivity and temporal resolution not possible with calcium indicators. However, one- and two-photon voltage imaging over prolonged periods with the same GEVI has not yet been demonstrated. Here, we report engineering of ASAP family GEVIs to enhance photostability by inversion of the fluorescence-voltage relationship. Two of the resulting GEVIs, ASAP4b and ASAP4e, respond to 100-mV depolarizations with ≥180% fluorescence increases, compared with the 50% fluorescence decrease of the parental ASAP3. With standard microscopy equipment, ASAP4e enables single-trial detection of spikes in mice over the course of minutes. Unlike GEVIs previously used for one-photon voltage recordings, ASAP4b and ASAP4e also perform well under two-photon illumination. By imaging voltage and calcium simultaneously, we show that ASAP4b and ASAP4e can identify place cells and detect voltage spikes with better temporal resolution than commonly used calcium indicators. Thus, ASAP4b and ASAP4e extend the capabilities of voltage imaging to standard one- and two-photon microscopes while improving the duration of voltage recordings.

Improving the quality of co-evolution intermolecular contact prediction with DisVis.

The steep rise in protein sequences and structures has paved the way for bioinformatics approaches to predict residue-residue interactions in protein complexes. Multiple sequence alignments are commonly used in contact predictions to identify co-evolving residues. These contacts, however, often include false positives (FPs), which may impair their use to predict three dimensional structures of biomolecular complexes and affect the accuracy of the generated models. Previously, we have developed DisVis to identify FP in mass spectrometry cross-linking data. DisVis allows to assess the accessible interaction space between two proteins consistent with a set of distance restraints. Here, we investigate if a similar approach could be applied to co-evolution predicted contacts in order to improve their precision prior to using them for modeling. We analyze co-evolution contact predictions with DisVis for a set of 26 protein-protein complexes. The DisVis-reranked and the original co-evolution contacts are then used to model the complexes with our integrative docking software HADDOCK using different filtering scenarios. Our results show that HADDOCK is robust with respect to the precision of the predicted contacts due to the 50% random contact removal during docking and can enhance the quality of docking predictions when combined with DisVis filtering for low precision contact data. DisVis can thus have a beneficial effect on low quality data, but overall HADDOCK can accommodate FP restraints without negatively impacting the quality of the resulting models. Other more precision-sensitive docking protocols might, however, benefit from the increased precision of the predicted contacts after DisVis filtering.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

Regulation of adenylyl cyclase 5 in striatal neurons confers the ability to detect coincident neuromodulatory signals.

Long-term potentiation and depression of synaptic activity in response to stimuli is a key factor in reinforcement learning. Strengthening of the corticostriatal synapses depends on the second messenger cAMP, whose synthesis is catalysed by the enzyme adenylyl cyclase 5 (AC5), which is itself regulated by the stimulatory Gαolf and inhibitory Gαi proteins. AC isoforms have been suggested to act as coincidence detectors, promoting cellular responses only when convergent regulatory signals occur close in time. However, the mechanism for this is currently unclear, and seems to lie in their diverse regulation patterns. Despite attempts to isolate the ternary complex, it is not known if Gαolf and Gαi can bind to AC5 simultaneously, nor what activity the complex would have. Using protein structure-based molecular dynamics simulations, we show that this complex is stable and inactive. These simulations, along with Brownian dynamics simulations to estimate protein association rates constants, constrain a kinetic model that shows that the presence of this ternary inactive complex is crucial for AC5's ability to detect coincident signals, producing a synergistic increase in cAMP. These results reveal some of the prerequisites for corticostriatal synaptic plasticity, and explain recent experimental data on cAMP concentrations following receptor activation. Moreover, they provide insights into the regulatory mechanisms that control signal processing by different AC isoforms.

Association of Both Inhibitory and Stimulatory Gα Subunits Implies Adenylyl Cyclase 5 Deactivation.

Adenylyl cyclase (AC) generates cyclic AMP required for a variety of cellular functions, and its regulation plays a major role in cellular signal transduction in eukaryotes and prokaryotes. All membrane-bound AC isoforms in eukaryotes can be activated by stimulatory G-proteins, but only AC1, AC5, and AC6 can be both stimulated and inhibited by active Gα subunits, Gαs and Gαi, respectively. In principle, these Gαi-sensitive AC isoforms could form both binary and ternary complexes with Gα subunits due to the noncompetitive association of inhibitory and stimulatory Gα. However, the formation and possible catalytic activity of a putative ternary complex have not yet been experimentally confirmed due to its proposed short-lived nature. Here, the catalytic activity of such a ternary complex consisting of apo AC5, stimulatory Gαolf, and inhibitory Gαi1 is investigated via classical molecular dynamics simulations. Trajectories of inhibited and stimulated binary complexes, AC5:Gαi1 and AC5:Gαolf, respectively, as well as Gα-free AC5 were also obtained to compare the sampled AC5 conformation in the ternary complex to those sampled under different Gα conditions. This comparison suggests that association of both Gα subunits results in an AC5 conformation similar to that sampled by the AC5:Gαi1 complex, indicating that the ternary complex mainly samples an inactive conformation.

Exploring the inhibition mechanism of adenylyl cyclase type 5 by n-terminal myristoylated Gαi1.

Adenylyl cyclase (AC) is an important messenger involved in G-protein-coupled-receptor signal transduction pathways, which is a well-known target for drug development. AC is regulated by activated stimulatory (Gαs) and inhibitory (Gαi) G proteins in the cytosol. Although experimental studies have shown that these Gα subunits can stimulate or inhibit AC's function in a non-competitive way, it is not well understood what the difference is in their mode of action as both Gα subunits appear structurally very similar in a non-lipidated state. However, a significant difference between Gαs and Gαi is that while Gαs does not require any lipidation in order to stimulate AC, N-terminal myristoylation is crucial for Gαi's inhibitory function as AC is not inhibited by non-myristoylated Gαi. At present, only the conformation of the complex including Gαs and AC has been resolved via X-ray crystallography. Therefore, understanding the interaction between Gαi and AC is important as it will provide more insight into the unknown mechanism of AC regulation. This study demonstrates via classical molecular dynamics simulations that the myristoylated Gαi1 structure is able to interact with apo adenylyl cyclase type 5 in a way that causes inhibition of the catalytic function of the enzyme, suggesting that Gα lipidation could play a crucial role in AC regulation and in regulating G protein function by affecting Gαi's active conformation.

Effect of N-Terminal Myristoylation on the Active Conformation of Gαi1-GTP.

G proteins are part of the G-protein-coupled receptor (GPCR) signal transduction cascade in which they transfer a signal from the membrane-embedded GPCR to other proteins in the cell. In the case of the inhibitory G-protein heterotrimer, permanent N-terminal myristoylation can transiently localize the Gαi subunit at the membrane as well as crucially influence Gαi's function in the GTP-bound conformation. The attachment of lipids to proteins is known to be essential for membrane trafficking; however, our results suggest that lipidation is also important for protein-protein interactions during signal transduction. Here we investigate the effect of myristoylation on the structure and dynamics of soluble Gαi1 and its possible implication for signal transduction. A 2 µs classical molecular dynamics simulation of a myristoylated Gαi1-GTP complex suggests that the myristoyl-induced conformational changes of the switch II and alpha helical domains create new possibilities for protein-protein interactions and emphasize the importance of permanent lipid attachment for the conformation and functional tunability of signaling proteins.

Interface refinement of low- to medium-resolution Cryo-EM complexes using HADDOCK2.4.

A wide range of cellular processes requires the formation of multimeric protein complexes. The rise of cryo-electron microscopy (cryo-EM) has enabled the structural characterization of these protein assemblies. The density maps produced can, however, still suffer from limited resolution, impeding the process of resolving structures at atomic resolution. In order to solve this issue, monomers can be fitted into low- to medium-resolution maps. Unfortunately, the models produced frequently contain atomic clashes at the protein-protein interfaces (PPIs), as intermolecular interactions are typically not considered during monomer fitting. Here, we present a refinement approach based on HADDOCK2.4 to remove intermolecular clashes and optimize PPIs. A dataset of 14 cryo-EM complexes was used to test eight protocols. The best-performing protocol, consisting of a semi-flexible simulated annealing refinement with centroid restraints on the monomers, was able to decrease intermolecular atomic clashes by 98% without significantly deteriorating the quality of the cryo-EM density fit.

Cyclization and Docking Protocol for Cyclic Peptide-Protein Modeling Using HADDOCK2.4.

An emerging class of therapeutic molecules are cyclic peptides with over 40 cyclic peptide drugs currently in clinical use. Their mode of action is, however, not fully understood, impeding rational drug design. Computational techniques could positively impact their design, but modeling them and their interactions remains challenging due to their cyclic nature and their flexibility. This study presents a step-by-step protocol for generating cyclic peptide conformations and docking them to their protein target using HADDOCK2.4. A dataset of 30 cyclic peptide-protein complexes was used to optimize both cyclization and docking protocols. It supports peptides cyclized via an N- and C-terminus peptide bond and/or a disulfide bond. An ensemble of cyclic peptide conformations is then used in HADDOCK to dock them onto their target protein using knowledge of the binding site on the protein side to drive the modeling. The presented protocol predicts at least one acceptable model according to the critical assessment of prediction of interaction criteria for each complex of the dataset when the top 10 HADDOCK-ranked single structures are considered (100% success rate top 10) both in the bound and unbound docking scenarios. Moreover, its performance in both bound and fully unbound docking is similar to the state-of-the-art software in the field, Autodock CrankPep. The presented cyclization and docking protocol should make HADDOCK a valuable tool for rational cyclic peptide-based drug design and high-throughput screening.

Gαi1 inhibition mechanism of ATP-bound adenylyl cyclase type 5.

Conversion of adenosine triphosphate (ATP) to the second messenger cyclic adenosine monophosphate (cAMP) is an essential reaction mechanism that takes place in eukaryotes, triggering a variety of signal transduction pathways. ATP conversion is catalyzed by the enzyme adenylyl cyclase (AC), which can be regulated by binding inhibitory, Gαi, and stimulatory, Gαs subunits. In the past twenty years, several crystal structures of AC in isolated form and complexed to Gαs subunits have been resolved. Nevertheless, the molecular basis of the inhibition mechanism of AC, induced by Gαi, is still far from being fully understood. Here, classical molecular dynamics simulations of the isolated holo AC protein type 5 and the holo binary complex AC5:Gαi have been analyzed to investigate the conformational impact of Gαi association on ATP-bound AC5. The results show that Gαi appears to inhibit the activity of AC5 by preventing the formation of a reactive ATP conformation.

Erratum: Publisher's Note: "Implications of short time scale dynamics on long time processes" (Struct. Dyn. 4, 061507 (2017)].

[This corrects the article DOI: 10.1063/1.4996448.].

How Rhodopsin Tunes the Equilibrium between Protonated and Deprotonated Forms of the Retinal Chromophore.

Rhodopsin is a photoactive G-protein-coupled receptor (GPCR) that converts dim light into a signal for the brain, leading to eyesight. Full activation of this GPCR is achieved after passing through several steps of the protein's photoactivation pathway. Key events of rhodopsin activation are the initial cis-trans photoisomerization of the covalently bound retinal moiety followed by conformational rearrangements and deprotonation of the chromophore's protonated Schiff base (PSB), which ultimately lead to full activation in the meta II state. PSB deprotonation is crucial for achieving full activation of rhodopsin; however, the specific structural rearrangements that have to take place to induce this pKa shift are not well understood. Classical molecular dynamics (MD) simulations were employed to identify intermediate states after the cis-trans isomerization of rhodopsin's retinal moiety. In order to select the intermediate state in which PSB deprotonation is experimentally known to occur, the validity of the intermediate configurations was checked through an evaluation of the optical properties in comparison with experiment. Subsequently, the selected state was used to investigate the molecular factors that enable PSB deprotonation at body temperature to obtain a better understanding of the difference between the protonated and the deprotonated state of the chromophore. To this end, the deprotonation reaction has been investigated by applying QM/MM MD simulations in combination with thermodynamic integration. The study shows that, compared to the inactive 11-cis-retinal case, trans-retinal rhodopsin is able to undergo PSB deprotonation due to a change in the conformation of the retinal and a consequent alteration in the hydrogen-bond (HB) network in which PSB and the counterion Glu113 are embedded. Besides the retinal moiety and Glu113, also two water molecules as well as Thr94 and Gly90 that are related to congenital night blindness are part of this essential HB network.

Does Proton Conduction in the Voltage-Gated H+ Channel hHv1 Involve Grotthuss-Like Hopping via Acidic Residues?

Hv1s are ubiquitous highly selective voltage-gated proton channels involved in male fertility, immunology, and the invasiveness of certain forms of breast cancer. The mechanism of proton extrusion in Hv1 is not yet understood, while it constitutes the first step toward the design of high-affinity drugs aimed at this important pharmacological target. In this contribution, we explore the details of the mechanism via an integrative approach, using classical and QM/MM molecular dynamics simulations of a monomeric hHv1 model. We propose that protons localize in three binding sites along the channel lumen, formed by three pairs of conserved negatively charged residues lining the pore: D174/E153, D112/D185, and E119/D123. Local rearrangements, involving notably a dihedral transition of F150, a conserved phenylalanine lining the permeation pathway, appear to allow protons to hop from one acidic residue to the next through a bridging water molecule. These results constitute a first attempt at rationalizing hHv1 selectivity for H+ and the role played by D112 in this process. They pave the way for further quantitative characterization of H+ transport in hHv1.

Implications of short time scale dynamics on long time processes.

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.