Proteins are key components in many processes in living cells, and physical interactions with other proteins and nucleic acids often form key parts of their functions. In many cases, large flexibility of proteins as they interact is key to their function. To understand the mechanisms of these processes, it is necessary to consider the 3D structures of such protein complexes. When such structures are not yet experimentally determined, protein docking has long been present to computationally generate useful structure models. However, protein docking has long had the limitation that the consideration of flexibility is usually limited to very small movements or very small structures. Methods have been developed which handle minor flexibility via normal mode or other structure sampling, but new methods are required to model ordered proteins which undergo large-scale conformational changes to elucidate their function at the molecular level. Here, we present Flex-LZerD, a framework for docking such complexes. Via partial assembly multidomain docking and an iterative normal mode analysis admitting curvilinear motions, we demonstrate the ability to model the assembly of a variety of protein-protein and protein-nucleic acid complexes.

Eukaryotic post-translational arginylation, mediated by the family of enzymes known as the arginyltransferases (ATE1s), is an important post-translational modification that can alter protein function and even dictate cellular protein half-life. Multiple major biological pathways are linked to the fidelity of this process, including neural and cardiovascular developments, cell division, and even the stress response. Despite this significance, the structural, mechanistic, and regulatory mechanisms that govern ATE1 function remain enigmatic. To that end, we have used X-ray crystallography to solve the crystal structure of ATE1 from the model organism Saccharomyces cerevisiae ATE1 (ScATE1) in the apo form. The three-dimensional structure of ScATE1 reveals a bilobed protein containing a GCN5-related N-acetyltransferase (GNAT) fold, and this crystalline behavior is faithfully recapitulated in solution based on size-exclusion chromatography-coupled small angle X-ray scattering (SEC-SAXS) analyses and cryo-EM 2D class averaging. Structural superpositions and electrostatic analyses point to this domain and its domain-domain interface as the location of catalytic activity and tRNA binding, and these comparisons strongly suggest a mechanism for post-translational arginylation. Additionally, our structure reveals that the N-terminal domain, which we have previously shown to bind a regulatory [Fe-S] cluster, is dynamic and disordered in the absence of metal bound in this location, hinting at the regulatory influence of this region. When taken together, these insights bring us closer to answering pressing questions regarding the molecular-level mechanism of eukaryotic post-translational arginylation.

Bacterial cell division begins with the formation of the Z-ring via polymerization of FtsZ and the localization of Z-ring beneath the inner membrane through membrane anchors. In Mycobacterium tuberculosis (Mtb), SepF is one such membrane anchor, but our understanding of the underlying mechanism is very limited. Here we used molecular dynamics simulations to characterize how SepF itself, a water-soluble protein, tethers to acidic membranes that mimic the Mtb inner membrane. In addition to an amphipathic helix (residues 1-12) at the N-terminus, membrane binding also occurs through two stretches of positively charged residues (Arg27-Arg37 and Arg95-Arg107) in the long linker preceding the FtsZ-binding core domain (residues 128-218). The additional interactions via the disordered linker stabilize the membrane tethering of SepF, and keep the core domain of SepF and hence the attached Z-ring close to the membrane. The resulting membrane proximity of the Z-ring in turn enables its interactions with and thus recruitment of two membrane proteins, FtsW and CrgA, at the late stage of cell division.

Many vertebrate genes generate multiple transcript variants that may encode functionally distinct protein isoforms, but the transcriptomes of various developmental stages are poorly defined. Identifying the transcriptome and its regulation during the normal developmental process is the key to deciphering the developmental stage-specific functions of genes. Here we presented a systematic assessment of the temporal alternative splicing (AS) events during the critical development stages to capture the dynamic gene expression changes and AS in zebrafish. An unexpected transcriptome complexity generated by AS was observed during zebrafish development. The patterns of AS events varied substantially among developmental stages despite the similarities in the total proportion of AS genes. We further found that AS afforded substantial functional diversification of genes through the generation of stage-specific AS events from broadly protein-coding genes as an essential developmental regulatory mechanism. Skipped exon (SE) showed the strongest signals among developmental AS (devAS), suggesting that devAS events generated by SE may be necessary for the normal development of zebrafish. Most developmental genes regulated by AS mechanisms were not modulated in terms of their overall expression levels, indicating that AS shaped the transcriptome independently from transcriptional regulation during development. 128-cell stage was a critical stage for gene transcription during embryonic development. Splicing factors as an essential developmental regulator underwent AS in the potential autoregulatory feedback loop and expressed multiple isoforms. Thus, zebrafish development was shaped by an interplay of programs controlling gene expression levels and AS. Overall, we provided a global view of developmental patterns of AS during zebrafish development and revealed that AS transitions were the crucial regulatory component of zebrafish embryonic development.

Recent advances in interrogating RNA folding dynamics have shown the classical model of RNA folding to be incomplete. Here, we pose three prominent questions for the field that are at the forefront of our understanding of the importance of RNA folding dynamics for RNA function. The first centers on the most appropriate biophysical framework to describe changes to the RNA folding energy landscape that a growing RNA chain encounters during transcriptional elongation. The second focuses on the potential ubiquity of strand displacement - a process by which RNA can rapidly change conformations - and how this process may be generally present in broad classes of seemingly different RNAs. The third raises questions about the potential importance and roles of cellular protein factors in RNA conformational switching. Answers to these questions will greatly improve our fundamental knowledge of RNA folding and function, drive biotechnological advances that utilize engineered RNAs, and potentially point to new areas of biology yet to be discovered.

Conformational dynamics are essential to macromolecular function. This is certainly true of RNA, whose ability to undergo programmed conformational dynamics is essential to create and regulate complex biological processes. However, methods to easily and simultaneously interrogate both the structure and conformational dynamics of fully functional RNAs in isolation and in complex with proteins have not historically been available. Due to its ability to image and classify single particles, cryogenic electron microscopy (cryo-EM) has the potential to address this gap and may be particularly amenable to exploring structural dynamics within the three-dimensional folds of biologically active RNAs. We discuss the possibilities and current limitations of applying cryo-EM to simultaneously study RNA structure and conformational dynamics, and present one example that illustrates this (as of yet) not fully realized potential.

RNA switches are versatile tools in synthetic biology for sensing and regulation applications. The discoveries of RNA-mediated translational and transcriptional control have facilitated the development of complex de novo designs of RNA switches. Specifically, RNA toehold-mediated switches, in which binding to the toehold sensing domain controls the transition between switch states via strand displacement, have been extensively adapted for coupling systems responses to specific trans-RNA inputs. This review highlights some of the challenges associated with applying these switches for native RNA detection in vivo, including transferability between organisms. The applicability and design considerations of toehold-mediated switches are discussed by highlighting twelve recently developed switch designs. This review finishes with future perspectives to address current gaps in the field, particularly regarding the power of structural prediction algorithms for improved in vivo functionality of RNA switches.

RNA structure probing experiments have emerged over the last decade as a straightforward way to determine the structure of RNA molecules in a number of different contexts. Although powerful, the ability of RNA to dynamically interconvert between, and to simultaneously populate, alternative structural configurations, poses a nontrivial challenge to the interpretation of data derived from these experiments. Recent efforts aimed at developing computational methods for the reconstruction of coexisting alternative RNA conformations from structure probing data are paving the way to the study of RNA structure ensembles, even in the context of living cells. In this review, we critically discuss these methods, their limitations and possible future improvements.

Synthetic riboswitches are powerful tools in synthetic biology in which sensing and execution are consolidated in a single RNA molecule. By using SELEX to select aptamers in vitro, synthetic riboswitches can in theory be engineered against any ligand of choice. Surprisingly, very few in vitro selected aptamers have been used for the engineering of synthetic riboswitches. In-depth studies of these aptamers suggest that the key characteristics of such regulatory active RNAs are their structural switching abilities and their binding dynamics. Conventional SELEX approaches seem to be inadequate to select for these characteristics, which may explain the lack of in vitro selected aptamers suited for engineering of synthetic riboswitches. In this review, we explore the functional principles of synthetic riboswitches, identify key characteristics of regulatory active in vitro selected aptamers and integrate these findings in context with available in vitro selection methods. Based on these insights, we propose to use a combination of capture-SELEX and subsequent functional screening for a more successful in vitro selection of aptamers that can be applied for the engineering of synthetic riboswitches.

The type III secretion system (T3SS) is indispensable for successful host cell infection by many Gram-negative pathogens. The molecular syringe delivers effector proteins that suppress the host immune response. Synthesis of T3SS components in Yersinia pseudotuberculosis relies on host body temperature, which induces the RNA thermometer (RNAT)-controlled translation of lcrF coding for a virulence master regulator that activates transcription of the T3SS regulon. The assembly of the secretion machinery follows a strict coordinated succession referred to as outside-in assembly, in which the membrane ring complex and the export apparatus represent the nucleation points. Two components essential for the initial assembly are YscJ and YscT. While YscJ connects the membrane ring complex with the export apparatus in the inner membrane, YscT is required for a functional export apparatus. Previous transcriptome-wide RNA structuromics data suggested the presence of unique intercistronic RNATs upstream of yscJ and yscT. Here, we show by reporter gene fusions that both upstream regions confer translational control. Moreover, we demonstrate the temperature-induced opening of the Shine-Dalgarno region, which facilitates ribosome binding, by in vitro structure probing and toeprinting methods. Rationally designed thermostable RNAT variants of the yscJ and yscT thermometers confirmed their physiological relevance with respect to T3SS assembly and host infection. Since we have shown in a recent study that YopN, the gatekeeper of type III secretion, also is under RNAT control, it appears that the synthesis, assembly and functionality of the Yersinia T3S machinery is coordinated by RNA-based temperature sensors at multiple levels.

Riboswitches are an outstanding example of genetic regulation mediated by RNA conformational switching. In these non-coding RNA elements, the occupancy status of a ligand-binding domain governs the mRNA's decision to form one of two mutually exclusive structures in the downstream expression platform. Temporal constraints upon the function of many riboswitches, requiring folding of complex architectures and conformational switching in a limited co-transcriptional timeframe, make them ideal model systems for studying these processes. In this review, we focus on the mechanism of ligand-directed conformational changes in one of the most widely distributed riboswitches in bacteria: the cobalamin family. We describe the architectural features of cobalamin riboswitches whose structures have been determined by x-ray crystallography, which suggest a direct physical role of cobalamin in effecting the regulatory switch. Next, we discuss a series of experimental approaches applied to several model cobalamin riboswitches that interrogate these structural models. As folding is central to riboswitch function, we consider the differences in folding landscapes experienced by RNAs that are produced in vitro and those that are allowed to fold co-transcriptionally. Finally, we highlight a set of studies that reveal the difficulties of studying cobalamin riboswitches outside the context of transcription and that co-transcriptional approaches are essential for developing a more accurate picture of their structure-function relationships in these switches. This understanding will be essential for future advancements in the use of small-molecule guided RNA switches in a range of applications such as biosensors, RNA imaging tools, and nucleic acid-based therapies.

MicroRNAs (miRNAs) are small, noncoding RNAs that mediate post-transcriptional downregulation of specific target genes. These transcripts are the products of a two-step processing pathway; primary miRNAs (pri-miRNAs) are processed by Drosha into individual precursor miRNA (pre-miRNA) hairpins, which are subsequently processed by Dicer into mature miRNAs. Single nucleotide polymorphisms (SNPs) that occur in pri-miRNAs, pre-miRNAs and mature miRNAs have been shown to affect the processing of specific target genes by modulating Drosha and Dicer processing or interactions with RNA binding proteins (RBPs). Using NMR and single-molecule optical tweezer experiments, we have investigated the conformational effects of a cancer-linked G/A mutation in the terminal loop of pri-miR-30c RNA, and how this influences binding by the SRSF3 and hnRNP A1 RBPs, which are implicated in its processing. Our results reveal that the wildtype and G/A variant pri-miR-30c RNAs adopt very similar elongated stem-loop structures, both of which are bound by SRSF3. However, while both wildtype and G/A pri-miR-30c RNAs can form dimeric kissing hairpin structures, the G to A mutation results in partial destabilization of the dimer in the variant transcript. This promotes recognition and binding by hnRNP A1, an RBP that enhances pri-miR-30c processing. Our data provide structural insight into the conformational effects of a G/A mutation in pri-miR-30c RNA and how this could affect processing and promote cancer.

Many heat shock genes in bacteria are regulated through a class of temperature-sensitive stem-loop (SL) RNAs called RNA thermometers (RNATs). One of the most widely studied RNATs is the Repression Of heat Shock Expression (ROSE) element associated with expression of heat shock proteins. Located in the 5'UTR, the RNAT contains one to three auxiliary hairpins upstream of it. Herein, we address roles of these upstream SLs in the folding and function of an RNAT. Bradyrhizobium japonicum is a nitrogen-fixing bacterium that experiences a wide range of temperatures in the soil and contains ROSE elements, each having multiple upstream SLs. The 5'UTR of the messenger (mRNA) for heat shock protein A (hspA) in B. japonicum has an intricate secondary structure containing three SLs upstream of the RNAT SL. While structure-function studies of the hspA RNAT itself have been reported, it has been unclear if these auxiliary SLs contribute to the temperature-sensing function of the ROSE elements. Herein, we show that the full length (FL) sequence has several melting transitions indicating that the ROSE element unfolds in a non-two-state manner. The upstream SLs are more stable than the RNAT itself, and a variant with disrupted base pairing in the SL immediately upstream of the RNAT has little influence on the melting of the RNAT. On the basis of these results and modeling of the co-transcriptional folding of the ROSE element, we propose that the upstream SLs function to stabilize the transcript and aid proper folding and dynamics of the RNAT.

The Sm protein Hfq chaperones small non-coding RNAs (sRNAs) in bacteria, facilitating sRNA regulation of target mRNAs. Hfq acts in part by remodeling the sRNA and mRNA structures, yet the basis for this remodeling activity is not understood. To understand how Hfq remodels RNA, we used single-molecule Förster resonance energy transfer (smFRET) to monitor conformational changes in OxyS sRNA upon Hfq binding. The results show that E. coli Hfq first compacts OxyS, bringing its 5' and 3 ends together. Next, Hfq destabilizes an internal stem-loop in OxyS, allowing the RNA to adopt a more open conformation that is stabilized by a conserved arginine on the rim of Hfq. The frequency of transitions between compact and open conformations depend on interactions with Hfqs flexible C-terminal domain (CTD), being more rapid when the CTD is deleted, and slower when OxyS is bound to Caulobacter crescentus Hfq, which has a shorter and more stable CTD than E. coli Hfq. We propose that the CTDs gate transitions between OxyS conformations that are stabilized by interaction with one or more arginines. These results suggest a general model for how basic residues and intrinsically disordered regions of RNA chaperones act together to refold RNA.

The miR-17 ∼ 92a polycistron, also known as oncomiR-1, is commonly overexpressed in multiple cancers and has several oncogenic properties. OncomiR-1 encodes six constituent microRNAs (miRs), each enzymatically processed with different efficiencies. However, the structural mechanism that regulates this differential processing remains unclear. Chemical probing of oncomiR-1 revealed that the Drosha cleavage sites of pri-miR-92a are sequestered in a four-way junction. NPSL2, an independent stem loop element, is positioned just upstream of pri-miR-92a and sequesters a crucial part of the sequence that constitutes the basal helix of pri-miR-92a. Disruption of the NPSL2 hairpin structure could promote the formation of a pri-miR-92a structure that is primed for processing by Drosha. Thus, NPSL2 is predicted to function as a structural switch, regulating pri-miR-92a processing. Here, we determined the solution structure of NPSL2 using solution NMR spectroscopy. This is the first high-resolution structure of an oncomiR-1 element. NPSL2 adopts a hairpin structure with a large, but highly structured, apical and internal loops. The 10-bp apical loop contains a pH-sensitive A+·C mismatch. Additionally, several adenosines within the apical and internal loops have elevated pKa values. The protonation of these adenosines can stabilize the NPSL2 structure through electrostatic interactions. Our study provides fundamental insights into the secondary and tertiary structure of an important RNA hairpin proposed to regulate miR biogenesis.

Messenger RNA regulatory elements, such as riboswitches, can display a high degree of flexibility. By characterizing their energy landscapes, and corresponding distributions of 3D configurations, structure-function relationships can be elucidated. Molecular dynamics simulation with enhanced sampling is an important strategy used to computationally access free energy landscapes characterizing the accessible 3D conformations of RNAs. While tertiary contacts are thought to play important roles in RNA dynamics, it is difficult, in explicit solvent, to sample the formation and breakage of tertiary contacts, such as helix-helix interactions, pseudoknot interactions, and junction interactions, while maintaining intact secondary structure elements. To this end, we extend previously developed collective variables and metadynamics efforts, to establish a simple metadynamics protocol, which utilizes only one collective variable, based on multiple tertiary contacts, to characterize the underlying free energy landscape of any RNA molecule. We develop a modified collective variable, the tertiary contacts distance (QTC), which can probe the formation and breakage of all or selectively chosen tertiary contacts of the RNA. The SAM-I riboswitch in the presence of three ionic and substrate conditions was investigated and validated against the structure ensemble previously generated using SAXS experiments. This efficient and easy to implement all-atom MD simulation based approach incorporating metadynamics to study RNA conformational dynamics can also be transferred to any other type of biomolecule.

Human hepatitis B virus (HBV) replication is initiated by the binding of the viral polymerase (P) to epsilon (ε), an ≈85-nucleotide (nt) cis-acting regulatory stem-loop RNA located at the 5'-end of the pre-genomic RNA (pgRNA). This interaction triggers P and pgRNA packaging and protein-primed reverse transcription and is therefore an attractive therapeutic target. Our recent nuclear magnetic resonance (NMR) structure of ε provides a useful starting point toward a detailed understanding of HBV replication, and hints at the functional importance of ε dynamics. Here, we present a detailed description of ε motions on the ps to ns and μs to ms time scales by NMR spin relaxation and relaxation dispersion, respectively. We also carried out molecular dynamics simulations to provide additional insight into ε conformational dynamics. These data outline a series of complex motions on multiple time scales within ε. Moreover, these motions occur in mostly conserved nucleotides from structural regions (i.e., priming loop, pseudo-triloop, and U43 bulge) that biochemical and mutational studies have shown to be essential for P binding, P-pgRNA packaging, protein-priming, and DNA synthesis. Taken together, our work implicates RNA dynamics as an integral feature that governs HBV replication.

RNA folding free energy change parameters are widely used to predict RNA secondary structure and to design RNA sequences. These parameters include terms for the folding free energies of helices and loops. Although the full set of parameters has only been traditionally available for the four common bases and backbone, it is well known that covalent modifications of nucleotides are widespread in natural RNAs. Covalent modifications are also widely used in engineered sequences. We recently derived a full set of nearest neighbor terms for RNA that includes N6-methyladenosine (m6A). In this work, we test the model using 98 optical melting experiments, matching duplexes with or without N6-methylation of A. Most experiments place RRACH, the consensus site of N6-methylation, in a variety of contexts, including helices, bulge loops, internal loops, dangling ends, and terminal mismatches. For matched sets of experiments that include either A or m6A in the same context, we find that the parameters for m6A are as accurate as those for A. Across all experiments, the root mean squared deviation between estimated and experimental free energy changes is 0.67 kcal/mol. We used the new experimental data to refine the set of nearest neighbor parameter terms for m6A. These parameters enable prediction of RNA secondary structures including m6A, which can be used to model how N6-methylation of A affects RNA structure.

Translational riboswitches are bacterial gene regulatory elements found in the 5'-untranslated region of mRNAs. They operate through a conformational refolding reaction that is triggered by a concentration change of a modulating small molecular ligand. The translation initiation region (TIR) is either released from or incorporated into base pairing interactions through the conformational switch. Hence, initiation of translation is regulated by the accessibility of the Shine-Dalgarno sequence and start codon. Interaction with the 30S ribosome is indispensable for the structural switch between functional OFF and ON states. However, on a molecular level it is still not fully resolved how the ribosome is accommodated near or at the translation initiation region in the context of translational riboswitches. The standby model of translation initiation postulates a binding site where the mRNA enters the ribosome and where it resides until the initiation site becomes unstructured and accessible. We here investigated the adenine-sensing riboswitch from Vibrio vulnificus. By application of a 19F labelling strategy for NMR spectroscopy that utilizes ligation techniques to synthesize differentially 19F labelled riboswitch molecules we show that nucleotides directly downstream of the riboswitch domain are first involved in productive interaction with the 30S ribosomal subunit. Upon the concerted action of ligand and the ribosomal protein rS1 the TIR becomes available and subsequently the 30S ribosome can slide towards the TIR. It will be interesting to see whether this is a general feature in translational riboswitches or if riboswitches exist where this region is structured and represent yet another layer of regulation.

Alternative splicing of the HIV transcriptome is controlled through cis regulatory elements functioning as enhancers or silencers depending on their context and the type of host RNA binding proteins they recruit. Splice site acceptor A3 (ssA3) is one of the least used acceptor sites in the HIV transcriptome and its activity determines the levels of tat mRNA. Splice acceptor 3 is regulated by a combination of cis regulatory sequences, auxiliary splicing factors, and presumably RNA structure. The mechanisms by which these multiple regulatory components coordinate to determine the frequency in which ssA3 is utilized is poorly understood. By NMR spectroscopy and phylogenetic analysis, we show that the ssA3 regulatory locus is conformationally heterogeneous and that the sequences that encompass the locus are conserved across most HIV isolates. Despite the conformational heterogeneity, the major stem loop (A3SL1) observed in vitro folds to base pair the Polypyrimdine Tract (PPyT) to the Exon Splicing Silencer 2p (ESS2p) element and to a conserved downstream linker. The 3D structure as determined by NMR spectroscopy further reveals that the A3 consensus cleavage site is embedded within a unique stereochemical environment within the apical loop, where it is surrounded by alternating base-base interactions. Despite being described as a receptor for hnRNP H, the ESS2p element is sequestered by base pairing to the 3' end of the PPyT and within this context it cannot form a stable complex with hnRNP H. By comparison, hnRNP A1 directly binds to the A3 consensus cleavage site located within the apical loop, suggesting that it can directly modulate U2AF assembly. Sequence mutations designed to destabilize the PPyT:ESS2p helix results in an increase usage of ssA3 within HIV-infected cells, consistent with the PPyT becoming more accessible for U2AF recognition. Additional mutations introduced into the downstream ESS2 element synergize with ESS2p to cause further increases in ssA3 usage. When taken together, our work provides a unifying picture by which cis regulatory sequences, splicing auxiliary factors and RNA structure cooperate to provide stringent control over ssA3. We describe this as the pair-and-lock mechanism to restrict access of the PPyT, and posit that it operates to regulate a subset of the heterogenous structures encompassing the ssA3 regulatory locus.

Allosteric regulation of proteins has been utilized to study various aspects of cell signaling, from unicellular events to organism-wide phenotypes. However, traditional methods of allosteric regulation, such as constitutively active mutants and inhibitors, lack tight spatiotemporal control. This often leads to unintended signaling consequences that interfere with data interpretation. To overcome these obstacles, researchers employed protein engineering approaches that enable tight control of protein function through allosteric mechanisms. These methods provide high specificity as well as spatial and temporal precision in regulation of protein activity in vitro and in vivo. In this review, we focus on the recent advancements in engineered allosteric regulation and discuss the various bioengineered allosteric techniques available now, from chimeric GPCRs to chemogenetic and optogenetic switches. We highlight the benefits and pitfalls of each of these techniques as well as areas in which future improvements can be made. Additionally, we provide a brief discussion on implementation of engineered allosteric regulation approaches, demonstrating that these tools can shed light on elusive biological events and have the potential to be utilized in precision medicine.

The cAMP- and cGMP-dependent protein kinases (PKA and PKG) are canonically activated by the corresponding cyclic nucleotides. However, both systems are also sensitive to a wide range of non-canonical allosteric effectors, such as reactive oxygen species, which induce the formation of regulatory inter- and intra-molecular disulfide bridges, and disease-related mutations (DRMs). Here, we present a combined analysis of representative non-canonical allosteric effectors for PKA and PKG, and we identify common molecular mechanisms underlying non-canonical allostery in these kinases, from shifts in dynamical regulatory equilibria to modulation of inter-protomer interactions. In addition, mutations may also drive oligomerization beyond dimerization, and possibly phase transitions, causing loss of kinase inhibitory function and amplifying the allosteric effects of DRMs. Hence non-canonical allosteric stimuli often result in constitutive kinase activation underlying either physiological control of downstream signaling pathways or pathological outcomes, from aortic aneurisms to cancer predisposition. Overall, PKA and PKG emerge as "pan-sensors" going well beyond canonical cyclic nucleotide activation, revealing their versatile roles as central signaling hubs.

Here, we discuss the principles of allosteric activating mutations, propagation downstream of the signals that they prompt, and allosteric drugs, with examples from the Ras signaling network. We focus on Abl kinase where mutations shift the landscape toward the active, imatinib binding-incompetent conformation, likely resulting in the high affinity ATP outcompeting drug binding. Recent pharmacological innovation extends to allosteric inhibitor (GNF-5)-linked PROTAC, targeting Bcr-Abl1 myristoylation site, and broadly, allosteric heterobifunctional degraders that destroy targets, rather than inhibiting them. Designed chemical linkers in bifunctional degraders can connect the allosteric ligand that binds the target protein and the E3 ubiquitin ligase warhead anchor. The physical properties and favored conformational state of the engineered linker can precisely coordinate the distance and orientation between the target and the recruited E3. Allosteric PROTACs, noncompetitive molecular glues, and bitopic ligands, with covalent links of allosteric ligands and orthosteric warheads, increase the effective local concentration of productively oriented and placed ligands. Through covalent chemical or peptide linkers, allosteric drugs can collaborate with competitive drugs, degrader anchors, or other molecules of choice, driving innovative drug discovery.

The protein MCL-1 is a crucial factor in regulating apoptosis, the programmed cell death, and thus plays a major role in numerous cancer types. The allosteric protein MCL-1 is naturally moderated by the BH3-only peptide BIM, which binds at its canonical binding groove. In its isolated form, BIM is disordered but assumes an α-helical shape when bound by MCL-1. The underlying binding mechanism (i.e., induced fit vs conformational selection), as well as the time scales of the signal cascade subsequent to binding, are not understood. Here, an artificially photoswitchable variant of the MCL-1/BIM complex was designed and investigated by transient infrared spectroscopy. By destabilizing the α-helix of BIM with a covalently linked azobenzene photoswitch, the dynamical response of the whole complex upon an ultrafast photo-perturbation was characterized. While the destabilized and partially unfolded BIM still binds to MCL-1, a step-like cascade of structural rearrangements of both, MCL-1 and BIM was detected, spanning a wide range of time scales from pico- to microseconds. The results indicate that BIM binds according to an induced fit mechanism, while the structural adaptations of MCL-1 may constitute an allosteric signal.

Allostery is a phenomenon that the protein activity is regulated when a non-functional site on it is bounded. This phenomenon is important in life process and disease therapy. However, it is difficult to study allostery due to the lack of knowledge. Facing this demand, we have created Allosteric Database (ASD) 10 years before to collect numerous kinds of allosteric data. In this review, we will introduce the 4 categories of data in ASD. For each category, we further reviewed how researchers applied ASD data to conduct studies. We focused on their research topics, analytical methods and conclusions. Several discoveries of new drug targets and allosteric modulators driven by ASD are also summarized. We hope this review could inspire researchers with new utilities of ASD data.

Allosteric mechanisms are pervasive in nature, but human-designed allosteric perturbagens are rare. The history of KRASG12C inhibitor development suggests that covalent chemistry may be a key to expanding the armamentarium of allosteric inhibitors. In that effort, irreversible targeting of a cysteine converted a non-deal allosteric binding pocket and low affinity ligands into a tractable drugging strategy. Here we examine the feasibility of expanding this approach to other allosteric pockets of RAS and kinase family members, given that both protein families are regulators of vital cellular processes that are often dysregulated in cancer and other human diseases. Moreover, these heavily studied families are the subject of numerous drug development campaigns that have resulted, sometimes serendipitously, in the discovery of allosteric inhibitors. We consequently conducted a comprehensive search for cysteines, a commonly targeted amino acid for covalent drugs, using AlphaFold-generated structures of those families. This new analysis presents potential opportunities for allosteric targeting of validated and understudied drug targets, with an emphasis on cancer therapy.

Allosteric proteins transition between 'inactive' and 'active' states. In general, such proteins assume distinct conformational states at the level of secondary, tertiary and/or quaternary structure. Different conformers of an allosteric protein can be antigenically dissimilar and induce antibodies with a highly distinctive specificities and neutralizing functional effects. Here we summarize studies on various functional types of monoclonal antibodies obtained against different allosteric conformers of the mannose-specific bacterial adhesin FimH - the most common cell attachment protein of Escherichia coli and other enterobacterial pathogens. Included are types of antibodies that activate the FimH function via interaction with ligand-induced binding sites or by wedging between domains as well as antibodies that inhibit FimH through orthosteric, parasteric, or novel dynasteric mechanisms. Understanding the molecular mechanism of antibody action against allosteric proteins provides insights on how to design antibodies with a desired functional effect, including those with neutralizing activity against bacterial and viral cell attachment proteins.

Inhibiting the main protease of SARS-CoV-2 is of great interest in tackling the COVID-19 pandemic caused by the virus. Most efforts have been centred on inhibiting the binding site of the enzyme. However, considering allosteric sites, distant from the active or orthosteric site, broadens the search space for drug candidates and confers the advantages of allosteric drug targeting. Here, we report the allosteric communication pathways in the main protease dimer by using two novel fully atomistic graph-theoretical methods: Bond-to-bond propensity, which has been previously successful in identifying allosteric sites in extensive benchmark data sets without a priori knowledge, and Markov transient analysis, which has previously aided in finding novel drug targets in catalytic protein families. Using statistical bootstrapping, we score the highest ranking sites against random sites at similar distances, and we identify four statistically significant putative allosteric sites as good candidates for alternative drug targeting.

The transporter BetP in C. glutamicum is essential in maintaining bacterial cell viability during hyperosmotic stress and functions by co-transporting betaine and Na+ into bacterial cells. Hyperosmotic stress leads to increased intracellular K+ concentrations which in turn promotes betaine binding. While structural details of multiple end state conformations of BetP have provided high resolution snapshots, how K+ sensing by the C-terminal domain is allosterically relayed to the betaine binding site is not well understood. In this study, we describe conformational dynamics in solution of BetP using amide hydrogen/deuterium exchange mass spectrometry. These reveal how K+ alters conformation of the disordered C- and N-terminal domains to allosterically reconfigure transmembrane helices 3, 8, and 10 to enhance betaine interactions. A map of the betaine binding site, at near single amino acid resolution, reveals a critical extrahelical H-bond mediated by TM3 with betaine.

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

Accurate development of allosteric modulators of GPCRs require a thorough assessment of their sequence, structure, and dynamics, toward gaining insights into their mechanisms of actions shared by family members, as well as dynamic features that distinguish subfamilies. Building on recent progress in the characterization of the signature dynamics of proteins, we analyzed here a dataset of 160 Class A GPCRs to determine their sequence similarities, structural landscape, and dynamic features across different species (human, bovine, mouse, squid, and rat), different activation states (active/inactive), and different subfamilies. The two dominant directions of variability across experimentally resolved structures, identified by principal component analysis of the dataset, shed light to cooperative mechanisms of activation, subfamily differentiation, and speciation of Class A GPCRs. The analysis reveals the functional significance of the conformational flexibilities of specific structural elements, including: the dominant role of the intracellular loop 3 (ICL3) together with the cytoplasmic ends of the adjoining helices TM5 and TM6 in enabling allosteric activation; the role of particular structural motifs at the extracellular loop 2 (ECL2) connecting TM4 and TM5 in binding ligands specific to different subfamilies; or even the differentiation of the N-terminal conformation across different species. Detailed analyses of the modes of motions accessible to the members of the dataset and their variations across members demonstrate how the active and inactive states of GPCRs obey distinct conformational dynamics. The collective fluctuations of the GPCRs are robustly defined in the active state, while the inactive conformers exhibit broad variance among members.

Drug research and development is a multidisciplinary field with its own successes. Yet, given the complexity of the process, it also faces challenges over the long development stages and even includes those that develop once a drug is marketed, i.e. drug toxicity and drug resistance. Better success can be achieved via well designed criteria in the early drug development stages. Here, we introduce the concepts of allostery and missense mutations, and argue that incorporation of these two intermittently linked biological phenomena into the early computational drug discovery stages would help to reduce the attrition risk in later stages of the process. We discuss the individual or in concert mechanisms of actions of mutations in allostery. Design of allosteric drugs is challenging compared to orthosteric drugs, yet they have been gaining popularity in recent years as alternative systems for the therapeutic regulation of proteins with an action-at-a-distance mode and non-invasive mechanisms. We propose an easy-to-apply computational allosteric drug discovery protocol which considers the mutation effect, and detail it with three case studies focusing on (1) analysis of effect of an allosteric mutation related to isoniazid drug resistance in tuberculosis; (2) identification of a cryptic pocket in the presence of an allosteric mutation of falcipain-2 as a malarial drug target; and (3) deciphering the effects of SARS-CoV-2 evolutionary mutations on a potential allosteric modulator with changes to allosteric communication paths.

Biological information processing networks rely on allosteric protein switches that dynamically interconvert biological signals. Construction of their artificial analogues is a central goal of synthetic biology and bioengineering. Receptor domain insertion is one of the leading methods for constructing chimeric protein switches. Here we present an in vitro expression-based platform for the analysis of chimeric protein libraries for which traditional cell survival or cytometric high throughput assays are not applicable. We utilise this platform to screen a focused library of chimeras between PQQ-glucose dehydrogenase and calmodulin. Using this approach, we identified 50 chimeras (approximately 23% of the library) that were activated by calmodulin-binding peptides. We analysed performance parameters of the active chimeras and demonstrated that their dynamic range and response times are anticorrelated, pointing to the existence of an inherent thermodynamic trade-off. We show that the structure of the ligand peptide affects both the response and activation kinetics of the biosensors suggesting that the structure of a ligand:receptor complex can influence the chimera's activation pathway. In order to understand the extent of structural changes in the reporter protein induced by the receptor domains, we have analysed one of the chimeric molecules by CD spectroscopy and hydrogen-deuterium exchange mass spectrometry. We concluded that subtle ligand-induced changes in the receptor domain propagated into the GDH domain and affected residues important for substrate and cofactor binding. Finally, we used one of the identified chimeras to construct a two-component rapamycin biosensor and demonstrated that core switch optimisation translated into improved biosensor performance.

While being a thoroughly studied model of dynamic allostery in a small protein, the pathway of signal transduction in the PDZ3 domain has not been fully determined. Here, we investigate peptide binding to the PDZ3 domain by conventional and fully data-driven analyses of molecular dynamics simulations. First, we identify isoleucine 37 as a key residue by widely used computational procedures such as cross-correlation and community network analysis. Simulations of the Ile37Ala mutant show disruption of the coordinated movements of spatially close regular elements of secondary structure. Then, we employ a recently developed unsupervised, data-driven procedure to determine an optimized reaction coordinate (slowest-relaxation eigenvector) of peptide binding. We use this reaction coordinate to improve sampling by restarting additional simulations from the transition state region. Significant differences in the distributions of some of the pairwise residue distances in the bound and unbound states emerge from the projection onto the optimized reaction coordinate. The unsupervised analysis shows that allosteric signaling is transduced from the β2 strand, which forms part of the peptide binding site, to the spatially adjacent β3 and β4 strands, and from there to the α3 helix. The domino-like transmission of a (peptide binding) signal along β strands and α helices that are close in three-dimensional space is likely to be a general mechanism of allostery in single-domain proteins.

While the accelerating quest for precision medicine requires new individually targeting and selective drugs, and the ability to work with so-called undruggable targets, the realm of allosteric drugs meeting this need remains largely uncharted. Generalizing the observations on two major drug targets with widely observed inherent allostery, GPCRs and kinases, we describe and discuss basic allosteric modes of action that are universally applicable in all types of structures and functions. Using examples of Class A GPCRs and CMGC protein kinases, we show how Allosteric Signalling and Probing Fingerprints can be used to identify potential allosteric sites and reveal effector-leads that may serve as a starting point for the development of allosteric drugs targeting these regulatory sites. A set of distinct characteristics of allosteric ligands was established, which highlights the versatility of their design and make them advantageous before their orthosteric counterparts in personalized medicine. We argue that rational design of allosteric drugs should begin with the search for latent sites or design of non-natural binding sites followed by fragment-based design of allosteric ligands and by the mutual adjustment of the site-ligand pair in order to achieve required drug efficacy. On the basis of the perturbative nature and reversibility of allosteric communication, we propose a generic protocol for computational design of allosteric effectors, enabling also the allosteric tuning of biologics, in obtaining allosteric control over protein functions.

Allosteric communication between distant protein sites represents a key mechanism of biomolecular regulation and signal transduction. Compared to other processes such as protein folding, however, the dynamical evolution of allosteric transitions is still not well understood. As an example of allosteric coupling between distant protein regions, we consider the global open-closed motion of the two domains of T4 lysozyme, which is triggered by local motion in the hinge region. Combining extensive molecular dynamics simulations with a correlation analysis of interresidue contacts, we identify a network of interresidue distances that move in a concerted manner. The cooperative process originates from a cogwheel-like motion of the hydrophobic core in the hinge region, which constitutes an evolutionary conserved and flexible transmission network. Through rigid contacts and the protein backbone, the small local changes of the hydrophobic core are passed on to the distant terminal domains and lead to the emergence of a rare global conformational transition. As in an Ising-type model, the cooperativity of the allosteric transition can be explained via the interaction of local fluctuations.

Understanding allostery in the Mycobacterium tuberculosis low molecular weight protein tyrosine phosphatase (MptpA) is a subject of great interest since MptpA is one of two protein tyrosine phosphatases (PTPs) from the pathogenic organism Mycobacterium tuberculosis expressed during host cell infection. Here, we combine computational modeling with solution NMR spectroscopy and we find that Q75 is an allosteric site. Removal of the polar side chain of Q75 by mutation to leucine results in a cascade of events that reposition the acid loop over the active site and relocates the catalytic aspartic acid (D126) at an optimal position for proton donation to the leaving aryl group of the substrate and for subsequent hydrolysis of the thiophosphoryl intermediate. The computational analysis is consistent with kinetic data, and NMR spectroscopy, showing that the Q75L mutant exhibits enhanced reaction kinetics with similar substrate binding affinity. We anticipate that our findings will motivate further studies on the possibility that MptpA remains passivated during the chronic state of infection and increases its activity as part of the pathogenic life cycle of M. tuberculosis possibly via allosteric means.

Allostery is a key biological control mechanism, and dynamic information flow provides a perspective to describe allosteric interactions in causal relationships. Here, as a novel implementation of the Gaussian Network Model (GNM) based Transfer Entropy (TE) calculations, we show that the dissection of dynamic information into subsets of slow dynamic modes discloses different layers of multi-directional allosteric pathways inherent in a given protein structure. In these subsets of slow modes, the degree of collectivity (Col) in the information transfer of residues with their TE values (TECol score) identifies distinct residues as powerful effectors, global information sources; showing themselves with a high dynamic capacity to collectively disseminate information to others. As exemplified on aspartate transcarbamoylase (ATCase), Na+/K+-adenosine triphosphatase (Na+/K+-ATPase), and human transient receptor potential melastatin 2 (TRPM2) along with a dataset of 20 proteins, these specific residues are associated with known active and allosteric sites. These information source residues, which collectively control others and lead allosteric communication pathways, hint at plausible binding sites for structure-based rational drug design.

Mutations in K-Ras GTPase replacing Gly12 with either Asp or Val are common in cancer. These mutations decelerate intrinsic and catalyzed GTP hydrolysis, leading to accumulation of K-Ras-GTP in cells. Signaling cascades initiated by K-Ras-GTP promote cell proliferation, survival, and invasion. Despite functional differences between the most frequent G12D mutation and the most aggressive and chemotherapy resistant G12V mutation, their long-suspected distinct structural features remain elusive. Using NMR, X-ray structures, and computational methods, we found that oncogenic mutants of K-Ras4B, the predominant splice variant of K-Ras, exhibit distinct conformational dynamics when GDP-bound, visiting the "active-like" conformational state similar to the one observed in GTP-bound K-Ras. This behavior distinguishes G12V from wild type and G12D K-Ras4B-GDP. The likely reason is interactions between the aliphatic sidechain of V12 and the Switch II region of K-Ras4BG12V-GDP, which are distinct in K-Ras4BG12D-GDP. In the X-ray structures, crystal contacts reduce the dynamics of the sidechain at position 12 by stabilizing the Switch I region of the protein. This explains why structural differences between G12V and G12D K-Ras have yet not been reported. Together, our results suggest a previously unknown mechanism of K-Ras activation. This mechanism relies on conformational dynamics caused by specific oncogenic mutations in the GDP-bound state. Our findings also imply that the therapeutic strategies decreasing the level of K-Ras-GTP by interfering with nucleotide exchange or by expediting GTP hydrolysis may work differently in different oncogenic mutants.

Hand-in-hand work of physics and evolution delivered protein universe with diversity of forms, sizes, and functions. Pervasiveness and advantageous traits of allostery made it an important component of the protein function regulation, calling for thorough investigation of its structural determinants and evolution. Learning directly from nature, we explored here allosteric communication in several major folds and repeat proteins, including α/β and β-barrels, β-propellers, Ig-like fold, ankyrin and α/β leucine-rich repeat proteins, which provide structural platforms for many different enzymatic and signalling functions. We obtained a picture of conserved allosteric communication characteristic in different fold types, modifications of the structure-driven signalling patterns via sequence-determined divergence to specific functions, as well as emergence and potential diversification of allosteric regulation in multi-domain proteins and oligomeric assemblies. Our observations will be instrumental in facilitating the engineering and de novo design of proteins with allosterically regulated functions, including development of therapeutic biologics. In particular, results described here may guide the identification of the optimal structural platforms (e.g. fold type, size, and oligomerization states) and the types of diversifications/perturbations, such as mutations, effector binding, and order-disorder transition. The tunable allosteric linkage across distant regions can be used as a pivotal component in the design/engineering of modular biological systems beyond the traditional scaffolding function.

The tumor protein 53 (p53) is involved in transcription-dependent and independent processes. Several p53 variants related to cancer have been found to impact protein stability. Other variants, on the contrary, might have little impact on structural stability and have local or long-range effects on the p53 interactome. Our group previously identified a loop in the DNA binding domain (DBD) of p53 (residues 207-213) which can recruit different interactors. Experimental structures of p53 in complex with other proteins strengthen the importance of this interface for protein-protein interactions. We here characterized with structure-based approaches somatic and germline variants of p53 which could have a marginal effect in terms of stability and act locally or allosterically on the region 207-213 with consequences on the cytosolic functions of this protein. To this goal, we studied 1132 variants in the p53 DBD with structure-based approaches, accounting also for protein dynamics. We focused on variants predicted with marginal effects on structural stability. We then investigated each of these variants for their impact on DNA binding, dimerization of the p53 DBD, and intramolecular contacts with the 207-213 region. Furthermore, we identified variants that could modulate long-range the conformation of the region 207-213 using a coarse-grain model for allostery and all-atom molecular dynamics simulations. Our predictions have been further validated using enhanced sampling methods for 15 variants. The methodologies used in this study could be more broadly applied to other p53 variants or cases where conformational changes of loop regions are essential in the function of disease-related proteins.

Allostery commonly refers to the mechanism that regulates protein activity through the binding of a molecule at a different, usually distal, site from the orthosteric site. The omnipresence of allosteric regulation in nature and its potential for drug design and screening render the study of allostery invaluable. Nevertheless, challenges remain as few computational methods are available to effectively predict allosteric sites, identify signalling pathways involved in allostery, or to aid with the design of suitable molecules targeting such sites. Recently, bond-to-bond propensity analysis has been shown successful at identifying allosteric sites for a large and diverse group of proteins from knowledge of the orthosteric sites and its ligands alone by using network analysis applied to energy-weighted atomistic protein graphs. To address the identification of signalling pathways, we propose here a method to compute and score paths of optimised propensity that link the orthosteric site with the identified allosteric sites, and identifies crucial residues that contribute to those paths. We showcase the approach with three well-studied allosteric proteins: h-Ras, caspase-1, and 3-phosphoinositide-dependent kinase-1 (PDK1). Key residues in both orthosteric and allosteric sites were identified and showed agreement with experimental results, and pivotal signalling residues along the pathway were also revealed, thus providing alternative targets for drug design. By using the computed path scores, we were also able to differentiate the activity of different allosteric modulators.

Herein we examine the determinants of the allosteric inhibition of the mitochondrial chaperone TRAP1 by a small molecule ligand. The knowledge generated is harnessed into the design of novel derivatives with interesting biological properties. TRAP1 is a member of the Hsp90 family of proteins, which work through sequential steps of ATP processing coupled to client-protein remodeling. Isoform selective inhibition of TRAP1 can provide novel information on the biomolecular mechanisms of molecular chaperones, as well as new insights into the development of small molecules with therapeutic potential. Our analysis of the interactions between an active first-generation allosteric ligand and TRAP1 shows how the small molecule induces long-range perturbations that influence the attainment of reactive poses in the active site. At the same time, the dynamic adaptation of the allosteric binding pocket to the presence of the first-generation compound sets the stage for the design of a set of second-generation ligands: the characterization of the formation/disappearance of pockets around the allosteric site that is used to guide optimize the ligands' fit for the allosteric site and improve inhibitory activities. The effects of the newly designed molecules are validated experimentally in vitro and in vivo. We discuss the implications of our approach as a promising strategy towards understanding the molecular determinants of allosteric regulation in chemical and molecular biology, and towards speeding up the design of allosteric small molecule modulators.

The effect of perturbation at the allosteric site was investigated through several replicas of molecular dynamics (MD) simulations conducted on bacterial phosphofructokinase (SaPFK). In our previous work, an alternative binding site was estimated to be allosteric in addition to the experimentally reported one. To highlight the effect of both allosteric sites on receptor's dynamics, MD runs were carried out on apo forms with and without perturbation. Perturbation was achieved via incorporating multiple bond restraints for residue pairs located at the allosteric site. Restraints applied to the predicted site caused one dimer to stiffen, whereas an increase in mobility was detected in the same dimer when the experimentally resolved site was restrained. Fluctuations in Cα-Cα distances which is used to disclose residues with high potential of communication indicated a marked increase in signal transmission within each dimer as the receptor switched to a restrained state. Cross-correlation of positional fluctuations indicated an overall decrease in the magnitude of both positive and negative correlations when restraints were employed on the predicted allosteric site whereas an exact opposite effect was observed for the reported site. Finally, mutual correspondence between positional fluctuations noticeably increased with restraints on predicted allosteric site, whereas an opposite effect was observed for restraints applied on experimentally reported one. In view of these findings, it is clear that the perturbation of either one of two allosteric sites effected the dynamics of the receptor with a distinct and contrasting character.

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often "cryptic" without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein-ligand systems.

Proteins fluctuate between different conformations in solution, and these conformational fluctuations can be important for protein function and allosteric regulation. The chorismate mutase from Saccharomyces cerevisiae (ScCM), a key enzyme in the biosynthesis of aromatic amino acids, is allosterically activated and inhibited by tryptophan and tyrosine, respectively. It was initially proposed that in the absence of effector, ScCM fluctuates between activated R and inhibited T conformations according to the Monod-Wyman-Changeux (MWC) model, although a more complex regulation pattern was later suggested by mutagenesis and kinetic data. Here we used NMR relaxation dispersion experiments to understand the conformational fluctuations on the microsecond-to-millisecond timescale that occur in ScCM. In the absence of allosteric effectors, ScCM did not exclusively exchange between T and R conformations, suggesting that the two-state MWC model is insufficient to explain conformational dynamics. Addition of tyrosine led to the quenching of much of the motion on this timescale, while new motions were identified in the presence of tryptophan. These new motions are consistent with conformational fluctuations into an alternative conformation that may be important for enzyme activity.

The family of coarse-grained models for protein dynamics known as Elastic Network Models (ENMs) require careful choice of parameters to represent well experimental measurements or fully-atomistic simulations. The most basic ENM that represents each protein residue by a node at the position of its C-alpha atom, all connected by springs of equal stiffness, up to a cut-off in distance. Even at this level a choice is required of the optimum cut-off distance and the upper limit of elastic normal modes taken in any sum for physical properties, such as dynamic correlation or allosteric effects on binding. Additionally, backbone-enhanced ENM (BENM) may improve the model by allocating a higher stiffness to springs that connect along the protein backbone. This work reports on the effect of varying these three parameters (distance and mode cutoffs, backbone stiffness) on the dynamical structure of three proteins, Catabolite Activator Protein (CAP), Glutathione S-transferase (GST), and the SARS-CoV-2 Main Protease (M pro ). Our main results are: (1) balancing B-factor and dispersion-relation predictions, a near-universal optimal value of 8.5 Å is advisable for ENMs; (2) inhomogeneity in elasticity brings the first mode containing spatial structure not well-resolved by the ENM typically within the first 20; (3) the BENM only affects modes in the upper third of the distribution, and, additionally to the ENM, is only able to model the dispersion curve better in this vicinity; (4) BENM does not typically affect fluctuation-allostery, which also requires careful treatment of the effector binding to the host protein to capture.

Intrinsically disordered proteins/regions (IDPRs) are a very large and functionally important class of proteins that participate in weak multivalent interactions in protein complexes. They are recalcitrant for interrogations using X-ray crystallography and cryo-EM. The IDPRs observed at the interface of the photosynthetic pigment protein complexes (PPCs) remain much less clear, e.g., the major cyanobacterial light-harvesting complex (PBS) contains an unstructured PB-loop insertion in the phycocyanobilin domain (PB domain) of ApcE (the largest polypeptide in PBS). Here, a joint platform is built to probe such structural domains. This platform is characterized by two-round progressive justifications of in silico models by using the structural mass spectrometry data. First, the AlphaFold-generated 3D structure of the PB domain (containing PB-loop) was justified in the context of PBS. Second, docking the AlphaFold-generated ApcG (a ligand) into the first-step justified structure (a receptor). The final ligand-receptor complex was then subjected to a second-round justification, again, by using unequivocal isotopically-encoded cross-links identified in LC-MS/MS. This work reveals a full-length PB-loop structure modelled in the PBS basal cylinder, free from any spatial conflicts against the other subunits in PBS. The structure of PB domain highlights the close associations of the intrinsically disordered PB-loop with its binding partners in PBS, including ApcG, another IDPR. The PB-loop region involved in the binding of photosystem II (PSII) is also discussed in the context of excitation energy transfer regulation. This work calls attention to the highly disordered, yet interrogatable interface between the light-harvesting antenna complexes and the reaction centers.

Allosteric regulation is the most direct and efficient way of regulating protein function, wherein proteins transmit the perturbations at one site to another distinct functional site. Deciphering the mechanism of allosteric regulation is of vital importance for the comprehension of both physiological and pathological events in vivo as well as the rational allosteric drug design. However, it remains challenging to elucidate dominant allosteric signal transduction pathways, especially for large and multi-component protein machineries where long-range allosteric regulation exits. One of the quintessential examples having long-range allosteric regulation is the ternary complex, SPRED1-RAS-neurofibromin type 1 (NF1, a RAS GTPase-activating protein), in which SPRED1 facilitates RAS-GTP hydrolysis by interacting with NF1 at a distal, allosteric site from the RAS binding site. To address the underlying mechanism, we performed extensive Gaussian accelerated molecular dynamics simulations and Markov state model analysis of KRAS-NF1 complex in the presence and absence of SPRED1. Our findings suggested that SPRED1 loading allosterically enhanced KRAS-NF1 binding, but hindered conformational transformation of the NF1 catalytic center for RAS hydrolysis. Moreover, we unveiled the possible allosteric pathways upon SPRED1 binding through difference contact network analysis. This study not only provided an in-depth mechanistic insight into the allosteric regulation of KRAS by SPRED1, but also shed light on the investigation of long-range allosteric regulation among complex macromolecular systems.

The human soluble epoxide hydrolase (hsEH) is a key regulator of epoxy fatty acid (EpFA) metabolism. Inhibition of sEH can maintain endogenous levels of beneficial EpFAs and reduce the levels of their corresponding diol products, thus ameliorating a variety of pathological conditions including cardiovascular, central nervous system and metabolic diseases. The quest for orthosteric drugs that bind directly to the catalytic crevice of hsEH has been prolonged and sustained over the past decades, but the disappointing outcome of clinical trials to date warrants alternative pharmacological approaches. Previously, we have shown that hsEH can be allosterically inhibited by the endogenous electrophilic lipid 15-deoxy-Δ12,14-Prostaglandin-J2, via covalent adduction to two cysteines, C423 and C522. In this study, we explore the properties and behaviour of three electrophilic lipids belonging to the class of the nitro fatty acids, namely 9- and 10-nitrooleate and 10-nitrolinoleate. Biochemical and biophysical investigations revealed that, in addition to C423 and C522, nitro fatty acids can covalently bind to additional nucleophilic residues in hsEH C-terminal domain (CTD), two of which predicted in this study to be latent allosteric sites. Systematic mapping of the protein mutational space and evaluation of possible propagation pathways delineated selected residues, both in the allosteric patches and in other regions of the enzyme, envisaged to play a role in allosteric signalling. The responses elicited by the ligands on the covalent adduction sites supports future fragment-based design studies of new allosteric effectors for hsEH with increased efficacy and selectivity.

Skp1(S-phase kinase-associated protein 1 - Homo sapiens) is an adapter protein of the SCF(Skp1-Cullin1-Fbox) complex, which links the constant components (Cul1-RBX) and the variable receptor (F-box proteins) in Ubiquitin E3 ligase. It is intriguing how Skp1 can recognise and bind to a variety of structurally different F-box proteins. For practical reasons, previous efforts have used truncated Skp1, and thus it has not been possible to track the crucial aspects of the substrate recognition process. In this background, we report the solution structure of the full-length Skp1 protein determined by NMR spectroscopy for the first time and investigate the sequence-dependent dynamics in the protein. The solution structure reveals that Skp1 has an architecture: β1-β2-H1-H2-L1-H3-L2-H4-H5-H6-H7(partially formed) and a long tail-like disordered C-terminus. Structural analysis using DALI (Distance Matrix Alignment) reveals conserved domain structure across species for Skp1. Backbone dynamics investigated using NMR relaxation suggest substantial variation in the motional timescales along the length of the protein. The loops and the C-terminal residues are highly flexible, and the (R2/R1) data suggests μs-ms timescale motions in the helices as well. Further, the dependence of amide proton chemical shift on temperature and curved profiles of their residuals indicate that the residues undergo transitions between native state and excited state. The curved profiles for several residues across the length of the protein suggest that there are native-like low-lying excited states, particularly for several C-terminal residues. Our results provide a rationale for how the protein can adapt itself, bind, and get functionally associated with other proteins in the SCF complex by utilising its flexibility and conformational sub-states.

The FimH protein of Escherichia coli is a model two-domain adhesin that is able to mediate an allosteric catch bond mechanism of bacterial cell attachment, where the mannose-binding lectin domain switches from an 'inactive' conformation with fast binding to mannose to an 'active' conformation with slow detachment from mannose. Because mechanical tensile force favors separation of the domains and, thus, FimH activation, it has been thought that the catch bonds can only be manifested in a fluidic shear-dependent mode of adhesion. Here, we used recombinant FimH variants with a weakened inter-domain interaction and show that a fast and sustained allosteric activation of FimH can also occur under static, non-shear conditions. Moreover, it appears that lectin domain conformational activation happens intrinsically at a constant rate, independently from its ability to interact with the pilin domain or mannose. However, the latter two factors control the rate of FimH deactivation. Thus, the allosteric catch bond mechanism can be a much broader phenomenon involved in both fast and strong cell-pathogen attachments under a broad range of hydrodynamic conditions. This concept that allostery can enable more effective receptor-ligand interactions is fundamentally different from the conventional wisdom that allostery provides a mechanism to turn binding off under specific conditions.

Enterobacteria phage P1 expresses two types of tail fibre, S and S'. Despite the wide usage of phage P1 for transduction, the host range and the receptor for its alternative S' tail fibre was never determined. Here, a ΔS-cin Δpac E. coli P1 lysogenic strain was generated to allow packaging of phagemid DNA into P1 phage having either S or S' tail fibre. P1(S') could transduce phagemid DNA into Shigella flexneri 2a 2457O, Shigella flexneri 5a M90T and Escherichia coli O3 efficiently. Mutational analysis of the O-antigen assembly genes and LPS inhibition assays indicated that P1(S') transduction requires at least one O-antigen unit. E. coli O111:B4 LPS produced a high neutralising effect against P1(S') transduction, indicating that this E. coli strain could be susceptible to P1(S')-mediated transduction. Mutations in the O-antigen modification genes of S. flexneri 2a 2457O and S. flexneri 5a M90T did not cause significant changes to P1(S') transduction efficiency. A higher transduction efficiency of P1(S') improved the delivery of a cas9 antimicrobial phagemid into both S. flexneri 2457O and M90T. These findings provide novel insights into P1 tropism-switching, by identifying the bacterial strains which are susceptible to P1(S')-mediated transduction, as well as demonstrating its potential for delivering a DNA sequence-specific Cas9 antimicrobial into clinically relevant S. flexneri.

Understanding molecular principles underlying Hsp90 chaperone functions and modulation of client activity is fundamental to dissect activation mechanisms of many proteins. In this work, we performed a computational investigation of the Hsp90-Hsp70-Hop-CR client complex to examine allosteric regulatory mechanisms underlying dynamic chaperone interactions and principles of chaperone-dependent client recognition and remodeling. Conformational dynamics analysis using high-resolution coarse-grained simulations and ensemble-based local frustration analysis suggest that the Hsp90 chaperone could recognize and recruit the GR client by invoking reciprocal dynamic exchanges near the intermolecular interfaces with the client. Using mutational scanning of the intermolecular residues in the Hsp90-Hsp70-Hop-GR complex, we identified binding energy hotspots in the regulatory complex. Perturbation-based network analysis and dynamic fluctuations-based modeling of allosteric residue potentials are employed for a detailed analysis of allosteric interaction networks and identification of conformational communication switches. We found that allosteric interactions between the Hsp90, the client-bound Hsp70 and Hop cochaperone can define two allosteric residue clusters that control client recruitment in which the intrinsic Hsp70 allostery is exploited to mediate integration of the Hsp70-bound client into the Hsp90 chaperone system. The results suggest a model of dynamics-driven allostery that enables efficient client recruitment and loading through allosteric couplings between intermolecular interfaces and communication switch centers. This study showed that the Hsp90 interactions with client proteins may operate under dynamic-based allostery in which ensembles of preexisting conformational states and intrinsic allosteric pathways present in the Hsp90 and Hsp70 chaperones can be exploited for recognition and integration of substrate proteins.

Accurate quantitative estimates of protein-membrane interactions are critical to studies of membrane proteins. Here, we demonstrate that thermodynamic analyses based on current hydropathy scales do not account for the significant and experimentally determined effects that Ca2+ or Mg2+ have on protein-membrane interactions. We examined distinct modes of interaction (interfacial partitioning and folding and transmembrane insertion) by studying three highly divergent peptides: Bid-BH3 (derived from apoptotic regulator Bid), peripherin-2-derived prph2-CTER, and the cancer-targeting pH-Low-Insertion-Peptide (pHLIP). Fluorescence experiments demonstrate that adding 1-2 mM of divalent cations led to a substantially more favorable bilayer partitioning and insertion, with free energy differences of 5-15 kcal/mol.

RecG, belonging to the category of Superfamily-2 plays a vital role in rescuing different kinds of stalled fork. The elemental mechanism of the helicase activity of RecG with several non-homologous stalled fork structures resembling intermediates formed during the process of DNA repair has been investigated in the present study to capture the dynamic stages of genetic rearrangement. The functional characterization has been exemplified through quantifying the response of the substrate in terms of their molecular heterogeneity and dynamical response by employing single-molecule fluorescence methods. An elevated processivity of RecG is observed for the stalled fork where progression of lagging daughter strand is ahead as compared to that of the leading strand. Through precise alteration of its function in terms of unwinding, depending upon the substrate DNA, RecG catalyzes the formation of Holliday junction from a stalled fork DNA. RecG is found to adopt an asymmetric mode of locomotion to unwind the lagging daughter strand for facilitating formation of Holliday junction that acts as a suitable intermediate for recombinational repair pathway. Our results emphasize the mechanism adopted by RecG during its 'sliding back' mode along the lagging daughter strand to be 'active translocation and passive unwinding'. This also provide clues as to how this helicase decides and controls the mode of translocation along the DNA to unwind.

As a subgroup of sorting nexins (SNXs) that contain regulator of G protein signaling homology (RH) domain, SNX-RH proteins, including SNX13, SNX14 and SNX25, were proposed to play bifunctional roles in protein sorting and GPCR signaling regulation. However, mechanistic details of SNX-RH proteins functioning via RH domain remain to be illustrated. Here, we delineate crystal structures of the RH domains of SNX13 and SNX25, revealing a homodimer of SNX13 RH domain mediated by unique extended α4 and α5 helices, and a thiol modulated homodimer of SNX25-RH triggered by a unique cysteine on α6 helix. Further studies showed that RH domains of SNX-RH do not possess binding capacity toward Gα subunits, owing to the lack of critical residues for interaction. Thus, this study identifies a group of novel non-canonical RH domains that can act as a dimerization module in sorting nexins, which provides structural basis for mechanism studies on SNX-RH protein functions.