A proposed molecular model for the interaction of calcineurin with the cyclosporin A-cyclophilin A complex.

Cyclosporin A (CsA) and FK506 are potent natural product immunosuppressants that induce their biological effects by forming an initial complex with cytosolic proteins termed immunophilins. These drug immunophilin complexes then bind to and inhibit the serine/threonine protein phosphatase calcineurin (CN). Two classes of immunophilin have been identified with cyclophilins (CyP's) being proteins specifically binding CsA and FKBPs specifically binding FK506. Solution and crystal structures of various CsA-CyP and FK506-FKBP complexes have been determined and show no apparent structural similarity between the two classes of drug protein complexes. These findings raise the question as to how, given their structural differences, these two complexes can both inhibit CN. While the crystal structure of the FK506-FKBP12-CN complex has been reported, no structure for a CsA-CyP CN complex has been determined. Here are reported studies that use various modelling strategies to construct a model for the interaction of the cyclosporin A- cyclophilin A complex with calcineurin. The first stage of constructing this model consisted of using conformational comparison of CsA and FK506, GRID and GROUP analysis and restrained molecular dynamics to dock CsA into the FK506 binding site of the FK506-FKBP12-CN structure. An initial model for the CsA-CyPA-CN complex was then constructed by superimposing the structure of the CsA-CyPA complex onto the docked CsA molecule. This model was then optimised with molecular dynamics simulations run on sterically clashing regions. The validity of the model for the CsA-CyPA-CN complex was then examined with respect to the effect of chemical modifications to CsA and amino acid substitutions within CyPA on the ability of the drug-immunophilin complex to inhibit calcineurin.     

Kinetic studies of RNA-protein interactions using surface plasmon resonance

Although structural, biochemical, and genetic studies have provided much insight into the determinants of specificity and affinity of proteins for RNA, little is currently known about the kinetics that underlie RNA-protein interactions. Protein-RNA complexes are dynamic, and the kinetics of binding and release could influence many processes, such as the ability of RNA-binding proteins to compete for binding sites, the sequential assembly of ribonucleoprotein complexes, and the ability of bound RNA to move between cellular compartments. Therefore, to attain a complete and biologically relevant understanding of RNA-protein interactions, complex formation must be studied not only in equilibrated reactions, but also as a dynamic process. BIACORE, a surface plasmon resonance-based biosensor technology, allows intermolecular interactions to be measured in real time, and can provide both equilibrium and kinetic information about complex formation. This technology is a powerful tool with which to study the dynamics of RNA-protein interactions. We have used BIACORE extensively to obtain detailed insight into the interaction between RNA and proteins carrying RNA recognition motif domains. Here we discuss the physical principles on which BIACORE is based, and the required instrumentation. We describe how to design well-controlled RNA-protein interaction experiments aimed at yielding high-quality data, and outline the steps required for data analysis. In addition, we present examples to illustrate how kinetic studies have provided us with unique insights into the interaction of the spliceosomal U1A protein and the neuronal HuD protein with their respective RNA targets.     

Structure-based prediction of DNA target sites by regulatory proteins

Regulatory proteins play a critical role in controlling complex spatial and temporal patterns of gene expression in higher organism, by recognizing multiple DNA sequences and regulating multiple target genes. Increasing amounts of structural data on the protein-DNA complex provides clues for the mechanism of target recognition by regulatory proteins. The analyses of the propensities of base-amino acid interactions observed in those structural data show that there is no one-to-one correspondence in the interaction, but clear preferences exist. On the other hand, the analysis of spatial distribution of amino acids around bases shows that even those amino acids with strong base preference such as Arg with G are distributed in a wide space around bases. Thus, amino acids with many different geometries can form a similar type of interaction with bases. The redundancy and structural flexibility in the interaction suggest that there are no simple rules in the sequence recognition, and its prediction is not straightforward. However, the spatial distributions of amino acids around bases indicate a possibility that the structural data can be used to derive empirical interaction potentials between amino acids and bases. Such information extracted from structural databases has been successfully used to predict amino acid sequences that fold into particular protein structures. We surmised that the structures of protein-DNA complexes could be used to predict DNA target sites for regulatory proteins, because determining DNA sequences that bind to a particular protein structure should be similar to finding amino acid sequences that fold into a particular structure. Here we demonstrate that the structural data can be used to predict DNA target sequences for regulatory proteins. Pairwise potentials that determine the interaction between bases and amino acids were empirically derived from the structural data. These potentials were then used to examine the compatibility between DNA sequences and the protein-DNA complex structure in a combinatorial "threading" procedure. We applied this strategy to the structures of protein-DNA complexes to predict DNA binding sites recognized by regulatory proteins. To test the applicability of this method in target-site prediction, we examined the effects of cognate and noncognate binding, cooperative binding, and DNA deformation on the binding specificity, and predicted binding sites in real promoters and compared with experimental data. These results show that target binding sites for several regulatory proteins are successfully predicted, and our data suggest that this method can serve as a powerful tool for predicting multiple target sites and target genes for regulatory proteins.     

Continuous association of Escherichia coli single-stranded DNA binding protein with stable complexes of recA protein and single-stranded DNA

The single-stranded DNA binding protein of Escherichia coli (SSB) stimulates recA protein promoted DNA strand exchange reactions by promoting and stabilizing the interaction between recA protein and single-stranded DNA (ssDNA). Utilizing the intrinsic tryptophan fluorescence of SSB, an ATP-dependent interaction has been detected between SSB and recA-ssDNA complexes. This interaction is continuous for periods exceeding 1 h under conditions that are optimal for DNA strand exchange. Our data suggest that this interaction does not involve significant displacement of recA protein in the complex by SSB when ATP is present. The properties of this interaction are consistent with the properties of SSB-stabilized recA-ssDNA complexes determined by other methods. The data are incompatible with models in which SSB is displaced after functioning transiently in the formation of recA-ssDNA complexes. A continuous association of SSB with recA-ssDNA complexes may therefore be an important feature of the mechanism by which SSB stimulates recA protein promoted reactions.     

A Bir1p Sli15p kinetochore passenger complex regulates septin organization during anaphase

Kinetochore-passenger complexes in metazoans have been proposed to coordinate the segregation of chromosomes in anaphase with the induction of cytokinesis. Passenger protein homologues in the budding yeast Saccharomyces cerevisiae play a critical role early in mitosis, ensuring proper biorientation of kinetochore-microtubule attachments. Our recent work has implicated the passenger protein Bir1p (Survivin) and the inner kinetochore complex centromere binding factor 3 (CBF3) in the regulation of septin dynamics during anaphase. Here, we present data that is consistent with there being multiple passenger protein complexes. Our data show that Bir1p links together a large passenger complex containing Ndc10p, Sli15p (INCENP), and Ipl1p (Aurora B) and that the interaction between Bir1p and Sli15p is specifically involved in regulating septin dynamics during anaphase. Neither conditional alleles nor mutants of BIR1 that disrupt the interaction between Bir1p and Sli15p resulted in mono-attached kinetochores, suggesting that the Bir1p-Sli15p complex functions in anaphase and independently from Sli15p-Ipl1p complexes. We present a model for how discrete passenger complexes coordinate distinct aspects of mitosis.     

Differential Dynamic Engagement within 24 SH3 Domain: Peptide Complexes Revealed by Co-Linear Chemical Shift Perturbation Analysis

There is increasing evidence for the functional importance of multiple dynamically populated states within single proteins. However, peptide binding by protein-protein interaction domains, such as the SH3 domain, has generally been considered to involve the full engagement of peptide to the binding surface with minimal dynamics and simple methods to determine dynamics at the binding surface for multiple related complexes have not been described. We have used NMR spectroscopy combined with isothermal titration calorimetry to comprehensively examine the extent of engagement to the yeast Abp1p SH3 domain for 24 different peptides. Over one quarter of the domain residues display co-linear chemical shift perturbation (CCSP) behavior, in which the position of a given chemical shift in a complex is co-linear with the same chemical shift in the other complexes, providing evidence that each complex exists as a unique dynamic rapidly inter-converting ensemble. The extent the specificity determining sub-surface of AbpSH3 is engaged as judged by CCSP analysis correlates with structural and thermodynamic measurements as well as with functional data, revealing the basis for significant structural and functional diversity amongst the related complexes. Thus, CCSP analysis can distinguish peptide complexes that may appear identical in terms of general structure and percent peptide occupancy but have significant local binding differences across the interface, affecting their ability to transmit conformational change across the domain and resulting in functional differences.     

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.     

Kinetic analysis of interaction of BRCA1 tandem breast cancer c-terminal domains with phosphorylated peptides reveals two binding conformations

Tandem breast cancer C-terminal (BRCT) domains, present in many DNA repair and cell cycle checkpoint signaling proteins, are phosphoprotein binding modules. The best-characterized tandem BRCT domains to date are from the protein BRCA1 (BRCA1-BRCT), an E3 ubiquitin ligase that has been linked to breast and ovarian cancer. While X-ray crystallography and NMR spectroscopy studies have uncovered the structural determinants of specificity of BRCA1-BRCT for phosphorylated peptides, a detailed kinetic and thermodynamic characterization of the interaction is also required to understand how structure and dynamics are connected and therefore better probe the mechanism of phosphopeptide recognition by BRCT domains. Through a global analysis of binding kinetics data obtained from surface plasmon resonance (SPR) and stopped-flow fluorescence spectroscopy, we show that the recognition mechanism is complex and best modeled by two equilibrium conformations of BRCA1-BRCT in the free state that both interact with a phosphopeptide, with dissociation constants ( K d) in the micromolar range. We show that the apparent global dissociation constant derived from this kinetic analysis is similar to the K d values measured using steady-state SPR, isothermal titration calorimetry, and fluorescence anisotropy. The dynamic nature of BRCA1-BRCT may facilitate the binding of BRCA1 to different phosphorylated protein targets.     

Interaction between bacteriophage PM2 protein IV and DNA.

The interaction between DNA and the structural protein IV of bacteriophage PM2 was studied by co-sedimentation, filter binding and electron microscopy. The co-sedimentation data and the sigmoid-shaped filter binding curve were interpreted in terms of co-operative binding. At a given DNA/protein input ratio, some DNA molecules were associated with a large amount of protein IV while others had no detectable protein bound to them. Electron microscopic examination of DNA-protein IV mixtures showed highly condensed DNA molecules alongside uncomplexed native DNA. Dissociation experiments revealed the presence of two types of complexes. Type I dissociated rapidly while type II had a long half-life. Dissociation of complexes obtained with increasing protein/DNA ratios suggested that the type I complex was a precursor of type II complex. Protein IV binds equally well to superhelical, relaxed or linear DNA as well as to single-stranded DNA. These observations lead to a model for the interaction and for the consequent alterations in the DNA structure.     

DNA strand arrangement within the SfiI-DNA complex: atomic force microscopy analysis

The SfiI restriction enzyme binds to DNA as a tetramer holding two usually distant DNA recognition sites together before cleavage of the four DNA strands. To elucidate structural properties of the SfiI-DNA complex, atomic force microscopy (AFM) imaging of the complexes under noncleaving conditions (Ca2+ instead of Mg2+ in the reaction buffer) was performed. Intramolecular complexes formed by protein interaction between two binding sites in one DNA molecule (cis interaction) as well as complexes formed by the interaction of two sites in different molecules (trans interaction) were analyzed. Complexes were identified unambiguously by the presence of a tall spherical blob at the DNA intersections. To characterize the path of DNA within the complex, the angles between the DNA helices in the proximity of the complex were systematically analyzed. All the data show clear-cut bimodal distributions centered around peak values corresponding to 60 degrees and 120 degrees. To unambiguously distinguish between the crossed and bent models for the DNA orientation within the complex, DNA molecules with different arm lengths flanking the SfiI binding site were designed. The analysis of the AFM images for complexes of this type led to the conclusion that the DNA recognition sites within the complex are crossed. The angles of 60 degrees or 120 degrees between the DNA helices correspond to a complex in which one of the helices is flipped with respect to the orientation of the other. Complexes formed by five different recognition sequences (5'-GGCCNNNNNGGCC-3'), with different central base pairs, were also analyzed. Our results showed that complexes containing the two possible orientations of the helices were formed almost equally. This suggests no preferential orientation of the DNA cognate site within the complex, suggesting that the central part of the DNA binding site does not form strong sequence specific contacts with the protein.     

A pseudo-particle approach for studying protein-ligand models truncated to their active sites

A molecular dynamics method has been developed to describe the structural and dynamic properties of protein-ligand complexes that are truncated to their active sites. The active site is comprised of the ligand and discontinuous, positionally unrestrained peptide chains. This truncated active-site complex is surrounded by big unspecific pseudo-particles representing the complete protein and the solvent. Thus, knowledge of the folding of the outer parts of the protein is not required, and the method can be applied to protein models, derived from homology modeling. The method has been tested using ligand complexes of adenylate kinase, retinol binding protein, HIV-1 protease, and human leucocyte antigen. Comparisons with their crystal structures and with results from time-demanding simulations of the whole complexes in explicit water solvent show that the ligand binding properties are conserved. Most of the hydrogen bonds between the ligand and the active-site residues are reproduced and, furthermore, the simulation time is reduced. © 1996 John Wiley & Sons, Inc.     

Characterization of molecular recognition features, MoRFs, and their binding partners

Molecular Recognition Features (MoRFs) are short, interaction-prone segments of protein disorder that undergo disorder-to-order transitions upon specific binding, representing a specific class of intrinsically disordered regions that exhibit molecular recognition and binding functions. MoRFs are common in various proteomes and occupy a unique structural and functional niche in which function is a direct consequence of intrinsic disorder. Example MoRFs collected from the Protein Data Bank (PDB) have been divided into three subtypes according to their structures in the bound state: alpha-MoRFs form alpha-helices, beta-MoRFs form beta-strands, and iota-MoRFs form structures without a regular pattern of backbone hydrogen bonds. These example MoRFs were indicated to be intrinsically disordered in the absence of their binding partners by several criteria. In this study, we used several geometric and physiochemical criteria to examine the properties of 62 alpha-, 20 beta-, and 176 iota-MoRF complex structures. Interface residues were examined by calculating differences in accessible surface area between the complex and isolated monomers. The compositions and physiochemical properties of MoRF and MoRF partner interface residues were compared to the interface residues of homodimers, heterodimers, and antigen-antibody complexes. Our analysis indicates that there are significant differences in residue composition and several geometric and physicochemical properties that can be used to discriminate, with a high degree of accuracy, between various interfaces in protein interaction data sets. Implications of these findings for the development of MoRF-partner interaction predictors are discussed. In addition, structural changes upon MoRF-to-partner complex formation were examined for several illustrative examples.     

Allosteric communication in cysteinyl tRNA synthetase: a network of direct and indirect readout

Protein structure networks are constructed for the identification of long-range signaling pathways in cysteinyl tRNA synthetase (CysRS). Molecular dynamics simulation trajectory of CysRS-ligand complexes were used to determine conformational ensembles in order to gain insight into the allosteric signaling paths. Communication paths between the anticodon binding region and the aminoacylation region have been identified. Extensive interaction between the helix bundle domain and the anticodon binding domain, resulting in structural rigidity in the presence of tRNA, has been detected. Based on the predicted model, six residues along the communication paths have been examined by mutations (single and double) and shown to mediate a coordinated coupling between anticodon recognition and activation of amino acid at the active site. This study on CysRS clearly shows that specific key residues, which are involved in communication between distal sites in allosteric proteins but may be elusive in direct structure analysis, can be identified from dynamics of protein structure networks.     

Anticancer activity of new imidazole derivative of 1R,2R-diaminocyclohexane palladium and platinum complexes as DNA fluorescent probes

The aim of this study was synthesis of two new water-soluble fluorescent palladium and platinum complexes with formulas of [Pt(DACH)(FIP)](NO₃)₂ and [Pd(DACH)(FIP)](NO₃)₂, respectively, where FIP is 2-(furan-2-yl)-1H-imidazo[4,5-f][1,10] phenanthroline and DACH is 1R,2R-diaminocyclohexane. Fluorescence spectroscopy, circular dichroism (CD), thermal denaturation measurement, ionic strength, and kinetic study displayed groove binding of Pt complex on DNA, while due to binding of Pd complex, B form of DNA convert to Z form. Due to electrostatic interaction of Pd complex with DNA, the DNA form is converted and it provides enough space for Pd complex to insert between base stacking of DNA. UV–vis study shows two complexes could denature the DNA at low concentrations in exothermic process and Pt complex is more active than Pd complex. Finally, the anticancer and growth inhibitory activities of synthesized complexes were investigated against human colon cancer cell line HCT116 after incubation time of 24 h using MTT assay and higher activity was observed for the platinum complex. Interaction of the two metal derivative complexes was studied by molecular docking and molecular dynamics simulation. The results showed that Pt complexes have higher negative docking energy and higher tendency for interaction with DNA, and exert more structural change on DNA.     

Simulation of the interaction between ScyTx and small conductance calcium-activated potassium channel by docking and MM-PBSA

Computational methods are employed to simulate interaction of scorpion toxin ScyTx in complex with the small conductance calcium-activated potassium channel rsk2. All of available 25 structures of ScyTx in the Protein Data Bank determined by NMR were considered for improving performance of rigid protein docking of ZDOCK. Four main binding modes were found among a large number of predicted complexes by using clustering analysis, screening with expert knowledge, energy minimization, and molecular dynamics simulations. The quality and validity of the resulting complexes were further evaluated by molecular dynamics simulations with the generalized Born solvation model and by calculation of relative binding free energies with the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) in the AMBER 7 suit of programs. The complex formed by the 22nd structure of the ScyTx and rsk2 channel was identified as the most favorable complex by using a combination of computational methods, which contain further introduction of flexibility without restraining residue side chain. From the resulted spatial structure of the ScyTx and rsk2 channel, ScyTx associates the mouth of the rsk2 channel with alpha-helix rather than beta-sheet. Structural analysis first revealed that Arg(13) played a novel and vital role of blocking the pore of the rsk2 channel, whose role is remarkably different from that of highly homologous scorpion toxin P05. Between the interfaces in the ScyTx-rsk2 complex, strong electrostatic interaction and hydrogen bonds exist between Arg(13) of ScyTx and Gly-Tyr-Gly-Asp sequential residues located in the four symmetrical chains of the pore region. Simultaneously, five hydrogen bonds between Arg(6) of ScyTx and Asp(341)(C), Val(366)(C), and Pro(367)(C), and electrostatic interaction between Arg(6) of ScyTx and Asp(364)(B) and Asp(341)(C) are also found by structural analysis. In addition, His(31) located at the C-terminal of ScyTx is surrounded by Val(342)(A), Asp(364)(A), Met(365)(A), Pro(367)(B), and Asn(366)(B) within a contact distance of 4.0 A. These simulation results are in good agreement with experimental data and can effectively explain the binding phenomena between ScyTx and the potassium channel at the level of molecular spatial structure. The consistency between results of molecular modeling and experimental data strongly suggests that our spatial structure model of the ScyTx-rsk2 complex is reasonable. Therefore, molecular docking combined with molecular dynamics simulations followed by molecular mechanics Poisson-Boltzmann surface area analysis is an attractive approach for modeling scorpion toxin-potassium channel complexes a priori for further biological studies.     

The structural basis for partitioning of the XRCC1/DNA ligase III-α BRCT-mediated dimer complexes

The ultimate step common to almost all DNA repair pathways is the ligation of the nicked intermediate to form contiguous double-stranded DNA. In the mammalian nucleotide and base excision repair pathways, the ligation step is carried out by ligase III-α. For efficient ligation, ligase III-α is constitutively bound to the scaffolding protein XRCC1 through interactions between the C-terminal BRCT domains of each protein. Although structural data for the individual domains has been available, no structure of the complex has been determined and several alternative proposals for this interaction have been advanced. Interpretation of the models is complicated by the formation of homodimers that, depending on the model, may either contribute to, or compete with heterodimer formation. We report here the structures of both homodimer complexes as well as the heterodimer complex. Structural characterization of the heterodimer formed from a longer XRCC1 BRCT domain construct, including residues comprising the interdomain linker region, revealed an expanded heterodimer interface with the ligase III-α BRCT domain. This enhanced linker-mediated binding interface plays a significant role in the determination of heterodimer/homodimer selectivity. These data provide fundamental insights into the structural basis of BRCT-mediated dimerization, and resolve questions related to the organization of this important repair complex.     

Toward the understanding of the structure and dynamics of protein-carbohydrate interactions: molecular dynamics studies of the complexes between hevein and oligosaccharidic ligands

Herein we study, through all atom molecular dynamics simulations, the complex between hevein and two N-acetylated chitin oligomers, namely N,N(')-diacetylchitobiose and N,N('),N(")-triacetylchitotriose. The results of the simulations for two disaccharide complexes and one trisaccharide complex show that a carbohydrate oligomer is able to move on the surface of the relatively flat binding pocket of hevein, therefore occupying different binding subpockets. Statistical analysis methods were also applied in order to define the principal overall motions in the complexes, showing how the different ligands in the simulations modulate the protein motions. The oligosaccharide binding can be considered as defined by a subtle balance between enthalpic (formation of intermolecular interactions between the ligand and the receptor) and entropic (due mainly to the possibility for the sugar to move on the surface of the protein domain) effects, determining multiple binding conformations. This structural and dynamical view could parallel the results obtained by regularly used restrained MD simulations based on NOE NMR data that provide a well defined structure for both the disaccharide and trisaccharide complexes, and agrees with the observations for longer oligosaccharide chains.