Remote loop evolution reveals a complex biological function for chitinase enzymes beyond the active site.

Loops are small secondary structural elements that play a crucial role in the emergence of new enzyme functions. However, the evolutionary molecular mechanisms how proteins acquire these loop elements and obtain new function is poorly understood. To address this question, we study glycoside hydrolase family 19 (GH19) chitinase-an essential enzyme family for pathogen degradation in plants. By revealing the evolutionary history and loops appearance of GH19 chitinase, we discover that one loop which is remote from the catalytic site, is necessary to acquire the new antifungal activity. We demonstrate that this remote loop directly accesses the fungal cell wall, and surprisingly, it needs to adopt a defined structure supported by long-range intramolecular interactions to perform its function. Our findings prove that nature applies this strategy at the molecular level to achieve a complex biological function while maintaining the original activity in the catalytic pocket, suggesting an alternative way to design new enzyme function.

Mapping the conformational landscape of the stimulatory heterotrimeric G protein.

Heterotrimeric G proteins serve as membrane-associated signaling hubs, in concert with their cognate G-protein-coupled receptors. Fluorine nuclear magnetic resonance spectroscopy was employed to monitor the conformational equilibria of the human stimulatory G-protein α subunit (Gsα) alone, in the intact Gsαß1γ2 heterotrimer or in complex with membrane-embedded human adenosine A2A receptor (A2AR). The results reveal a concerted equilibrium that is strongly affected by nucleotide and interactions with the ßγ subunit, the lipid bilayer and A2AR. The α1 helix of Gsα exhibits significant intermediate timescale dynamics. The α4ß6 loop and α5 helix undergo membrane/receptor interactions and order-disorder transitions respectively, associated with G-protein activation. The αN helix adopts a key functional state that serves as an allosteric conduit between the ßγ subunit and receptor, while a significant fraction of the ensemble remains tethered to the membrane and receptor upon activation.

Mechanistic insights into the nickel-dependent allosteric response of the Helicobacter pylori NikR transcription factor.

In Helicobacter pylori, the nickel-responsive NikR transcription factor plays a key role in regulating intracellular nickel concentrations, which is an essential process for survival of this pathogen in the acidic human stomach. Nickel binding to H. pylori NikR (HpNikR) allosterically activates DNA binding to target promoters encoding genes involved in nickel homeostasis and acid adaptation, to either activate or repress their transcription. We previously showed that HpNikR adopts an equilibrium between an open conformation and DNA-binding competent cis and trans states. Nickel binding slows down conformational exchange between these states and shifts the equilibrium toward the binding-competent states. The protein then becomes stabilized in a cis conformation upon binding the ureA promoter. Here, we investigate how nickel binding creates this response and how it is transmitted to the DNA-binding domains. Through mutagenesis, DNA-binding studies, and computational methods, the allosteric response to nickel was found to be propagated from the nickel-binding sites to the DNA-binding domains via the ß-sheets of the metal-binding domain and a network of residues at the inter-domain interface. Our computational results suggest that nickel binding increases protein rigidity to slow down the conformational exchange. A thymine base in the ureA promoter sequence, known to be critical for high affinity DNA binding by HpNikR, was also found to be important for the allosteric response, while a modified version of this promoter further highlighted the importance of the DNA sequence in modulating the response. Collectively, our results provide insights into regulation of a key protein for H. pylori survival.

Probing Conformational Dynamics of Antibodies with Geometric Simulations.

This chapter describes the application of constrained geometric simulations for prediction of antibody structural dynamics. We utilize constrained geometric simulations method FRODAN, which is a low computational complexity alternative to molecular dynamics (MD) simulations that can rapidly explore flexible motions in protein structures. FRODAN is highly suited for conformational dynamics analysis of large proteins, complexes, intrinsically disordered proteins, and dynamics that occurs on longer biologically relevant time scales that are normally inaccessible to classical MD simulations. This approach predicts protein dynamics at an all-atom scale while retaining realistic covalent bonding, maintaining dihedral angles in energetically good conformations while avoiding steric clashes in addition to performing other geometric and stereochemical criteria checks. In this chapter, we apply FRODAN to showcase its applicability for probing functionally relevant dynamics of IgG2a, including large-amplitude domain-domain motions and motions of complementarity determining region (CDR) loops. As was suggested in previous experimental studies, our simulations show that antibodies can explore a large range of conformational space.

Low pH structure of heliorhodopsin reveals chloride binding site and intramolecular signaling pathway.

Within the microbial rhodopsin family, heliorhodopsins (HeRs) form a phylogenetically distinct group of light-harvesting retinal proteins with largely unknown functions. We have determined the 1.97 Å resolution X-ray crystal structure of Thermoplasmatales archaeon SG8-52-1 heliorhodopsin (TaHeR) in the presence of NaCl under acidic conditions (pH 4.5), which complements the known 2.4 Å TaHeR structure acquired at pH 8.0. The low pH structure revealed that the hydrophilic Schiff base cavity (SBC) accommodates a chloride anion to stabilize the protonated retinal Schiff base when its primary counterion (Glu-108) is neutralized. Comparison of the two structures at different pH revealed conformational changes connecting the SBC and the extracellular loop linking helices A-B. We corroborated this intramolecular signaling transduction pathway with computational studies, which revealed allosteric network changes propagating from the perturbed SBC to the intracellular and extracellular space, suggesting TaHeR may function as a sensory rhodopsin. This intramolecular signaling mechanism may be conserved among HeRs, as similar changes were observed for HeR 48C12 between its pH 8.8 and pH 4.3 structures. We additionally performed DEER experiments, which suggests that TaHeR forms possible dimer-of-dimer associations which may be integral to its putative functionality as a light sensor in binding a transducer protein.

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex.

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

Allosteric modulation of the adenosine A2A receptor by cholesterol.

Cholesterol is a major component of the cell membrane and commonly regulates membrane protein function. Here, we investigate how cholesterol modulates the conformational equilibria and signaling of the adenosine A2A receptor (A2AR) in reconstituted phospholipid nanodiscs. This model system conveniently excludes possible effects arising from cholesterol-induced phase separation or receptor oligomerization and focuses on the question of allostery. GTP hydrolysis assays show that cholesterol weakly enhances the basal signaling of A2AR while decreasing the agonist EC50. Fluorine nuclear magnetic resonance (19F NMR) spectroscopy shows that this enhancement arises from an increase in the receptor's active state population and a G-protein-bound precoupled state. 19F NMR of fluorinated cholesterol analogs reveals transient interactions with A2AR, indicating a lack of high-affinity binding or direct allosteric modulation. The combined results suggest that the observed allosteric effects are largely indirect and originate from cholesterol-mediated changes in membrane properties, as shown by membrane fluidity measurements and high-pressure NMR.

Mutation-Induced Long-Range Allosteric Interactions in the Spike Protein Determine the Infectivity of SARS-CoV-2 Emerging Variants.

The emergence of a variety of highly transmissible SARS-CoV-2 variants, the causative agent of COVID-19, with multiple spike mutations poses serious challenges in overcoming the ongoing deadly pandemic. It is, therefore, essential to understand how these variants gain enhanced ability to evade immune responses with a higher rate of spreading infection. To address this question, here we have individually assessed the effects of SARS-CoV-2 variant-specific spike (S) protein receptor-binding domain (RBD) mutations E484K, K417N, L452Q, L452R, N501Y, and T478K that characterize and differentiate several emerging variants. Despite the hundreds of apparently neutral mutations observed in the domains other than the RBD, we have shown that each RBD mutation site is differentially engaged in an interdomain allosteric network involving mutation sites from a distant domain, affecting interactions with the human receptor angiotensin-converting enzyme-2 (ACE2). This allosteric network couples the residues of the N-terminal domain (NTD) and the RBD, which are modulated by the RBD-specific mutations and are capable of propagating mutation-induced perturbations between these domains through a combination of structural changes and effector-dependent modulations of dynamics. One key feature of this network is the inclusion of compensatory mutations segregated into three characteristically different clusters, where each cluster residue site is allosterically coupled with specific RBD mutation sites. Notably, each RBD mutation acted like a positive allosteric modulator; nevertheless, K417N was shown to have the largest effects among all of the mutations on the allostery and thereby holds the highest binding affinity with ACE2. This result will be useful for designing the targeted control measure and therapeutic efforts aiming at allosteric modulators.

The accuracy of NMR protein structures in the Protein Data Bank.

The program ANSURR measures the accuracy of NMR structures by comparing rigidity obtained from experimental backbone chemical shifts and from structures. We report on ANSURR analysis of 7,000 PDB NMR ensembles within the Protein Data Bank, which can be found at ansurr.com. The accuracy of NMR structures progressively improved up until 2005, but since then, it has plateaued. Most structures have accurate secondary structure, but are generally too floppy, particularly in loops. Thus, there is a need for more experimental restraints in loops. Currently, the best predictors of accuracy are Ramachandran distribution and the number of NOE restraints per residue. The precision of structures within the ensemble correlates well with accuracy, as does the number of hydrogen bond restraints per residue. Structure accuracy is improved when other components (such as additional polypeptide chains or ligands) are included.

Delineating the conformational landscape of the adenosine A2A receptor during G protein coupling.

G-protein-coupled receptors (GPCRs) represent a ubiquitous membrane protein family and are important drug targets. Their diverse signaling pathways are driven by complex pharmacology arising from a conformational ensemble rarely captured by structural methods. Here, fluorine nuclear magnetic resonance spectroscopy (19F NMR) is used to delineate key functional states of the adenosine A2A receptor (A2AR) complexed with heterotrimeric G protein (Gαsß1γ2) in a phospholipid membrane milieu. Analysis of A2AR spectra as a function of ligand, G protein, and nucleotide identifies an ensemble represented by inactive states, a G-protein-bound activation intermediate, and distinct nucleotide-free states associated with either partial- or full-agonist-driven activation. The Gßγ subunit is found to be critical in facilitating ligand-dependent allosteric transmission, as shown by 19F NMR, biochemical, and computational studies. The results provide a mechanistic basis for understanding basal signaling, efficacy, precoupling, and allostery in GPCRs.

Probing Allosteric Mechanism with Long-Range Rigidity Transmission Across Protein Networks.

Allosteric transmission refers to regulation of protein function at a distance. "Allostery" involves regulation and/or signal transduction induced by a perturbation event. Allostery, which has been coined the "second secret of life," is a fundamental property of most dynamics proteins. Most of critical questions surrounding allostery are largely unresolved. One of the key puzzles is to describe the physical mechanism of distant coupled conformational change. Another hot research area surrounding allostery is detection of allosteric pathways or regions (residues) in the protein that are the most critical for transmission of allosteric information. Using techniques inspired by mathematical rigidity theory and mechanical linkages, we have previously proposed a mechanistic model and description of allosteric transmission and an accompanying computational method, the Rigidity Transmission Allostery (RTA) algorithm. The RTA algorithm and method are designed to predict if mechanical perturbation of rigidity, for example, due to ligand binding, at one site of the protein can transmit and propagate across a protein structure and in turn cause a change in available conformational degrees of freedom and a change in conformation at a second distant site, equivalently resulting in allosteric transmission. The RTA algorithm is computationally very fast and can rapidly scan many unknown sites for allosteric transmission, identifying potential novel allosteric sites and quantify their allosteric effect. In this chapter we will discuss the rigidity-based mechanistic model of allosteric communication. As a case illustrative study, we will demonstrate RTA analysis on a G protein coupled receptor (GPCR) human adenosine A2A receptor. Our method gives important implications and a novel prospective for general mechanistic description of allosteric communication.

A method for validating the accuracy of NMR protein structures.

We present a method that measures the accuracy of NMR protein structures. It compares random coil index [RCI] against local rigidity predicted by mathematical rigidity theory, calculated from NMR structures [FIRST], using a correlation score (which assesses secondary structure), and an RMSD score (which measures overall rigidity). We test its performance using: structures refined in explicit solvent, which are much better than unrefined structures; decoy structures generated for 89 NMR structures; and conventional predictors of accuracy such as number of restraints per residue, restraint violations, energy of structure, ensemble RMSD, Ramachandran distribution, and clashscore. Restraint violations and RMSD are poor measures of accuracy. Comparisons of NMR to crystal structures show that secondary structure is equally accurate, but crystal structures are typically too rigid in loops, whereas NMR structures are typically too floppy overall. We show that the method is a useful addition to existing measures of accuracy.

Allosteric regulation of glutamate dehydrogenase deamination activity.

Glutamate dehydrogenase (GDH) is a key enzyme interlinking carbon and nitrogen metabolism. Recent discoveries of the GDH specific role in breast cancer, hyperinsulinism/hyperammonemia (HI/HA) syndrome, and neurodegenerative diseases have reinvigorated interest on GDH regulation, which remains poorly understood despite extensive and long standing studies. Notwithstanding the growing evidence of the complexity of allosteric network behind GDH regulation, identifications of allosteric factors and associated mechanisms are paramount to deepen our understanding of the complex dynamics that regulate GDH enzymatic activity. Combining structural analyses of cryo-electron microscopy data with molecular dynamic simulations, here we show that the cofactor NADH is a key player in the GDH regulation process. Our structural analysis indicates that, binding to the regulatory sites in proximity of the antenna region, NADH acts as a positive allosteric modulator by enhancing both the affinity of the inhibitor GTP binding and inhibition of GDH catalytic activity. We further show that the binding of GTP to the NADH-bound GDH activates a triangular allosteric network, interlinking the inhibitor with regulatory and catalytic sites. This allostery produces a local conformational rearrangement that triggers an anticlockwise rotational motion of interlinked alpha-helices with specific tilted helical extension. This structural transition is a fundamental switch in the GDH enzymatic activity. It introduces a torsional stress, and the associated rotational shift in the Rossmann fold closes the catalytic cleft with consequent inhibition of the deamination process. In silico mutagenesis examinations further underpin the molecular basis of HI/HA dominant mutations and consequent over-activity of GDH through alteration of this allosteric communication network. These results shed new light on GDH regulation and may lay new foundation in the design of allosteric agents.

Feedforward Control of Plant Nitrate Transporter NRT1.1 Biphasic Adaptive Activity.

Defective nitrate signaling in plants causes disorder in nitrogen metabolism, and it negatively affects nitrate transport systems, which toggle between high- and low-affinity modes in variable soil nitrate conditions. Recent discovery of a plasma membrane nitrate transceptor protein NRT1.1-a transporter cum sensor-provides a clue on this toggling mechanism. However, the general mechanistic description still remains poorly understood. Here, we illustrate adaptive responses and regulation of NRT1.1-mediated nitrate signaling in a wide range of extracellular nitrate concentrations. The results show that the homodimeric structure of NRT1.1 and its dimeric switch play an important role in eliciting specific cytosolic calcium waves sensed by the calcineurin-B-like calcium sensor CBL9, which activates the kinase CIPK23, in low nitrate concentration that is, however, impeded in high nitrate concentration. Nitrate binding at the high-affinity unit initiates NRT1.1 dimer decoupling and priming of the Thr101 site for phosphorylation by CIPK23. This phosphorylation stabilizes the NRT1.1 monomeric state, acting as a high-affinity nitrate transceptor. However, nitrate binding in both monomers, retaining the unmodified NRT1.1 state through dimerization, attenuates CIPK23 activity and thereby maintains the low-affinity mode of nitrate signaling and transport. This phosphorylation-led modulation of NRT1.1 activity shows bistable behavior controlled by an incoherent feedforward loop, which integrates nitrate-induced positive and negative regulatory effects on CIPK23. These results, therefore, advance our molecular understanding of adaptation in fluctuating nutrient availability and are a way forward for improving plant nitrogen use efficiency.

Substrate-Based Allosteric Regulation of a Homodimeric Enzyme.

Many enzymes operate through half-of-the sites reactivity wherein a single protomer is catalytically engaged at one time. In the case of the homodimeric enzyme, fluoroacetate dehalogenase, substrate binding triggers closing of a regulatory cap domain in the empty protomer, preventing substrate access to the remaining active site. However, the empty protomer serves a critical role by acquiring more disorder upon substrate binding, thereby entropically favoring the forward reaction. Empty protomer dynamics are also allosterically coupled to the bound protomer, driving conformational exchange at the active site and progress along the reaction coordinate. Here, we show that at high concentrations, a second substrate binds along the substrate-access channel of the occupied protomer, thereby dampening interprotomer dynamics and inhibiting catalysis. While a mutation (K152I) abrogates second site binding and removes inhibitory effects, it also precipitously lowers the maximum catalytic rate, implying a role for the allosteric pocket at low substrate concentrations, where only a single substrate engages the enzyme at one time. We show that this outer pocket first desolvates the substrate, whereupon it is deposited in the active site. Substrate binding to the active site then triggers the empty outer pocket to serve as an interprotomer allosteric conduit, enabling enhanced dynamics and sampling of activation states needed for catalysis. These allosteric networks and the ensuing changes resulting from second substrate binding are delineated using rigidity-based allosteric transmission theory and validated by nuclear magnetic resonance and functional studies. The results illustrate the role of dynamics along allosteric networks in facilitating function.

Mechanistic insights into allosteric regulation of the A2A adenosine G protein-coupled receptor by physiological cations.

Cations play key roles in regulating G-protein-coupled receptors (GPCRs), although their mechanisms are poorly understood. Here, 19F NMR is used to delineate the effects of cations on functional states of the adenosine A2A GPCR. While Na+ reinforces an inactive ensemble and a partial-agonist stabilized state, Ca2+ and Mg2+ shift the equilibrium toward active states. Positive allosteric effects of divalent cations are more pronounced with agonist and a G-protein-derived peptide. In cell membranes, divalent cations enhance both the affinity and fraction of the high affinity agonist-bound state. Molecular dynamics simulations suggest high concentrations of divalent cations bridge specific extracellular acidic residues, bringing TM5 and TM6 together at the extracellular surface and allosterically driving open the G-protein-binding cleft as shown by rigidity-transmission allostery theory. An understanding of cation allostery should enable the design of allosteric agents and enhance our understanding of GPCR regulation in the cellular milieu.

Repertoire Analysis of Antibody CDR-H3 Loops Suggests Affinity Maturation Does Not Typically Result in Rigidification.

Antibodies can rapidly evolve in specific response to antigens. Affinity maturation drives this evolution through cycles of mutation and selection leading to enhanced antibody specificity and affinity. Elucidating the biophysical mechanisms that underlie affinity maturation is fundamental to understanding B-cell immunity. An emergent hypothesis is that affinity maturation reduces the conformational flexibility of the antibody's antigen-binding paratope to minimize entropic losses incurred upon binding. In recent years, computational and experimental approaches have tested this hypothesis on a small number of antibodies, often observing a decrease in the flexibility of the complementarity determining region (CDR) loops that typically comprise the paratope and in particular the CDR-H3 loop, which contributes a plurality of antigen contacts. However, there were a few exceptions and previous studies were limited to a small handful of cases. Here, we determined the structural flexibility of the CDR-H3 loop for thousands of recent homology models of the human peripheral blood cell antibody repertoire using rigidity theory. We found no clear delineation in the flexibility of naïve and antigen-experienced antibodies. To account for possible sources of error, we additionally analyzed hundreds of human and mouse antibodies in the Protein Data Bank through both rigidity theory and B-factor analysis. By both metrics, we observed only a slight decrease in the CDR-H3 loop flexibility when comparing affinity matured antibodies to naïve antibodies, and the decrease was not as drastic as previously reported. Further analysis, incorporating molecular dynamics simulations, revealed a spectrum of changes in flexibility. Our results suggest that rigidification may be just one of many biophysical mechanisms for increasing affinity.

Suppressing allostery in epitope mapping experiments using millisecond hydrogen / deuterium exchange mass spectrometry.

Localization of the interface between the candidate antibody and its antigen target, commonly known as epitope mapping, is a critical component of the development of therapeutic monoclonal antibodies. With the recent availability of commercial automated systems, hydrogen / deuterium eXchange (HDX) is rapidly becoming the tool for mapping epitopes preferred by researchers in both industry and academia. However, this approach has a significant drawback in that it can be confounded by 'allosteric' structural and dynamic changes that result from the interaction, but occur far from the point(s) of contact. Here, we introduce a 'kinetic' millisecond HDX workflow that suppresses allosteric effects in epitope mapping experiments. The approach employs a previously introduced microfluidic apparatus that enables millisecond HDX labeling times with on-chip pepsin digestion and electrospray ionization. The 'kinetic' workflow also differs from conventional HDX-based epitope mapping in that the antibody is introduced to the antigen at the onset of HDX labeling. Using myoglobin / anti-myoglobin as a model system, we demonstrate that at short 'kinetic' workflow labeling times (i.e., 200 ms), the HDX signal is already fully developed at the 'true' epitope, but is still largely below the significance threshold at allosteric sites. Identification of the 'true' epitope is supported by computational docking predictions and allostery modeling using the rigidity transmission allostery algorithm.

The role of dimer asymmetry and protomer dynamics in enzyme catalysis.

Freeze-trapping x-ray crystallography, nuclear magnetic resonance, and computational techniques reveal the distribution of states and their interconversion rates along the reaction pathway of a bacterial homodimeric enzyme, fluoroacetate dehalogenase (FAcD). The crystal structure of apo-FAcD exhibits asymmetry around the dimer interface and cap domain, priming one protomer for substrate binding. This asymmetry is dynamically averaged through conformational exchange on a millisecond time scale. During catalysis, the protomer conformational exchange rate becomes enhanced, the empty protomer exhibits increased local disorder, and water egresses. Computational studies identify allosteric pathways between protomers. Water release and enhanced dynamics associated with catalysis compensate for entropic losses from substrate binding while facilitating sampling of the transition state. The studies provide insights into how substrate-coupled allosteric modulation of structure and dynamics facilitates catalysis in a homodimeric enzyme.

Hyperphosphorylation of intrinsically disordered tau protein induces an amyloidogenic shift in its conformational ensemble.

Tau is an intrinsically disordered protein (IDP) whose primary physiological role is to stabilize microtubules in neuronal axons at all stages of development. In Alzheimer's and other tauopathies, tau forms intracellular insoluble amyloid aggregates known as neurofibrillary tangles, a process that appears in many cases to be preceded by hyperphosphorylation of tau monomers. Understanding the shift in conformational bias induced by hyperphosphorylation is key to elucidating the structural factors that drive tau pathology, however, as an IDP, tau is not amenable to conventional structural characterization. In this work, we employ a straightforward technique based on Time-Resolved ElectroSpray Ionization Mass Spectrometry (TRESI-MS) and Hydrogen/Deuterium Exchange (HDX) to provide a detailed picture of residual structure in tau, and the shifts in conformational bias induced by hyperphosphorylation. By comparing the native and hyperphosphorylated ensembles, we are able to define specific conformational biases that can easily be rationalized as enhancing amyloidogenic propensity. Representative structures for the native and hyperphosphorylated tau ensembles were generated by refinement of a broad sample of conformations generated by low-computational complexity modeling, based on agreement with the TRESI-HDX profiles.