Kinesin-5 Allosteric Inhibitors Uncouple the Dynamics of Nucleotide,Microtubule, and Neck-Linker Binding Sites

Kinesin motor domains couple cycles of ATP hydrolysis to cycles of microtubule binding and conformational changes that result in directional force and movement on microtubules. The general principles of this mechanochemical coupling have been established; however, fundamental atomistic details of the underlying allosteric mechanisms remain unknown. This lack of knowledge hampers the development of new inhibitors and limits our understanding of how disease-associated mutations in distal sites can interfere with the fidelity of motor domain function. Here, we combine unbiased molecular-dynamics simulations, bioinformatics analysis, and mutational studies to elucidate the structural dynamic effects of nucleotide turnover and allosteric inhibition of the kinesin-5 motor. Multiple replica simulations of ATP-, ADP-, and inhibitor-bound states together with network analysis of correlated motions were used to create a dynamic protein structure network depicting the internal dynamic coordination of functional regions in each state. This analysis revealed the intervening residues involved in the dynamic coupling of nucleotide, microtubule, neck-linker, and inhibitor binding sites. The regions identified include the nucleotide binding switch regions, loop 5, loop 7, α4-α5-loop 13, α1, and β4-β6-β7. Also evident were nucleotide- and inhibitor-dependent shifts in the dynamic coupling paths linking functional sites. In particular, inhibitor binding to the loop 5 region affected β-sheet residues and α1, leading to a dynamic decoupling of nucleotide, microtubule, and neck-linker binding sites. Additional analyses of point mutations, including P131 (loop 5), Q78/I79 (α1), E166 (loop 7), and K272/I273 (β7) G325/G326 (loop 13), support their predicted role in mediating the dynamic coupling of distal functional surfaces. Collectively, our results and approach, which we make freely available to the community, provide a framework for explaining how binding events and point mutations can alter dynamic couplings that are critical for kinesin motor domain function.     

Computational optimization of the SARS-CoV-2 receptor-binding-motif affinity for human ACE2

The coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is responsible for the coronavirus disease 2019 pandemic, and the closely related SARS-CoV coronavirus enter cells by binding at the human angiotensin converting enzyme 2 (hACE2). The stronger hACE2 affinity of SARS-CoV-2 has been connected with its higher infectivity. In this work, we study hACE2 complexes with the receptor-binding domains (RBDs) of the human SARS-CoV-2 and human SARS-CoV viruses, using all-atom molecular dynamics simulations and computational protein design with a physics-based energy function. The molecular dynamics simulations identify charge-modifying substitutions between the CoV-2 and CoV RBDs, which either increase or decrease the hACE2 affinity of the SARS-CoV-2 RBD. The combined effect of these mutations is small, and the relative affinity is mainly determined by substitutions at residues in contact with hACE2. Many of these findings are in line and interpret recent experiments. Our computational protein design calculations redesign positions 455, 493, 494, and 501 of the SARS-CoV-2 receptor binding motif, which contact hACE2 in the complex and are important for ACE2 recognition. Sampling is enhanced by an adaptive importance sampling Monte Carlo method. Sequences with increased affinity replace CoV-2 glutamine by a negative residue at position 493; serine by a nonpolar or aromatic residue or an asparagine at position 494; and asparagine by valine or threonine at position 501. Substitutions at positions 455 and 501 have a smaller effect on affinity. Substitutions suggested by our design are seen in viral sequences encountered in other species, including bat and pangolin. Our results might be used to identify potential virus strains with higher human infectivity and assist in the design of peptide-based or peptidomimetic compounds with the potential to inhibit SARS-CoV-2 binding at hACE2.     

Kinesin-5 Allosteric Inhibitors Uncouple the Dynamics of Nucleotide,Microtubule, and Neck-Linker Binding Sites

Kinesin motor domains couple cycles of ATP hydrolysis to cycles of microtubule binding and conformational changes that result in directional force and movement on microtubules. The general principles of this mechanochemical coupling have been established; however, fundamental atomistic details of the underlying allosteric mechanisms remain unknown. This lack of knowledge hampers the development of new inhibitors and limits our understanding of how disease-associated mutations in distal sites can interfere with the fidelity of motor domain function. Here, we combine unbiased molecular-dynamics simulations, bioinformatics analysis, and mutational studies to elucidate the structural dynamic effects of nucleotide turnover and allosteric inhibition of the kinesin-5 motor. Multiple replica simulations of ATP-, ADP-, and inhibitor-bound states together with network analysis of correlated motions were used to create a dynamic protein structure network depicting the internal dynamic coordination of functional regions in each state. This analysis revealed the intervening residues involved in the dynamic coupling of nucleotide, microtubule, neck-linker, and inhibitor binding sites. The regions identified include the nucleotide binding switch regions, loop 5, loop 7, α4-α5-loop 13, α1, and β4-β6-β7. Also evident were nucleotide- and inhibitor-dependent shifts in the dynamic coupling paths linking functional sites. In particular, inhibitor binding to the loop 5 region affected β-sheet residues and α1, leading to a dynamic decoupling of nucleotide, microtubule, and neck-linker binding sites. Additional analyses of point mutations, including P131 (loop 5), Q78/I79 (α1), E166 (loop 7), and K272/I273 (β7) G325/G326 (loop 13), support their predicted role in mediating the dynamic coupling of distal functional surfaces. Collectively, our results and approach, which we make freely available to the community, provide a framework for explaining how binding events and point mutations can alter dynamic couplings that are critical for kinesin motor domain function.     

Dynamic coupling of residues within proteins as a mechanistic foundation of many enigmatic pathogenic missense variants

Many pathogenic missense mutations are found in protein positions that are neither well-conserved nor fall in any known functional domains. Consequently, we lack any mechanistic underpinning of dysfunction caused by such mutations. We explored the disruption of allosteric dynamic coupling between these positions and the known functional sites as a possible mechanism for pathogenesis. In this study, we present an analysis of 591 pathogenic missense variants in 144 human enzymes that suggests that allosteric dynamic coupling of mutated positions with known active sites is a plausible biophysical mechanism and evidence of their functional importance. We illustrate this mechanism in a case study of β-Glucocerebrosidase (GCase) in which a vast majority of 94 sites harboring Gaucher disease-associated missense variants are located some distance away from the active site. An analysis of the conformational dynamics of GCase suggests that mutations on these distal sites cause changes in the flexibility of active site residues despite their distance, indicating a dynamic communication network throughout the protein. The disruption of the long-distance dynamic coupling caused by missense mutations may provide a plausible general mechanistic explanation for biological dysfunction and disease.     

The Inherent Dynamics and Interaction Sites of the SARS-CoV-2 Nucleocapsid N-Terminal Region

The nucleocapsid protein is one of four structural proteins encoded by SARS-CoV-2 and plays a central role in packaging viral RNA and manipulating the host cell machinery, yet its dynamic behavior and promiscuity in nucleotide binding has made standard structural methods to address its atomic-resolution details difficult. To begin addressing the SARS-CoV-2 nucleocapsid protein interactions with both RNA and the host cell along with its dynamic behavior, we have specifically focused on the folded N-terminal domain (NTD) and its flanking regions using nuclear magnetic resonance solution studies. Studies performed here reveal a large repertoire of interactions, which includes a temperature-dependent self-association mediated by the disordered flanking regions that also serve as binding sites for host cell cyclophilin-A while nucleotide binding is largely mediated by the central NTD core. NMR studies that include relaxation experiments have revealed the complicated dynamic nature of this viral protein. Specifically, while much of the N-terminal core domain exhibits micro-millisecond motions, a central β-hairpin shows elevated inherent flexibility on the pico-nanosecond timescale and the serine/arginine-rich region of residues 176–209 undergoes multiple exchange phenomena. Collectively, these studies have begun to reveal the complexities of the nucleocapsid protein dynamics and its preferred interaction sites with its biological targets.     

Mapping the Structural and Dynamical Features of Kinesin Motor Domains

Kinesin motor proteins drive intracellular transport by coupling ATP hydrolysis to conformational changes that mediate directed movement along microtubules. Characterizing these distinct conformations and their interconversion mechanism is essential to determining an atomic-level model of kinesin action. Here we report a comprehensive principal component analysis of 114 experimental structures along with the results of conventional and accelerated molecular dynamics simulations that together map the structural dynamics of the kinesin motor domain. All experimental structures were found to reside in one of three distinct conformational clusters (ATP-like, ADP-like and Eg5 inhibitor-bound). These groups differ in the orientation of key functional elements, most notably the microtubule binding α4–α5, loop8 subdomain and α2b-β4-β6-β7 motor domain tip. Group membership was found not to correlate with the nature of the bound nucleotide in a given structure. However, groupings were coincident with distinct neck-linker orientations. Accelerated molecular dynamics simulations of ATP, ADP and nucleotide free Eg5 indicate that all three nucleotide states could sample the major crystallographically observed conformations. Differences in the dynamic coupling of distal sites were also evident. In multiple ATP bound simulations, the neck-linker, loop8 and the α4–α5 subdomain display correlated motions that are absent in ADP bound simulations. Further dissection of these couplings provides evidence for a network of dynamic communication between the active site, microtubule-binding interface and neck-linker via loop7 and loop13. Additional simulations indicate that the mutations G325A and G326A in loop13 reduce the flexibility of these regions and disrupt their couplings. Our combined results indicate that the reported ATP and ADP-like conformations of kinesin are intrinsically accessible regardless of nucleotide state and support a model where neck-linker docking leads to a tighter coupling of the microtubule and nucleotide binding regions. Furthermore, simulations highlight sites critical for large-scale conformational changes and the allosteric coupling between distal functional sites.     

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.     

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.     

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE): quantitative RNA structure analysis at single nucleotide resolution

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) interrogates local backbone flexibility in RNA at single-nucleotide resolution under diverse solution environments. Flexible RNA nucleotides preferentially sample local conformations that enhance the nucleophilic reactivity of 2'-hydroxyl groups toward electrophiles, such as N-methylisatoic anhydride (NMIA). Modified sites are detected as stops in an optimized primer extension reaction, followed by electrophoretic fragment separation. SHAPE chemistry scores local nucleotide flexibility at all four ribonucleotides in a single experiment and discriminates between base-paired versus unconstrained or flexible residues with a dynamic range of 20-fold or greater. Quantitative SHAPE reactivity information can be used to establish the secondary structure of an RNA, to improve the accuracy of structure prediction algorithms, to monitor structural differences between related RNAs or a single RNA in different states, and to detect ligand binding sites. SHAPE chemistry rarely needs significant optimization and requires two days to complete for an RNA of 100-200 nucleotides.     

FASN inhibitor TVB-3166 prevents S-acylation of the spike protein of human coronaviruses

The spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other coronaviruses mediates host cell entry and is S-acylated on multiple phylogenetically conserved cysteine residues. Multiple protein acyltransferase enzymes have been reported to post-translationally modify spike proteins; however, strategies to exploit this modification are lacking. Using resin-assisted capture MS, we demonstrate that the spike protein is S-acylated in SARS-CoV-2-infected human and monkey epithelial cells. We further show that increased abundance of the acyltransferase ZDHHC5 associates with increased S-acylation of the spike protein, whereas ZDHHC5 knockout cells had a 40% reduction in the incorporation of an alkynyl-palmitate using click chemistry detection. We also found that the S-acylation of the spike protein is not limited to palmitate, as clickable versions of myristate and stearate were also labelled the protein. Yet, we observed that ZDHHC5 was only modified when incubated with alkyne-palmitate, suggesting it has specificity for this acyl-CoA, and that other ZDHHC enzymes may use additional fatty acids to modify the spike protein. Since multiple ZDHHC isoforms may modify the spike protein, we also examined the ability of the FASN inhibitor TVB-3166 to prevent S-acylation of the spike proteins of SARS-CoV-2 and human CoV-229E. We show that treating cells with TVB-3166 inhibited S-acylation of expressed spike proteins and attenuated the ability of SARS-CoV-2 and human CoV-229E to spread in vitro. Our findings further substantiate the necessity of CoV spike protein S-acylation and demonstrate that de novo fatty acid synthesis is critical for the proper S-acylation of the spike protein.     

Multisite phosphorylation and binding alter conformational dynamics of the 4E-BP2 protein

Intrinsically disordered proteins (IDPs) play critical roles in regulatory protein interactions, but detailed structural/dynamic characterization of their ensembles remain challenging, both in isolation and when they form dynamic “fuzzy” complexes. Such is the case for mRNA cap-dependent translation initiation, which is regulated by the interaction of the predominantly folded eukaryotic initiation factor 4E (eIF4E) with the intrinsically disordered eIF4E binding proteins (4E-BPs) in a phosphorylation-dependent manner. Single-molecule Förster resonance energy transfer showed that the conformational changes of 4E-BP2 induced by binding to eIF4E are non-uniform along the sequence; while a central region containing both motifs that bind to eIF4E expands and becomes stiffer, the C-terminal region is less affected. Fluorescence anisotropy decay revealed a non-uniform segmental flexibility around six different labeling sites along the chain. Dynamic quenching of these fluorescent probes by intrinsic aromatic residues measured via fluorescence correlation spectroscopy report on transient intra- and inter-molecular contacts on nanosecond-to-microsecond timescales. Upon hyperphosphorylation, which induces folding of ～40 residues in 4E-BP2, the quenching rates decreased at most labeling sites. The chain dynamics around sites in the C-terminal region far away from the two binding motifs significantly increased upon binding to eIF4E, suggesting that this region is also involved in the highly dynamic 4E-BP2:eIF4E complex. Our time-resolved fluorescence data paint a sequence-level rigidity map of three states of 4E-BP2 differing in phosphorylation or binding status and distinguish regions that form contacts with eIF4E. This study adds complementary structural and dynamics information to recent studies of 4E-BP2, and it constitutes an important step toward a mechanistic understanding of this important IDP via integrative modeling.     

D-dimer and CoV-2 spike-immune complexes contribute to the production of PGE2 and proinflammatory cytokines in monocytes

An overreactive inflammatory response and coagulopathy are observed in patients with severe form of COVID-19. Since increased levels of D-dimer (DD) are associated with coagulopathy in COVID-19, we explored whether DD contributes to the aberrant cytokine responses. Here we show that treatment of healthy human monocytes with DD induced a dose dependent increase in production of pyrogenic mediator, Prostaglandin E2 (PGE2) and inflammatory cytokines, IL-6 and IL-8. The DD-induced PGE2 and inflammatory cytokines were enhanced significantly by co-treatment with immune complexes (IC) of SARS CoV-2 recombinant S protein or of pseudovirus containing SARS CoV-2 S protein (PVCoV-2) coated with spike-specific chimeric monoclonal antibody (MAb) containing mouse variable and human Fc regions. The production of PGE2 and cytokines in monocytes activated with DD and ICs was sensitive to the inhibitors of β2 integrin and FcγRIIa, and to the inhibitors of calcium signaling, Mitogen-Activated Protein Kinase (MAPK) pathway, and tyrosine-protein kinase. Importantly, strong increase in PGE2 and in IL-6/IL-8/IL-1β cytokines was observed in monocytes activated with DD in the presence of IC of PVCoV-2 coated with plasma from hospitalized COVID-19 patients but not from healthy donors. The IC of PVCoV-2 with convalescent plasma induced much lower levels of PGE2 and cytokines compared with plasma from hospitalized COVID-19 patients. PGE2 and IL-6/IL-8 cytokines produced in monocytes activated with plasma-containing IC, correlated well with the levels of spike binding antibodies and not with neutralizing antibody titers. Our study suggests that a combination of high levels of DD and high titers of spike-binding antibodies that can form IC with SARS CoV-2 viral particles might accelerate the inflammatory status of lung infiltrating monocytes leading to increased lung pathology in patients with severe form of COVID-19.     

NMR studies of internal dynamics of serine proteinase protein inhibitors: Binding region mobilities of intact and reactive-site hydrolyzed Cucurbita maxima trypsin inhibitor (CMTI)-III of the squash family and comparison with those of counterparts of CMTI-V of the potato I family.

Serine proteinase protein inhibitors follow the standard mechanism of inhibition (Laskowski M Jr, Kato I, 1980, Annu Rev Biochem 49:593-626), whereby an enzyme-catalyzed equilibrium between intact (I) and reactive-site hydrolyzed inhibitor (I*) is reached. The hydrolysis constant, Khyd, is defined as [I*]/[I]. Here, we explore the role of internal dynamics in the resynthesis of the scissile bond by comparing the internal mobility data of intact and cleaved inhibitors belonging to two different families. The inhibitors studied are recombinant Cucurbita maxima trypsin inhibitor III (rCMTI-III; Mr 3 kDa) of the squash family and rCMTI-V (Mr approximately 7 kDa) of the potato I family. These two inhibitors have different binding loop-scaffold interactions and different Khyd values--2.4 (CMTI-III) and 9 (CMTI-V)--at 25 degrees C. The reactive-site peptide bond (P1-P1') is that between Arg5 and Ile6 in CMTI-III, and that between Lys44 and Asp45 in CMTI-V. The order parameters (S2) of backbone NHs of uniformly 15N-labeled rCMTI-III and rCMTI-III* were determined from measurements of 15N spin-lattice and spin-spin relaxation rates, and [1H]-15N steady-state heteronuclear Overhauser effects, using the model-free formalism, and compared with the data reported previously for rCMTI-V and rCMTI-V*. The backbones of rCMTI-III [(S2) = 0.71] and rCMTI-III* [(S2) = 0.63] are more flexible than those of rCMTI-V [(S2) = 0.83] and rCMTI-V* [(S2) = 0.85]. The binding loop residues, P4-P1, in the two proteins show the following average order parameters: 0.57 (rCMTI-III) and 0.44 (rCMTI-III*); 0.70 (rCMTI-V) and 0.40 (rCMTI-V*). The P1'-P4' residues, on the other hand, are associated with (S2) values of 0.56 (rCMTI-III) and 0.47 (rCMTI-III*); and 0.73 (rCMTI-V) and 0.83 (rCMTI-V*). The newly formed C-terminal (Pn residues) gains a smaller magnitude of flexibility in rCMTI-III* due to the Cys3-Cys20 crosslink. In contrast, the newly formed N-terminal (Pn' residues) becomes more flexible only in rCMTI-III*, most likely due to lack of an interaction between the P1' residue and the scaffold in rCMTI-III. Thus, diminished flexibility gain of the Pn residues and, surprisingly, increased flexibility of the Pn' residues seem to facilitate the resynthesis of the P1-P1' bond, leading to a lower Khyd value.     

Calcium-induced flexibility changes in the troponin C-troponin I complex

The contraction of vertebrate striated muscle is modulated by Ca(2+) binding to the regulatory protein troponin C (TnC). Ca(2+) binding causes conformational changes in TnC which alter its interaction with the inhibitory protein troponin I (TnI), initiating the regulatory process. We have used the frequency domain method of fluorescence resonance energy transfer (FRET) to measure distances and distance distributions between specific sites in the TnC-TnI complex in the presence and absence of Ca(2+) or Mg(2+). Using sequences based on rabbit skeletal muscle proteins, we prepared functional, binary complexes of wild-type TnC and a TnI mutant which contains no Cys residues and a single Trp residue at position 106 within the TnI inhibitory region. We used TnI Trp-106 as the FRET donor, and we introduced energy acceptor groups into TnC by labeling at Met-25 with dansyl aziridine or at Cys-98 with N-(iodoacetyl)-N'-(1-sulfo-5-naphthyl)ethylenediamine. Our distance distribution measurements indicate that the TnC-TnI complex is relatively rigid in the absence of Ca(2+), but becomes much more flexible when Ca(2+) binds to regulatory sites in TnC. This increased flexibility may be propagated to the whole thin filament, helping to release the inhibition of actomyosin ATPase activity and allowing the muscle to contract. This is the first report of distance distributions between TnC and TnI in their binary complex.     

The effects of Ca2+ binding on the dynamic properties of a designed Ca2+-binding protein

The effects of Ca(2+) binding on the dynamic properties of Ca(2+)-binding proteins are important in Ca(2+) signaling. To understand the role of Ca(2+) binding, we have successfully designed a Ca(2+)-binding site in the domain 1 of rat CD2 (denoted as Ca.CD2) with the desired structure and retained function. In this study, the backbone dynamic properties of Ca.CD2 have been investigated using (15)N spin relaxation NMR spectroscopy to reveal the effect of Ca(2+) binding on the global and local dynamic properties without the complications of multiple interactive Ca(2+) binding and global conformational change. Like rat CD2 (rCD2) and human CD2 (hCD2), residues involved in the recognition of the target molecule CD48 exhibit high flexibility. Mutations N15D and N17D that introduce the Ca(2+) ligands increase the flexibility of the neighboring residues. Ca(2+)-induced local dynamic changes occur mainly at the residues proximate to the Ca(2+)-binding pocket or the residues in loop regions. The beta-strand B of Ca.CD2 that provides two Asp for the Ca(2+) undergoes an S(2) decrease upon the Ca(2+) binding, while the DE-loop that provides one Asn and one Asp undergoes an S(2) increase. Our study suggests that Ca(2+) binding has a differential effect on the rigidity of the residues depending on their flexibility and location within the secondary structure.     

Dynamic Coupling and Allosteric Networks in the α Subunit of Heterotrimeric G Proteins

G protein α subunits cycle between active and inactive conformations to regulate a multitude of intracellular signaling cascades. Important structural transitions occurring during this cycle have been characterized from extensive crystallographic studies. However, the link between observed conformations and the allosteric regulation of binding events at distal sites critical for signaling through G proteins remain unclear. Here we describe molecular dynamics simulations, bioinformatics analysis, and experimental mutagenesis that identifies residues involved in mediating the allosteric coupling of receptor, nucleotide, and helical domain interfaces of Gα_i. Most notably, we predict and characterize novel allosteric decoupling mutants, which display enhanced helical domain opening, increased rates of nucleotide exchange, and constitutive activity in the absence of receptor activation. Collectively, our results provide a framework for explaining how binding events and mutations can alter internal dynamic couplings critical for G protein function.     

Oncogenic mutations on Rac1 affect global intrinsic dynamics underlying GTP and PAK1 binding

Rac1 is a small member of the Rho GTPase family. One of the most important downstream effectors of Rac1 is a serine/threonine kinase, p21-activated kinase 1 (PAK1). Mutational activation of PAK1 by Rac1 has oncogenic signaling effects. Here, although we focus on Rac1-PAK1 interaction by atomic-force-microscopy-based single-molecule force spectroscopy experiments, we explore the effect of active mutations on the intrinsic dynamics and binding interactions of Rac1 by Gaussian network model analysis and molecular dynamics simulations. We observe that Rac1 oncogenic mutations are at the hinges of three global modes of motion, suggesting the mechanical changes as potential markers of oncogenicity. Indeed, the dissociation of wild-type Rac1-PAK1 complex shows two distinct unbinding dynamic states that are reduced to one with constitutively active Q61L and oncogenic Y72C mutant Rac1, as revealed by single-molecule force spectroscopy experiments. Q61L and Y72C mutations change the mechanics of the Rac1-PAK1 complex by increasing the elasticity of the protein and slowing down the transition to the unbound state. On the other hand, Rac1’s intrinsic dynamics reveal more flexible GTP and PAK1-binding residues on switches I and II with Q61L, Y72C, oncogenic P29S and Q61R, and negative T17N mutations. The cooperativity in the fluctuations of GTP-binding sites around the p-loop and switch I decreases in all mutants, mostly in Q61L, whereas some PAK1-binding residues display enhanced coupling with GTP-binding sites in Q61L and Y72C and within each other in P29S. The predicted binding free energies of the modeled Rac1-PAK1 complexes show that the change in the dynamic behavior likely means a more favorable PAK1 interaction. Overall, these findings suggest that the active mutations affect intrinsic functional dynamic events and alter the mechanics underlying the binding of Rac1 to GTP and upstream and downstream partners including PAK1.     

Probing the ATP-induced conformational flexibility of the PcrA helicase protein using molecular dynamics simulation

Helicases are enzymes that unwind double-stranded DNA (dsDNA) into its single-stranded components. It is important to understand the binding and unbinding of ATP from the active sites of helicases, as this knowledge can be used to elucidate the functionality of helicases during the unwinding of dsDNA. In this work, we investigated the unbinding of ATP and its effect on the active-site residues of the helicase PcrA using molecular dynamic simulations. To mimic the unbinding process of ATP from the active site of the helicase, we simulated the application of an external force that pulls ATP from the active site and computed the free-energy change during this process. We estimated an energy cost of ~85 kJ/mol for the transformation of the helicase from the ATP-bound state (1QHH) to the ATP-free state (1PJR). Unbinding led to conformational changes in the residues of the protein at the active site. Some of the residues at the ATP-binding site were significantly reoriented when the ATP was pulled. We observed a clear competition between reorientation of the residues and energy stabilization by hydrogen bonds between the ATP and active-site residues. We also checked the flexibility of the PcrA protein using a principal component analysis of domain motion. We found that the ATP-free state of the helicase is more flexible than the ATP-bound state.     

Altered flexibility in the substrate-binding site of related native and engineered high-alkaline Bacillus subtilisins.

High-alkaline serine proteases have been successfully applied as protein degrading components of detergent formulations and are subject to extensive protein engineering efforts to improve their stability and performance. Dynamics has been suggested to play an important role in determining enzyme activity and specificity and it is therefore of interest to establish how local changes in internal mobility affect protein stability, specificity and performance. Here we present the dynamic properties of the 269 residue serine proteases subtilisin PB92 (Maxacal(TM)) and subtilisin BLS (Savinase(TM)), secreted by Bacillus lentus, and an engineered quadruple variant, DSAI, that has improved washing performance. T1, T2 and heteronuclear NOE measurements of the 15N nuclei indicate that for all three proteins the majority of the backbone is very rigid, with only a limited number of residues being involved in local mobility. Many of the residues that constitute the S1 and S4 pockets, determining substrate specificity, are flexible in solution. In contrast, the backbone amides of the residues that constitute the catalytic triad do not exhibit any motion. Subtilisins PB92, BLS and DSAI demonstrate similar but not identical NMR relaxation rates. A detailed analysis of local flexibility indicates that the motion of residues Thr143 and Ala194 becomes more restricted in subtilisin BLS and DSAI. Noteworthy, the loop regions involved in substrate binding become more structured in the engineered variant as compared with the two native proteases, suggesting a relation between altered dynamics and performance. Similar conclusions have been established by X-ray crystallograpic methods, as shown in the accompanying paper.