Zooming in on protons: Neutron structure of protein kinase A trapped in a product complex

The αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic subunit of protein kinase A (PKA-C). This long-range synergistic action is involved in substrate recognition and fidelity, and it may also regulate PKA's association with regulatory subunits and other binding partners. To date, a complete understanding of this intramolecular mechanism is still lacking. Here, we integrated NMR(Nuclear Magnetic Resonance)-restrained molecular dynamics simulations and a Markov State Model to characterize the free energy landscape and conformational transitions of PKA-C. We found that the apoenzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the αC-β4 loop. We validated the second excited state by analyzing the F100A mutant of PKA-C, assessing its structural response to ATP and substrate binding. While PKA-CF100A preserves its catalytic efficiency with Kemptide, this mutation rearranges the αC-β4 loop conformation, interrupting the coupling of the two lobes and abolishing the allosteric binding cooperativity. The highly conserved αC-β4 loop emerges as a pivotal element to control the synergistic binding of nucleotide and substrate, explaining how mutations or insertions near or within this motif affect the function and drug sensitivity in homologous kinases.

Neutron diffraction from a microgravity-grown crystal reveals the active site hydrogens of the internal aldimine form of tryptophan synthase

Pyridoxal 5'-phosphate (PLP), the biologically active form of vitamin B6, is an essential cofactor in many biosynthetic pathways. The emergence of PLP-dependent enzymes as drug targets and biocatalysts, such as tryptophan synthase (TS), has underlined the demand to understand PLP-dependent catalysis and reaction specificity. The ability of neutron diffraction to resolve the positions of hydrogen atoms makes it an ideal technique to understand how the electrostatic environment and selective protonation of PLP regulates PLP-dependent activities. Facilitated by microgravity crystallization of TS with the Toledo Crystallization Box, we report the 2.1 Å joint X-ray/neutron (XN) structure of TS with PLP in the internal aldimine form. Positions of hydrogens were directly determined in both the α- and β-active sites, including PLP cofactor. The joint XN structure thus provides insight into the selective protonation of the internal aldimine and the electrostatic environment of TS necessary to understand the overall catalytic mechanism.

Protein Kinase Structure and Dynamics: Role of the αC-β4 Loop

Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, Local Spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.

The αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic (C) subunit of protein kinase A (PKA). Not only this long-range synergistic action is involved in substrate recognition and fidelity, but it is likely to regulate PKA association with regulatory subunits and other binding partners. To date, a complete understanding of the molecular determinants for this intramolecular mechanism is still lacking. Here, we used an integrated NMR-restrained molecular dynamics simulations and a Markov Model to characterize the free energy landscape and conformational transitions of the catalytic subunit of protein kinase A (PKA-C). We found that the apo-enzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the tip of the αC-β4 loop. To experimentally validate the second excited state, we mutated F100 into alanine and used NMR spectroscopy to characterize the binding thermodynamics and structural response of ATP and a prototypical peptide substrate. While the activity of PKA-CF100A toward a prototypical peptide substrate is unaltered and the enzyme retains its affinity for ATP and substrate, this mutation rearranges the αC-β4 loop conformation interrupting the allosteric coupling between nucleotide and substrate. The highly conserved αC-β4 loop emerges as a pivotal element able to modulate the synergistic binding between nucleotide and substrate and may affect PKA signalosome. These results may explain how insertion mutations within this motif affect drug sensitivity in other homologous kinases.

Revealing protonation states and tracking substrate in serine hydroxymethyltransferase with room-temperature X-ray and neutron crystallography

Pyridoxal 5'-phosphate (PLP)-dependent enzymes utilize a vitamin B6-derived cofactor to perform a myriad of chemical transformations on amino acids and other small molecules. Some PLP-dependent enzymes, such as serine hydroxymethyltransferase (SHMT), are promising drug targets for the design of small-molecule antimicrobials and anticancer therapeutics, while others have been used to synthesize pharmaceutical building blocks. Understanding PLP-dependent catalysis and the reaction specificity is crucial to advance structure-assisted drug design and enzyme engineering. Here we report the direct determination of the protonation states in the active site of Thermus thermophilus SHMT (TthSHMT) in the internal aldimine state using room-temperature joint X-ray/neutron crystallography. Conserved active site architecture of the model enzyme TthSHMT and of human mitochondrial SHMT (hSHMT2) were compared by obtaining a room-temperature X-ray structure of hSHMT2, suggesting identical protonation states in the human enzyme. The amino acid substrate serine pathway through the TthSHMT active site cavity was tracked, revealing the peripheral and cationic binding sites that correspond to the pre-Michaelis and pseudo-Michaelis complexes, respectively. At the peripheral binding site, the substrate is bound in the zwitterionic form. By analyzing the observed protonation states, Glu53, but not His residues, is proposed as the general base catalyst, orchestrating the retro-aldol transformation of L-serine into glycine.

A critical evaluation of protein kinase regulation by activation loop autophosphorylation

Phosphorylation of proteins is a ubiquitous mechanism of regulating their function, localization, or activity. Protein kinases, enzymes that use ATP to phosphorylate protein substrates are, therefore, powerful signal transducers in eukaryotic cells. The mechanism of phosphoryl-transfer is universally conserved among protein kinases, which necessitates the tight regulation of kinase activity for the orchestration of cellular processes with high spatial and temporal fidelity. In response to a stimulus, many kinases enhance their own activity by autophosphorylating a conserved amino acid in their activation loop, but precisely how this reaction is performed is controversial. Classically, kinases that autophosphorylate their activation loop are thought to perform the reaction in trans, mediated by transient dimerization of their kinase domains. However, motivated by the recently discovered regulation mechanism of activation loop cis-autophosphorylation by a kinase that is autoinhibited in trans, we here review the various mechanisms of autoregulation that have been proposed. We provide a framework for critically evaluating biochemical, kinetic, and structural evidence for protein kinase dimerization and autophosphorylation, and share some thoughts on the implications of these mechanisms within physiological signaling networks.

An NMR portrait of functional and dysfunctional allosteric cooperativity in cAMP-dependent protein kinase A

The cAMP-dependent protein kinase A (PKA) is the archetypical eukaryotic kinase. The catalytic subunit (PKA-C) structure is highly conserved among the AGC-kinase family. PKA-C is a bilobal enzyme with a dynamic N-lobe, harbouring the Adenosine-5'-triphosphate (ATP) binding site and a more rigid helical C-lobe. The substrate-binding groove resides at the interface of the two lobes. A distinct feature of PKA-C is the positive binding cooperativity between nucleotide and substrate. Several PKA-C mutations lead to the development of adenocarcinomas, myxomas, and other rare forms of liver tumours. Nuclear magnetic resonance (NMR) spectroscopy shows that these mutations disrupt the allosteric communication between the two lobes, causing a drastic decrease in binding cooperativity. The loss of cooperativity correlates with changes in substrate fidelity and reduced kinase affinity for the endogenous protein kinase inhibitor (PKI). The similarity between PKI and the inhibitory sequence of the kinase regulatory subunits suggests that the overall mechanism of regulation of the kinase may be disrupted. We surmise that a reduced or obliterated cooperativity may constitute a common trait for both orthosteric and allosteric mutations of PKA-C that may lead to dysregulation and disease.

Molecular-dynamics simulation methods for macromolecular crystallography

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein-solvent interface. To address this limitation, a set of MD-MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

An N⋯H⋯N low-barrier hydrogen bond preorganizes the catalytic site of aspartate aminotransferase to facilitate the second half-reaction

Pyridoxal 5'-phosphate (PLP)-dependent enzymes have been extensively studied for their ability to fine-tune PLP cofactor electronics to promote a wide array of chemistries. Neutron crystallography offers a straightforward approach to studying the electronic states of PLP and the electrostatics of enzyme active sites, responsible for the reaction specificities, by enabling direct visualization of hydrogen atom positions. Here we report a room-temperature joint X-ray/neutron structure of aspartate aminotransferase (AAT) with pyridoxamine 5'-phosphate (PMP), the cofactor product of the first half reaction catalyzed by the enzyme. Between PMP NSB and catalytic Lys258 Nζ amino groups an equally shared deuterium is observed in an apparent low-barrier hydrogen bond (LBHB). Density functional theory calculations were performed to provide further evidence of this LBHB interaction. The structural arrangement and the juxtaposition of PMP and Lys258, facilitated by the LBHB, suggests active site preorganization for the incoming ketoacid substrate that initiates the second half-reaction.

From structure to the dynamic regulation of a molecular switch: A journey over 3 decades

It is difficult to imagine where the signaling community would be today without the Protein Data Bank. This visionary resource, established in the 1970s, has been an essential partner for sharing information between academics and industry for over 3 decades. We describe here the history of our journey with the protein kinases using cAMP-dependent protein kinase as a prototype. We summarize what we have learned since the first structure, published in 1991, why our journey is still ongoing, and why it has been essential to share our structural information. For regulation of kinase activity, we focus on the cAMP-binding protein kinase regulatory subunits. By exploring full-length macromolecular complexes, we discovered not only allostery but also an essential motif originally attributed to crystal packing. Massive genomic data on disease mutations allows us to now revisit crystal packing as a treasure chest of possible protein:protein interfaces where the biological significance and disease relevance can be validated. It provides a new window into exploring dynamic intrinsically disordered regions that previously were deleted, ignored, or attributed to crystal packing. Merging of crystallography with cryo-electron microscopy, cryo-electron tomography, NMR, and millisecond molecular dynamics simulations is opening a new world for the signaling community where those structure coordinates, deposited in the Protein Data Bank, are just a starting point!

Direct Observation of Protonation State Modulation in SARS-CoV-2 Main Protease upon Inhibitor Binding with Neutron Crystallography

The main protease (3CL Mpro) from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, is an essential enzyme for viral replication with no human counterpart, making it an attractive drug target. To date, no small-molecule clinical drugs are available that specifically inhibit SARS-CoV-2 Mpro. To aid rational drug design, we determined a neutron structure of Mpro in complex with the α-ketoamide inhibitor telaprevir at near-physiological (22 °C) temperature. We directly observed protonation states in the inhibitor complex and compared them with those in the ligand-free Mpro, revealing modulation of the active-site protonation states upon telaprevir binding. We suggest that binding of other α-ketoamide covalent inhibitors can lead to the same protonation state changes in the Mpro active site. Thus, by studying the protonation state changes induced by inhibitors, we provide crucial insights to help guide rational drug design, allowing precise tailoring of inhibitors to manipulate the electrostatic environment of SARS-CoV-2 Mpro.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of Mpro, a cysteine protease, have been determined, facilitating structure-based drug design. Mpro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, Mpro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV Mpro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 Mpro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of Mpro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N δ (HD) and N ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 Mpro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

Unusual zwitterionic catalytic site of SARS-CoV-2 main protease revealed by neutron crystallography

The main protease (3CL Mpro) from SARS-CoV-2, the etiological agent of COVID-19, is an essential enzyme for viral replication. 3CL Mpro possesses an unusual catalytic dyad composed of Cys145 and His41 residues. A critical question in the field has been what the protonation states of the ionizable residues in the substrate-binding active-site cavity are; resolving this point would help understand the catalytic details of the enzyme and inform rational drug development against this pernicious virus. Here, we present the room-temperature neutron structure of 3CL Mpro, which allowed direct determination of hydrogen atom positions and, hence, protonation states in the protease. We observe that the catalytic site natively adopts a zwitterionic reactive form in which Cys145 is in the negatively charged thiolate state and His41 is doubly protonated and positively charged, instead of the neutral unreactive state usually envisaged. The neutron structure also identified the protonation states, and thus electrical charges, of all other amino acid residues and revealed intricate hydrogen-bonding networks in the active-site cavity and at the dimer interface. The fine atomic details present in this structure were made possible by the unique scattering properties of the neutron, which is an ideal probe for locating hydrogen positions and experimentally determining protonation states at near-physiological temperature. Our observations provide critical information for structure-assisted and computational drug design, allowing precise tailoring of inhibitors to the enzyme's electrostatic environment.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M ^pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M ^pro , a cysteine protease, have been determined, facilitating structure-based drug design. M ^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, M ^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV M ^pro , but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M ^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M ^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α -ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N _δ (HD) and N _ϵ (HE) protonation of His41 and His164, respectively, the α -ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M ^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

β adrenergic receptor modulated signaling in glioma models: promoting β adrenergic receptor-β arrestin scaffold-mediated activation of extracellular-regulated kinase 1/2 may prove to be a panacea in the treatment of intracranial and spinal malignancy and extra-neuraxial carcinoma

Neoplastically transformed astrocytes express functionally active cell surface β adrenergic receptors (βARs). Treatment of glioma models in vitro and in vivo with β adrenergic agonists variably amplifies or attenuates cellular proliferation. In the majority of in vivo models, β adrenergic agonists generally reduce cellular proliferation. However, treatment with β adrenergic agonists consistently reduces tumor cell invasive potential, angiogenesis, and metastasis. β adrenergic agonists induced decreases of invasive potential are chiefly mediated through reductions in the expression of matrix metalloproteinases types 2 and 9. Treatment with β adrenergic agonists also clearly reduce tumoral neoangiogenesis, which may represent a putatively useful mechanism to adjuvantly amplify the effects of bevacizumab. Bevacizumab is a monoclonal antibody targeting the vascular endothelial growth factor receptor. We may accordingly designate βagonists to represent an enhancer of bevacizumab. The antiangiogenic effects of β adrenergic agonists may thus effectively render an otherwise borderline effective therapy to generate significant enhancement in clinical outcomes. β adrenergic agonists upregulate expression of the major histocompatibility class II DR alpha gene, effectively potentiating the immunogenicity of tumor cells to tumor surveillance mechanisms. Authors have also demonstrated crossmodal modulation of signaling events downstream from the β adrenergic cell surface receptor and microtubular polymerization and depolymerization. Complex effects and desensitization mechanisms of the β adrenergic signaling may putatively represent promising therapeutic targets. Constant stimulation of the β adrenergic receptor induces its phosphorylation by β adrenergic receptor kinase (βARK), rendering it a suitable substrate for alternate binding by β arrestins 1 or 2. The binding of a β arrestin to βARK phosphorylated βAR promotes receptor mediated internalization and downregulation of cell surface receptor and contemporaneously generates a cell surface scaffold at the βAR. The scaffold mediated activation of extracellular regulated kinase 1/2, compared with protein kinase A mediated activation, preferentially favors cytosolic retention of ERK1/2 and blunting of nuclear translocation and ensuant pro-transcriptional activity. Thus, βAR desensitization and consequent scaffold assembly effectively retains the cytosolic homeostatic functions of ERK1/2 while inhibiting its pro-proliferative effects. We suggest these mechanisms specifically will prove quite promising in developing primary and adjuvant therapies mitigating glioma growth, angiogenesis, invasive potential, and angiogenesis. We suggest generating compounds and targeted mutations of the β adrenergic receptor favoring β arrestin binding and scaffold facilitated activation of ERK1/2 may hold potential promise and therapeutic benefit in adjuvantly treating most or all cancers. We hope our discussion will generate fruitful research endeavors seeking to exploit these mechanisms.

Catalytic Subunit of PKA as a Prototype of the Eukaryotic Protein Kinase Family

The catalytic subunit of protein kinase A (PKAc) is conserved in all eukaryotic protein kinases. PKAc consists of two lobes that form the catalytic cleft containing the ATP-binding, peptide-binding site, and catalytic sites. During folding, PKAc secondary structures organize so that the non-polar regions form a globular core, while mobile loops and tails are exposed and can act as regulatory elements. De novo synthesized PKAc is phosphorylated at the T-loop, resulting in the formation of the active center capable of high-affinity binding of co-substrates. The ATP-molecule "sticks" the two lobes together, whereas the binding of peptide substrate completes the active center formation. The resulting catalytic triad (γ-phosphate of ATP, hydroxyl of Ser/Thr residue of the protein substrate, and Asp166 carboxyl) occupies a position optimal for catalysis. During the catalytic cycle, dynamic reorganization of polar and hydrophobic interactions ensures PKAc transition from the open to the closed conformation and vice versa. Understanding the structural basis of functioning of eukaryotic protein kinases (ePKs) is essential for successful design of ePK modulators.

Molecular Insights into the Trapping Effect of Ca2+ in Protein Kinase A: A Molecular Dynamics Study

Protein kinase A has become a model system for the study of kinases, and therefore, a comprehensive understanding of the underlying molecular mechanisms in its catalytic cycle is of crucial importance. One of the aspects that has received recent attention is the role that metal cofactors play in the catalytic cycle. Although Mg²⁺ is the well-known physiological ion used by protein kinases, Ca²⁺ ions can also assist the phosphoryl transfer reaction but with lower catalytic activities. This inhibitory effect has been attributed to the ability of Ca²⁺ to trap the reaction products at the active site, and it has been proposed as a possible regulatory mechanism of the enzyme. Thus, in order to get a clearer understanding of these molecular events, computational simulations in the product state of PKA, in the presence of Mg²⁺ and Ca²⁺ ions, were performed through molecular dynamics (MD). Different protonation states of the active site were considered in order to model the different mechanistic pathways that have been proposed. Our results show that different protonation states of the phosphorylated serine residue at the peptide substrate (pSer21), as well as the protonation state of residue Asp166, can have a marked influence on the flexibility of regions surrounding the active site. This is the case of the glycine-rich loop, a structural motif that is directly involved in the release of the products from the PKA active site. MD simulations were capable to reproduce the crystallographic conformations but also showed other conformations not previously reported in the crystal structures that may be involved in enhancing the affinity of pSP20 to PKA in the presence of Ca²⁺. Hydrogen bonding interactions at the PKA-pSP20 interface were influenced whether by the protonation state of the active site or by the metal cofactor used by the enzyme. Altogether, our results provide molecular aspects into the inhibitory mechanism of Ca²⁺ in PKA and suggest which is the most probable protonation state of the phosphorylated product at the active site.

What are the current limits on determination of protonation state using neutron macromolecular crystallography?

The rate of deposition of models determined by neutron diffraction, or a hybrid approach that combines X-ray and neutron diffraction, has increased in recent years. The benefit of neutron diffraction is that hydrogen atom (H) positions are detectable, allowing for the determination of protonation state and water molecule orientation. This study analyses all neutron models deposited in the Protein Data Bank to date, focusing on protonation state and properties of H (or deuterium, D) atoms as well as the details of water molecules. In particular, clashes and hydrogen bonds involving H or D atoms are investigated. As water molecules are typically the least reproducible part of a structural model, their positions in neutron models were compared to those in homologous high-resolution X-ray structures. For models determined by joint refinement against X-ray and neutron data, the water structure comparison was also carried out for models re-refined against the X-ray data alone. The homologues have generally fewer conserved water molecules where X-ray only was used and the positions of equivalent waters vary more than in the case of the hybrid X-ray model. As neutron diffraction data are generally less complete than X-ray data, the influence of neutron data completeness on nuclear density maps was also analyzed. We observe and discuss systematic map quality deterioration as result of data incompleteness.

Proton transfer and drug binding details revealed in neutron diffraction studies of wild-type and drug resistant HIV-1 protease

HIV-1 protease is an essential therapeutic target for the design and development of antiviral inhibitors to treat AIDS. We used room temperature neutron crystallography to accurately determine hydrogen atom positions in several protease complexes with clinical drugs, amprenavir and darunavir. Hydrogen bonding interactions were carefully mapped to provide an unprecedented picture of drug binding to the protease target. We demonstrate that hydrogen atom positions within the enzyme catalytic site can be altered by introducing drug resistant mutations and by protonating surface residues that trigger proton transfer reactions between the catalytic Asp residues and the hydroxyl group of darunavir. When protein perdeuteration is not feasible, we validate the use of initial H/D exchange with unfolded protein and partial deuteration in pure D2O with hydrogenous glycerol to maximize deuterium incorporation into the protein, with no detrimental effects on the growth of quality crystals suitable for neutron diffraction experiments.

A shared vision for macromolecular crystallography over the next five years

Macromolecular crystallography (MX) is the dominant means of determining the three-dimensional structures of biological macromolecules, but the method has reached a critical juncture. New diffraction-limited storage rings and upgrades to the existing sources will provide beamlines with higher flux and brilliance, and even the largest detectors can collect at rates of several hundred hertz. Electron cryomicroscopy is successfully competing for structural biologists' most exciting projects. As a result, formerly scarce beam time is becoming increasingly abundant, and beamlines must innovate to attract users and ensure continued funding. Here, we will show how data collection has changed over the preceding five years and how alternative methods have emerged. We then explore how MX at synchrotrons might develop over the next five years. We predict that, despite the continued dominance of rotation crystallography, applications previously considered niche or experimental, such as serial crystallography, pink-beam crystallography, and crystallography at energies above 25 keV and below 5 keV, will rise in prominence as beamlines specialize to offer users the best value. Most of these emerging methods will require new hardware and software. With these advances, MX will more efficiently provide the high-resolution structures needed for drug development. MX will also be able to address a broader range of questions than before and contribute to a deeper understanding of biological processes in the context of integrative structural biology.

How Electrostatic Coupling Enables Conformational Plasticity in a Tyrosine Kinase

Protein kinases are important cellular signaling molecules involved in cancer and a multitude of other diseases. It is well-known that inactive kinases display a remarkable conformational plasticity; however, the molecular mechanisms remain poorly understood. Conformational heterogeneity presents an opportunity but also a challenge in kinase drug discovery. The ability to predictively model various conformational states could accelerate selective inhibitor design. Here we performed a proton-coupled molecular dynamics study to explore the conformational landscape of a c-Src kinase. Starting from a completely inactive structure, the simulations captured all major types of conformational states without the use of a target structure, mutation, or bias. The simulations allowed us to test the experimental hypotheses regarding the mechanism of DFG flip, its coupling to the αC-helix movement, and the formation of regulatory spine. Perhaps the most significant finding is how key titratable residues, such as DFG-Asp, αC-Glu, and HRD-Asp, change protonation states dependent on the DFG, αC, and activation loop conformations. Our data offer direct evidence to support a long-standing hypothesis that protonation of Asp favors the DFG-out state and explain why DFG flip is also possible in simulations with deprotonated Asp. The simulations also revealed intermediate states, among which a unique DFG-out/α-C state formed as DFG-Asp is moved into a back pocket forming a salt bridge with catalytic Lys, which can be tested in selective inhibitor design. Our finding of how proton coupling enables the remarkable conformational plasticity may shift the paradigm of computational studies of kinases which assume fixed protonation states. Understanding proton-coupled conformational dynamics may hold a key to further innovation in kinase drug discovery.