Ensemble- and Rigidity Theory-Based Perturbation Approach To Analyze Dynamic Allostery

Molecular Mechanisms Underlying Single Nucleotide Polymorphism-Induced Reactivity Decrease in CYP2D6

Cytochrome P450 2D6 (CYP2D6) is one of the most important enzymes involved in drug metabolism. Genetic polymorphism can influence drug metabolism by CYP2D6 such that a therapy is seriously affected by under- or overdosing of drugs. However, a general explanation at the atomistic level for poor activity is missing so far. Here we show for the 20 most common single nucleotide polymorphisms (SNPs) of CYP2D6 that poor metabolism is driven by four mechanisms. We found in extensive all-atom molecular dynamics simulations that the rigidity of the I-helix (central helix), distance between central phenylalanines (stabilizing bound substrate), availability of basic residues on the surface of CYP2D6 (binding of cytochrome P450 reductase), and position of arginine 132 (electron transfer to heme) are essential for an extensive function of the enzyme. These results were applied to SNPs with unknown effects, and potential SNPs that may lead to poor drug metabolism were identified. The revealed molecular mechanisms might be important for other drug-metabolizing cytochrome P450 enzymes.

NRIMD, a Web Server for Analyzing Protein Allosteric Interactions Based on Molecular Dynamics Simulation

Long-range allosteric communication between distant sites and active sites in proteins is central to biological regulation but still poorly characterized, limiting the development of protein engineering and drug design. Addressing this gap, NRIMD is an open-access web server for analyzing long-range interactions in proteins from molecular dynamics (MD) simulations, such as the effect of mutations at distal sites or allosteric ligand binding at allosteric sites on the active center. Based on our recent works on neural relational inference using graph neural networks, this cloud-based web server accepts MD simulation data on any length of residues in the alpha-carbon skeleton format from mainstream MD software. The input trajectory data are validated at the frontend deployed on the cloud and then processed on the backend deployed on a high-performance computer system with a collection of complementary tools. The web server provides a one-stop-shop MD analysis platform to predict long-range interactions and their paths between distant sites and active sites. It provides a user-friendly interface for detailed analysis and visualization. To the best of our knowledge, NRIMD is the first-of-its-kind online service to provide comprehensive long-range interaction analysis on MD simulations, which significantly lowers the barrier of predictions on protein long-range interactions using deep learning. The NRIMD web server is publicly available at https://nrimd.luddy.indianapolis.iu.edu/.

Using deep learning and molecular dynamics simulations to unravel the regulation mechanism of peptides as noncompetitive inhibitor of xanthine oxidase

Xanthine oxidase (XO) is a crucial enzyme in the development of hyperuricemia and gout. This study focuses on LWM and ALPM, two food-derived inhibitors of XO. We used molecular docking to obtain three systems and then conducted 200 ns molecular dynamics simulations for the Apo, LWM, and ALPM systems. The results reveal a stronger binding affinity of the LWM peptide to XO, potentially due to increased hydrogen bond formation. Notable changes were observed in the XO tunnel upon inhibitor binding, particularly with LWM, which showed a thinner, longer, and more twisted configuration compared to ALPM. The study highlights the importance of residue F914 in the allosteric pathway. Methodologically, we utilized the perturbed response scan (PRS) based on Python, enhancing tools for MD analysis. These findings deepen our understanding of food-derived anti-XO inhibitors and could inform the development of food-based therapeutics for reducing uric acid levels with minimal side effects.

Regulation of STING activity in DNA sensing by ISG15 modification

Sensing of human immunodeficiency virus type 1 (HIV-1) DNA is mediated by the cyclic GMP-AMP synthase-stimulator of interferon genes (cGAS-STING) signaling axis. Signal transduction and regulation of this cascade is achieved by post-translational modifications. Here we show that cGAS-STING-dependent HIV-1 sensing requires interferon-stimulated gene 15 (ISG15). ISG15 deficiency inhibits STING-dependent sensing of HIV-1 and STING agonist-induced antiviral response. Upon external stimuli, STING undergoes ISGylation at residues K224, K236, K289, K347, K338, and K370. Inhibition of STING ISGylation at K289 suppresses STING-mediated type Ⅰ interferon induction by inhibiting its oligomerization. Of note, removal of STING ISGylation alleviates gain-of-function phenotype in STING-associated vasculopathy with onset in infancy (SAVI). Molecular modeling suggests that ISGylation of K289 is an important regulator of oligomerization. Taken together, our data demonstrate that ISGylation at K289 is crucial for STING activation and represents an important regulatory step in DNA sensing of viruses and autoimmune responses.

Combining structural and coevolution information to unveil allosteric sites

Understanding allosteric regulation in biomolecules is of great interest to pharmaceutical research and computational methods emerged during the last decades to characterize allosteric coupling. However, the prediction of allosteric sites in a protein structure remains a challenging task. Here, we integrate local binding site information, coevolutionary information, and information on dynamic allostery into a structure-based three-parameter model to identify potentially hidden allosteric sites in ensembles of protein structures with orthosteric ligands. When tested on five allosteric proteins (LFA-1, p38-α, GR, MAT2A, and BCKDK), the model successfully ranked all known allosteric pockets in the top three positions. Finally, we identified a novel druggable site in MAT2A confirmed by X-ray crystallography and SPR and a hitherto unknown druggable allosteric site in BCKDK validated by biochemical and X-ray crystallography analyses. Our model can be applied in drug discovery to identify allosteric pockets.

Uncoupling of Voltage- and Ligand-Induced Activation in HCN2 Channels by Glycine Inserts

Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are tetramers that generate electrical rhythmicity in special brain neurons and cardiomyocytes. The channels are activated by membrane hyperpolarization. The binding of cAMP to the four available cyclic nucleotide-binding domains (CNBD) enhances channel activation. We analyzed in the present study the mechanism of how the effect of cAMP binding is transmitted to the pore domain. Our strategy was to uncouple the C-linker (CL) from the channel core by inserting one to five glycine residues between the S6 gate and the A′-helix (constructs 1G to 5G). We quantified in full-length HCN2 channels the resulting functional effects of the inserted glycines by current activation as well as the structural dynamics and statics using molecular dynamics simulations and Constraint Network Analysis. We show functionally that already in 1G the cAMP effect on activation is lost and that with the exception of 3G and 5G the concentration-activation relationships are shifted to depolarized voltages with respect to HCN2. The strongest effect was found for 4G. Accordingly, the activation kinetics were accelerated by all constructs, again with the strongest effect in 4G. The simulations reveal that the average residue mobility of the CL and CNBD domains is increased in all constructs and that the junction between the S6 and A′-helix is turned into a flexible hinge, resulting in a destabilized gate in all constructs. Moreover, for 3G and 4G, there is a stronger downward displacement of the CL-CNBD than in HCN2 and the other constructs, resulting in an increased kink angle between S6 and A′-helix, which in turn loosens contacts between the S4-helix and the CL. This is suggested to promote a downward movement of the S4-helix, similar to the effect of hyperpolarization. In addition, exclusively in 4G, the selectivity filter in the upper pore region and parts of the S4-helix are destabilized. The results provide new insights into the intricate activation of HCN2 channels.

Neural relational inference to learn long-range allosteric interactions in proteins from molecular dynamics simulations

Protein allostery is a biological process facilitated by spatially long-range intra-protein communication, whereby ligand binding or amino acid change at a distant site affects the active site remotely. Molecular dynamics (MD) simulation provides a powerful computational approach to probe the allosteric effect. However, current MD simulations cannot reach the time scales of whole allosteric processes. The advent of deep learning made it possible to evaluate both spatially short and long-range communications for understanding allostery. For this purpose, we applied a neural relational inference model based on a graph neural network, which adopts an encoder-decoder architecture to simultaneously infer latent interactions for probing protein allosteric processes as dynamic networks of interacting residues. From the MD trajectories, this model successfully learned the long-range interactions and pathways that can mediate the allosteric communications between distant sites in the Pin1, SOD1, and MEK1 systems. Furthermore, the model can discover allostery-related interactions earlier in the MD simulation trajectories and predict relative free energy changes upon mutations more accurately than other methods.

Respiratory and C4-photosynthetic NAD-malic enzyme coexist in bundle sheath cell mitochondria and evolved via association of differentially adapted subunits

In plant mitochondria, nicotinamide adenine dinucleotide-malic enzyme (NAD-ME) has a housekeeping function in malate respiration. In different plant lineages, NAD-ME was independently co-opted in C4 photosynthesis. In the C4 Cleome species, Gynandropsis gynandra and Cleome angustifolia, all NAD-ME genes (NAD-MEα, NAD-MEβ1, and NAD-MEβ2) were affected by C4 evolution and are expressed at higher levels than their orthologs in the C3 species Tarenaya hassleriana. In T. hassleriana, the NAD-ME housekeeping function is performed by two heteromers, NAD-MEα/β1 and NAD-MEα/β2, with similar biochemical properties. In both C4 species, this role is restricted to NAD-MEα/β2. In the C4 species, NAD-MEα/β1 is exclusively present in the leaves, where it accounts for most of the enzymatic activity. Gynandropsis gynandra NAD-MEα/β1 (GgNAD-MEα/β1) exhibits high catalytic efficiency and is differentially activated by the C4 intermediate aspartate, confirming its role as the C4-decarboxylase. During C4 evolution, NAD-MEβ1 lost its catalytic activity; its contribution to the enzymatic activity results from a stabilizing effect on the associated α-subunit and the acquisition of regulatory properties. We conclude that in bundle sheath cell mitochondria of C4 species, the functions of NAD-ME as C4 photosynthetic decarboxylase and as a housekeeping enzyme coexist and are performed by isoforms that combine the same α-subunit with differentially adapted β-subunits.

Mutual information analysis of the dynamic correlation between side chains in proteins

Protein dynamics play an essential role in function regulation. In recent years, many experimental and theoretical studies have shown that changes in protein fluctuations in the backbone and side chains fulfill a pivotal role associated with amino acid mutations, chemical modifications, and ligand binding. The dynamic correlations between protein side chains have not been sufficiently studied, and no reliable analysis method has been available so far. Therefore, we developed a method to evaluate the dynamic correlation between protein side chains using mutual information and molecular dynamics simulations. To eliminate the structural superposition errors dealing with conventional analysis methods, and to accurately extract the intrinsic fluctuation properties of the side chains, we employed distance principal component analysis (distPCA). The motion of the side chain was then projected onto the eigenvector space obtained by distPCA, and the mutual information between the projected motions was calculated. The proposed method was then applied to a small protein "eglin c" and the mutants. The results show that even a single mutation significantly changed the dynamic correlations and also suggest that the dynamic change is deeply related to the stability. Those results indicate that our developed method could be useful for analyzing the molecular mechanism of allosteric communication in proteins.

F/G Region Rigidity is Inversely Correlated to Substrate Promiscuity of Human CYP Isoforms Involved in Metabolism

Of 57 human cytochrome P450 (CYP) enzymes, 12 metabolize 90% of xenobiotics. To our knowledge, no study has addressed the relation between enzyme dynamics and substrate promiscuity for more than three CYPs. Here, we show by constraint dilution simulations with the Constraint Network Analysis for the 12 isoforms that structural rigidity of the F/G region is significantly inversely correlated to the enzymes' substrate promiscuity. This highlights the functional importance of structural dynamics of the substrate tunnel.

Aqueous ionic liquids redistribute local enzyme stability via long-range perturbation pathways

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for biocatalysis due to their unique properties. On the other hand, the incubation of enzymes in IL or aIL often reduces enzyme activity. Recent studies proposed various aIL-induced effects to explain the reduction, classified as direct effects, e.g., local dehydration or competitive inhibition, and indirect effects, e.g., structural perturbations or disturbed catalytic site integrity. However, the molecular origin of indirect effects has largely remained elusive. Here we show by multi-μs long molecular dynamics simulations, free energy computations, and rigidity analyses that aIL favorably interact with specific residues of Bacillus subtilis Lipase A (BsLipA) and modify the local structural stability of this model enzyme by inducing long-range perturbations of noncovalent interactions. The perturbations percolate over neighboring residues and eventually affect the catalytic site and the buried protein core. Validation against a complete experimental site saturation mutagenesis library of BsLipA (3620 variants) reveals that the residues of the perturbation pathways are distinguished sequence positions where substitutions highly likely yield significantly improved residual activity. Our results demonstrate that identifying these perturbation pathways and specific IL ion-residue interactions there effectively predicts focused variant libraries with improved aIL tolerance.

Allosteric Type and Pathways Are Governed by the Forces of Protein-Ligand Binding

Allostery is central to many cellular processes, by up- or down-regulating target function. However, what determines the allosteric type remains elusive and currently it is impossible to predict whether the allosteric compounds would activate or inhibit target function before experimental studies. We demonstrated that the allosteric type and allosteric pathways are governed by the forces imposed by ligand binding to target protein using the anisotropic network model and developed an allosteric type prediction method (AlloType). AlloType correctly predicted 13 of the 16 allosteric systems in the data set with experimentally determined protein and complex structures as well as verified allosteric types, which was also used to identify allosteric pathways. When applied to glutathione peroxidase 4, a protein with no complex structure information, AlloType could still be able to predict the allosteric type of the recently reported allosteric activators, demonstrating its potential application in designing specific allosteric drugs and uncovering allosteric mechanisms.

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.

Analyzing In Silico the Relationship Between the Activation of the Edema Factor and Its Interaction With Calmodulin

Molecular dynamics (MD) simulations have been recorded on the complex between the edema factor (EF) of Bacilllus anthracis and calmodulin (CaM), starting from a structure with the orthosteric inhibitor adefovir bound in the EF catalytic site. The starting structure has been destabilized by alternately suppressing different co-factors, such as adefovir ligand or ions, revealing several long-distance correlations between the conformation of CaM, the geometry of the CaM/EF interface, the enzymatic site and the overall organization of the complex. An allosteric communication between CaM/EF interface and the EF catalytic site, highlighted by these correlations, was confirmed by several bioinformatics approaches from the literature. A network of hydrogen bonds and stacking interactions extending from the helix V of of CaM, and the residues of the switches A, B and C, and connecting to catalytic site residues, is a plausible candidate for the mediation of allosteric communication. The greatest variability in volume between the different MD conditions was also found for cavities present at the EF/CaM interface and in the EF catalytic site. The similarity between the predictions from literature and the volume variability might introduce the volume variability as new descriptor of allostery.

Structural and dynamic insights revealing how lipase binding domain MD1 of Pseudomonas aeruginosa foldase affects lipase activation

Folding and cellular localization of many proteins of Gram-negative bacteria rely on a network of chaperones and secretion systems. Among them is the lipase-specific foldase Lif, a membrane-bound steric chaperone that tightly binds (K_D = 29 nM) and mediates folding of the lipase LipA, a virulence factor of the pathogenic bacterium P. aeruginosa. Lif consists of five-domains, including a mini domain MD1 essential for LipA folding. However, the molecular mechanism of Lif-assisted LipA folding remains elusive. Here, we show in in vitro experiments using a soluble form of Lif (sLif) that isolated MD1 inhibits sLif-assisted LipA activation. Furthermore, the ability to activate LipA is lost in the variant sLif_Y99A, in which the evolutionary conserved amino acid Y99 from helix α1 of MD1 is mutated to alanine. This coincides with an approximately three-fold reduced affinity of the variant to LipA together with increased flexibility of sLif_Y99A in the complex as determined by polarization-resolved fluorescence spectroscopy. We have solved the NMR solution structures of P. aeruginosa MD1 and variant MD1_Y99A revealing a similar fold indicating that a structural modification is likely not the reason for the impaired activity of variant sLif_Y99A. Molecular dynamics simulations of the sLif:LipA complex in connection with rigidity analyses suggest a long-range network of interactions spanning from Y99 of sLif to the active site of LipA, which might be essential for LipA activation. These findings provide important details about the putative mechanism for LipA activation and point to a general mechanism of protein folding by multi-domain steric chaperones.

Basal Histamine H4 Receptor Activation: Agonist Mimicry by the Diphenylalanine Motif

Histamine H4 receptor (H4 R) orthologues are G-protein-coupled receptors (GPCRs) that exhibit species-dependent basal activity. In contrast to the basally inactive mouse H4 R (mH4 R), human H4 R (hH4 R) shows a high degree of basal activity. We have performed long-timescale molecular dynamics simulations and rigidity analyses on wild-type hH4 R, the experimentally characterized hH4 R variants S179M, F169V, F169V+S179M, F168A, and on mH4 R to investigate the molecular nature of the differential basal activity. H4 R variant-dependent differences between essential motifs of GPCR activation and structural stabilities correlate with experimentally determined basal activities and provide a molecular explanation for the differences in basal activation. Strikingly, during the MD simulations, F16945.55 dips into the orthosteric binding pocket only in the case of hH4 R, thus adopting the role of an agonist and contributing to the stabilization of the active state. The results shed new light on the molecular mechanism of basal H4 R activation that are of importance for other GPCRs.

Does Ras Activate Raf and PI3K Allosterically?

The mechanism through which oncogenic Ras activates its effectors is vastly important to resolve. If allostery is at play, then targeting allosteric pathways could help in quelling activation of MAPK (Raf/MEK/ERK) and PI3K (PI3K/Akt/mTOR) cell proliferation pathways. On the face of it, allosteric activation is reasonable: Ras binding perturbs the conformational ensembles of its effectors. Here, however, we suggest that at least for Raf, PI3K, and NORE1A (RASSF5), that is unlikely. Raf's long disordered linker dampens effective allosteric activation. Instead, we suggest that the high-affinity Ras–Raf binding relieves Raf's autoinhibition, shifting Raf's ensemble from the inactive to the nanocluster-mediated dimerized active state, as Ras also does for NORE1A. PI3K is recruited and allosterically activated by RTK (e.g., EGFR) at the membrane. Ras restrains PI3K's distribution and active site orientation. It stabilizes and facilitates PIP₂ binding at the active site and increases the PI3K residence time at the membrane. Thus, RTKs allosterically activate PI3Kα; however, merging their action with Ras accomplishes full activation. Here we review their activation mechanisms in this light and draw attention to implications for their pharmacology.

The Membrane-Integrated Steric Chaperone Lif Facilitates Active Site Opening of Pseudomonas aeruginosa Lipase A

Lipases are essential and widely used biocatalysts. Hence, the production of lipases requires a detailed understanding of the molecular mechanism of its folding and secretion. Lipase A from Pseudomonas aeruginosa, PaLipA, constitutes a prominent example that has additional relevance because of its role as a virulence factor in many diseases. PaLipA requires the assistance of a membrane-integrated steric chaperone, the lipase-specific foldase Lif, to achieve its enzymatically active state. However, the molecular mechanism of how Lif activates its cognate lipase has remained elusive. Here, we show by molecular dynamics simulations at the atomistic level and potential of mean force computations that Lif catalyzes the activation process of PaLipA by structurally stabilizing an intermediate PaLipA conformation, particularly a β-sheet in the region of residues 17-30, such that the opening of PaLipA's lid domain is facilitated. This opening allows substrate access to PaLipA's catalytic site. A surprising and so far not fully understood aspect of our study is that the open state of PaLipA is unstable compared to the closed one according to our computational and in vitro biochemical results. We thus speculate that further interactions of PaLipA with the Xcp secretion machinery and/or components of the extracellular matrix contribute to the remaining activity of secreted PaLipA. © 2019 Wiley Periodicals, Inc.

Biallelic mutation of human SLC6A6 encoding the taurine transporter TAUT is linked to early retinal degeneration

We previously reported that inactivation of the transmembrane taurine transporter (TauT or solute carrier 6a6) causes early retinal degeneration in mice. Compatible with taurine's indispensability for cell volume homeostasis, protein stabilization, cytoprotection, antioxidation, and immuno- and neuromodulation, mice develop multisystemic dysfunctions (hearing loss; liver fibrosis; and behavioral, heart, and skeletal muscle abnormalities) later on. Here, by genetic, cell biologic, in vivo ¹H-magnetic resonance spectroscopy and molecular dynamics simulation studies, we conducted in-depth characterization of a novel disorder: human TAUT deficiency. Loss of TAUT function due to a homozygous missense mutation caused panretinal degeneration in 2 brothers. TAUT_p.A78E still localized in the plasma membrane but is predicted to impact structural stabilization. ³H-taurine uptake by peripheral blood mononuclear cells was reduced by 95%, and taurine levels were severely reduced in plasma, skeletal muscle, and brain. Extraocular dysfunctions were not yet detected, but significantly increased urinary excretion of 8-oxo-7,8-dihydroguanosine indicated generally enhanced (yet clinically unapparent) oxidative stress and RNA oxidation, warranting continuous broad surveillance.-Preising, M. N., Görg, B., Friedburg, C., Qvartskhava, N., Budde, B. S., Bonus, M., Toliat, M. R., Pfleger, C., Altmüller, J., Herebian, D., Beyer, M., Zöllner, H. J., Wittsack, H.-J., Schaper, J., Klee, D., Zechner, U., Nürnberg, P., Schipper, J., Schnitzler, A., Gohlke, H., Lorenz, B., Häussinger, D., Bolz, H. J. Biallelic mutation of human SLC6A6 encoding the taurine transporter TAUT is linked to early retinal degeneration.

Cosolvent-Enhanced Sampling and Unbiased Identification of Cryptic Pockets Suitable for Structure-Based Drug Design

Modulating protein activity with small-molecules binding to cryptic pockets offers great opportunities to overcome hurdles in drug design. Cryptic sites are atypical binding sites in proteins that are closed in the absence of a stabilizing ligand and are thus inherently difficult to identify. Many studies have proposed methods to predict cryptic sites. However, a general approach to prospectively sample open conformations of these sites and to identify cryptic pockets in an unbiased manner suitable for structure-based drug design remains elusive. Here, we describe an all-atom, explicit cosolvent, molecular dynamics (MD) simulations-based workflow to sample the open states of cryptic sites and identify opened pockets, in a manner that does not require a priori knowledge about these sites. Furthermore, the workflow relies on a target-independent parametrization that only distinguishes between binding pockets for peptides or small molecules. We validated our approach on a diverse test set of seven proteins with crystallographically determined cryptic sites. The known cryptic sites were found among the three highest-ranked predicted cryptic sites, and an open site conformation was sampled and selected for most of the systems. Crystallographic ligand poses were well reproduced by docking into these identified open conformations for five of the systems. When the fully open state could not be reproduced, we were still able to predict the location of the cryptic site, or identify other cryptic sites that could be retrospectively validated with knowledge of the protein target. These characteristics render our approach valuable for investigating novel protein targets without any prior information.

Advances in the Computational Identification of Allosteric Sites and Pathways in Proteins

With the increasing difficulty to develop new drugs and the emergence of resistance to traditional orthosteric-site inhibitors, the search for alternatives is finally approaching the focus on allosteric sites. Allosteric sites offer opportunities to regulate many pharmacologically targeted pathways by inhibition or activation. In addition, allosteric sites tend to be less conserved than the functional site, which may facilitate the design of specific effectors in the protein families for which specific orthosteric inhibitors have proved difficult to design. Furthermore, recent evidence suggests that all proteins might be susceptible of allosteric regulation, increasing the space of druggable targets. Computational identification of allosteric sites has therefore become an active field of research. The problem can be approached from two sides: (1) the identification of allosteric-communication pathways between the functional site and potential allosteric sites and (2) the functional-site-independent identification of allosteric sites. While the first approach tends to be more laborious and thus restricted to a single protein, the second tends to be more amenable to larger-scale analysis, thus providing tools for the two drug discovery scenarios: the analysis of known targets and the screening for new potential targets. Here, I show some basic concepts and methods useful to the identification of allosteric sites and pathways, in line with these two approaches. I describe them in some detail to build a clear framework, at the risk of losing the interest of experts. Examples of recent studies involving these methods are also illustrated, focusing on the techniques rather than on their findings on allosterism.

Recognition motif and mechanism of ripening inhibitory peptides in plant hormone receptor ETR1

Synthetic peptides derived from ethylene-insensitive protein 2 (EIN2), a central regulator of ethylene signalling, were recently shown to delay fruit ripening by interrupting protein-protein interactions in the ethylene signalling pathway. Here, we show that the inhibitory peptide NOP-1 binds to the GAF domain of ETR1 - the prototype of the plant ethylene receptor family. Site-directed mutagenesis and computational studies reveal the peptide interaction site and a plausible molecular mechanism for the ripening inhibition.