T cell receptor-ligand interactions: a conformational preequilibrium or an induced fit.

Kinetic parameters of T cell receptor (TCR) interactions with its ligand have been proposed to control T cell activation. Analysis of kinetic data obtained has so far produced conflicting insights; here, we offer a consideration of this problem. As a model system, association and dissociation of a soluble TCR (sT1) and its specific ligand, an azidobenzoic acid derivative of the peptide SYIPSAEK-(ABA)I (residues 252-260 from Plasmodium berghei circumsporozoite protein), bound to class I MHC H-2K(d)-encoded molecule (MHCp) were studied by surface plasmon resonance. The association time courses exhibited biphasic patterns. The fast and dominant phase was assigned to ligand association with the major fraction of TCR molecules, whereas the slow component was attributed to the presence of traces of TCR dimers. The association rate constant derived for the fast phase, assuming a reversible, single-step reaction mechanism, was relatively slow and markedly temperature-dependent, decreasing from 7.0 x 10(3) at 25 degrees C to 1.8 x 10(2) M(-1).s(-1) at 4 degrees C. Hence, it is suggested that these observed slow rate constants are the result of unresolved elementary steps of the process. Indeed, our analysis of the kinetic data shows that the time courses of TCR-MHCp interaction fit well to two different, yet closely related mechanisms, where an induced fit or a preequilibrium of two unbound TCR conformers are operational. These mechanisms may provide a rationale for the reported conformational flexibility of the TCR and its unusual ligand recognition properties, which combine high specificity with considerable crossreactivity.

Induced-fit binding of π-electron donor substrates to macrocyclic aromatic ether-imide-sulfones: a versatile approach to molecular assembly

Novel macrocyclic receptors which bind electron-donor aromatic substrates via π-stacking donor- acceptor interactions are obtained by cyclo-imidization of an amine-functionalized arylether-sulfone with pyromellitic- and 1,4,5,8-naphthalene-tetracarboxylic dianhydrides. These macrocycles complex with a wide variety of π-donor substrates including tetrathiafulvalene, naphthalene, anthracene, pyrene, perylene, and functional derivatives of these polycyclic hydrocarbons. The resulting supramolecular assemblies range from simple 1:1 complexes, to [2]- and [3]-pseudorotaxanes, and even (as a result of crystallographic disorder) an apparent polyrotaxane. Direct, five-component self-assembly of a metal-centred [3]pseudorotaxane is also observed, on complexation of a macrocyclic ether-imide with 8-hydroxyquinoline in the presence of palladium(II) ions. Binding studies in solution were carried out by 1H NMR and UV-visible spectroscopy, and the stoichiometries of binding were confirmed by Job plots based on charge-transfer absorption bands. The highest association constants are found for strong π-donor guests with large surface-areas, notably perylene and 1-hydroxypyrene, for which Ka values of 1.4 x 103 and 2.3 x 103 M-1 respectively are found. Single crystal X-ray analyses of the receptors and their derived complexes reveal large, induced-fit distortions of the macrocyclic frameworks as a result of complexation. These structures provide compelling evidence for the existence of strong, attractive forces between the electronically-complementary aromatic π-systems of host and guest.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

T-Cell Receptors Binding Orientation over Peptide/MHC Class I Is Driven by Long-Range Interactions.

Crystallographic data about T-Cell Receptor - peptide - major histocompatibility complex class I (TCRpMHC) interaction have revealed extremely diverse TCR binding modes triggering antigen recognition. Understanding the molecular basis that governs TCR orientation over pMHC is still a considerable challenge. We present a simplified rigid approach applied on all non-redundant TCRpMHC crystal structures available. The CHARMM force field in combination with the FACTS implicit solvation model is used to study the role of long-distance interactions between the TCR and pMHC. We demonstrate that the sum of the coulomb interactions and the electrostatic solvation energies is sufficient to identify two orientations corresponding to energetic minima at 0° and 180° from the native orientation. Interestingly, these results are shown to be robust upon small structural variations of the TCR such as changes induced by Molecular Dynamics simulations, suggesting that shape complementarity is not required to obtain a reliable signal. Accurate energy minima are also identified by confronting unbound TCR crystal structures to pMHC. Furthermore, we decompose the electrostatic energy into residue contributions to estimate their role in the overall orientation. Results show that most of the driving force leading to the formation of the complex is defined by CDR1,2/MHC interactions. This long-distance contribution appears to be independent from the binding process itself, since it is reliably identified without considering neither short-range energy terms nor CDR induced fit upon binding. Ultimately, we present an attempt to predict the TCR/pMHC binding mode for a TCR structure obtained by homology modeling. The simplicity of the approach and the absence of any fitted parameters make it also easily applicable to other types of macromolecular protein complexes.

Lack of Enantioselectivity in the SULT1A3-catalyzed Sulfoconjugation of Normetanephrine Enantiomers: An In Vitro and Computational Study.

(1R)-Normetanephrine is the natural stereoisomeric substrate for sulfotransferase 1A3 (SULT1A3)-catalyzed sulfonation. Nothing appears known on the enantioselectivity of the reaction despite its potential significance in the metabolism of adrenergic amines and in clinical biochemistry. We confronted the kinetic parameters of the sulfoconjugation of synthetic (1R)-normetanephrine and (1S)-normetanephrine by recombinant human SULT1A3 to a docking model of each normetanephrine enantiomer with SULT1A3 and the 3'-phosphoadenosine-5'-phosphosulfate cofactor on the basis of molecular modeling and molecular dynamics simulations of the stability of the complexes. The K(M) , V(max) , and k(cat) values for the sulfonation of (1R)-normetanephrine, (1S)-normetanephrine, and racemic normetanephrine were similar. In silico models were consistent with these findings as they showed that the binding modes of the two enantiomers were almost identical. In conclusion, SULT1A3 is not substrate-enantioselective toward normetanephrine, an unexpected finding explainable by a mutual adaptability between the ligands and SULT1A3 through an "induced-fit model" in the catalytic pocket. Chirality, 00:000-000, 2012.© 2012 Wiley Periodicals, Inc.

Um paradigma da química medicinal: a flexibilidade dos ligantes e receptores

In general, molecular modeling techniques applied in medicinal chemistry have been static and drug based. However the active site geometry and the intrinsic flexibility of both receptor and ligand are fundamental properties for molecular recognition and drug action. As a consequence, the use of dynamic models to describe the ligand-receptor complex is becoming a more common procedure. In this work we discuss the relevance of considering the receptor structure in medicinal chemistry studies as well as the flexibility of the ligand-receptor complex.

Molecular Interactions of Prodiginines with the BH3 Domain of Anti-Apoptotic Bcl-2 Family Members.

Prodigiosin and obatoclax, members of the prodiginines family, are small molecules with anti-cancer properties that are currently under preclinical and clinical trials. The molecular target(s) of these agents, however, is an open question. Combining experimental and computational techniques we find that prodigiosin binds to the BH3 domain in some BCL-2 protein families, which play an important role in the apoptotic programmed cell death. In particular, our results indicate a large affinity of prodigiosin for MCL-1, an anti-apoptotic member of the BCL-2 family. In melanoma cells, we demonstrate that prodigiosin activates the mitochondrial apoptotic pathway by disrupting MCL-1/BAK complexes. Computer simulations with the PELE software allow the description of the induced fit process, obtaining a detailed atomic view of the molecular interactions. These results provide new data to understand the mechanism of action of these molecules, and assist in the development of more specific inhibitors of anti-apoptotic BCL-2 proteins.

Structural studies on enzymes of biotechnological and biomedical interest

Structural studies of proteins aim at elucidating the atomic details of molecular interactions in biological processes of living organisms. These studies are particularly important in understanding structure, function and evolution of proteins and in defining their roles in complex biological settings. Furthermore, structural studies can be used for the development of novel properties in biomolecules of environmental, industrial and medical importance. X-ray crystallography is an invaluable tool to obtain accurate and precise information about the structure of proteins at the atomic level. Glutathione transferases (GSTs) are amongst the most versatile enzymes in nature. They are able to catalyze a wide variety of conjugation reactions between glutathione (GSH) and non-polar components containing an electrophilic carbon, nitrogen or sulphur atom. Plant GSTs from the Tau class (a poorly characterized class) play an important role in the detoxification of xenobiotics and stress tolerance. Structural studies were performed on a Tau class fluorodifen-inducible glutathione transferase from Glycine max (GmGSTU4-4) complexed with GSH (2.7 Å) and a product analogue Nb-GSH (1.7 Å). The three-dimensional structure of the GmGSTU4-4-GSH complex revealed that GSH binds in different conformations in the two subunits of the dimer: in an ionized form in one subunit and a non-ionized form in the second subunit. Only the ionized form of the substrate may lead to the formation of a catalytically competent complex. Structural comparison between the GSH and Nb-GSH bound complexes revealed significant differences with respect to the hydrogen-bonding, electrostatic interaction pattern, the upper part of -helix H4 and the C-terminus of the enzyme. These differences indicate an intrasubunit modulation between the G-and Hsites suggesting an induced-fit mechanism of xenobiotic substrate binding. A novel binding site on the surface of the enzyme was also revealed. Bacterial type-II L-asparaginases are used in the treatment of haematopoietic diseases such as acute lymphoblastic leukaemia (ALL) and lymphomas due to their ability to catalyze the conversion of L-asparagine to L-aspartate and ammonia. Escherichia coli and Erwinia chrysanthemi asparaginases are employed for the treatment of ALL for over 30 years. However, serious side-effects affecting the liver and pancreas have been observed due to the intrinsic glutaminase activity of the administered enzymes. Structural studies on Helicobacter pylori L-asparaginase (HpA) were carried out in an effort to discover novel L-asparaginases with potential chemotherapeutic utility in ALL treatment. Detailed analysis of the active site geometry revealed structurally significant differences between HpA and other Lasparaginases that may be important for the biological activities of the enzyme and could be further exploited in protein engineering efforts.

Statistical coupling analysis of aspartic proteinases based on crystal structures of the Trichoderma reesei enzyme and its complex with pepstatin A

The crystal structures of an aspartic proteinase from Trichoderma reesei (TrAsP) and of its complex with a competitive inhibitor, pepstatin A, were solved and refined to crystallographic R-factors of 17.9% (R(free)=21.2%) at 1.70 angstrom resolution and 15.81% (R(free) = 19.2%) at 1.85 angstrom resolution, respectively. The three-dimensional structure of TrAsP is similar to structures of other members of the pepsin-like family of aspartic proteinases. Each molecule is folded in a predominantly beta-sheet bilobal structure with the N-terminal and C-terminal domains of about the same size. Structural comparison of the native structure and the TrAsP-pepstatin complex reveals that the enzyme undergoes an induced-fit, rigid-body movement upon inhibitor binding, with the N-terminal and C-terminal lobes tightly enclosing the inhibitor. Upon recognition and binding of pepstatin A, amino acid residues of the enzyme active site form a number of short hydrogen bonds to the inhibitor that may play an important role in the mechanism of catalysis and inhibition. The structures of TrAsP were used as a template for performing statistical coupling analysis of the aspartic protease family. This approach permitted, for the first time, the identification of a network of structurally linked residues putatively mediating conformational changes relevant to the function of this family of enzymes. Statistical coupling analysis reveals coevolved continuous clusters of amino acid residues that extend from the active site into the hydrophobic cores of each of the two domains and include amino acid residues from the flap regions, highlighting the importance of these parts of the protein for its enzymatic activity. (C) 2008 Elsevier Ltd. All rights reserved.

Inhibition of myotoxic activity of Bothrops asper myotoxin II by the anti-trypanosomal drug surarnin

Suramin, a synthetic polysulfonated compound, developed initially for the treatment of African trypanosomiasis and onchocerciasis, is currently used for the treatment of several medically relevant disorders. Suramin, heparin, and other polyanions inhibit the myotoxic activity of Lys49 phospholipase A(2) analogues both in vitro and in vivo, and are thus of potential importance as therapeutic agents in the treatment of viperid snake bites. Due to its conformational flexibility around the single bonds that link the central phenyl rings to the secondary amide backbone, the symmetrical suramin molecule binds by an induced-fit mechanism complementing the hydrophobic surfaces of the dimer and adopts a novel conformation that lacks C2 symmetry in the dimeric crystal structure of the suramin-Bothrops asper myotoxin II complex. The simultaneous binding of suramin at the surfaces of the two monomers partially restricts access to the nominal active sites and significantly changes the overall charge of the interfacial recognition face of the protein, resulting in the inhibition of myotoxicity. (c) 2005 Elsevier Ltd. All rights reserved.

Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

Novel, potent, and radio-iodinatable somatostatin receptor 1 (sst1) selective analogues

The proposed sst(1) pharmacophore (J. Med. Chem. 2005, 48, 523-533) derived from the NMR structures of a family of mono- and dicyclic undecamers was used to design octa-, hepta-, and hexamers with high affinity and selectivity for the somatostatin sst(1) receptor. These compounds were tested for their in vitro binding properties to all five somatostatin (SRIF) receptors using receptor autoradiography; those with high SRIF receptor subtype 1 (sst(1)) affinity and selectivity were shown to be agonists when tested functionally in a luciferase reporter gene assay. Des-AA(1,4-6,10,12,13)-[DTyr(2),DAgl(NMe,2naphthoyl)(8),IAmp(9)]-SRIF-Thr-NH(2) (25) was radio-iodinated ((125)I-25) and specifically labeled sst(1)-expressing cells and tissues. 3D NMR structures were calculated for des-AA(1,4-6,10,12,13)-[DPhe(2),DTrp(8),IAmp(9)]-SRIF-Thr-NH(2) (16), des-AA(1,2,4-6,10,12,13)-[DAgl(NMe,2naphthoyl)(8),IAmp(9)]-SRIF-Thr-NH(2) (23), and des-AA(1,2,4-6,10,12,13)-[DAgl(NMe,2naphthoyl)(8),IAmp(9),Tyr(11)]-SRIF-NH(2) (27) in DMSO. Though the analogues have the sst(1) pharmacophore residues at the previously determined distances from each other, the positioning of the aromatic residues in 16, 23, and 27 is different from that described earlier, suggesting an induced fit mechanism for sst(1) binding of these novel, less constrained sst(1)-selective family members.

Anti-peptide aptamers recognize amino acid sequence and bind a protein epitope.

In vitro selection of nucleic acid binding species (aptamers) is superficially similar to the immune response. Both processes produce biopolymers that can recognize targets with high affinity and specificity. While antibodies are known to recognize the sequence and conformation of protein surface features (epitopes), very little is known about the precise interactions between aptamers and their epitopes. Therefore, aptamers that could recognize a particular epitope, a peptide fragment of human immunodeficiency virus type I Rev, were selected from a random sequence RNA pool. Several of the selected RNAs could bind the free peptide more tightly than a natural RNA ligand, the Rev-binding element. In accord with the hypothesis that protein and nucleic acid binding cusps are functionally similar, interactions between aptamers and the peptide target could be disrupted by sequence substitutions. Moreover, the aptamers appeared to be able to bind peptides with different solution conformations, implying an induced fit mechanism for binding. Just as anti-peptide antibodies can sometimes recognize the corresponding epitope when presented in a protein, the anti-peptide aptamers were found to specifically bind to Rev.

Countering cooperative effects in protease inhibitors using constrained b-strand-mimicking templates in focused combinatorial libraries

A major problem in de novo design of enzyme inhibitors is the unpredictability of the induced fit, with the shape of both ligand and enzyme changing cooperatively and unpredictably in response to subtle structural changes within a ligand. We have investigated the possibility of dampening the induced fit by using a constrained template as a replacement for adjoining segments of a ligand. The template preorganizes the ligand structure, thereby organizing the local enzyme environment. To test this approach, we used templates consisting of constrained cyclic tripeptides, formed through side chain to main chain linkages, as structural mimics of the protease-bound extended beta-strand conformation of three adjoining amino acid residues at the N- or C-terminal sides of the scissile bond of substrates. The macrocyclic templates were derivatized to a range of 30 structurally diverse molecules via focused combinatorial variation of nonpeptidic appendages incorporating a hydroxyethylamine transition-state isostere. Most compounds in the library were potent inhibitors of the test protease (HIV-1 protease). Comparison of crystal structures for five protease-inhibitor complexes containing an N-terminal macrocycle and three protease-inhibitor complexes containing a C-terminal macrocycle establishes that the macrocycles fix their surrounding enzyme environment, thereby permitting independent variation of acyclic inhibitor components with only local disturbances to the protease. In this way, the location in the protease of various acyclic fragments on either side of the macrocyclic template can be accurately predicted. This type of templating strategy minimizes the problem of induced fit, reducing unpredictable cooperative effects in one inhibitor region caused by changes to adjacent enzyme-inhibitor interactions. This idea might be exploited in template-based approaches to inhibitors of other proteases, where a beta-strand mimetic is also required for recognition, and also other protein-binding ligands where different templates may be more appropriate.