Scoring function based approach for locating binding sites and understanding recognition mechanism of protein-DNA complexes

Protein-DNA recognition plays an essential role in the regulation of gene expression. Understanding the recognition mechanism of protein-DNA complexes is a challenging task in molecular and computational biology. In this work, a scoring function based approach has been developed for identifying the binding sites and delineating the important residues for binding in protein-DNA complexes. This approach considers both the repulsive interactions and the effect of distance between atoms in protein and DNA. The results showed that positively charged, polar, and aromatic residues are important for binding. These residues influence the formation of electrostatic, hydrogen bonding, and stacking interactions. Our observation has been verified with experimental binding specificity of protein-DNA complexes and found to be in good agreement with experiments. The comparison of protein-RNA and protein-DNA complexes reveals that the contribution of phosphate atoms in DNA is twice as large as in protein-RNA complexes. Furthermore, we observed that the positively charged, polar, and aromatic residues serve as hotspot residues in protein-RNA complexes, whereas other residues also altered the binding specificity in protein-DNA complexes. Based on the results obtained in the present study and related reports, a plausible mechanism has been proposed for the recognition of protein-DNA complexes.     

Molecular recognition. Design of new receptors for complexation and catalysis

Abstract

In recent years there has been intense activity in the design of synthetic molecules capable of enzyme-like recognition and binding of small substrates.¹ Two fundamental approaches have been taken. The first has generally involved non-directional binding forces (such as solvophobic, π-stacking and dispersion interactions) in water-soluble cyclophane frameworks.² This approach led to extremely important quantitative insights into the hydrophobic effect and the enthalpic and entropic contributions of solvent reorganization to binding.³ However, the weakly oriented nature of the binding interactions has resulted in only moderate substrate selectivity beyond the shape recognition permitted by the cavity. In nature such selectivity is a prerequisite for the chiral recognition and catalytic activity of enzymes and is achieved by hydrogen bonding and electrostatic interactions. The second major approach to artificial receptors makes use of these more directional interactions by incorporating several hydrogen bonding groups into a cleft or cavity of defined geometry.⁴ The resulting hosts form strong and selective complexes to those substrates with complementary shape and hydrogen bonding characteristics.⁵ In these cases, however, the binding free energy is solvent dependent, diminishing to zero as the polarity of the medium increases, due to the strong solvation of the hydrogen bonding sites. A central goal in contemporary molecular recognition research must be to develop receptors that effectively use directed hydrogen bonding interactions in competitive solvents. Success will probably require combining strong (possibly charged) hydrogen bonding groups with hydrophobic sites capable not only of effective apolar association with the substrate but also of protecting the polar sites from full solvation.     

Electrostatics in protein binding and function

Protein electrostatic properties stem from the proportion and distribution of polar and charged residues. Polar and charged residues regulate the electrostatic properties by forming short-range interactions, like salt-bridges and hydrogen-bonds, and by defining the over-all electrostatic environment in the protein. Electrostatics play a major role in defining the mechanisms of protein-protein complex formation, molecular recognitions, thermal stabilities, conformational adaptabilities and protein movements. For example:- Functional hinges, or flexible regions of the protein, lack short-range electrostatic interactions; Thermophilic proteins have higher electrostatic interactions than their mesophilic counter parts; Increase in binding specificity and affinity involve optimization of electrostatics; High affinity antibodies have higher, and stronger, electrostatic interactions with their antigens; Rigid parts of proteins have higher and stronger electrostatic interactions. In this review we address the significance of electrostatics in protein folding, binding and function. We discuss that the electrostatic properties are evolutionally selected by a protein to perform an specific function. We also provide bona fide examples to illustrate this. Additionally, using continuum electrostatic and molecular dynamics approaches we show that the "hot-spot" inter-molecular interactions in a very specific antibody-antigen binding are mainly established through charged residues. These "hot-spot" molecular interactions stay intact even during high temperature molecular dynamics simulations, while the other inter-molecular interactions, of lesser functional significance, disappear. This further corroborates the significance of charge-charge interactions in defining binding mechanisms. High affinity binding frequently involves "electrostatic steering". The forces emerge from over-all electrostatic complementarities and by the formation of charged and polar interactions. We demonstrate that although the high affinity binding of barnase-barstar and anti-hen egg white lysozyme (HEL) antibody-HEL complexes involve different molecular mechanisms, it is electrostatically regulated in both the cases. These observations, and several other studies, suggest that a fine tuning of local and global electrostatic properties are essential for protein binding and function.     

蛋白质-RNA相互作用界面预测与设计

蛋白质-RNA之间的相互作用是蛋白质在细胞里面行使功能的重要方式之一. 结构生物学家利用实验手段可以得到蛋白质-RNA复合物的三维结构, 通过原子水平的晶体结构来解释蛋白质与RNA的识别过程. 但实验取得蛋白质-RNA的复合物结构非常困难, 耗钱、耗时, 同时受限于其相互作用强度. 因而利用理论的方法对蛋白质-RNA相互作用界面进行预测与设计在生物医学研究中十分重要. 本文主要综述了近期蛋白质-RNA相互作用界面预测与设计方面的进展, 包括以下几个方面: (1) 蛋白质-RNA分子对接算法以及对接前后存在的构象变化的处理; (2) 蛋白质-RNA 识别机制的研究; (3) 基于蛋白质-RNA 相互作用界面的分子设计. 蛋白质-RNA分子对接算法逐步完善将有助于我们对大量未知功能的蛋白质与RNA进行功能注释, 而基于生物大分子相互作用界面的分子设计将在药物设计领域中有广阔的应用前景.     

Interactions between hydrophobic and ionic solutes in aqueous guanidinium chloride and urea solutions: lessons for protein denaturation mechanism

In order to clarify the mechanism of denaturant-induced unfolding of proteins we have calculated the interactions between hydrophobic and ionic species in aqueous guanidinium chloride and urea solutions using molecular dynamics simulations. Hydrophobic association is not significantly changed in urea or guanidinium chloride solutions. The strength of interaction between ion pairs is greatly diminished by the guanidinium ion. Although the changes in electrostatic interactions in urea are small, examination of structures, using appropriate pair functions, of urea and water around the solutes show strong hydrogen bonding between urea's carbonyl oxygen and the positively charged solute. Our results strongly suggest protein denaturation occurs by the direct interaction model according to which the most commonly used denaturants unfold proteins by altering electrostatic interactions either by solvating the charged residues or by engaging in hydrogen bonds with the protein backbone. To further validate the direct interaction model we show that, in urea and guanidinium chloride solutions, unfolding of an unusually stable helix (H1) from mouse PrPC (residues 144-153) occurs by hydrogen bonding of denaturants to charged side chains and backbone carbonyl groups.     

Aromatic Rings as Molecular Determinants for the Molecular Recognition of Protein Kinase Inhibitors

Protein kinases are key enzymes in many signal transduction pathways, and play a crucial role in cellular proliferation, differentiation, and various cell regulatory processes. However, aberrant function of kinases has been associated with cancers and many other diseases. Consequently, competitive inhibition of the ATP binding site of protein kinases has emerged as an effective means of curing these diseases. Over the past three decades, thousands of protein kinase inhibitors (PKIs) with varying molecular frames have been developed. Large-scale data mining of the Protein Data Bank resulted in a database of 2139 non-redundant high-resolution X-ray crystal structures of PKIs bound to protein kinases. This provided us with a unique opportunity to study molecular determinants for the molecular recognition of PKIs. A chemoinformatic analysis of 2139 PKIs resulted in findings that PKIs are “flat” molecules with high aromatic ring counts and low fractions of sp³ carbon. All but one PKI possessed one or more aromatic rings. More importantly, it was found that the average weighted hydrogen bond count is inversely proportional to the number of aromatic rings. Based on this linear relationship, we put forward the exchange rule of hydrogen bonding interactions and non-bonded π-interactions. Specifically, a loss of binding affinity caused by a decrease in hydrogen bonding interactions is compensated by a gain in binding affinity acquired by an increase in aromatic ring-originated non-bonded interactions (i.e., π–π stacking interactions, CH–π interactions, cation–π interactions, etc.), and vice versa. The very existence of this inverse relationship strongly suggests that both hydrogen bonding and aromatic ring-originated non-bonded interactions are responsible for the molecular recognition of PKIs. As an illustration, two representative PKI–kinase complexes were employed to examine the relative importance of different modes of non-bonded interactions for the molecular recognition of PKIs. For this purpose, two FDA-approved PKI drugs, ibrutinib and lenvatinib, were chosen. The binding pockets of both PKIs were thoroughly examined to identify all non-bonded intermolecular interactions. Subsequently, the strengths of interaction energies between ibrutinib and its interacting residues in tyrosine kinase BTK were quantified by means of the double hybrid DFT method B2PLYP. The resulting energetics for the binding of ibrutinib in tyrosine kinase BTK showed that CH–π interactions and π–π stacking interactions between aromatic rings of the drug and hydrophobic residues in its binding pocket dominate the binding interactions. Thus, this work establishes that, in addition to hydrogen bonding, aromatic rings function as important molecular determinants for the molecular recognition of PKIs. In conclusion, our findings support the following pharmacophore model for ATP-competitive kinase inhibitors: a small molecule features a scaffold of one or more aromatic rings which is linked with one or more hydrophilic functional groups. The former has the structural role of acting as a scaffold and the functional role of participating in aromatic ring-originated non-bonded interactions with multiple hydrophobic regions in the ATP binding pocket of kinases. The latter ensure water solubility and form hydrogen bonds with the hinge region and other hydrophilic residues of the ATP binding pocket.     

硫脲类阴离子受体的研究进展

阴离子识别是超分子化学研究的重要内容之一,其关键环节是构筑可识别阴离子的结合受体,后者以非共价键力如静电作用、疏水作用、氢键等与阴离子结合.本文详细评述了近5年来硫脲类阴离子识别受体的设计、结构及其阴离子识别作用的研究进展.     

(1)H NMR investigation of solvent effects in aromatic stacking interactions

One of the marquis challenges in modern Organic Chemistry concerns the design and synthesis of abiotic compounds that emulate the exquisite complex structures and/or functions of biological macromolecules. Oligomers possessing the propensity to adopt well-defined compact conformations, or foldamers, have been attained utilizing hydrogen bonding, torsional restriction, and solvophobic interactions.(1) In this laboratory, aromatic electron donor--acceptor interactions have been exploited in the design of aedamers--foldamers that adopt a novel, pleated secondary structure in aqueous solution. Herein is reported detailed (1)H NMR binding studies of aedamer monomers that were carried out in solvents and solvent mixtures covering a broad polarity range. Curve-fitting analysis of the binding data using a model that incorporated the formation of higher order and self-associated complexes yielded a linear free energy relationship between the free energy of complexation and the empirical solvent polarity parameter, E(T)(30). From these studies, the association of electron-rich and electron-deficient aedamer monomers was seen to be driven primarily by hydrophobic interactions in polar solvents. However, the magnitude of these interactions is modulated to a significant extent by the geometry of the donor--acceptor complex, which, in turn, is dictated by the electrostatic complementarity between the electron-deficient and electron-rich aromatic faces of the monomers.     

Composite Aromatic Boxes for Enzymatic Transformations of Quaternary Ammonium Substrates

Cation–π interactions to cognate ligands in enzymes have key roles in ligand binding and enzymatic catalysis. We have deciphered the key functional role of both charged and aromatic residues within the choline binding subsite of CTP:phosphocholine cytidylyltransferase and choline kinase from Plasmodium falciparum. Comparison of quaternary ammonium binding site structures revealed a general composite aromatic box pattern of enzyme recognition sites, well distinguished from the aromatic box recognition site of receptors.     

Analysis and prediction of RNA-binding residues using sequence, evolutionary conservation, and predicted secondary structure and solvent accessibility

Identification and prediction of RNA-binding residues (RBRs) provides valuable insights into the mechanisms of protein-RNA interactions. We analyzed the contributions of a wide range of factors including amino acid sequence, evolutionary conservation, secondary structure and solvent accessibility, to the prediction/characterization of RBRs. Five feature sets were designed and feature selection was performed to find and investigate relevant features. We demonstrate that (1) interactions with positively charged amino acids Arg and Lys are preferred by the egatively charged nucleotides; (2) Gly provides flexibility for the RNA binding sites; (3) Glu with negatively charged side chain and several hydrophobic residues such as Leu, Val, Ala and Phe are disfavored in the RNA-binding sites; (4) coil residues, especially in long segments, are more flexible (than other secondary structures) and more likely to interact with RNA; (5) helical residues are more rigid and consequently they are less likely to bind RNA; and (6) residues partially exposed to the solvent are more likely to form RNA-binding sites. We introduce a novel sequence-based predictor of RBRs, RBRpred, which utilizes the selected features. RBRpred is comprehensively tested on three datasets with varied atom distance cutoffs by performing both five-fold cross validation and jackknife tests and achieves Matthew's correlation coefficient (MCC) of 0.51, 0.48 and 0.42, respectively. The quality is comparable to or better than that for state-of-the-art predictors that apply the distancebased cutoff definition. We show that the most important factor for RBRs prediction is evolutionary conservation, followed by the amino acid sequence, predicted secondary structure and predicted solvent accessibility. We also investigate the impact of using native vs. predicted secondary structure and solvent accessibility. The predictions are sufficient for the RBR prediction and the knowledge of the actual solvent accessibility helps in predictions for lower distance cutoffs.     

Individually double minimum-distance definition of protein–RNA binding residues and application to structure-based prediction

Identifying protein–RNA binding residues is essential for understanding the mechanism of protein–RNA interactions. So far, rigid distance thresholds are commonly used to define protein–RNA binding residues. However, after investigating 182 non-redundant protein–RNA complexes, we find that it would be unsuitable for a certain amount of complexes since the distances between proteins and RNAs vary widely. In this work, a novel definition method was proposed based on a flexible distance cutoff. This method can fully consider the individual differences among complexes by setting a variable tolerance limit of protein–RNA interactions, i.e. the double minimum-distance by which different distance thresholds are achieved for different complexes. In order to validate our method, a comprehensive comparison between our flexible method and traditional rigid methods was implemented in terms of interface structure, amino acid composition, interface area and interaction force, etc. The results indicate that this method is more reasonable because it incorporates the specificity of different complexes by extracting the important residues lost by rigid distance methods and discarding some redundant residues. Finally, to further test our double minimum-distance definition strategy, we developed a classifier to predict those binding sites derived from our new method by using structural features and a random forest machine learning algorithm. The model achieved a satisfactory prediction performance and the accuracy on independent data sets reaches to 85.0%. To the best of our knowledge, it is the first prediction model to define positive and negative samples using a flexible cutoff. So the comparison analysis and modeling results have demonstrated that our method would be a very promising strategy for more precisely defining protein–RNA binding sites.     

Electronic structure, molecular electrostatic potential and hydrogen bonding in DMSO-X complexes (X = ethanol, methanol and water)

In the present work, we have studied the electronic structure, molecular electrostatic potential (MEP) and hydrogen bonding in DMSO-ethanol, DMSO-methanol and DMSO-water complexes by employing the MP2 method. Different conformers were simulated on the basis of possible binding sites guided by molecular electrostatic potential topology. The stronger hydrogen bonded interaction lowers the energy of the conformer. Molecular electron density topology and natural bond orbital analysis were used to explain the strength of interactions. Experimental vibrations are also compared with the calculated normal vibrations. Blue shift is predicted for SC vibration in experimental and theoretical spectra as well. Molecular electrostatic potential and topology are used to understand the interaction strength of the conformer.     

Theoretical investigation of the hydrogen bond interactions of methanol and dimethylamine with hydrazone and its derivatives

The interactions between hydrogen bond donors (dimethylamine (DMA) and methanol (MeOH)) and acceptors (formaldehyde dimethylhydrazone, acetaldehyde N,N-dimethylhydrazone and N-nitrosodimethylamine) were theoretically investigated by DFT. The hydrogen bonding interactions were found on several bonding sites of the acceptors. The important properties of structure, binding energy, enthalpy of formation, Gibbs free energy of formation and equilibrium constant were investigated. Compared to the monomer, the DMA complexes show a small red shift of the NH-stretching vibrational transition but a significantly intensity enhancement. On the other hand, the MeOH complexes have a large red shift but a relatively small intensity enhancement of the OH-stretching transition. Atoms-in-molecules analysis revealed that several types of hydrogen bond interaction were present in the complexes. Since natural bond orbital analysis overestimated the effect of charge transfer, the more reliable localized molecular orbital energy decomposition analysis was performed and it shows that the major contribution to the total interaction energy is the electrostatic interaction. All these parameters suggest that the hydrogen bond donor strength of MeOH is substantially greater than DMA.

     

The CH–π Interaction in Protein–Carbohydrate Binding: Bioinformatics and In Vitro Quantification

The molecular recognition of carbohydrates by proteins plays a key role in many biological processes including immune response, pathogen entry into a cell, and cell–cell adhesion (e.g., in cancer metastasis). Carbohydrates interact with proteins mainly through hydrogen bonding, metal-ion-mediated interaction, and non-polar dispersion interactions. The role of dispersion-driven CH–π interactions (stacking) in protein–carbohydrate recognition has been underestimated for a long time considering the polar interactions to be the main forces for saccharide interactions. However, over the last few years it turns out that non-polar interactions are equally important. In this study, we analyzed the CH–π interactions employing bioinformatics (data mining, structural analysis), several experimental (isothermal titration calorimetry (ITC), X-ray crystallography), and computational techniques. The Protein Data Bank (PDB) has been used as a source of structural data. The PDB contains over 12 000 protein complexes with carbohydrates. Stacking interactions are very frequently present in such complexes (about 39 % of identified structures). The calculations and the ITC measurement results suggest that the CH–π stacking contribution to the overall binding energy ranges from 4 up to 8 kcal mol⁻¹. All the results show that the stacking CH–π interactions in protein–carbohydrate complexes can be considered to be a driving force of the binding in such complexes.     

Pseudoreceptor model for ryanodine derivatives at calcium release channels

This paper describes the generation of a pseudoreceptor model for ryanodine receptor (RyR) modulating ryanoids in rabbit skeletal muscle. For this purpose, the molecular modelling software PrGen was applied to correlate experimentally determined and calculated free energies of binding for a set of 15 ryanodine derivatives. The final model indicates a narrow cleft with hydrogen bond donor and acceptor capacities (represented by an Asn) as most crucial for binding the pyrrole carboxylate substituent at C3 of ryanodine. In addition, hydrophobic residues flank the aromatic pyrrole ring (Tyr, Phe, and Ile). Two of those residues (Tyr and Ile) interact with the 2-isopropyl moiety, which seems to contribute to binding. Opposite to the pyrrole locus, a second hydrophobic region (represented by a Leu) restricts ryanodine derivatives in their longitudinal axis and leads to the discrimination of equatorial and axial positioned methyl groups and of polar substituents at C9. Finally, a charged glutamate residue generates strong hydrogen bonding and electrostatic interactions with the hydroxyl groups at C10 and C15. For this binding-site model – composed of six amino acid residues – a correlation for the training set ligands of R = 0.99 (Q² = 0.975) and a root mean square (rms) deviation of 0.568 kcal/mol for the prediction of the binding energies of four test set ligands was obtained. Based on this pseudoreceptor model the putative topology of the real binding site of ryanoids will be discussed.     

Tuning the physisorption of molecular hydrogen: binding to aromatic,hetero-aromatic and metal-organic framework materials

We present a study on the binding properties of molecular hydrogen to several polar aromatic molecules and to a model for the metal-oxide corner of the metal organic framework materials recently investigated as promising supports for hydrogen storage. Density functional theory employing the Perdew Wang exchange-correlation functional and second order Møller-Plesset calculations are used to determine the equilibrium structures of complexes with molecular hydrogen and their stability. It is found that for most hetero-aromatics the edge sites for molecular hydrogen physisorption have stabilities comparable to the top sites. The DFT predicted binding energies compare favorably with those estimated at MP2 level, and get closer to the MP2 results for increased electrostatic contributions (induced by the polar aromatics) to the intermolecular interaction. Vibrational frequencies are also computed at the DFT level, and infrared activities of the H₂ stretching frequency are compared for the various complexes. Pyrrole, pyridine and n-oxide pyridine are predicted to form the more stable complexes among one-ring aromatics. The computed binding energies to metal-organic framework materials are in good agreement with experimental observations. It is suggested that replacement of the organic linker in MOF materials with some of the more efficient aromatics investigated here might contribute to enhance the H₂ storage properties of mixed inorganic–organic materials.     

Coevolution of protein and RNA structures within a highly conserved ribosomal domain

The X-ray crystal structure of a ribosomal L11-rRNA complex with chloroplast-like mutations in both protein and rRNA is presented. The global structure is almost identical to that of the wild-type (bacterial) complex, with only a small movement of the protein alpha helix away from the surface of the RNA required to accommodate the altered protein residue. In contrast, the specific hydrogen bonding pattern of the mutated residues is substantially different, and now includes a direct interaction between the protein side chain and an RNA base edge and a water-mediated contact. Comparison of the two structures allows the observations of sequence variation and relative affinities of wild-type and mutant complexes to be clearly rationalized, but reinforces the concept that there is no single simple code for protein-RNA recognition.     

Do collective atomic fluctuations account for cooperative effects? Molecular dynamics studies of the U1A-RNA complex

A complete understanding of gene expression relies on a comprehensive understanding of the protein-RNA recognition process. However, the study of protein-RNA recognition is complicated by many factors that contribute to both binding affinity and specificity, including structure, energetics, dynamical motions, and cooperative interactions. Several recent studies have suggested that energetic coupling between residues contributes to formation of the complex between the U1A protein and stem loop 2 of U1 snRNA as a consequence of a cooperative network of interactions. We have performed molecular dynamics simulations on the U1A-RNA complex, including explicit water and counterions, and analyzed the results based on the calculated positional cross-correlations of atomic fluctuations. The results indicate that cross-correlations calculated on a per residue basis agree well with the observed inter-residue cooperativity and predict that the networks identified to date may also be coupled into an extensive hyper-network that reflects the intrinsic rigidity of the RNA recognition motif. In addition, we report a comparison of the MD calculated correlations with the results of a positional covariance analysis based on the sequences of 330 RNA recognition motifs, including U1A. The calculated inter-residue cross-correlations agree very well with the results of the sites exhibiting positional covariance. Collectively, these results strongly support the hypothesis that collective fluctuations contribute to cooperativity and the corresponding observed thermodynamic coupling. Predictions of additional sites in U1A that may be involved in cooperative networks are advanced.     

RNA architecture dictates the conformations of a bound peptide.

BACKGROUND: The biological function of several viral and bacteriophage proteins, and their arginine-rich subdomains, involves RNA-mediated interactions. It has been shown recently that bound peptides adopt either beta-hairpin or alpha-helical conformations in viral and phage peptide-RNA complexes. We have compared the structures of the arginine-rich peptide domain of HIV-1 Rev bound to two RNA aptamers to determine whether RNA architecture can dictate the conformations of a bound peptide. RESULTS: The core-binding segment of the HIV-1 Rev peptide class II RNA aptamer complex spans the two-base bulge and hairpin loop of the bound RNA and the carboxy-terminal segment of the bound peptide. The bound peptide is anchored in place by backbone and sidechain intermolecular hydrogen bonding and van der Waals stacking interactions. One of the bulge bases participates in U*(A*U) base triple formation, whereas the other is looped out and flaps over the bound peptide in the complex. The seven-residue hairpin loop is closed by a sheared G*A mismatch pair with several pyrimidines looped out of the hairpin fold. CONCLUSIONS: Our structural studies establish that RNA architecture dictates whether the same HIV-1 Rev peptide folds into an extended or alpha-helical conformation on complex formation. Arginine-rich peptides can therefore adapt distinct secondary folds to complement the tertiary folds of their RNA targets. This contrasts with protein-RNA complexes in which elements of RNA secondary structure adapt to fit within the tertiary folds of their protein targets.