Allosteric Signal Transmission Research Articles

β-sheets mediate the conformational change and allosteric signal transmission between the AsLOV2 termini.

Avena sativa phototropin 1 light-oxygen-voltage 2 domain (AsLOV2) is a model protein of Per-Arnt-Sim (PAS) superfamily, characterized by conformational changes in response to external environmental stimuli. This conformational change begins with the unfolding of the N-terminal A'α helix in the dark state followed by the unfolding of the C-terminal Jα helix. The light state is characterized by the unfolded termini and the subsequent modifications in hydrogen bond patterns. In this photoreceptor, β-sheets are identified as crucial components for mediating allosteric signal transmission between the two termini. Through combined experimental and computational investigations, the Hβ and Iβ strands are recognized as the most critical and influential β-sheets in AsLOV2's allosteric mechanism. To elucidate the role of these β-sheets, we introduced 13 distinct mutations (F490L, N492A, L493A, F494L, H495L, L496F, Q497A, R500A, F509L, Q513A, L514A, D515V, and T517V) and conducted comprehensive molecular dynamics simulations. In-depth hydrogen bond analyses emphasized the role of two hydrogen bonds, Asn482-Leu453 and Gln479-Val520, in the observed distinct behaviors of L493A, L496F, Q497A, and D515V mutants. This illustrates the role of β-sheets in the transmission of the allosteric signal upon the photoactivation of the light state.

Allosteric pathways of SARS and SARS-CoV-2 spike protein identified by neural relational inference.

The receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein must undergo a crucial conformational transition to invade human cells. It is intriguing that this transition is accompanied by a synchronized movement of the entire spike protein. Therefore, it is possible to design allosteric regulators targeting non-RBD but hindering the conformational transition of RBD. To understand the allosteric mechanism in detail, we establish a computational framework by integrating coarse-grained molecular dynamic simulations and a state-of-the-art neural network model called neural relational inference. Leveraging this framework, we have elucidated the allosteric pathway of the SARS-CoV-2 spike protein at the residue level and identified the molecular mechanisms involved in the transmission of allosteric signals. The movement of D614 is coupled with that of Q321. This interaction subsequently influences the movement of K528/K529, ultimately coupling with the movement of RBD during conformational changes. Mutations that weaken the interactions within this pathway naturally block the allosteric signal transmission, thereby modulating the conformational transitions. This observation also offers a rationale for the distinct allosteric patterns observed in the SARS-CoV spike protein. Our result provides a useful method for analyzing the dynamics of potential viral variants in the future.

Investigating the Activation Mechanism Differences between Human and Mouse cGAS by Molecular Dynamics Simulations.

Cyclic GMP-AMP synthase (cGAS) has been widely investigated as a drug target for its crucial role in innate immunity. However, the inhibitors designed using mouse model were often shown to be ineffective for humans. This outcome indicates that the activation mechanisms of human and mouse cGAS (mcGAS) are different. The cGAS activation is achieved by dimerization via binding to DNA, the detailed mechanism of which, however, is not entirely clear. To investigate these mechanisms, molecular dynamics (MD) simulations were performed on several states of four types of cGAS, namely, the mcGAS, the wild-type and A- and C-type mutations of human cGAS (hcGAS). We find that sequence differences between hcGAS and mcGAS can directly affect the protein structure stability, especially that of the siteB domain. The sequence and structural differences also contribute to DNA-binding differences. In addition, the conformational fluctuations of cGAS are found to correlate with the regulation of catalytic capacity. More importantly, we illustrate that dimerization enhances the correlation among distant residues and significantly reinforces the allosteric signal transmission among the DNA-binding interfaces and the catalytic pocket, which facilitates rapid immune response to cytosolic DNA. We conclude that siteB domain plays a prominent role in mcGAS activation, while siteA domain is key to hcGAS activation.

Improving tyrosol production efficiency through shortening the allosteric signal transmission distance of pyruvate decarboxylase.

Tyrosol is an important chemical in medicine and chemical industries, which can be synthesized by a four-enzyme cascade pathway constructed in our previous study. However, the low catalytic efficiency of pyruvate decarboxylase from Candida tropicalis (CtPDC) in this cascade is a rate-limiting step. In this study, we resolved the crystal structure of CtPDC and investigated the mechanism of allosteric substrate activation and decarboxylation of this enzyme toward 4-hydroxyphenylpyruvate (4-HPP). In addition, based on the molecular mechanism and structural dynamic changes, we conducted protein engineering of CtPDC to improve decarboxylation efficiency. The conversion of the best mutant, CtPDCQ112G/Q162H/G415S/I417V (CtPDCMu5), had over two-fold improvement compared to the wild-type. Molecular dynamic (MD) simulation revealed that the key catalytic distances and allosteric transmission pathways were shorter in CtPDCMu5 than in the wild type. Furthermore, when CtPDC in the tyrosol production cascade was replaced with CtPDCMu5, the tyrosol yield reached 38 g·L-1 with 99.6% conversion and 1.58 g·L-1·h-1 space-time yield in 24 h through further optimization of the conditions. Our study demonstrates that protein engineering of the rate-limiting enzyme in the tyrosol synthesis cascade provides an industrial-scale platform for the biocatalytic production of tyrosol. KEY POINTS: • Protein engineering of CtPDC based on allosteric regulation improved the catalytic efficiency of decarboxylation. • The application of the optimum mutant of CtPDC removed the rate-limiting bottleneck in the cascade. • The final titer of tyrosol reached 38 g·L-1 in 24 h in 3 L bioreactor.

Identification of residues involved in allosteric signal transmission from amino acid binding site of pyruvate kinase muscle isoform 2.

PKM2 is a rate-limiting enzyme in the glycolytic process and is involved in regulating tumor proliferation. Several amino acids (AAs) such as Asn, Asp, Val, and Cys have been shown to bind to the AA binding pocket of PKM2 and modulate its oligomeric state, substrate binding affinity, and activity. Although previous studies have attributed that the main chain and side chain of bound AAs are responsible for initiating signal to regulate PKM2, the signal transduction pathway remains elusive. To identify the residues involved in signal transfer process, N70 and N75 located at two ends of a β strand connecting the active site and AA binding pocket were altered. Biochemical studies of these variants with various AA ligands (Asn, Asp, Val, and Cys), illustrate that N70 and N75, along with β1 connecting these residues are part of the signal transduction pathway between the AA binding pocket and the active site. The results demonstrate that mutation of N70 to D prevents the transfer of the inhibitory signal mediated by Val and Cys, whereas N75 to L alteration blocks the activating signal initiated by Asn and Asp. Taken together, this study confirms that N70 is one of the residues responsible for transmitting the inhibitory signal and N75 is involved in the activation signal flow.

Dynamics and Allosteric Information Pathways of Unphosphorylated c-Cbl.

Human c-Cbl is a RING-type ligase and plays a central role in the protein degradation cascade. To elucidate its conformational changes related to substrate binding, we performed molecular dynamics simulations of different variants/states of c-Cbl for a cumulative time of 68 μs. Our simulations demonstrate that before the substrate binds, the RING domain samples a broad set of conformational states at a biologically relevant salt concentration, including the closed, partially open, and fully open states, whereas substrate binding leads to a restricted conformational sampling. Phe378 and the C-terminal region play an essential role in stabilizing the partially open state. To visualize the allosteric signal transmission pathways from the substrate-binding site to the 40 Å apart RING domain and identify the critical residues for allostery, we have created a subgraph from the optimal and suboptimal paths. Redundant paths are seen in the SH2 domain where the substrate binds, while the major bottlenecks are found at the junction between the SH2 domain and the linker helix region as well as that between the SH2 domain and the 4H bundle. These bottlenecks separate the paths into two overall routes. The nodes/residues at the bottlenecks on the subgraph are considered allosteric hot spots. This subgraph approach provides a general tool for network visualization and determination of critical residues for allostery. The structurally and allosterically critical residues identified in our work are testable and would provide valuable insights into the emerging strategies for drug discovery, such as targeted protein degradation.

Allosteric mechanism for SL RNA recognition by polypyrimidine tract binding protein RRM1: An atomistic MD simulation and network-based study

Polypyrimidine tract-binding protein (PTB), an RNA-binding protein, is involved in the regulation of diverse processes in mRNA metabolism. However, the allosteric modulation of its binding with RNA remains unclear. We explore the dynamic characteristics of PTB RNA recognition motif 1 (RRM1) in its RNA-free and wild-type/mutant RNA-bound states to understand the issues using molecular dynamics (MD) simulation, perturbation response scanning (PRS) and protein structure network (PSN) models. It is found that RNA binding strengthens RRM1 stability, while L151G mutation in α3 helix far away from the interface makes the complex unstable. The latter is caused by long-distance dynamic couplings, which makes intermolecular electrostatic and entropy energies unfavorable. The weakened couplings between interface β sheets and C-terminal parts upon mutation reveal RNA recognition is co-regulated by these regions. Interestingly, PRS analysis reveals the allostery caused by the perturbation on α3 helix has already been pre-encoded in the equilibrium dynamics of the protein structure. PSN analysis shows the details of the allosteric signal transmission, revealing the necessity of strong couplings between α3 helix and interface for maintaining the high binding affinity. This study sheds light on the mechanisms of PTB allostery and RNA recognition and can provide important information for drug design.

Programmable allosteric DNA regulations for molecular networks and nanomachines.

Structure-based molecular regulations have been widely adopted to modulate protein networks in cells and recently developed to control allosteric DNA operations in vitro. However, current examples of programmable allosteric signal transmission through integrated DNA networks are stringently constrained by specific design requirements. Developing a new, more general, and programmable scheme for establishing allosteric DNA networks remains challenging. Here, we developed a general strategy for programmable allosteric DNA regulations that can be finely tuned by varying the dimensions, positions, and number of conformational signals. By programming the allosteric signals, we realized fan-out/fan-in DNA gates and multiple-layer DNA cascading networks, as well as expanding the approach to long-range allosteric signal transmission through tunable DNA origami nanomachines ~100 nm in size. This strategy will enable programmable and complex allosteric DNA networks and nanodevices for nanoengineering, chemical, and biomedical applications displaying sense-compute-actuate molecular functionalities.

Frustration-driven allosteric regulation and signal transmission in the SARS-CoV-2 spike omicron trimer structures: a crosstalk of the omicron mutation sites allosterically regulates tradeoffs of protein stability and conformational adaptability.

Dissecting the regulatory principles underlying function and activity of the SARS-CoV-2 spike protein at the atomic level is of paramount importance for understanding the mechanisms of virus transmissibility and immune escape. In this work, we introduce a hierarchical computational approach for atomistic modeling of allosteric mechanisms in the SARS-CoV-2 Omicron spike proteins and present evidence of a frustration-based allostery as an important energetic driver of the conformational changes and spike activation. By examining conformational landscapes and the residue interaction networks in the SARS-CoV-2 Omicron spike protein structures, we have shown that the Omicron mutational sites are dynamically coupled and form a central engine of the allosterically regulated spike machinery that regulates the balance and tradeoffs between conformational plasticity, protein stability, and functional adaptability. We have found that the Omicron mutational sites at the inter-protomer regions form regulatory hotspot clusters that control functional transitions between the closed and open states. Through perturbation-based modeling of allosteric interaction networks and diffusion analysis of communications in the closed and open spike states, we have quantified the allosterically regulated activation mechanism and uncover specific regulatory roles of the Omicron mutations. Atomistic reconstruction of allosteric communication pathways and kinetic modeling using Markov transient analysis reveal that the Omicron mutations form the inter-protomer electrostatic bridges that operate as a network of coupled regulatory switches that could control global conformational changes and signal transmission in the spike protein. The results of this study have revealed distinct and yet complementary roles of the Omicron mutation sites as a network of hotspots that enable allosteric modulation of structural stability and conformational changes which are central for spike activation and virus transmissibility.

Allosteric Effect of Nanobody Binding on Ligand-Specific Active States of the β2 Adrenergic Receptor.

Nanobody binding stabilizes G-protein-coupled receptors (GPCR) in a fully active state and modulates their affinity for bound ligands. However, the atomic-level basis for this allosteric regulation remains elusive. Here, we investigate the conformational changes induced by the binding of a nanobody (Nb80) on the active-like β2 adrenergic receptor (β2AR) via enhanced sampling molecular dynamics simulations. Dimensionality reduction analysis shows that Nb80 stabilizes structural features of the β2AR with an ∼14 Å outward movement of transmembrane helix 6 and a close proximity of transmembrane (TM) helices 5 and 7, and favors the fully active-like conformation of the receptor, independent of ligand binding, in contrast to the conditions under which no intracellular binding partner is bound, in which case the receptor is only stabilized in an intermediate-active state. This activation is supported by the residues located at hotspots located on TMs 5, 6, and 7, as shown by supervised machine learning methods. Besides, ligand-specific subtle differences in the conformations assumed by intracellular loop 2 and extracellular loop 2 are captured from the trajectories of various ligand-bound receptors in the presence of Nb80. Dynamic network analysis further reveals that Nb80 binding triggers tighter and stronger local communication networks between the Nb80 and the ligand-binding sites, primarily involving residues around ICL2 and the intracellular end of TM3, TM5, TM6, as well as ECL2, ECL3, and the extracellular ends of TM6 and TM7. In particular, we identify unique allosteric signal transmission mechanisms between the Nb80-binding site and the extracellular domains in conformations modulated by a full agonist, BI167107, and a G-protein-biased partial agonist, salmeterol, involving mainly TM1 and TM2, and TM5, respectively. Altogether, our results provide insights into the effect of intracellular binding partners on the GPCR activation mechanism, which should be taken into account in structure-based drug discovery.

Light-Induced Change of Arginine Conformation Modulates the Rate of Adenosine Triphosphate to Cyclic Adenosine Monophosphate Conversion in the Optogenetic System Containing Photoactivated Adenylyl Cyclase.

We report the first computational characterization of an optogenetic system composed of two photosensing BLUF (blue light sensor using flavin adenine dinucleotide) domains and two catalytic adenylyl cyclase (AC) domains. Conversion of adenosine triphosphate (ATP) to the reaction products, cyclic adenosine monophosphate (cAMP) and pyrophosphate (PPi), catalyzed by ACs initiated by excitation in photosensing domains has emerged in the focus of modern optogenetic applications because of the request in photoregulated enzymes that modulate cellular concentrations of signaling messengers. The photoactivated AC from the soil bacterium Beggiatoa sp. (bPAC) is an important model showing a considerable increase in the ATP to cAMP conversion rate in the catalytic domain after the illumination of the BLUF domain. The 1 μs classical molecular dynamics simulations reveal that the activation of the BLUF domain leading to tautomerization of Gln49 in the chromophore-binding pocket results in switching of the position of the side chain of Arg278 in the active site of AC. Allosteric signal transmission pathways between Gln49 from BLUF and Arg278 from AC were revealed by the dynamical network analysis. The Gibbs energy profiles of the ATP → cAMP + PPi reaction computed using QM(DFT(ωB97X-D3/6-31G**))/MM(CHARMM) molecular dynamics simulations for both Arg278 conformations in AC clarify the reaction mechanism. In the light-activated system, the corresponding arginine conformation stabilizes the pentacoordinated phosphorus of the α-phosphate group in the transition state, thus lowering the activation energy. Simulations of the bPAC system with the Tyr7Phe replacement in the BLUF demonstrate occurrence of both arginine conformations in an equal ratio, explaining the experimentally observed intermediate catalytic activity of the bPAC-Y7F variant as compared with the dark and light states of the wild-type bPAC.

Allosteric signaling in C-linker and cyclic nucleotide-binding domain of HCN2 channels

Opening of hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels is controlled by membrane hyperpolarization and binding of cyclic nucleotides to the tetrameric cyclic nucleotide-binding domain (CNBD), attached to the C-linker (CL) disk. Confocal patch-clamp fluorometry revealed pronounced cooperativity of ligand binding among protomers. However, by which pathways allosteric signal transmission occurs remained elusive. Here, we investigate how changes in the structural dynamics of the CL-CNBD of mouse HCN2 upon cAMP binding relate to inter- and intrasubunit signal transmission. Applying a rigidity-theory-based approach, we identify two intersubunit and one intrasubunit pathways that differ in allosteric coupling strength between cAMP-binding sites or toward the CL. These predictions agree with results from electrophysiological and patch-clamp fluorometry experiments. Our results map out distinct routes within the CL-CNBD that modulate different cAMP-binding responses in HCN2 channels. They signify that functionally relevant submodules may exist within and across structurally discernable subunits in HCN channels.

In silico study of intrinsic dynamics of full-length apo-ACE2 and RBD-ACE2 complex

The key step for SARS-CoV-2 to infect human cells is the membrane fusion triggered by the binding of the viral extracellular Spike protein to the human extracellular receptor, the angiotensin-converting enzyme 2 (ACE2). Although the Cryo-electron microscopy (Cryo-EM) uncovered the static atomic details of ACE2 homodimers, there is still a lack of research on the kinetic and thermodynamic properties of these full-length structures. This information is helpful to understand and interpret the role of ACE2 in the cell entry of SARS-CoV-2. In order to obtain this information, we performed microsecond-scale conventional and accelerated molecular dynamics (MD) simulations of full-length all-atomic systems of the RBD-ACE2 complex, the normal and torsional conformations of the apo-ACE2 homodimer. The comparative analysis of these systems showed that there were differences in their allosteric signal pathways and motion trends. These results may be helpful to further explore the cell entry mechanism of SARS-CoV-2. Moreover, the binding free energy and hydrogen bond distribution analysis of RBD-ACE2 binding interface provided the binding motifs that may be critical to allosteric signal transmission and RBD binding. These multi-conformational binding motifs can be used as targets or templates for the inhibitor design of the cell entry of SARS-CoV-2.

An investigation into the allosteric mechanism of GPCR A2A adenosine receptor with trajectory-based information theory and complex network model

G protein-coupled receptors (GPCRs), a large superfamily of transmembrane (TM) proteins, allosterically transduce the signal of ligand binding in the extracellular (EC) domain to couple to effector proteins in the intracellular (IC) domain, therefore forming the largest class of drug targets. The A2A adenosine receptor (A2AAR), a class-A GPCR, has been extensively studied as it offers numerous possibilities for therapeutic applications. However, the mechanism of allosteric communication between EC and IC domains is not completely clear. In this work, we utilize torsional mutual information to quantify the correlated motions of residue pairs from its molecular dynamics (MD) simulation trajectories, and further use the complex network model to obtain allosteric pipelines and hubs. The identified allosteric communication pipelines mainly transmit the signal from EC domain to the cytoplasmic ends of TM helix 5 (TM5), TM6 and TM7. The allosteric hubs, mostly located at TM5, TM6 and TM7, play an important role in mediating allosteric signal transmission to keep the receptor rigid and prevent G protein from binding to IC domain, which can explain the reason why their mutations distant from ligand-binding site do not affect the ligand binding affinity but affect the ligand efficacy. Additionally, we identify the key residues located in antagonist ZM241385 binding pocket which mediate multiple allosteric pathways and have been experimentally proven to play a critical role in affecting the ligand potency. This study is helpful for understanding the allosteric communication mechanism of A2AAR, and can provide valuable information for the structure-based drug design of GPCRs. Communicated by Ramaswamy H. Sarma

Predicted Hotspot Residues Involved in Allosteric Signal Transmission in Pro-Apoptotic Peptide-Mcl1 Complexes.

Mcl1 is a primary member of the Bcl–2 family—anti–apoptotic proteins (AAP)—that is overexpressed in several cancer pathologies. The apoptotic regulation is mediated through the binding of pro-apoptotic peptides (PAPs) (e.g., Bak and Bid) at the canonical hydrophobic binding groove (CBG) of Mcl1. Although all PAPs form amphipathic α-helices, their amino acid sequences vary to different degree. This sequence variation exhibits a central role in the binding partner selectivity towards different AAPs. Thus, constructing a novel peptide or small organic molecule with the ability to mimic the natural regulatory process of PAP is essential to inhibit various AAPs. Previously reported experimental binding free energies (BFEs) were utilized in the current investigation aimed to understand the mechanistic basis of different PAPs targeted to mMcl1. Molecular dynamics (MD) simulations used to estimate BFEs between mMcl1—PAP complexes using Molecular Mechanics-Generalized Born Solvent Accessible (MMGBSA) approach with multiple parameters. Predicted BFE values showed an excellent agreement with the experiment (R2 = 0.92). The van–der Waals (ΔGvdw) and electrostatic (ΔGele) energy terms found to be the main energy components that drive heterodimerization of mMcl1—PAP complexes. Finally, the dynamic network analysis predicted the allosteric signal transmission pathway involves more favorable energy contributing residues. In total, the results obtained from the current investigation may provide valuable insights for the synthesis of a novel peptide or small organic inhibitor targeting Mcl1.

Exploring Mechanisms of Communication Switching in the Hsp90-Cdc37 Regulatory Complexes with Client Kinases through Allosteric Coupling of Phosphorylation Sites: Perturbation-Based Modeling and Hierarchical Community Analysis of Residue Interaction Networks.

Understanding molecular principles underlying chaperone-based modulation of kinase client activity is critically important to dissect functions and activation mechanisms of many oncogenic proteins. The recent experimental studies have suggested that phosphorylation sites in the Hsp90 and Cdc37 proteins can serve as conformational communication switches of chaperone regulation and kinase interactions. However, a mechanism of allosteric coupling between phosphorylation sites in the Hsp90 and Cdc37 during client binding is poorly understood, and the molecular signatures underpinning specific roles of phosphorylation sites in the Hsp90 regulation remain unknown. In this work, we employed a combination of evolutionary analysis, coarse-grained molecular simulations together with perturbation-based network modeling and scanning of the unbound and bound Hsp90 and Cdc37 structures to quantify allosteric effects of phosphorylation sites and identify unique signatures that are characteristic for communication switches of kinase-specific client binding. By using network-based metrics of the dynamic intercommunity bridgeness and community centrality, we characterize specific signatures of phosphorylation switches involved in allosteric regulation. Through perturbation-based analysis of the dynamic residue interaction networks, we show that mutations of kinase-specific phosphorylation switches can induce long-range effects and lead to a global rewiring of the allosteric network and signal transmission in the Hsp90-Cdc37-kinase complex. We determine a specific group of phosphorylation sites in the Hsp90 where mutations may have a strong detrimental effect on allosteric interaction network, providing insight into the mechanism of phosphorylation-induced communication switching. The results demonstrate that kinase-specific phosphorylation switches of communications in the Hsp90 may be partly predisposed for their regulatory role based on preexisting allosteric propensities.

A study on allosteric communication in U1A-snRNA binding interactions: network analysis combined with molecular dynamics data

The allosteric regulation during the binding interactions between small nuclear RNAs (snRNAs) and the associated protein factors is critical to the function of spliceosomes in alternative RNA splicing. Although network models combined with molecular dynamics simulations have shown to be powerful tools for the analysis of protein allostery, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems. In this work, we use a residual network model combined with a coarse-grained Gaussian network model (GNM) to investigate the binding interactions between the snRNA and the human U1A protein which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle, and to identify the residues that play an important role in the allosteric communication in U1A during this process. We also utilize the Girvan-Newman method to detect the structural organization in U1A-snRNA recognition and interactions. Our results reveal that: (Ι) not only the residues at the binding sites that are traditionally considered to play a major role in U1A-snRNA association, but those residues that are far away from the RNA binding interface participate in the U1A's allosteric signal transmission induced by the RNA binding; (Π) the structure of U1A protein is well organized with different communities acting different roles for its RNA binding and allosteric regulation. The study demonstrates that the combination of the residual network and elastic network models is an effective and efficient method which can be readily extended to the investigation of the allosteric communication for other macromolecular interaction systems.

Local Ordering at the N-H Sites of the Rho GTPase Binding Domain of Plexin-B1: Impact of Dimerization.

We have developed a new molecular dynamics (MD) based method for describing analytically local potentials at mobile N-H sites in proteins. Here we apply it to the monomer and dimer of the Rho GTPase binding domain (RBD) of the transmembrane receptor plexin-B1 to gain insight into dimerization, which can compete with Rho GTPase binding. In our method, the local potential is given by linear combinations, u(DL,K), of the real combinations of the Wigner rotation matrix elements, DL,K, with L = 1-4 and appropriate symmetry. The combination that "fits best" the corresponding MD potential of mean force, u(MD), is the potential we are seeking, u(DL,K-BEST). For practical reasons the fitting process involves probability distributions, Peq ∝ exp(-u), instead of potentials, u. The symmetry of the potential, u(DL,K), may be related to the irreducible representations of the D2h point group. The monomer (dimer) potentials have mostly Ag and B2u (B1u and B2u) symmetry. For the monomer, the associated probability distributions are generally dispersed in space, shallow, and centered at the "reference N-H orientation" (defined in section 3.1. below); for the dimer many are more concentrated, deep and centered away from the "reference N-H orientation". The u(DL,K) functions provide a consistent description of the potential energy landscape at protein N-H sites. The L1-loop of the plexin-B1 RBD is not seen in the crystal structure, and many resonances of the L4 loop are missing in the NMR 15N-1H HSQC spectrum of the dimer; we suggest reasons for these features. An allosteric signal transmission pathway was reported previously for the monomer. We find that it has shallow N-H potentials at its ends, which become deeper as one proceeds toward the middle, complementing structurally the previously derived dynamic picture. Prospects of this study include correlating u(DL,K-BEST) with MD force-fields, and using them without further adjustment in NMR relaxation analysis schemes.

Hinge Twists and Population Shifts Deliver Regulated Catalysis for ATP-PRT in Histidine Biosynthesis

Allosteric regulation plays an important role in the control of metabolic flux in biosynthetic pathways. In microorganisms, many enzymes in these pathways adopt different strategies of allostery to allow the tuning of their activities in response to metabolic demand. Thus, it is important to uncover the mechanism of allosteric signal transmission to fully comprehend the complex control of enzyme function and its evolution. ATP-phosphoribosyltransferase (ATP-PRT), as the first enzyme in the histidine biosynthetic pathway, is allosterically regulated by histidine and offers a good platform for the study of allostery. Two forms of ATP-PRT, namely long and short forms, were discovered that show different arrangements of their regulatory machinery. Crystal structures of the long-form ATP-PRT have revealed overall conformational changes in the inhibited state, but the observed changes in the active state are quite subtle, making the elucidation of its allosteric mechanism difficult. Here, we combine computational methods (ligand docking, quantum mechanics/molecular mechanics optimization, and molecular dynamic simulations) with experimental studies to probe the signal transmission between remote allosteric and active sites. Our results reveal that distinct conformational ensembles of the catalytic domain with different dynamic properties exist in the ligand-free and histidine-bound enzymes. These ensembles display different capabilities in supporting the catalytic and allosteric function of ATP-PRT. The findings give insight into the underlying mechanism of allostery and allow us to propose that the hinge twisting within the catalytic domain is the key for both enhancement of catalysis and provision of regulation in ATP-PRT enzymes.

Interpreting the Dynamics of Binding Interactions of snRNA and U1A Using a Coarse-Grained Model

The binding interactions of small nuclear RNAs (snRNA) and the associated protein factors are critical to the function of spliceosomes in alternatively splicing primary RNA transcripts. Although molecular dynamics simulations are a powerful tool to interpret the mechanism of biological processes, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems, such as snRNA-protein complexes. We extend the coarse-grained Gaussian network model, which models the RNA-protein complexes as a harmonic chain of Cα, P, and O4′ atoms, to investigating the impact of the snRNA-binding interaction on the dynamic stability of the human U1A protein, which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle. The results reveal that the first and third loops and the C-terminal helix regions of the U1A domain undergo a significant loss of flexibility upon the RNA binding due to the forming of mostly electrostatic and hydrogen bond interactions with RNA 5′ stem and loop. By examining the residues whose mutations significantly change the binding free energy between U1A and snRNA, the Gaussian network model-based calculations show that not only the residues at the binding sites that are traditionally considered to play a major role in U1A-RNA association but also those residues that are far away from the RNA-binding interface can participate in the long-range allosteric signal transmission; these calculations are quantitatively consistent with the data observed in the recent snRNA binding experiments. The study demonstrates a useful avenue to utilize the simplified elastic network model to investigate the dynamics characteristics of the biologically important macromolecular interactions.