Mutagenesis Data Research Articles

Binding stoichiometry and structural model of the HIV-1 Rev/importin β complex.

HIV-1 Rev mediates the nuclear export of intron-containing viral RNA transcripts and is essential for viral replication. Rev is imported into the nucleus by the host protein importin β (Impβ), but how Rev associates with Impβ is poorly understood. Here, we report biochemical, mutational, and biophysical studies of the Impβ/Rev complex. We show that Impβ binds two Rev monomers through independent binding sites, in contrast to the 1:1 binding stoichiometry observed for most Impβ cargos. Peptide scanning data and charge-reversal mutations identify the N-terminal tip of Rev helix α2 within Rev's arginine-rich motif (ARM) as a primary Impβ-binding epitope. Cross-linking mass spectrometry and compensatory mutagenesis data combined with molecular docking simulations suggest a structural model in which one Rev monomer binds to the C-terminal half of Impβ with Rev helix α2 roughly parallel to the HEAT-repeat superhelical axis, whereas the other monomer binds to the N-terminal half. These findings shed light on the molecular basis of Rev recognition by Impβ and highlight an atypical binding behavior that distinguishes Rev from canonical cellular Impβ cargos.

QM/MM Well-Tempered Metadynamics Study of the Mechanism of XBP1 mRNA Cleavage by Inositol Requiring Enzyme 1α RNase.

A range of in silico methodologies were herein employed to study the unconventional XBP1 mRNA cleavage mechanism performed by the unfolded protein response (UPR) mediator Inositol Requiring Enzyme 1α (IRE1). Using Protein–RNA molecular docking along with a series of extensive restrained/unrestrained atomistic molecular dynamics (MD) simulations, the dynamical behavior of the system was evaluated and a reliable model of the IRE1/XBP1 mRNA complex was constructed. From a series of well-converged quantum mechanics molecular mechanics well-tempered metadynamics (QM/MM WT-MetaD) simulations using the Grimme dispersion interaction corrected semiempirical parametrization method 6 level of theory (PM6-D3) and the AMBER14SB-OL3 force field, the free energy profile of the cleavage mechanism was determined, along with intermediates and transition state structures. The results show two distinct reaction paths based on general acid–general base type mechanisms, with different activation energies that perfectly match observations from experimental mutagenesis data. The study brings unique atomistic insights into the cleavage mechanism of XBP1 mRNA by IRE1 and clarifies the roles of the catalytic residues H910 and Y892. Increased understanding of the details in UPR signaling can assist in the development of new therapeutic agents for its modulation.

Hotspot Identification and Drug Design of Protein-Protein Interaction Modulators Using the Fragment Molecular Orbital Method.

Protein-protein interactions (PPIs) are essential for the function of many proteins. Aberrant PPIs have the potential to lead to disease, making PPIs promising targets for drug discovery. There are over 64,000 PPIs in the human interactome reference database; however, to date, very few PPI modulators have been approved for clinical use. Further development of PPI-specific therapeutics is highly dependent on the availability of structural data and the existence of reliable computational tools to explore the interface between two interacting proteins. The fragment molecular orbital (FMO) quantum mechanics method offers comprehensive and computationally inexpensive means of identifying the strength (in kcal/mol) and the chemical nature (electrostatic or hydrophobic) of the molecular interactions taking place at the protein-protein interface. We have integrated FMO and PPI exploration (FMO-PPI) to identify the residues that are critical for protein-protein binding (hotspots). To validate this approach, we have applied FMO-PPI to a dataset of protein-protein complexes representing several different protein subfamilies and obtained FMO-PPI results that are in agreement with published mutagenesis data. We observed that critical PPIs can be divided into three major categories: interactions between residues of two proteins (intermolecular), interactions between residues within the same protein (intramolecular), and interactions between residues of two proteins that are mediated by water molecules (water bridges). We extended our findings by demonstrating how this information obtained by FMO-PPI can be utilized to support the structure-based drug design of PPI modulators (SBDD-PPI).

Structural basis for tunable affinity and specificity of LxCxE-dependent protein interactions with the retinoblastoma protein family

The retinoblastoma protein (Rb) and its homologs p107 and p130 are critical regulators of gene expression during the cell cycle and are commonly inactivated in cancer. Rb proteins use their "pocket domain" to bind an LxCxE sequence motif in other proteins, many of which function with Rb proteins to co-regulate transcription. Here, we present binding data and crystal structures of the p107 pocket domain in complex with LxCxE peptides from the transcriptional co-repressor proteins HDAC1, ARID4A, and EID1. Our results explain why Rb and p107 have weaker affinity for cellular LxCxE proteins compared with the E7 protein from human papillomavirus, which has been used as the primary model for understanding LxCxE motif interactions. Our structural and mutagenesis data also identify and explain differences in Rb and p107 affinities for some LxCxE-containing sequences. Our study provides new insights into how Rb proteins bind their cell partners with varying affinity and specificity.

The Bacillus subtilis open reading frame ysgA encodes the SPOUT methyltransferase RlmP forming 2'-O-methylguanosine at position 2553 in the A-loop of 23S rRNA.

A previous bioinformatic analysis predicted that the ysgA open reading frame of Bacillus subtilis encodes an RNA methyltransferase of the SPOUT superfamily. Here we show that YsgA is the 2'-O-methyltransferase that targets position G2553 (Escherichia coli numbering) of the A-loop of 23S rRNA. This was shown by a combination of biochemical and mass spectrometry approaches using both rRNA extracted from B. subtilis wild-type or ΔysgA cells and in vitro synthesized rRNA. When the target G2553 is mutated, YsgA is able to methylate the ribose of adenosine. However, it cannot methylate cytidine nor uridine. The enzyme modifies free 23S rRNA but not the fully assembled ribosome nor the 50S subunit, suggesting that the modification occurs early during ribosome biogenesis. Nevertheless, ribosome subunits assembly is unaffected in a B. subtilis ΔysgA mutant strain. The crystal structure of the recombinant YsgA protein, combined with mutagenesis data, outlined in this article highlights a typical SPOUT fold preceded by an L7Ae/L30 (eL8/eL30 in a new nomenclature) amino-terminal domain.

Alchemical Free Energy Calculations to Investigate Protein-Protein Interactions: the Case of the CDC42/PAK1 Complex.

Here, we show that alchemical free energy calculations can quantitatively compute the effect of mutations at the protein–protein interface. As a test case, we have used the protein complex formed by the small Rho-GTPase CDC42 and its downstream effector PAK1, a serine/threonine kinase. Notably, the CDC42/PAK1 complex offers a wealth of structural, mutagenesis, and binding affinity data because of its central role in cellular signaling and cancer progression. In this context, we have considered 16 mutations in the CDC42/PAK1 complex and obtained excellent agreement between computed and experimental data on binding affinity. Importantly, we also show that a careful analysis of the side-chain conformations in the mutated amino acids can considerably improve the computed estimates, solving issues related to sampling limitations. Overall, this study demonstrates that alchemical free energy calculations can conveniently be integrated into the design of experimental mutagenesis studies.

Molecular mechanism of the wake-promoting agent TAK-925

The OX2 orexin receptor (OX2R) is a highly expressed G protein-coupled receptor (GPCR) in the brain that regulates wakefulness and circadian rhythms in humans. Antagonism of OX2R is a proven therapeutic strategy for insomnia drugs, and agonism of OX2R is a potentially powerful approach for narcolepsy type 1, which is characterized by the death of orexinergic neurons. Until recently, agonism of OX2R had been considered ‘undruggable.’ We harness cryo-electron microscopy of OX2R-G protein complexes to determine how the first clinically tested OX2R agonist TAK-925 can activate OX2R in a highly selective manner. Two structures of TAK-925-bound OX2R with either a Gq mimetic or Gi reveal that TAK-925 binds at the same site occupied by antagonists, yet interacts with the transmembrane helices to trigger activating microswitches. Our structural and mutagenesis data show that TAK-925’s selectivity is mediated by subtle differences between OX1 and OX2 receptor subtypes at the orthosteric pocket. Finally, differences in the polarity of interactions at the G protein binding interfaces help to rationalize OX2R’s coupling selectivity for Gq signaling. The mechanisms of TAK-925’s binding, activation, and selectivity presented herein will aid in understanding the efficacy of small molecule OX2R agonists for narcolepsy and other circadian disorders.

Structure and drug development of the human formylpeptide receptors FPR1 and FPR2

The formylpeptide receptors (FPRs) mediate pattern recognition of formylated peptides derived from invading pathogens or mitochondria from dead host cells. They can also sense other structurally distinct native peptides and even lipid mediators to either promote or resolve inflammation. Pharmacological targeting of FPRs represents a novel therapeutic approach in treating inflammatory diseases. However, the molecular mechanisms underlying FPR ligand recognition are elusive. We report cryo‐EM structures of Gi‐coupled FPR1 and FPR2 bound to a formylpeptide and Gi‐coupled FPR2 bound to two synthetic peptide and small‐molecule agonists. Together with mutagenesis data, our structures reveal the molecular mechanism of formylpeptide recognition by FPRs and structural variations of FPR1 and FPR2 leading to their different ligand preferences. Structural analysis also suggests that diverse FPR agonists sample a conserved activation chamber at the bottom of ligand‐binding pockets to activate FPRs. Based on our structures, we further screen several new chemicals as potential drugs. Our results provide a basis for structure‐based drug design on FPRs.

How clustered protocadherin binding specificity is tuned for neuronal self-/nonself-recognition.

The stochastic expression of fewer than 60 clustered protocadherin (cPcdh) isoforms provides diverse identities to individual vertebrate neurons and a molecular basis for self-/nonself-discrimination. cPcdhs form chains mediated by alternating cis and trans interactions between apposed membranes, which has been suggested to signal self-recognition. Such a mechanism requires that cPcdh cis dimers form promiscuously to generate diverse recognition units, and that trans interactions have precise specificity so that isoform mismatches terminate chain growth. However, the extent to which cPcdh interactions fulfill these requirements has not been definitively demonstrated. Here, we report biophysical experiments showing that cPcdh cis interactions are promiscuous, but with preferences favoring formation of heterologous cis dimers. Trans homophilic interactions are remarkably precise, with no evidence for heterophilic interactions between different isoforms. A new C-type cPcdh crystal structure and mutagenesis data help to explain these observations. Overall, the interaction characteristics we report for cPcdhs help explain their function in neuronal self-/nonself-discrimination.

Conformational Transitions in Yeast Chorismate Mutase Important for Allosteric Regulation as Identified by Nuclear Magnetic Resonance Spectroscopy

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.

Structural insights into ligand recognition, activation, and signaling of the α2A adrenergic receptor.

The α2A adrenergic receptor (α2AAR) is a G protein (heterotrimeric guanine nucleotide–binding protein)–coupled receptor that mediates important physiological functions in response to the endogenous neurotransmitters norepinephrine and epinephrine, as well as numerous chemically distinct drugs. However, the molecular mechanisms of drug actions remain poorly understood. Here, we report the cryo–electron microscopy structures of the human α2AAR-GoA complex bound to norepinephrine and three imidazoline derivatives (brimonidine, dexmedetomidine, and oxymetazoline). Together with mutagenesis and functional data, these structures provide important insights into the molecular basis of ligand recognition, activation, and signaling at the α2AAR. Further structural analyses uncover different molecular determinants between α2AAR and βARs for recognition of norepinephrine and key regions that determine the G protein coupling selectivity. Overall, our studies provide a framework for understanding the signal transduction of the adrenergic system at the atomic level, which will facilitate rational structure-based discovery of safer and more effective medications for α2AAR.

Force-Generation by the Trans-Envelope Tol-Pal System.

The Tol-Pal system spans the cell envelope of Gram-negative bacteria, transducing the potential energy of the proton motive force (PMF) into dissociation of the TolB-Pal complex at the outer membrane (OM), freeing the lipoprotein Pal to bind the cell wall. The primary physiological role of Tol-Pal is to maintain OM integrity during cell division through accumulation of Pal molecules at division septa. How the protein complex couples the PMF at the inner membrane into work at the OM is unknown. The effectiveness of this trans-envelope energy transduction system is underscored by the fact that bacteriocins and bacteriophages co-opt Tol-Pal as part of their import/infection mechanisms. Mechanistic understanding of this process has been hindered by a lack of structural data for the inner membrane TolQ-TolR stator, of its complexes with peptidoglycan (PG) and TolA, and of how these elements combined power events at the OM. Recent studies on the homologous stators of Ton and Mot provide a starting point for understanding how Tol-Pal works. Here, we combine ab initio protein modeling with previous structural data on sub-complexes of Tol-Pal as well as mutagenesis, crosslinking, co-conservation analysis and functional data. Through this composite pooling of in silico, in vitro, and in vivo data, we propose a mechanism for force generation in which PMF-driven rotary motion within the stator drives conformational transitions within a long TolA helical hairpin domain, enabling it to reach the TolB-Pal complex at the OM.

Molecular recognition of formylpeptides and diverse agonists by the formylpeptide receptors FPR1 and FPR2

The formylpeptide receptors (FPRs) mediate pattern recognition of formylated peptides derived from invading pathogens or mitochondria from dead host cells. They can also sense other structurally distinct native peptides and even lipid mediators to either promote or resolve inflammation. Pharmacological targeting of FPRs represents a novel therapeutic approach in treating inflammatory diseases. However, the molecular mechanisms underlying FPR ligand recognition are elusive. We report cryo-EM structures of Gi-coupled FPR1 and FPR2 bound to a formylpeptide and Gi-coupled FPR2 bound to two synthetic peptide and small-molecule agonists. Together with mutagenesis data, our structures reveal the molecular mechanism of formylpeptide recognition by FPRs and structural variations of FPR1 and FPR2 leading to their different ligand preferences. Structural analysis also suggests that diverse FPR agonists sample a conserved activation chamber at the bottom of ligand-binding pockets to activate FPRs. Our results provide a basis for rational drug design on FPRs.

Directed mutational scanning reveals a balance between acidic and hydrophobic residues in strong human activation domains.

SUMMARYAcidic activation domains are intrinsically disordered regions of the transcription factors that bind coactivators. The intrinsic disorder and low evolutionary conservation of activation domains have made it difficult to identify the sequence features that control activity. To address this problem, we designed thousands of variants in seven acidic activation domains and measured their activities with a high-throughput assay in human cell culture. We found that strong activation domain activity requires a balance between the number of acidic residues and aromatic and leucine residues. These findings motivated a predictor of acidic activation domains that scans the human proteome for clusters of aromatic and leucine residues embedded in regions of high acidity. This predictor identifies known activation domains and accurately predicts previously unidentified ones. Our results support a flexible acidic exposure model of activation domains in which the acidic residues solubilize hydrophobic motifs so that they can interact with coactivators. A record of this paper’s transparent peer review process is included in the supplemental information.

Predicting and interpreting large-scale mutagenesis data using analyses of protein stability and conservation.

Understanding and predicting the functional consequences of single amino acid changes is central in many areas of protein science. Here, we collect and analyze experimental measurements of effects of >150,000 variants in 29 proteins. We use biophysical calculations to predict changes in stability for each variant and assess them in light of sequence conservation. We find that the sequence analyses give more accurate prediction of variant effects than predictions of stability and that about half of the variants that show loss of function do so due to stability effects. We construct a machine learning model to predict variant effects from protein structure and sequence alignments and show how the two sources of information support one another and enable mechanistic interpretations. Together, our results show how one can leverage large-scale experimental assessments of variant effects to gain deeper and general insights into the mechanisms that cause loss of function.

Unveiling the Gating Mechanism of CRAC Channel: A Computational Study.

CRAC channel is ubiquitous and its importance in the regulation of the immune system is testified by the severe immunodeficiencies caused by its mutations. In this work we took advantage of the availability of open and closed structures of this channel to run for the first time simulations of the whole gating process reaching the relevant time-scale with an enhanced sampling technique, Targeted Molecular Dynamics. Our simulations highlighted a complex allosteric propagation of the conformational change from peripheral helices, where the activator STIM1 binds, to the central pore helices. In agreement with mutagenesis data, our simulations revealed the key role of residue H206 whose displacement creates an empty space behind the hydrophobic region of the pore, thus releasing a steric brake and allowing the opening of the channel. Conversely, the process of pore closing culminates with the formation of a bubble that occludes the pore even in the absence of steric block. This mechanism, known as “hydrophobic gating”, has been observed in an increasing number of biological ion channels and also in artificial nanopores. Our study therefore shows promise not only to better understand the molecular origin of diseases caused by disrupted calcium signaling, but also to clarify the mode of action of hydrophobically gated ion channels, possibly even suggesting strategies for the biomimetic design of synthetic nanopores.

The active site region plays a critical role in Na+ binding to thrombin

The catalytic activity of thrombin and other enzymes of the blood coagulation and complement cascades is enhanced significantly by binding of Na+ to a site >15 Å away from the catalytic residue S195, buried within the 180 and 220 loops that also contribute to the primary specificity of the enzyme. Rapid kinetics support a binding mechanism of conformational selection where the Na+-binding site is in equilibrium between open (N) and closed (N∗) forms and the cation binds selectively to the N form. Allosteric transduction of this binding step produces enhanced catalytic activity. Molecular details on how Na+ gains access to this site and communicates allosterically with the active site remain poorly defined. In this study, we show that the rate of the transition is strongly correlated with the analogous transition that governs the interaction of synthetic and physiologic substrates with the active site. This correlation supports the active site as the likely point of entry for Na+ to its binding site. Mutagenesis and structural data rule out an alternative path through the pore defined by the 180 and 220 loops. We suggest that the active site communicates allosterically with the Na+ site through a network of H-bonded water molecules that embeds the primary specificity pocket. Perturbation of the mobility of S195 and its H-bonding capabilities alters interaction with this network and influences the kinetics of Na+ binding and allosteric transduction. These findings have general mechanistic relevance for Na+-activated proteases and allosteric enzymes.

Structural Basis of Arrestin Selectivity for Active Phosphorylated G Protein-Coupled Receptors.

Arrestins are a small family of proteins that bind G protein-coupled receptors (GPCRs). Arrestin binds to active phosphorylated GPCRs with higher affinity than to all other functional forms of the receptor, including inactive phosphorylated and active unphosphorylated. The selectivity of arrestins suggests that they must have two sensors, which detect receptor-attached phosphates and the active receptor conformation independently. Simultaneous engagement of both sensors enables arrestin transition into a high-affinity receptor-binding state. This transition involves a global conformational rearrangement that brings additional elements of the arrestin molecule, including the middle loop, in contact with a GPCR, thereby stabilizing the complex. Here, we review structural and mutagenesis data that identify these two sensors and additional receptor-binding elements within the arrestin molecule. While most data were obtained with the arrestin-1-rhodopsin pair, the evidence suggests that all arrestins use similar mechanisms to achieve preferential binding to active phosphorylated GPCRs.

Structural basis of human ghrelin receptor signaling by ghrelin and the synthetic agonist ibutamoren

The hunger hormone ghrelin activates the ghrelin receptor GHSR to stimulate food intake and growth hormone secretion and regulate reward signaling. Acylation of ghrelin at Ser3 is required for its agonistic action on GHSR. Synthetic agonists of GHSR are under clinical evaluation for disorders related to appetite and growth hormone dysregulation. Here, we report high-resolution cryo-EM structures of the GHSR-Gi signaling complex with ghrelin and the non-peptide agonist ibutamoren as an investigational new drug. Our structures together with mutagenesis data reveal the molecular basis for the binding of ghrelin and ibutamoren. Structural comparison suggests a salt bridge and an aromatic cluster near the agonist-binding pocket as important structural motifs in receptor activation. Notable structural variations of the Gi and GHSR coupling are observed in our cryo-EM analysis. Our results provide a framework for understanding GHSR signaling and developing new GHSR agonist drugs.

Customized exogenous ferredoxin functions as an efficient electron carrier

Ferredoxin (Fdx) is regarded as the main electron carrier in biological electron transfer and acts as an electron donor in metabolic pathways of many organisms. Here, we screened a self-sufficient P450-derived reductase PRF with promising production yield of 9OHAD (9α-hydroxy4-androstene-3,17-dione) from AD, and further proved the importance of [2Fe–2S] clusters of ferredoxin-oxidoreductase in transferring electrons in steroidal conversion. The results of truncated Fdx domain in all oxidoreductases and mutagenesis data elucidated the indispensable role of [2Fe–2S] clusters in the electron transfer process. By adding the independent plant-type Fdx to the reaction system, the AD (4-androstene-3,17-dione) conversion rate have been significantly improved. A novel efficient electron transfer pathway of PRF + Fdx + KshA (KshA, Rieske-type oxygenase of 3-ketosteroid-9-hydroxylase) in the reaction system rather than KshAB complex system was proposed based on analysis of protein–protein interactions and redox potential measurement. Adding free Fdx created a new conduit for electrons to travel from reductase to oxygenase. This electron transfer pathway provides new insight for the development of efficient exogenous Fdx as an electron carrier.Graphical