Molecular Simulations Reveal an Unresolved Conformation of the Type IA Protein Kinase A Regulatory Subunit and Suggest Its Role in the cAMP Regulatory Mechanism.

We identify a previously unresolved, unrecognized, and highly stable conformation of the protein kinase A (PKA) regulatory subunit RIα. This conformation, which we term the "Flipback" structure, bridges conflicting characteristics in crystallographic structures and solution experiments of the PKA RIα heterotetramer. Our simulations reveal a hinge residue, G235, in the B/C helix that is conserved through all isoforms of RI. Brownian dynamics simulations suggest that the Flipback conformation plays a role in cAMP association to the A domain of the R subunit.

Development of Allosteric Small Molecule APOBEC3B Inhibitors from In Silico Screening.

APOBEC3B cytosine deaminase contributes to the mutational burdens of tumors, resulting in tumor progression and therapy resistance. Small molecule APOBEC3B inhibitors have potential to slow or mitigate these detrimental outcomes. Through molecular dynamics (MD) simulations and computational solvent mapping analysis, we identified a novel putative allosteric pocket on the C-terminal domain of APOBEC3B (A3Bctd), and virtually screened the ChemBridge Diversity Set (N~110,000) against both the active and potential allosteric sites. Selected high-scoring compounds were subsequently purchased, characterized for purity and composition, and tested in biochemical assays, which yielded 13 hit compounds. Orthogonal NMR assays verified binding to the target protein. Initial selectivity studies suggest these compounds preferentially target A3Bctd over related deaminase APOBEC3A (A3A), and MD simulations indicate this selectivity may be due to the steric repulsion from H56 that is unique to A3A. Taken together, our studies represent the first virtual screening effort against A3Bctd that has yielded candidate inhibitors suitable for further development.

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.

Multiscale computational modeling of the effects of 2'-deoxy-ATP on cardiac muscle calcium handling.

2'-Deoxy-ATP (dATP), a naturally occurring near analog of ATP, is a well-documented myosin activator that has been shown to increase contractile force, improve pump function, and enhance lusitropy in the heart. Calcium transients in cardiomyocytes with elevated levels of dATP show faster calcium decay compared with cardiomyocytes with basal levels of dATP, but the mechanisms behind this are unknown. Here, we design and utilize a multiscale computational modeling framework to test the hypothesis that dATP acts on the sarcoendoplasmic reticulum calcium-ATPase (SERCA) pump to accelerate calcium re-uptake into the sarcoplasmic reticulum during cardiac relaxation. Gaussian accelerated molecular dynamics simulations of human cardiac SERCA2A in the E1 apo, ATP-bound and dATP-bound states showed that dATP forms more stable contacts in the nucleotide binding pocket of SERCA and leads to increased closure of cytosolic domains. These structural changes ultimately lead to changes in calcium binding, which we assessed using Brownian dynamics simulations. We found that dATP increases calcium association rate constants to SERCA and that dATP binds to apo SERCA more rapidly than ATP. Using a compartmental ordinary differential equation model of human cardiomyocyte excitation-contraction coupling, we found that these increased association rate constants contributed to the accelerated rates of calcium transient decay observed experimentally. This study provides clear mechanistic evidence of enhancements in cardiac SERCA2A pump function due to interactions with dATP.

Erratum: Conserved patterns hidden within group A Streptococcus M protein hypervariability recognize human C4b-binding protein.

This corrects the article DOI: 10.1038/nmicrobiol.2016.155.

Conserved patterns hidden within group A Streptococcus M protein hypervariability recognize human C4b-binding protein.

No vaccine exists against group A Streptococcus (GAS), a leading cause of worldwide morbidity and mortality. A severe hurdle is the hypervariability of its major antigen, the M protein, with >200 different M types known. Neutralizing antibodies typically recognize M protein hypervariable regions (HVRs) and confer narrow protection. In stark contrast, human C4b-binding protein (C4BP), which is recruited to the GAS surface to block phagocytic killing, interacts with a remarkably large number of M protein HVRs (apparently ∼90%). Such broad recognition is rare, and we discovered a unique mechanism for this through the structure determination of four sequence-diverse M proteins in complexes with C4BP. The structures revealed a uniform and tolerant 'reading head' in C4BP, which detected conserved sequence patterns hidden within hypervariability. Our results open up possibilities for rational therapies that target the M-C4BP interaction, and also inform a path towards vaccine design.

Bridging scales through multiscale modeling: a case study on protein kinase A.

The goal of multiscale modeling in biology is to use structurally based physico-chemical models to integrate across temporal and spatial scales of biology and thereby improve mechanistic understanding of, for example, how a single mutation can alter organism-scale phenotypes. This approach may also inform therapeutic strategies or identify candidate drug targets that might otherwise have been overlooked. However, in many cases, it remains unclear how best to synthesize information obtained from various scales and analysis approaches, such as atomistic molecular models, Markov state models (MSM), subcellular network models, and whole cell models. In this paper, we use protein kinase A (PKA) activation as a case study to explore how computational methods that model different physical scales can complement each other and integrate into an improved multiscale representation of the biological mechanisms. Using measured crystal structures, we show how molecular dynamics (MD) simulations coupled with atomic-scale MSMs can provide conformations for Brownian dynamics (BD) simulations to feed transitional states and kinetic parameters into protein-scale MSMs. We discuss how milestoning can give reaction probabilities and forward-rate constants of cAMP association events by seamlessly integrating MD and BD simulation scales. These rate constants coupled with MSMs provide a robust representation of the free energy landscape, enabling access to kinetic, and thermodynamic parameters unavailable from current experimental data. These approaches have helped to illuminate the cooperative nature of PKA activation in response to distinct cAMP binding events. Collectively, this approach exemplifies a general strategy for multiscale model development that is applicable to a wide range of biological problems.