Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

Structural Differences at Quadruplex-Duplex Interfaces Enable Ligand-Induced Topological Transitions

Quadruplex-duplex (QD) junctions, which represent unique structural motifs of both biological and technological significance, have been shown to constitute high-affinity binding sites for various ligands. A QD hybrid construct based on a human telomeric sequence, which harbors a duplex stem-loop in place of a short lateral loop, is structurally characterized by NMR. It folds into two major species with a (3+1) hybrid and a chair-type (2+2) antiparallel quadruplex domain coexisting in a K+ buffer solution. The antiparallel species is stabilized by an unusual capping structure involving a thymine and protonated adenine base AH+ of the lateral loop facing the hairpin duplex to form a T·AH+ ·G·C quartet with the interfacial G·C base pair at neutral pH. Addition and binding of Phen-DC3 to the QD hybrid mixture by its partial intercalation at corresponding QD junctions leads to a topological transition with exclusive formation of the (3+1) hybrid fold. In agreement with the available experimental data, such an unprecedented discrimination of QD junctions by a ligand can be rationalized following an induced fit mechanism.

A Kinetic Transition Network Model Reveals the Diversity of Protein Dimer Formation Mechanisms

Protein homodimers have been classified as three-state or two-state dimers depending on whether a folded monomer forms before association, but the details of the folding-binding mechanisms are poorly understood. Kinetic transition networks of conformational states have provided insight into the folding mechanisms of monomeric proteins, but extending such a network to two protein chains is challenging as all the relative positions and orientations of the chains need to be included, greatly increasing the number of degrees of freedom. Here, we present a simplification of the problem by grouping all states of the two chains into two layers: a dissociated and an associated layer. We combined our two-layer approach with the Wako-Saito-Muñoz-Eaton method and used Transition Path Theory to investigate the dimer formation kinetics of eight homodimers. The analysis reveals a remarkable diversity of dimer formation mechanisms. Induced folding, conformational selection, and rigid docking are often simultaneously at work, and their contribution depends on the protein concentration. Pre-folded structural elements are always present at the moment of association, and asymmetric binding mechanisms are common. Our two-layer network approach can be combined with various methods that generate discrete states, yielding new insights into the kinetics and pathways of flexible binding processes.

Drastic alterations in the loop structure around colchicine upon complex formation with an engineered lipocalin indicate a conformational selection mechanism

Using Anticalin technology, a lipocalin protein dubbed Colchicalin, with the ability to bind the toxic plant alkaloid colchicine with picomolar affinity, has previously been engineered, thus offering a potential antidote in vivo and also allowing its sensitive detection in biological samples. To further analyze the mode of ligand recognition, the crystal structure of Colchicalin is now reported in its unliganded form and is compared with the colchicine complex. A superposition of the protein structures revealed major rearrangements in the four structurally variable loops of the engineered lipocalin. Notably, the binding pocket in the unbound protein is largely occupied by the inward-bent loop #3, in particular Ile97, as well as by the phenylalanine side chain at position 71 in loop #2. Upon binding of colchicine, a dramatic shift of loop #3 by up to 11.1 Å occurs, in combination with a side-chain flip of Phe71, thus liberating the necessary space within the ligand pocket. Interestingly, the proline residue at the neighboring position 72, which arose during the combinatorial engineering of Colchicalin, remained in a cis configuration in both structures. These findings provide a striking example of a conformational adaptation mechanism, which is a long-known phenomenon for antibodies in immunochemistry, during the recognition of a small ligand by an engineered lipocalin, thus illustrating the general similarity between the mode of antigen/ligand binding by immunoglobulins and lipocalins.

Transient conformations in the unliganded FK506 binding domain of FKBP51 correspond to two distinct inhibitor-bound states

Members of the FK506-binding protein (FKBP) family regulate a range of important physiological processes. Unfortunately, current therapeutics such as FK506 and rapamycin exhibit only modest selectivity among these functionally distinct proteins. Recent progress in developing selective inhibitors has been reported for FKBP51 and FKBP52, which act as mutual antagonists in the regulation of steroid hormone signaling. Two structurally similar inhibitors yield distinct protein conformations at the binding site. Localized conformational transition in the binding site of the unliganded FK1 domain of FKBP51 is suppressed by a K58T mutation that also suppresses the binding of these inhibitors. Here, it is shown that the changes in amide hydrogen exchange kinetics arising from this K58T substitution are largely localized to this structural region. Accurate determination of the hydroxide-catalyzed exchange rate constants in both the wildtype and K58T variant proteins impose strong constraints upon the pattern of amide exchange reactivities within either a single or a pair of transient conformations that could give rise to the differences between these two sets of measured rate constants. Poisson-Boltzmann continuum dielectric calculations provide moderately accurate predictions of the structure-dependent hydrogen exchange reactivity for solvent-exposed protein backbone amides. Applying such calculations to the local protein conformations observed in the two inhibitor-bound FKBP51 domains demonstrated that the experimentally determined exchange rate constants for the wildtype domain are robustly predicted by a population-weighted sum of the experimental hydrogen exchange reactivity of the K58T variant and the predicted exchange reactivities in model conformations derived from the two inhibitor-bound protein structures.

Agonist efficiency links binding and gating in a nicotinic receptor

Receptors signal by switching between resting (C) and active (O) shapes ('gating') under the influence of agonists. The receptor's maximum response depends on the difference in agonist binding energy, O minus C. In nicotinic receptors, efficiency (η) represents the fraction of agonist binding energy applied to a local rearrangement (an induced fit) that initiates gating. In this receptor, free energy changes in gating and binding can be interchanged by the conversion factor η. Efficiencies estimated from concentration-response curves (23 agonists, 53 mutations) sort into five discrete classes (%): 0.56 (17), 0.51(32), 0.45(13), 0.41(26), and 0.31(12), implying that there are 5 C versus O binding site structural pairs. Within each class efficacy and affinity are corelated linearly, but multiple classes hide this relationship. η unites agonist binding with receptor gating and calibrates one link in a chain of coupled domain rearrangements that comprises the allosteric transition of the protein.

Binding to the conserved and stably folded guide RNA pseudoknot induces Cas12a conformational changes during ribonucleoprotein assembly

Ribonucleoproteins (RNPs) comprise one or more RNA and protein molecules that interact to form a stable complex, which commonly involves conformational changes in the more flexible RNA components. Here, we propose that Cas12a RNP assembly with its cognate CRISPR RNA (crRNA) guide instead proceeds primarily through Cas12a conformational changes during binding to more stable, prefolded crRNA 5' pseudoknot handles. Phylogenetic reconstructions and sequence and structure alignments revealed that the Cas12a proteins are divergent in sequence and structure while the crRNA 5' repeat region, which folds into a pseudoknot and anchors binding to Cas12a, is highly conserved. Molecular dynamics simulations of three Cas12a proteins and their cognate guides revealed substantial flexibility for unbound apo-Cas12a. In contrast, crRNA 5' pseudoknots were predicted to be stable and independently folded. Limited trypsin hydrolysis, differential scanning fluorimetry, thermal denaturation, and CD analyses supported conformational changes of Cas12a during RNP assembly and an independently folded crRNA 5' pseudoknot. This RNP assembly mechanism may be rationalized by evolutionary pressure to conserve CRISPR loci repeat sequence, and therefore guide RNA structure, to maintain function across all phases of the CRISPR defense mechanism.

Inhibition of the Human Neuronal Sodium Channel Nav1.9 by Arachidonyl-2-Chloroethylamide, An Analogue of Anandamide in a hNav1.9/rNav1.4 Chimera, An Experimental and in Silico Study

Cannabinoids regulate analgesia, which has aroused much interest in identifying new pharmacological therapies in the management of refractory pain. Voltage-gated Na+ channels (Na_vs) play an important role in inflammatory and neuropathic pain. In particular, Na_v1.9 is involved in nociception and the understanding of its pharmacology has lagged behind because it is difficult to express in heterologous systems. Here, we utilized the chimeric channel hNa_v1.9_C4, that comprises the extracellular and transmembrane domains of hNa_v1.9, co-expressed with the ß1 subunit on CHO-K1 cells to characterize the electrophysiological effects of ACEA, a synthetic surrogate of the endogenous cannabinoid anandamide. ACEA induced a tonic block, decelerated the fast inactivation, markedly shifted steady-state inactivation in the hyperpolarized direction, decreasing the window current and showed use-dependent block, with a high affinity for the inactivated state (k_i = 0.84 µM). Thus, we argue that ACEA possess a local anaesthetic-like profile. To provide a mechanistic understanding of its mode of action at the molecular level, we combined induced fit docking with Monte Carlo simulations and electrostatic complementarity. In agreement with the experimental evidence, our computer simulations revealed that ACEA binds Tyr1599 of the local anaesthetics binding site of the hNa_v1.9, contacting residues that bind cannabinol (CBD) in the Na_vMs channel. ACEA adopted a conformation remarkably similar to the crystallographic conformation of anandamide on a non-homologous protein, obstructing the Na+ permeation pathway below the selectivity filter to occupy a highly conserved binding pocket at the intracellular side. These results describe a mechanism of action, possibly involved in cannabinoid analgesia.

Competition in drug binding and … the race to equilibrium

Binding kinetics has become a popular topic in pharmacology due to its potential contribution to the selectivity and duration of drug action. Yet, the overall kinetic aspects of complex binding mechanisms are still merely described in terms of elaborate algebraic equations. Interestingly, it has been recommended some 10 years ago to examine such mechanisms in terms of binding fluxes instead of the conventional rate constants. Alike the velocity of product formation in enzymology, those fluxes refer to the velocity by which one target species converts into another one. Novel binding flux-based approaches are utilized to get a better visual insight into the "competition" between two drugs/ligands for a single target as well as between induced fit- and conformational selection pathways for a single ligand within a thermodynamic cycle. The present data were obtained by differential equation-based simulations. Early on, the ligand-binding steps "race" to equilibrium (i.e., when their forward and reverse fluxes are equal) at their individual pace. The overall/global equilibrium is only reached later on. For the competition association assays, this parting might produce a transient "overshoot" of one of the bound target species. A similar overshoot may also show up within a thermodynamic cycle and, at first glance, suggest that the induced fit pathway dominates. Yet, present findings show that under certain circumstances, it could rather be the other way round. Novel binding flux-based approaches offer visually attractive insights into crucial aspects of "complex" binding mechanisms under non-equilibrium conditions.

In silico identification of potential Hsp90 inhibitors via ensemble docking, DFT and molecular dynamics simulations

The molecular chaperone heat shock protein 90 (Hsp90) has emerged as one of the most exciting targets for anticancer drug development and Hsp90 inhibitors are potentially useful chemotherapeutic agents in cancer. Within the current study, Hsp90 inhibitors that entered different phases of clinical trials were subjected to Zinc15 structure query to find similar compounds (≥ 78%). Obtained small molecules (1-29) with defined similarity cut-off were docked into ensemble of Hsp90-α NTDs. Docked complexes were ranked on the basis of binding modes and Gibbs free energies as Hsp90 binders (cut-off point; ΔG_b ≤ -12 kcal/mol). Top-ranked compounds were subjected to energy decomposition analysis per residue of binding pocket via density functional theory (DFT) calculations in B3LYP level of theory. Subsequent MD simulations of the top-ranked complexes were performed for 100 ns to explore the stable binding modes during a reasonable period in explicit water. Results of molecular docking and intermolecular binding analysis indicated that H-bond, hydrophobic and salt bridge interactions were determinant forces in complex formation. Compounds 19 and 20 were well accommodated in binding pocket of Hsp90 via relatively varied conformations. It was revealed that Asn51 and Phe138 were key residues that interacted stably to 19 and 20. Although primary mechanism of action for proposed molecules are unknown and yet to be explored, results of the present study revealed key structural features for future structure-guided optimization toward potent inhibitors of Hsp90-α NTD. HighlightsHsp90 inhibitors that entered different phases of clinical trials were subjected to Zinc15 based structure query to afford potential enzyme inhibitors 19 and 20.Quantum chemical calculations confirmed docking results and verified pivotal role of a conserved residues (Asn51, Leu103, Phe138 and Tyr139) in making effective hydrogen bonds.MD simulations of top-ranked docked derivatives revealed the achievement of stable binding modes with less conformational variation of 20 than 19 in the active site of Hsp90-α NTD.H-bond, hydrophobic contacts and salt bridge interactions were determinant forces in binding interactions of in silico hits.Resorcinol and isoxazole were important structural motifs of in silico hits in binding to the active site of Hsp90-α NTD.Communicated by Ramaswamy H. Sarma.

Antigen binding by conformational selection in near-germline antibodies

Conformational flexibility in antibody-combining sites has been hypothesized to facilitate polyspecificity toward multiple unique epitopes and enable the limited germline repertoire to match an overwhelming diversity of potential antigens; however, elucidating the mechanisms of antigen recognition by flexible antibodies has been understandably challenging. Here, multiple liganded and unliganded crystal structures of the near-germline anticarbohydrate antibodies S25–2 and S25–39 are reported, which reveal an unprecedented diversity of complementarity-determining region H3 conformations in apparent equilibrium. These structures demonstrate that at least some germline or near-germline antibodies are flexible entities sensitive to their chemical environments, with conformational selection available as an evolved mechanism that preserves the inherited ability to recognize common pathogens while remaining adaptable to new threats.

Crystal structures of PigF, an O-methyltransferase involved in the prodigiosin synthetic pathway, reveal an induced-fit substrate-recognition mechanism

Prodigiosin, a red linear tripyrrole pigment, is a typical secondary metabolite with numerous biological functions, such as anticancer, antibacterial and immunosuppressant activities, and is synthesized through a bifurcated biosynthesis pathway from 4-methoxy-2,2'-bipyrrole-5-carbaldehyde (MBC) and 2-methyl-3-n-amylpyrrole (MAP). The last step in the biosynthetic pathway of MBC is catalysed by PigF, which transfers a methyl group to 4-hydroxy-2,20-bipyrrole-5-carbaldehyde (HBC) to form the final product MBC. However, the catalytic mechanism of PigF is still elusive. In this study, crystal structures of apo PigF and S-adenosylhomocysteine (SAH)-bound PigF were determined. PigF forms a homodimer and each monomer consists of two domains: a C-terminal catalytic domain and an N-terminal dimerization domain. Apo PigF adopts an open conformation, while the structure of the complex with the product SAH adopts a closed conformation. The binding of SAH induces dramatic conformational changes of PigF, suggesting an induced-fit substrate-binding mechanism. Further structural comparison suggests that this induced-fit substrate-recognition mechanism may generally exist in O-methyltransferases. Docking and mutation studies identified three key residues (His98, His247 and Asp248) that are crucial for enzyme activity. The essential function of His247 and Asp248 and structure analysis suggests that both residues are involved in activation of the HBC substrate of PigF. The invariance of Asp248 in PigF further confirmed its essential role. The invariance and essential role of His98 in PigF suggests that it is involved in correctly positioning the substrate. This study provides new insight into the catalytic mechanism of PigF, reveals an induced-fit substrate-recognition model for PigF and broadens the understanding of O-methyltransferases.

Computational Molecular Characterization of the Interaction of Acetylcholine and the NMDA Receptor to Explain the Direct Glycine-Competitive Potentiation of NMDA-Mediated Neuronal Currents

The activation of N-methyl-d-aspartate receptor (NMDAR) is triggered by the closure of bilobed (D1 and D2) clamshell-like clefts upon binding glycine (Gly) and glutamate. There is evidence that cholinergic compounds modulate NMDAR-mediated currents via direct receptor-ligand interactions; however, molecular bases are unknown. Here, we first propose a mechanistic structure-based explanation for the observed ACh-induced submaximal potentiation of NMDA-elicited currents in striatal neurons by predicting competitive inhibition with Gly. Then, the model was validated, in principle, by confirming that the coapplication of Gly and ACh significantly reduces these neuronal currents. Finally, we delineate the interplay of ACh with the NMDAR by a combination of computational strategies. Crystallographic ACh-bound complexes were studied, revealing a similar ACh binding environment on the GluN1 subunit of the NMDAR. We illustrate how ACh can occupy X-ray monomeric open, dimeric "semiopen" cleft conformations obtained by molecular dynamics and a full-active cryo-EM NMDAR structure, explaining the suboptimal NMDAR electrophysiological activity under the "Venus Flytrap model". At an evolutionary biology level, the binding mode of ACh coincides with that of the homologous ornithine-bound periplasmic LAO binding protein complex. Our computed results indicate an analogous mechanism of action, inasmuch as ACh may stabilize the GluN1 subunit "semiclosed" conformations by inducing direct and indirect D1-to-D2 interdomain bonds. Additionally, an alternative binding site was detected, shared by the known NMDAR allosteric modulators. Experimental and computed results strongly suggest that ACh acts as a Gly-competitive, submaximal potentiating agent of the NMDAR, possibly constituting a novel chemotype for multitarget-directed drug development, e.g., to treat Alzheimer's, and it may lead to a new understanding of glutamatergic neurotransmission.

Investigation of sugar binding kinetics of the E. coli sugar/H+ symporter XylE using solid supported membrane-based electrophysiology

Bacterial transporters are difficult to study using conventional electrophysiology because of their low transport rates and the small size of bacterial cells. Here, we applied solid-supported membrane–based electrophysiology to derive kinetic parameters of sugar translocation by the Escherichia coli xylose permease (XylE), including functionally relevant mutants. Many aspects of the fucose permease (FucP) and lactose permease (LacY) have also been investigated, which allow for more comprehensive conclusions regarding the mechanism of sugar translocation by transporters of the major facilitator superfamily. In all three of these symporters, we observed sugar binding and transport in real time to determine K_M, V_max, K_D, and k_obs values for different sugar substrates. K_D and k_obs values were attainable because of a conserved sugar-induced electrogenic conformational transition within these transporters. We also analyzed interactions between the residues in the available X-ray sugar/H⁺ symporter structures obtained with different bound sugars. We found that different sugars induce different conformational states, possibly correlating with different charge displacements in the electrophysiological assay upon sugar binding. Finally, we found that mutations in XylE altered the kinetics of glucose binding and transport, as Q175 and L297 are necessary for uncoupling H⁺ and d-glucose translocation. Based on the rates for the electrogenic conformational transition upon sugar binding (>300 s⁻¹) and for sugar translocation (2 s⁻¹ − 30 s⁻¹ for different substrates), we propose a multiple-step mechanism and postulate an energy profile for sugar translocation. We also suggest a mechanism by which d-glucose can act as an inhibitor for XylE.

Induced fit versus conformational selection: From rate constants to fluxes… and back to rate constants

Induced fit- (IF) and conformational selection (CS) binding mechanisms have long been regarded to be mutually exclusive. Yet, they are now increasingly considered to produce the final ligand-target complex alongside within a thermodynamic cycle. This viewpoint benefited from the introduction of binding fluxes as a tool for analyzing the overall behavior of such cycle. This study aims to provide more vivid and applicable insights into this emerging field. In this respect, combining differential equation- based simulations and hitherto little explored alternative modes of calculation provide concordant information about the intricate workings of such cycle. In line with previous reports, we observe that the relative contribution of IF increases with the ligand concentration at equilibrium. Yet the baseline contribution may vary from one case to another and simulations as well as calculations show that this parameter is essentially regulated by the dissociation rate of both pathways. Closer attention should be paid to how the contributions of IF and CS compare at physiologically relevant drug/ligand concentrations. To this end, a simple equation discloses how changing a limited set of "microscopic" rate constants can extend the concentration range at which CS contributes most effectively. Finally, it could also be beneficial to extend the utilization of flux- based approaches to more physiologically relevant time scales and alternative binding models.

Kinetic Constraints in the Specific Interaction between Phosphorylated Ubiquitin and Proteasomal Shuttle Factors

Ubiquitin (Ub) specifically interacts with the Ub-associating domain (UBA) in a proteasomal shuttle factor, while the latter is involved in either proteasomal targeting or self-assembly coacervation. PINK1 phosphorylates Ub at S65 and makes Ub alternate between C-terminally relaxed (pUb_RL) and retracted conformations (pUb_RT). Using NMR spectroscopy, we show that pUb_RL but not pUb_RT preferentially interacts with the UBA from two proteasomal shuttle factors Ubqln2 and Rad23A. Yet discriminatorily, Ubqln2-UBA binds to pUb more tightly than Rad23A does and selectively enriches pUb_RL upon complex formation. Further, we determine the solution structure of the complex between Ubqln2-UBA and pUb_RL and uncover the thermodynamic basis for the stronger interaction. NMR kinetics analysis at different timescales further suggests an indued-fit binding mechanism for pUb-UBA interaction. Notably, at a relatively low saturation level, the dissociation rate of the UBA-pUb_RL complex is comparable with the exchange rate between pUb_RL and pUb_RT. Thus, a kinetic constraint would dictate the interaction between Ub and UBA, thus fine-tuning the functional state of the proteasomal shuttle factors.

Structural mechanism underlying ligand binding and activation of PPARγ

Ligands bind to an occluded orthosteric ligand-binding pocket within the nuclear receptor ligand-binding domain. Molecular simulations have revealed theoretical ligand entry/exit pathways to the orthosteric pocket; however, it remains unclear whether ligand binding proceeds through induced fit or conformational selection mechanisms. Here, using nuclear magnetic resonance spectroscopy, isothermal titration calorimetry, and surface plasmon resonance analysis, we provide evidence that structurally distinct agonists bind peroxisome proliferator-activated receptor γ (PPARγ) via a two-step induced fit mechanism involving an initial fast kinetic step followed by a slow conformational change. The agonist encounter complex binding pose is suggested in crystal structures where ligands bind to a surface pore suggested as a ligand entry site in molecular simulations. Our findings suggest an activation mechanism for PPARγ whereby agonist binding occurs through an initial encounter complex followed by a transition of the ligand into the final binding pose within the orthosteric pocket, inducing a transcriptionally active conformation.

In Silico Drug Repurposing by Structural Alteration after Induced Fit: Discovery of a Candidate Agent for Recovery of Nucleotide Excision Repair in Xeroderma Pigmentosum Group D Mutant (R683W)

Xeroderma pigmentosum complementation group D (XPD) is a UV-sensitive syndrome and a rare incurable genetic disease which is caused by the genetic mutation of the excision repair cross-complementation group 2 gene (ERCC2). Patients who harbor only XPD R683W mutant protein develop severe photosensitivity and progressive neurological symptoms. Cultured cells derived from patients with XPD (XPD R683W cells) demonstrate a reduced nucleotide excision repair (NER) ability. We hope to ameliorate clinical symptoms if we can identify candidate agents that would aid recovery of the cells' NER ability. To investigate such candidates, we created in silico methods of drug repurposing (in silico DR), a strategy that utilizes the recovery of ATP-binding in the XPD R683W protein after the induced fit. We chose 4E1RCat and aprepitant as the candidates for our in silico DR, and evaluated them by using the UV-induced unscheduled DNA synthesis (UDS) assay to verify the recovery of NER in XPD R683W cells. UDS values of the cells improved about 1.4-1.7 times after 4E1RCat treatment compared with solvent-only controls; aprepitant showed no positive effect. In this study, therefore, we succeeded in finding the candidate agent 4E1RCat for XPD R683W. We also demonstrated that our in silico DR method is a cost-effective approach for drug candidate discovery.

Kinetic Evidence for an Induced Fit Mechanism in the Binding of the Substrate Camphor by Cytochrome P450cam

Bacterial cytochrome P450 (P450) 101A1 (P450cam) has served as a prototype among the P450 enzymes and has high catalytic activity towards its cognate substrate, camphor. X-ray crystallography and NMR and IR spectroscopy have demonstrated the existence of multiple conformations of many P450s, including P450cam. Kinetic studies have indicated that substrate binding to several P450s is dominated by a conformational selection process, in which the substrate binds an individual conformer(s) of the unliganded enzyme. P450cam was found to differ in that binding of the substrate camphor is dominated by an induced fit mechanism, in which the enzyme binds camphor and then changes conformation, as evidenced by the equivalence of binding eigenvalues observed when varying both camphor and P450cam concentrations. The accessory protein putidaredoxin had no effect on substrate binding. Estimation of the rate of dissociation of the P450cam·camphor complex (15 s-1) and fitting of the data yield a minimal kinetic mechanism in which camphor binds (1.5 × 107 M-1 s-1) and the initial P450cam•camphor complex undergoes a reversible equilibrium (k forward 112 s-1, k reverse 28 s-1) to a final complex. This induced fit mechanism differs from those reported for several mammalian P450s and bacterial P450BM-3, indicative of the diversity of how P450s recognize multiple substrates. However, similar behavior was not observed with the alternate substrates (+)-α-pinene and 2-adamantanone, which probably utilize a conformational selection process.

Guest-Reaction Driven Cage to Conjoined Twin-Cage Mitosis-Like Host Transformation

We report here a guest-reaction-induced mitosis-like host transformation from a known Pd₄ L₂ cage 1 to a conjoined Pd₆ L₃ twin-cage 2 featuring two separate cavities. The encapsulation of 1-hydroxymethyl-2-naphthol (G1), a known ortho-quinone methide (o-QMs) precursor, within the hydrophobic cavity of cage 1 is found crucial to realize the cage to twin-cage conversion. Confined G1 molecules within the nanocavity undergo self-coupling dimerization reaction to form 2,2'-dihydroxy-1,1'-dinaphthylmethane (G2) which then triggers the cage to twin-cage mitosis. The same conversion also proceeds, in a much faster rate, via the direct templation of G2, confirming the induced-fit transformation mechanism. The structure of the (G2)₂ ⊂2 host-guest complex has been established by X-ray crystallographic study, where cis- to trans- conformational switch on one bridging ligand is revealed.

Snapshots during the catalytic cycle of a histidine acid phytase reveal an induced-fit structural mechanism

Highly engineered phytases, which sequentially hydrolyze the hexakisphosphate ester of inositol known as phytic acid, are routinely added to the feeds of monogastric animals to improve phosphate bioavailability. New phytases are sought as starting points to further optimize the rate and extent of dephosphorylation of phytate in the animal digestive tract. Multiple inositol polyphosphate phosphatases (MINPPs) are clade 2 histidine phosphatases (HP2P) able to carry out the stepwise hydrolysis of phytate. MINPPs are not restricted by a strong positional specificity making them attractive targets for development as feed enzymes. Here, we describe the characterization of a MINPP from the Gram-positive bacterium Bifidobacterium longum (BlMINPP). BlMINPP has a typical HP2P-fold but, unusually, possesses a large α-domain polypeptide insertion relative to other MINPPs. This insertion, termed the U-loop, spans the active site and contributes to substrate specificity pockets underpopulated in other HP2Ps. Mutagenesis of U-loop residues reveals its contribution to enzyme kinetics and thermostability. Moreover, four crystal structures of the protein along the catalytic cycle capture, for the first time in an HP2P, a large ligand-driven α-domain motion essential to allow substrate access to the active site. This motion recruits residues both downstream of a molecular hinge and on the U-loop to participate in specificity subsites, and mutagenesis identified a mobile lysine residue as a key determinant of positional specificity of the enzyme. Taken together, these data provide important new insights to the factors determining stability, substrate recognition, and the structural mechanism of hydrolysis in this industrially important group of enzymes.

Kinetic and thermodynamic analysis defines roles for two metal ions in DNA polymerase specificity and catalysis

Magnesium ions play a critical role in catalysis by many enzymes and contribute to the fidelity of DNA polymerases through a two-metal ion mechanism. However, specificity is a kinetic phenomenon and the roles of Mg²⁺ ions in each step in the catalysis have not been resolved. We first examined the roles of Mg²⁺ by kinetic analysis of single nucleotide incorporation catalyzed by HIV reverse transcriptase. We show that Mg.dNTP binding induces an enzyme conformational change at a rate that is independent of free Mg²⁺ concentration. Subsequently, the second Mg²⁺ binds to the closed state of the enzyme-DNA-Mg.dNTP complex (K_d = 3.7 mM) to facilitate catalysis. Weak binding of the catalytic Mg²⁺ contributes to fidelity by sampling the correctly aligned substrate without perturbing the equilibrium for nucleotide binding at physiological Mg²⁺ concentrations. An increase of the Mg²⁺ concentration from 0.25 to 10 mM increases nucleotide specificity (k_cat/K_m) 12-fold largely by increasing the rate of the chemistry relative to the rate of nucleotide release. Mg²⁺ binds very weakly (K_d ≤ 37 mM) to the open state of the enzyme. Analysis of published crystal structures showed that HIV reverse transcriptase binds only two metal ions prior to incorporation of a correct base pair. Molecular dynamics simulations support the two-metal ion mechanism and the kinetic data indicating weak binding of the catalytic Mg²⁺. Molecular dynamics simulations also revealed the importance of the divalent cation cloud surrounding exposed phosphates on the DNA. These results enlighten the roles of the two metal ions in the specificity of DNA polymerases.

Structural analysis of a novel substrate-free form of the aminoglycoside 6'-N-acetyltransferase from Enterococcus faecium

Aminoglycoside acetyltransferases (AACs) catalyze the transfer of an acetyl group between acetyl-CoA and an aminoglycoside, producing CoA and an acetylated aminoglycoside. AAC(6')-Ii enzymes target the amino group linked to the 6' C atom in an aminoglycoside. Several structures of the AAC(6')-Ii from Enterococcus faecium [Ef-AAC(6')-Ii] have been reported to date. However, the detailed mechanism of its enzymatic function remains elusive. In this study, the crystal structure of Ef-AAC(6')-Ii was determined in a novel substrate-free form. Based on structural analysis, it is proposed that Ef-AAC(6')-Ii sequentially undergoes conformational selection and induced fit for substrate binding. These results therefore provide a novel viewpoint on the mechanism of action of Ef-AAC(6')-Ii.

Enhanced sampling of protein conformational states for dynamic cross-docking within the protein-protein docking server SwarmDock

The formation of specific protein-protein interactions is often a key to a protein's function. During complex formation, each protein component will undergo a change in the conformational state, for some these changes are relatively small and reside primarily at the sidechain level; however, others may display notable backbone adjustments. One of the classic problems in the protein-docking field is to be able to a priori predict the extent of such conformational changes. In this work, we investigated three protocols to find the most suitable input structure conformations for cross-docking, including a robust sampling approach in normal mode space. Counterintuitively, knowledge of the theoretically best combination of normal modes for unbound-bound transitions does not always lead to the best results. We used a novel spatial partitioning library, Aether Engine (see Supplementary Materials), to efficiently search the conformational states of 56 receptor/ligand pairs, including a recent CAPRI target, in a systematic manner and selected diverse conformations as input to our automated docking server, SwarmDock, a server that allows moderate conformational adjustments during the docking process. In essence, here we present a dynamic cross-docking protocol, which when benchmarked against the simpler approach of just docking the unbound components shows a 10% uplift in the quality of the top docking pose.

Mechanisms of promiscuity among drug metabolizing enzymes and drug transporters

Detoxication, or 'drug-metabolizing', enzymes and drug transporters exhibit remarkable substrate promiscuity and catalytic promiscuity. In contrast to substrate-specific enzymes that participate in defined metabolic pathways, individual detoxication enzymes must cope with substrates of vast structural diversity, including previously unencountered environmental toxins. Presumably, evolution selects for a balance of 'adequate' kcat /KM values for a wide range of substrates, rather than optimizing kcat /KM for any individual substrate. However, the structural, energetic, and metabolic properties that achieve this balance, and hence optimize detoxication, are not well understood. Two features of detoxication enzymes that are frequently cited as contributions to promiscuity include the exploitation of highly reactive versatile cofactors, or cosubstrates, and a high degree of flexibility within the protein structure. This review examines these intuitive mechanisms in detail and clarifies the contributions of the classic ligand binding models 'induced fit' (IF) and 'conformational selection' (CS) to substrate promiscuity. The available literature data for drug metabolizing enzymes and transporters suggest that IF is exploited by these promiscuous detoxication enzymes, as it is with substrate-specific enzymes, but the detoxication enzymes uniquely exploit 'IFs' to retain a wide range of substrates at their active sites. In contrast, whereas CS provides no catalytic advantage to substrate-specific enzymes, promiscuous enzymes may uniquely exploit it to recruit a wide range of substrates. The combination of CS and IF, for recruitment and retention of substrates, can potentially optimize the promiscuity of drug metabolizing enzymes and drug transporters.

The Arabidopsis (ASHH2) CW domain binds monomethylated K4 of the histone H3 tail through conformational selection

Chromatin post-translational modifications are thought to be important for epigenetic effects on gene expression. Methylation of histone N-terminal tail lysine residues constitutes one of many such modifications, executed by families of histone lysine methyltransferase (HKMTase). One such protein is ASHH2 from the flowering plant Arabidopsis thaliana, equipped with the interaction domain, CW, and the HKMTase domain, SET. The CW domain of ASHH2 is a selective binder of monomethylation at lysine 4 on histone H3 (H3K4me1) and likely helps the enzyme dock correctly onto chromatin sites. The study of CW and related interaction domains has so far been emphasizing lock-key models, missing important aspects of histone-tail CW interactions. We here present an analysis of the ASHH2 CW-H3K4me1 complex using NMR and molecular dynamics, as well as mutation and affinity studies of flexible coils. β-augmentation and rearrangement of coils coincide with changes in the flexibility of the complex, in particular the η1, η3 and C-terminal coils, but also in the β1 and β2 strands and the C-terminal part of the ligand. Furthermore, we show that mutating residues with outlier dynamic behaviour affect the complex binding affinity despite these not being in direct contact with the ligand. Overall, the binding process is consistent with conformational selection. We propose that this binding mechanism presents an advantage when searching for the correct post-translational modification state among the highly modified and flexible histone tails, and also that the binding shifts the catalytic SET domain towards the nucleosome. DATABASES: Structural data are available in the PDB database under the accession code 6QXZ. Resonance assignments for CW42 in its apo- and holo-forms are available in the BMRB database under the accession code 27251.

Large-Scale Virtual Screening Against the MET Kinase Domain Identifies a New Putative Inhibitor Type

By using an ensemble-docking strategy, we undertook a large-scale virtual screening campaign in order to identify new putative hits against the MET kinase target. Following a large molecular dynamics sampling of its conformational space, a set of 45 conformers of the kinase was retained as docking targets to take into account the flexibility of the binding site moieties. Our screening funnel started from about 80,000 chemical compounds to be tested in silico for their potential affinities towards the kinase binding site. The top 100 molecules selected—thanks to the molecular docking results—were further analyzed for their interactions, and 25 of the most promising ligands were tested for their ability to inhibit MET activity in cells. F0514-4011 compound was the most efficient and impaired this scattering response to HGF (Hepatocyte Growth Factor) with an IC50 of 7.2 μM. Interestingly, careful docking analysis of this molecule with MET suggests a possible conformation halfway between classical type-I and type-II MET inhibitors, with an additional region of interaction. This compound could therefore be an innovative seed to be repositioned from its initial antiviral purpose towards the field of MET inhibitors. Altogether, these results validate our ensemble docking strategy as a cost-effective functional method for drug development.

Inhibitory Activity of Flavonoids, Chrysoeriol and Luteolin-7-O-Glucopyranoside, on Soluble Epoxide Hydrolase from Capsicum chinense

: Three flavonoids derived from the leaves of Capsicum chinense Jacq. were identified as chrysoeriol (1), luteolin-7-O-glucopyranoside (2), and isorhamnetin-7-O-glucopyranoside (3). They had IC50 values of 11.6±2.9, 14.4±1.5, and 42.7±3.5 µg/mL against soluble epoxide hydrolase (sEH), respectively. The three inhibitors (1-3) were found to non-competitively bind into the allosteric site of the enzyme with Ki values of 10.5±3.2, 11.9 ±2.8 and 38.0±4.1 µg/mL, respectively. The potential inhibitors 1 and 2 were located at the left edge ofa U-tube shape that contained the enzyme active site. Additionally, we observed changes in several factors involved in the binding of these complexes under 300 K and 1 bar. Finally, it was confirmed that each inhibitor, 1 and 2, could be complexed with sEH by the "induced fit" and "lock-and-key" models.

Coarse-grained and atomic resolution biomolecular docking with the ATTRACT approach

The ATTRACT protein-protein docking program has been employed to predict protein-protein complex structures in CAPRI rounds 38-45. For 11 out of 16 targets acceptable or better quality solutions have been submitted (~70%). It includes also several cases of peptide-protein docking and the successful prediction of the geometry of carbohydrate-protein interactions. The option of combining rigid body minimization and simultaneous optimization in collective degrees of freedom based on elastic network modes was employed and systematically evaluated. Application to a large benchmark set indicates a modest improvement in docking performance compared to rigid docking. Possible further improvements of the docking approach in particular at the scoring and the flexible refinement steps are discussed.

Structure-Based Analysis of Cryptic-Site Opening

Many proteins in their unbound structures have cryptic sites that are not appropriately sized for drug binding. We consider here 32 proteins from the recently published CryptoSite set with validated cryptic sites, and study whether the sites remain cryptic in all available X-ray structures of the proteins solved without any ligand bound near the sites. It was shown that only few of these proteins have binding pockets that never form without ligand binding. Sites that are cryptic in some structures but spontaneously form in others are also rare. In most proteins the forming of pockets is affected by mutations or ligand binding at locations far from the cryptic site. To further explore these mechanisms, we applied adiabatic biased molecular dynamics simulations to guide the proteins from their ligand-free structures to ligand-bound conformations, and studied the distribution of druggability scores of the pockets located at the cryptic sites.

D Amino Acids Highlight the Catalytic Power of the Ribosome

The possible mechanism(s) by which ribosomes make peptide bonds during protein synthesis have been explored for decades. Yet, there is no agreement on how the catalytic site, the peptidyl transferase center (PTC), promotes this reaction. Here, we discuss the results of recent investigations of translation with D amino acids that provide fresh insights into that longstanding question.

Binding-induced folding under unfolding conditions: Switching between induced fit and conformational selection mechanisms

The chemistry of protein-ligand binding is the basis of virtually every biological process. Ligand binding can be essential for a protein to function in the cell by stabilizing or altering the conformation of a protein, particularly for partially or completely unstructured proteins. However, the mechanisms by which ligand binding impacts disordered proteins or influences the role of disorder in protein folding is not clear. To gain insight into this question, the mechanism of folding induced by the binding of a Pro-rich peptide ligand to the SH3 domain of phosphatidylinositol 3-kinase unfolded in the presence of urea has been studied using kinetic methods. Under strongly denaturing conditions, folding was found to follow a conformational selection (CS) mechanism. However, under mildly denaturing conditions, a ligand concentration-dependent switch in the mechanism was observed. The folding mechanism switched from being predominantly a CS mechanism at low ligand concentrations to being predominantly an induced fit (IF) mechanism at high ligand concentrations. The switch in the mechanism manifests itself as an increase in the reaction flux along the IF pathway at high ligand concentrations. The results indicate that, in the case of intrinsically disordered proteins too, the folding mechanism is determined by the concentration of the ligand that induces structure formation.

Mismatch Recognition by Saccharomyces cerevisiae Msh2-Msh6: Role of Structure and Dynamics

The mismatch repair (MMR) pathway maintains genome integrity by correcting errors such as mismatched base pairs formed during DNA replication. In MMR, Msh2–Msh6, a heterodimeric protein, targets single base mismatches and small insertion/deletion loops for repair. By incorporating the fluorescent nucleoside base analog 6-methylisoxanthopterin (6-MI) at or adjacent to a mismatch site to probe the structural and dynamic elements of the mismatch, we address how Msh2–Msh6 recognizes these mismatches for repair within the context of matched DNA. Fluorescence quantum yield and rotational correlation time measurements indicate that local base dynamics linearly correlate with Saccharomyces cerevisiae Msh2–Msh6 binding affinity where the protein exhibits a higher affinity (K_D ≤ 25 nM) for mismatches that have a significant amount of dynamic motion. Energy transfer measurements measuring global DNA bending find that mismatches that are both well and poorly recognized by Msh2–Msh6 experience the same amount of protein-induced bending. Finally, base-specific dynamics coupled with protein-induced blue shifts in peak emission strongly support the crystallographic model of directional binding, in which Phe 432 of Msh6 intercalates 3′ of the mismatch. These results imply an important role for local base dynamics in the initial recognition step of MMR.

Profiling the interaction of 1-phenylbenzimidazoles to cyclooxygenases

In the design of 1-phenylbenzimidazoles as model cyclooxygenase (COX) inhibitors, docking to a series of crystallographic COX structures was performed to evaluate their potential for high-affinity binding and to reproduce the interaction profile of well-known COX inhibitors. The effect of ligand-specific induced fit on the calculations was also studied. To quantitatively compare the pattern of interactions of model compounds to the profile of several cocrystallized COX inhibitors, a geometric parameter, denominated ligand-receptor contact distance (LRCD), was developed. The interaction profile of several model complexes showed similarity to the profile of COX complexes with inhibitors such as iodosuprofen, iodoindomethacin, diclofenac, and flurbiprofen. Shaping of high-affinity binding sites upon ligand-specific induced fit mostly determined both the affinity and the binding mode of the ligands in the docking calculations. The results suggest potential of 1-phenylbenzimidazole derivatives as COX inhibitors on the basis of their predicted affinity and interaction profile to COX enzymes. The analyses also provided insights into the role of induced fit in COX enzymes. While inhibitors produce different local structural changes at the COX ligand binding site, induced fit allows inhibitors in diverse chemical classes to share characteristic interaction patterns that ensure key contacts to be achieved. Different interaction patterns may also be associated with different inhibitory mechanisms.

Human cytochrome P450 enzymes bind drugs and other substrates mainly through conformational-selection modes

Cytochrome P450 (P450) enzymes are major catalysts involved in the oxidations of most drugs, steroids, carcinogens, fat-soluble vitamins, and natural products. The binding of substrates to some of the 57 human P450s and other mammalian P450s is more complex than a two-state system and has been proposed to involve mechanisms such as multiple ligand occupancy, induced-fit, and conformational-selection. Here, we used kinetic analysis of binding with multiple concentrations of substrates and computational modeling of these data to discern possible binding modes of several human P450s. We observed that P450 2D6 binds its ligand rolapitant in a mechanism involving conformational-selection. P450 4A11 bound the substrate lauric acid via conformational-selection, as did P450 2C8 with palmitic acid. Binding of the steroid progesterone to P450 21A2 was also best described by a conformational-selection model. Hexyl isonicotinate binding to P450 2E1 could be described by either a conformational-selection or an induced-fit model. Simulation of the binding of the ligands midazolam, bromocriptine, testosterone, and ketoconazole to P450 3A4 was consistent with an induced-fit or a conformational-selection model, but the concentration dependence of binding rates for varying both P450 3A4 and midazolam concentrations revealed discordance in the parameters, indicative of conformational-selection. Binding of the P450s 2C8, 2D6, 3A4, 4A11, and 21A2 was best described by conformational-selection, and P450 2E1 appeared to fit either mode. These findings highlight the complexity of human P450-substrate interactions and that conformational-selection is a dominant feature of many of these interactions.

Substrate-Induced Formation of Ribosomal Decoding Center for Accurate and Rapid Genetic Code Translation

Accurate translation of genetic information is crucial for synthesis of functional proteins in all organisms. We use recent experimental data to discuss how induced fit affects accuracy of initial codon selection on the ribosome by aminoacyl transfer RNA in ternary complex ( T₃) with elongation factor Tu (EF-Tu) and guanosine-5'-triphosphate (GTP). We define actual accuracy ([Formula: see text]) of a particular protein synthesis system as its current accuracy and the effective selectivity ([Formula: see text]) as [Formula: see text] in the limit of zero ribosomal binding affinity for T₃. Intrinsic selectivity ([Formula: see text]), defined as the upper thermodynamic limit of [Formula: see text], is determined by the free energy difference between near-cognate and cognate T₃ in the pre-GTP hydrolysis state on the ribosome. [Formula: see text] is much larger than [Formula: see text], suggesting the possibility of a considerable increase in [Formula: see text] and [Formula: see text] at negligible kinetic cost. Induced fit increases [Formula: see text] and [Formula: see text] without affecting [Formula: see text], and aminoglycoside antibiotics reduce [Formula: see text] and [Formula: see text] at unaltered [Formula: see text].

The Molecular Mechanism of Transport by the Mitochondrial ADP/ATP Carrier

Mitochondrial ADP/ATP carriers transport ADP into the mitochondrial matrix for ATP synthesis, and ATP out to fuel the cell, by cycling between cytoplasmic-open and matrix-open states. The structure of the cytoplasmic-open state is known, but it has proved difficult to understand the transport mechanism in the absence of a structure in the matrix-open state. Here, we describe the structure of the matrix-open state locked by bongkrekic acid bound in the ADP/ATP-binding site at the bottom of the central cavity. The cytoplasmic side of the carrier is closed by conserved hydrophobic residues, and a salt bridge network, braced by tyrosines. Glycine and small amino acid residues allow close-packing of helices on the matrix side. Uniquely, the carrier switches between states by rotation of its three domains about a fulcrum provided by the substrate-binding site. Because these features are highly conserved, this mechanism is likely to apply to the whole mitochondrial carrier family.

Video Abstract

•

Structure of the matrix-open state of the mitochondrial ADP/ATP carrier solved

•

The inhibitor bongkrekic acid locks the state by occupying the substrate-binding site

•

Conformational changes during transport are highly dynamic, using six mobile elements

•

Roles of all conserved sequence features in mitochondrial carriers are now explained

Binding of a Telomestatin Derivative Changes the Mechanical Anisotropy of a Human Telomeric G-Quadruplex

Mechanical anisotropy is an essential property for biomolecules to assume structural and functional roles in mechanobiology. However, there is insufficient information on the mechanical anisotropy of ligand-biomolecule complexes. Herein, we investigated the mechanical property of individual human telomeric G-quadruplexes bound to telomestatin, using optical tweezers. Stacking of the ligand to the G-tetrad planes changes the conformation of the G-quadruplex, which resembles a balloon squeezed in certain directions. Such a squeezed balloon effect strengthens the G-tetrad planes, but dislocates and weakens the loops in the G-quadruplex upon ligand binding. These dynamic interactions indicate that the binding between the ligand and G-quadruplex follows the induced-fit model. We anticipate that the altered mechanical anisotropy of the ligand-G-quadruplex complex can add additional level of regulations on the motor enzymes that process DNA or RNA molecules.

Structural characterization of free-state and product-state Mycobacterium tuberculosis methionyl-tRNA synthetase reveals an induced-fit ligand-recognition mechanism

Mycobacterium tuberculosis (MTB) caused 10.4 million cases of tuberculosis and 1.7 million deaths in 2016. The incidence of multidrug-resistant and extensively drug-resistant MTB is becoming an increasing threat to public health and the development of novel anti-MTB drugs is urgently needed. Methionyl-tRNA synthetase (MetRS) is considered to be a valuable drug target. However, structural characterization of M. tuberculosis MetRS (MtMetRS) was lacking for decades, thus hampering drug design. Here, two high-resolution crystal structures of MtMetRS are reported: the free-state structure (apo form; 1.9 Å resolution) and a structure with the intermediate product methionyl-adenylate (Met-AMP) bound (2.4 Å resolution). It was found that free-state MtMetRS adopts a previously unseen conformation that has never been observed in other MetRS homologues. The pockets for methionine and AMP are not formed in free-state MtMetRS, suggesting that it is in a nonproductive conformation. Combining these findings suggests that MtMetRS employs an induced-fit mechanism in ligand binding. By comparison with the structure of human cytosolic MetRS, additional pockets specific to MtMetRS that could be used for anti-MTB drug design were located.

Magnesium Activates Microsecond Dynamics to Regulate Integrin-Collagen Recognition

Integrin receptors bind collagen via metal-mediated interactions that are modulated by magnesium (Mg²⁺) levels in the extracellular matrix. Nuclear magnetic resonance-based relaxation experiments, isothermal titration calorimetry, and adhesion assays reveal that Mg²⁺ functions as both a structural anchor and dynamic switch of the α₁β₁ integrin I domain (α₁I). Specifically, Mg²⁺ binding activates micro- to millisecond timescale motions of residues distal to the binding site, particularly those surrounding the salt bridge at helix 7 and near the metal ion-dependent adhesion site. Mutagenesis of these residues impacts α₁I functional activity, thereby suggesting that Mg-bound α₁I dynamics are important for collagen binding and consequent allosteric rearrangement of the low-affinity closed to high-affinity open conformation. We propose a multistep recognition mechanism for α₁I-Mg-collagen interactions involving both conformational selection and induced-fit processes. Our findings unravel the multifaceted role of Mg²⁺ in integrin-collagen recognition and assist in elucidating the molecular mechanisms by which metals regulate protein-protein interactions.

Confirmation for conformational selection

NMR studies settle part of a long-standing debate about the mechanism used by the Hsp70 chaperone to recognize substrates.

Structural insights into humanization of anti-tissue factor antibody 10H10

Murine antibody 10H10 raised against human tissue factor is unique in that it blocks the signaling pathway, and thus inhibits angiogenesis and tumor growth without interfering with coagulation. As a potential therapeutic, the antibody was humanized in a two-step procedure. Antigen-binding loops were grafted onto selected human frameworks and the resulting chimeric antibody was subjected to affinity maturation by using phage display libraries. The results of humanization were analyzed from the structural perspective through comparison of the structure of a humanized variant with the parental mouse antibody. This analysis revealed several hot spots in the framework region that appear to affect antigen binding, and therefore should be considered in human germline selection. In addition, some positions in the Vernier zone, e.g., residue 71 in the heavy chain, that are traditionally thought to be crucial appear to tolerate amino acid substitutions without any effect on binding. Several humanized variants were produced using both short and long forms of complementarity-determining region (CDR) H2 following the difference in the Kabat and Martin definitions. Comparison of such pairs indicated consistently higher thermostability of the variants with short CDR H2. Analysis of the binding data in relation to the structures singled out the ImMunoGeneTics information system® germline IGHV1-2*01 as dubious owing to two potentially destabilizing mutations as compared to the other alleles of the same germline and to other human germlines.

Understanding the mechanism of binding between Gab2 and the C terminal SH3 domain from Grb2

Gab2 is a large disordered protein that regulates several cellular signalling pathways and is overexpressed in different forms of cancer. Because of its disordered nature, a detailed characterization of the mechanisms of recognition between Gab2 and its physiological partners is particularly difficult. Here we provide a detailed kinetic characterization of the binding reaction between Gab2 and the C-terminal SH3 domain of the growth factor receptor-bound protein 2 (Grb2). We demonstrate that Gab2 folds upon binding following an induced fit type mechanism, whereby recognition is characterized by the formation of an intermediate, in which Gab2 is primarily disordered. In this scenario, folding of Gab2 into the bound conformation occurs only after binding. However, an alanine scanning of the proline residues of Gab2 suggests that the intermediate contains some degree of native-like structure, which might play a role for the recognition event to take place. The results, which represent a fundamental step forward in the understanding of this functional protein-protein interaction, are discussed on the light of previous structural works on these proteins.

Probing protein flexibility reveals a mechanism for selective promiscuity

Many eukaryotic regulatory proteins adopt distinct bound and unbound conformations, and use this structural flexibility to bind specifically to multiple partners. However, we lack an understanding of how an interface can select some ligands, but not others. Here, we present a molecular dynamics approach to identify and quantitatively evaluate the interactions responsible for this selective promiscuity. We apply this approach to the anticancer target PD-1 and its ligands PD-L1 and PD-L2. We discover that while unbound PD-1 exhibits a hard-to-drug hydrophilic interface, conserved specific triggers encoded in the cognate ligands activate a promiscuous binding pathway that reveals a flexible hydrophobic binding cavity. Specificity is then established by additional contacts that stabilize the PD-1 cavity into distinct bound-like modes. Collectively, our studies provide insight into the structural basis and evolution of multiple binding partners, and also suggest a biophysical approach to exploit innate binding pathways to drug seemingly undruggable targets.

DOI:http://dx.doi.org/10.7554/eLife.22889.001

Many proteins need to interact with other proteins to carry out their various tasks in cells. Such interactions are essential for almost all biological processes and are often disrupted in disease. Cells have thousands of different types of proteins and each has a unique shape that determines which other proteins it can bind to.

It was previously thought that two proteins bind to each other in a manner similar to that of a lock and a key, in which the rigid shape of one protein meshes perfectly with the rigid shape of its partner. However, many proteins are flexible and adopt different shapes depending on whether they are attached to their partner, or not. Moreover, an individual protein may also bind to several different partners, each requiring that protein to adopt several different shapes. These observations have challenged the lock and key model and suggest that flexibility in the structure of a protein plays a key role in its binding to other proteins. However, it is not clear how structural flexibility enables a protein to bind to several different partners while being selective enough to prevent the protein from binding to the wrong ones.

A protein called PD-1 is involved in immune responses in humans and is an emerging target for drugs to treat cancer. Pabon and Camacho used computer simulations to model PD-1’s structural flexibility and to find out how this enables the protein to form different shapes when it binds to different partners. The experiments show that the region of PD-1 that binds to other proteins adopts a different shape in the absence and presence of its partners. The binding partners make initial contact with PD-1 via specific features that they share in common. This causes PD-1 to change shape, uncovering a surface of PD-1 that is flexible and is able to accommodate a variety of partners. After this, the binding partners form additional contacts with PD-1 that are specific to each partner.

These findings suggest that the ability of a protein to bind to several different partners is unlocked by certain structures that are present in the binding partners. These structures are found in proteins produced by many different organisms, suggesting that this mechanism is likely to be widespread in nature. This work may open up new avenues for designing drugs to target PD-1 and other proteins that contribute to disease but have so far been impossible to target with drugs.

DOI:http://dx.doi.org/10.7554/eLife.22889.002

Are induced fit protein conformational changes caused by ligand-binding predictable? A molecular dynamics investigation

In this work, the ability of molecular dynamics simulations (MD) to prospectively predict regions of ligand binding sites that could undergo induced fit effects was investigated. Conventional MD was run on 39 apo structures (no ligand), and the resulting trajectories were compared to a set of 147 holo X-ray structures (ligand-bound). It was observed from the simulations, in the absence of the ligands, that structures exhibiting large residue conformational changes indicated higher likelihood of induced fit effects. Nevertheless, the simulation results did not perform better than using the normalized crystallographic structural factors as predictors of active-site rigid residues (87% predictive power) and mobile residues (47% predictive power). While the simulations could not produce full active sites conformations similar to holo-like states, it was found that the simulations could reproduce bound state conformations of individual residues. These results suggest potential issues in the use of unligated simulation frames directly for drug design applications such as ligand docking, and an overall caution in the use of protein flexibility in docking protocols should be emphasized. © 2017 Wiley Periodicals, Inc.

Rate Constants and Mechanisms of Protein-Ligand Binding

Whereas protein-ligand binding affinities have long-established prominence, binding rate constants and binding mechanisms have gained increasing attention in recent years. Both new computational methods and new experimental techniques have been developed to characterize the latter properties. It is now realized that binding mechanisms, like binding rate constants, can and should be quantitatively determined. In this review, we summarize studies and synthesize ideas on several topics in the hope of providing a coherent picture of and physical insight into binding kinetics. The topics include microscopic formulation of the kinetic problem and its reduction to simple rate equations; computation of binding rate constants; quantitative determination of binding mechanisms; and elucidation of physical factors that control binding rate constants and mechanisms.

Induced fit of the peptidyl-transferase center of the ribosome and conformational freedom of the esterified amino acids

The catalytic site of most enzymes can efficiently handle only one substrate. In contrast, the ribosome is capable of polymerizing at a similar rate at least 20 different kinds of amino acids from aminoacyl-tRNA carriers while using just one catalytic site, the peptidyl-transferase center (PTC). An induced-fit mechanism has been uncovered in the PTC, but a possible connection between this mechanism and the uniform handling of the substrates has not been investigated. We present an analysis of published ribosome structures supporting the hypothesis that the induced fit eliminates unreactive rotamers predominantly populated for some A-site aminoacyl esters before induction. We show that this hypothesis is fully consistent with the wealth of kinetic data obtained with these substrates. Our analysis reveals that induction constrains the amino acids into a reactive conformation in a side-chain independent manner. It allows us to highlight the rationale of the PTC structural organization, which confers to the ribosome the very unusual ability to handle large as well as small substrates.

Genome-scale characterization of RNA tertiary structures and their functional impact by RNA solvent accessibility prediction

As most RNA structures are elusive to structure determination, obtaining solvent accessible surface areas (ASAs) of nucleotides in an RNA structure is an important first step to characterize potential functional sites and core structural regions. Here, we developed RNAsnap, the first machine-learning method trained on protein-bound RNA structures for solvent accessibility prediction. Built on sequence profiles from multiple sequence alignment (RNAsnap-prof), the method provided robust prediction in fivefold cross-validation and an independent test (Pearson correlation coefficients, r, between predicted and actual ASA values are 0.66 and 0.63, respectively). Application of the method to 6178 mRNAs revealed its positive correlation to mRNA accessibility by dimethyl sulphate (DMS) experimentally measured in vivo (r = 0.37) but not in vitro (r = 0.07), despite the lack of training on mRNAs and the fact that DMS accessibility is only an approximation to solvent accessibility. We further found strong association across coding and noncoding regions between predicted solvent accessibility of the mutation site of a single nucleotide variant (SNV) and the frequency of that variant in the population for 2.2 million SNVs obtained in the 1000 Genomes Project. Moreover, mapping solvent accessibility of RNAs to the human genome indicated that introns, 5' cap of 5' and 3' cap of 3' untranslated regions, are more solvent accessible, consistent with their respective functional roles. These results support conformational selections as the mechanism for the formation of RNA-protein complexes and highlight the utility of genome-scale characterization of RNA tertiary structures by RNAsnap. The server and its stand-alone downloadable version are available at http://sparks-lab.org.

Predicting protein conformational changes for unbound and homology docking: learning from intrinsic and induced flexibility

Predicting protein conformational changes from unbound structures or even homology models to bound structures remains a critical challenge for protein docking. Here we present a study directly addressing the challenge by reducing the dimensionality and narrowing the range of the corresponding conformational space. The study builds on cNMA-our new framework of partner- and contact-specific normal mode analysis that exploits encounter complexes and considers both intrinsic and induced flexibility. First, we established over a CAPRI (Critical Assessment of PRedicted Interactions) target set that the direction of conformational changes from unbound structures and homology models can be reproduced to a great extent by a small set of cNMA modes. In particular, homology-to-bound interface root-mean-square deviation (iRMSD) can be reduced by 40% on average with the slowest 30 modes. Second, we developed novel and interpretable features from cNMA and used various machine learning approaches to predict the extent of conformational changes. The models learned from a set of unbound-to-bound conformational changes could predict the actual extent of iRMSD with errors around 0.6 Å for unbound proteins in a held-out benchmark subset, around 0.8 Å for unbound proteins in the CAPRI set, and around 1 Å even for homology models in the CAPRI set. Our results shed new insights into origins of conformational differences between homology models and bound structures and provide new support for the low-dimensionality of conformational adjustment during protein associations. The results also provide new tools for ensemble generation and conformational sampling in unbound and homology docking. Proteins 2017; 85:544-556. © 2016 Wiley Periodicals, Inc.

Protein-protein and peptide-protein docking and refinement using ATTRACT in CAPRI

The ATTRACT coarse-grained docking approach in combination with various types of atomistic, flexible refinement methods has been applied to predict protein-protein and peptide-protein complexes in CAPRI rounds 28-36. For a large fraction of CAPRI targets (12 out of 18), at least one model of acceptable or better quality was generated, corresponding to a success rate of 67%. In particular, for several peptide-protein complexes excellent predictions were achieved. In several cases, a combination of template-based modeling and extensive molecular dynamics-based refinement yielded medium and even high quality solutions. In one particularly challenging case, the structure of an ubiquitylation enzyme bound to the nucleosome was correctly predicted as a set of acceptable quality solutions. Based on the experience with the CAPRI targets, new interface refinement approaches and methods for ab-initio peptide-protein docking have been developed. Failures and possible improvements of the docking method with respect to scoring and protein flexibility will also be discussed. Proteins 2017; 85:391-398. © 2016 Wiley Periodicals, Inc.