Recognition and cooperation between the ATP-dependent RNA helicase RhlB and ribonuclease RNase E

The Escherichia coli protein RhlB is an ATP-dependent motor that unfolds structured RNA for destruction by partner ribonucleases. In E. coli, and probably many other related gamma-proteobacteria, RhlB associates with the essential endoribonuclease RNase E as part of the multi-enzyme RNA degradosome assembly. The interaction with RNase E boosts RhlB's ATPase activity by an order of magnitude. Here, we examine the origins and implications of this effect. The location of the interaction sites on both RNase E and RhlB are refined and analysed using limited protease digestion, domain cross-linking and homology modelling. These data indicate that RhlB's carboxy-terminal RecA-like domain engages a segment of RNase E that is no greater than 64 residues. The interaction between RhlB and RNase E has two important consequences: first, the interaction itself stimulates the unwinding and ATPase activities of RhlB; second, RhlB gains proximity to two RNA-binding sites on RNase E, with which it cooperates to unwind RNA. Our homology model identifies a pattern of residues in RhlB that may be key for recognition of RNase E and which may communicate the activating effects. Our data also suggest that the association with RNase E may partially repress the RNA-binding activity of RhlB. This repression may in fact permit the interplay of the helicase and adjacent RNA binding segments as part of a process that steers substrates to either processing or destruction, depending on context, within the RNA degradosome assembly.     

Flexible relaxation of rigid-body docking solutions

Molecular Dynamics (MD) simulations have been performed on a set of rigid-body docking poses, carried out over 25 protein-protein complexes. The results show that fully flexible relaxation increases the fraction of native contacts (NC) by up to 70% for certain docking poses. The largest increase in the fraction of NC is observed for docking poses where anchor residues are able to sample their bound conformation. For each MD simulation, structural snap-shots were clustered and the centre of each cluster used as the MD-relaxed docking pose. A comparison between two energy-based scoring schemes, the first calculated for the MD-relaxed poses, the second for energy minimized poses, shows that the former are better in ranking complexes with large hydrophobic interfaces. Furthermore, complexes with large interfaces are generally ranked well, regardless of the type of relaxation method chosen, whereas complexes with small hydrophobic interfaces remain difficult to rank. In general, the results indicate that current force-fields are able to correctly describe direct intermolecular interactions between receptor and ligand molecules. However, these force-fields still fail in cases where protein-protein complexes are stabilized by subtle energy contributions.     

Structural evidence for a functionally relevant second camphor binding site in P450cam: model for substrate entry into a P450 active site

P450cam has long served as a prototype for the cytochrome P450 (CYP) gene family. But, little is known about how substrate enters its active site pocket, and how access is achieved in a way that minimizes exposure of the reactive heme. We hypothesize that P450cam may first bind substrate transiently near the mobile F-G helix that covers the active site pocket. Such a two-step binding process is kinetically required if P450cam rarely populates an open conformation-as suggested by previous literature and the inability to obtain a crystal structure of P450cam in an open conformation. Such a mechanism would minimize exposure of the heme by allowing P450cam to stay in a closed conformation as long as possible, since only brief flexing into an open conformation would be required to allow substrate entry. To test this model, we have attempted to dock a second camphor molecule into the crystal structure of camphor-bound P450cam. The docking identified only one potential entry site pocket, a well-defined cavity on the F-helix side of the F-G flap, 16 A from the heme iron. Location of this entry site pocket is consistent with our NMR T1 relaxation-based measurements of distances for a camphor that binds in fast exchange (active site camphor is known to bind in slow exchange). Presence of a second camphor binding site is also confirmed with [(1)H-(13)C] HSQC titrations of (13)CH3-threonine labeled P450cam. To confirm that camphor can bind outside of the active site pocket, (13)CH3-S-pyridine was bound to the heme iron to physically block the active site, and to serve as an NMR chemical shift probe. Titration of this P450cam-pyridine complex confirms that camphor can bind to a site outside the active site pocket, with an estimated Kd of 43 microM. The two-site binding model that is proposed based on these data is analogous to that recently proposed for CYP3A4, and is consistent with recent crystal structures of P450cam bound to tethered-substrates, which force a partially opened conformation.     

Critical evaluation of methods to incorporate entropy loss upon binding in high-throughput docking

Proper accounting of the positional/orientational/conformational entropy loss associated with protein-ligand binding is important to obtain reliable predictions of binding affinity. Herein, we critically examine two simplified statistical mechanics-based approaches, namely a constant penalty per rotor method, and a more rigorous method, referred to here as the partition function-based scoring (PFS) method, to account for such entropy losses in high-throughput docking calculations. Our results on the estrogen receptor beta and dihydrofolate reductase proteins demonstrate that, while the constant penalty method over-penalizes molecules for their conformational flexibility, the PFS method behaves in a more "DeltaG-like" manner by penalizing different rotors differently depending on their residual entropy in the bound state. Furthermore, in contrast to no entropic penalty or the constant penalty approximation, the PFS method does not exhibit any bias towards either rigid or flexible molecules in the hit list. Preliminary enrichment studies using a lead-like random molecular database suggest that an accurate representation of the "true" energy landscape of the protein-ligand complex is critical for reliable predictions of relative binding affinities by the PFS method.     

Molecular definitions of fatty acid hydroxylases in Arabidopsis thaliana

Towards defining the function of Arabidopsis thaliana fatty acid hydroxylases, five members of the CYP86A subfamily have been heterologously expressed in baculovirus-infected Sf9 cells and tested for their ability to bind a range of fatty acids including unsubstituted (lauric acid (C12:0) and oleic acid (C18:1)) and oxygenated (9,10-epoxystearic acid and 9,10-dihydroxystearic acid). Comparison between these five P450s at constant P450 content over a range of concentrations for individual fatty acids indicates that binding of different fatty acids to CYP86A2 always results in a higher proportion of high spin state heme than binding titrations conducted with CYP86A1 or CYP86A4. In comparison to these three, CYP86A7 and CYP86A8 produce extremely low proportions of high spin state heme even with the most effectively bound fatty acids. In addition to their previously demonstrated lauric acid hydroxylase activities, all CYP86A proteins are capable of hydroxylating oleic acid but not oxygenated 9,10-epoxystearic acid. Homology models have been built for these five enzymes that metabolize unsubstituted fatty acids and sometimes bind oxygenated fatty acids. Comparison of the substrate binding modes and predicted substrate access channels indicate that all use channel pw2a consistent with the crystal structures and models of other fatty acid-metabolizing P450s in bacteria and mammals. Among these P450s, those that bind internally oxygenated fatty acids contain polar residues in their substrate binding cavity that help stabilize these charged/polar groups within their largely hydrophobic catalytic site.     

A novel approach to local similarity of protein binding sites substantially improves computational drug design results

We present a novel notion of binding site local similarity based on the analysis of complete protein environments of ligand fragments. Comparison of a query protein binding site (target) against the 3D structure of another protein (analog) in complex with a ligand enables ligand fragments from the analog complex to be transferred to positions in the target site, so that the complete protein environments of the fragment and its image are similar. The revealed environments are similarity regions and the fragments transferred to the target site are considered as binding patterns. The set of such binding patterns derived from a database of analog complexes forms a cloud-like structure (fragment cloud), which is a powerful tool for computational drug design. It has been shown on independent test sets that the combined use of a traditional energy-based score together with the cloud-based score responsible for the quality of embedding of a ligand into the fragment cloud improves the self-docking and screening results dramatically. The usage of a fragment cloud as a source of positioned molecular fragments fitting the binding protein environment has been validated by reproduction of experimental ligand optimization results.     

Strategies to search and design stabilizers of protein-protein interactions: a feasibility study

Since protein-protein interactions play a pivotal role in the communication on the molecular level in virtually every biological system and process, the search and design for modulators of such interactions is of utmost importance. In recent years many inhibitors for specific protein-protein interactions have been developed, however, in only a few cases, small and druglike molecules are able to interfere in the complex formation of proteins. On the other hand, there are several small molecules known to modulate protein-protein interactions by means of stabilizing an already assembled complex. To achieve this goal, a ligand is binding to a pocket, which is located rim-exposed at the interface of the interacting proteins, for example as the phytotoxin Fusicoccin, which stabilizes the interaction of plant H+-ATPase and 14-3-3 protein by nearly a factor of 100. To suggest alternative leads, we performed a virtual screening campaign to discover new molecules putatively stabilizing this complex. Furthermore, we screen a dataset of 198 transient recognition protein-protein complexes for cavities, which are located rim-exposed at their interfaces. We provide evidence for high similarity between such rim-exposed cavities and usual ligands accommodating active sites of enzymes. This analysis suggests that rim-exposed cavities at protein-protein interfaces are druggable binding sites. Therefore, the principle of stabilizing protein-protein interactions seems to be a promising alternative to the approach of the competitive inhibition of such interactions by small molecules.     

A protein-specifically adapted scoring function for the reranking of docking solutions

In this work, we developed a protein-specifically adapted scoring function and applied it to the reranking of protein-protein docking solutions generated with a conventional docking program. The approach was validated using experimentally determined structures of the bacterial HPr-protein in complex with four structurally nonhomologous binding partners as an example. A sufficiently large data basis for the generation of protein-specifically adapted pair potentials was generated by modeling all orthologous complexes for each type of interaction resulting in a total of 224 complexes. The parameters for potential generation were systematically varied and resulted in a total of 66,132 different scoring functions that were tested for their ability of successful reranking of 1000 docking solutions generated from modeled structures of the unbound binding partners. Parameters that proved critical for the generation of good scoring functions were the distance cutoff used for the generation of the pair potential, and an additional cutoff that allows a proper weighting of conserved and nonconserved contacts in the interface. Compared to the original scoring function, application of this novel type of scoring functions resulted in a significant accumulation of acceptable docking solutions within the first 10 ranks. Depending on the type of complex investigated one to five acceptable complex geometries are found among the 10 highest-ranked solutions and for three of the four systems tested, an acceptable solution was placed on the first rank.     

Invariant Ser211 is involved in the catalysis of PD-L4, type I RIP from Phytolacca dioica leaves

Multiple sequence alignment analysis of ribosome inactivating proteins (RIPs) has revealed the occurrence of an invariant seryl residue in proximity of the catalytic tryptophan. The involvement of this seryl residue in the catalytic mechanism of RIPs was investigated by site-directed mutagenesis in PD-L4, type 1 RIP isolated from Phytolacca dioica leaves. We show that the replacement of Ser211 with Ala apparently does not influence the N-beta-glycosidase activity on ribosomes (determined as IC(50) in a cell-free system), but it reduces the adenine polynucleotide glycosylase activity (APG), assayed spectrophotometrically on other substrates such as DNA, rRNA, and poly(A). The ability of PD-L4 to deadenylate polynucleotides appears more sensitive to the Ser211Ala replacement when poly(A) is used as substrate, as only 33% activity is retained by the mutant, while with more complex and heterogeneous substrates such as DNA and rRNA, its APG activity is 73% and 66%, respectively. While the mutated protein shows a conserved secondary structure by CD, it also exhibits a remarkably enhanced tryptophan fluorescence. This indicates that, although the overall protein tridimensional structure is maintained, removal of the hydroxyl group locally affects the environment of a Trp residue. Modelling and docking analyses confirm the interaction between Ser211 and Trp207, which is located within the active site, thus affecting RIP adenine polynucleotide glycosylase activity. Data accumulated so far confirm the potential involvement of Ser211 in the catalytic mechanism of type 1 RIP PD-L4 and a possible role in stabilizing the conformation of Trp207 side chain, which participates actively in the protein enzymatic activity.     

Ensemble docking of multiple protein structures: considering protein structural variations in molecular docking

One approach to incorporate protein flexibility in molecular docking is the use of an ensemble consisting of multiple protein structures. Sequentially docking each ligand into a large number of protein structures is computationally too expensive to allow large-scale database screening. It is challenging to achieve a good balance between docking accuracy and computational efficiency. In this work, we have developed a fast, novel docking algorithm utilizing multiple protein structures, referred to as ensemble docking, to account for protein structural variations. The algorithm can simultaneously dock a ligand into an ensemble of protein structures and automatically select an optimal protein structure that best fits the ligand by optimizing both ligand coordinates and the conformational variable m, where m represents the m-th structure in the protein ensemble. The docking algorithm was validated on 10 protein ensembles containing 105 crystal structures and 87 ligands in terms of binding mode and energy score predictions. A success rate of 93% was obtained with the criterion of root-mean-square deviation <2.5 A if the top five orientations for each ligand were considered, comparable to that of sequential docking in which scores for individual docking are merged into one list by re-ranking, and significantly better than that of single rigid-receptor docking (75% on average). Similar trends were also observed in binding score predictions and enrichment tests of virtual database screening. The ensemble docking algorithm is computationally efficient, with a computational time comparable to that for docking a ligand into a single protein structure. In contrast, the computational time for the sequential docking method increases linearly with the number of protein structures in the ensemble. The algorithm was further evaluated using a more realistic ensemble in which the corresponding bound protein structures of inhibitors were excluded. The results show that ensemble docking successfully predicts the binding modes of the inhibitors, and discriminates the inhibitors from a set of noninhibitors with similar chemical properties. Although multiple experimental structures were used in the present work, our algorithm can be easily applied to multiple protein conformations generated by computational methods, and helps improve the efficiency of other existing multiple protein structure(MPS)-based methods to accommodate protein flexibility.