Protein-protein docking with backbone flexibility

Computational protein-protein docking methods currently can create models with atomic accuracy for protein complexes provided that the conformational changes upon association are restricted to the side chains. However, it remains very challenging to account for backbone conformational changes during docking, and most current methods inherently keep monomer backbones rigid for algorithmic simplicity and computational efficiency. Here we present a reformulation of the Rosetta docking method that incorporates explicit backbone flexibility in protein-protein docking. The new method is based on a "fold-tree" representation of the molecular system, which seamlessly integrates internal torsional degrees of freedom and rigid-body degrees of freedom. Problems with internal flexible regions ranging from one or more loops or hinge regions to all of one or both partners can be readily treated using appropriately constructed fold trees. The explicit treatment of backbone flexibility improves both sampling in the vicinity of the native docked conformation and the energetic discrimination between near-native and incorrect models.     

Biologically enhanced sampling geometric docking and backbone flexibility treatment with multiconformational superposition

An efficient biologically enhanced sampling geometric docking method is presented based on the FTDock algorithm to predict the protein-protein binding modes. The active site data from different sources, such as biochemical and biophysical experiments or theoretical analyses of sequence data, can be incorporated in the rotation-translation scan. When discretizing a protein onto a 3-dimensional (3D) grid, a zero value is given to grid points outside a sphere centered on the geometric center of specified residues. In this way, docking solutions are biased toward modes where the interface region is inside the sphere. We also adopt a multiconformational superposition scheme to represent backbone flexibility in the proteins. When these procedures were applied to the targets of CAPRI, a larger number of hits and smaller ligand root-mean-square deviations (RMSDs) were obtained at the conformational search stage in all cases, and especially Target 19. With Target 18, only 1 near-native structure was retained by the biologically enhanced sampling geometric docking method, but this number increased to 53 and the least ligand RMSD decreased from 8.1 A to 2.9 A after performing multiconformational superposition. These results were obtained after the CAPRI prediction deadlines.     

Side-chain and backbone flexibility in protein core design.

We have developed a computational approach for the design and prediction of hydrophobic cores that includes explicit backbone flexibility. The program consists of a two-stage combination of a genetic algorithm and monte carlo sampling using a torsional model of the protein. Backbone structures are evaluated either by a canonical force-field or a constraining potential that emphasizes the preservation of local geometry. The utility of the method for protein design and engineering is explored by designing three novel hydrophobic core variants of the protein 434 cro. We use the new method to evaluate these and previously designed 434 cro variants, as well as a series of phage T4 lysozyme variants. In order to properly evaluate the influence of backbone flexibility, we have also analyzed the effects of varying amounts of side-chain flexibility on the performance of fixed backbone methods. Comparison of results using a fixed versus flexible backbone reveals that, surprisingly, the two methods are almost equivalent in their abilities to predict relative experimental stabilities, but only when full side-chain flexibility is allowed. The prediction of core side-chain structure can vary dramatically between methods. In some, but not all, cases the flexible backbone method is a better predictor of structure. The development of a flexible backbone approach to core design is particularly important for attempts at de novo protein design, where there is no prior knowledge of a precise backbone structure.     

Increased backbone flexibility in threonine45-phosphorylated hirudin upon pH change.

Protein phosphorylation on serine/threonine side chains represents a major regulatory event in the posttranslational control of protein functionality, where it is thought to operate at the level of structural changes in the polypeptide chain. However, key questions about molecular aspects of phosphate ester induced conformational alterations remain open. Among these concerns are the radius of action of the phosphate ester group, its effective ionic state, and its interplay with distinct bonds of the polypeptide chain. Primarily to define short-range effects upon threonine phosphorylation, the native 65 amino acid protein hirudin, conformationally restrained by a proline flanking the pThr(45) site and three intramolecular disulfide bonds, was structurally characterized in both the phosphorylated and the unphosphorylated state in solution. Circular dichroism and hydrogen exchange experiments (MALDI-TOF) showed that structural changes were caused by Thr(45)-Pro(46) phosphorylation only when the phosphate ester group was in its dianionic state. The spatial arrangement of the amino acids, monitored by 1H NMR spectroscopy, appears to be affected within a radius of about 10 A around the pThr(45)-OgammaH, with phosphorylation resulting in a loss of structure and increased flexibility within a segment of at least seven amino acid residues. Thus, the transition from the monoanionic to the dianionic phosphate group over the pH range 5.2-8.5 represents a general phosphorylation-dependent conformational switch operating at physiological pH values.     

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.     

Repeat motions and backbone flexibility in designed proteins with different numbers of identical consensus tetratricopeptide repeats

The tetratricopeptide repeat (TPR) is a 34-residue helix-turn-helix motif that occurs as three or more tandem repeats in a wide variety of proteins. We have determined the repeat motions and backbone fluctuations of proteins containing two or three consensus TPR repeats (CTPR2 and CPTR3, respectively) using 15N NMR relaxation measurements. Rotational diffusion tensors calculated from these data for each repeat within each TPR protein indicate that there is a high degree of motional correlation between different repeats in the same protein. This is consistent with the prevailing view that repeat proteins, such as CTPR2 and CTPR3, behave as single cooperatively folded domains. The internal motions of backbone NH groups were determined using the Lipari-Szabo model-free formalism. For most residues, there was a clear separation between the influence of internal motion and the influence of global rotational tumbling on the observed magnetic relaxation. The local internal motions are highly restricted in most of the helical elements, with slightly greater flexibility in the linker elements. Comparisons between CTPR2 and CTPR3 indicate that an addition of a TPR repeat to the C-terminus (before the solvation helix) of CTPR2 slightly reduces the flexibility of the preceding helix.     

Incorporating biochemical information and backbone flexibility in RosettaDock for CAPRI rounds 6-12

In CAPRI rounds 6-12, RosettaDock successfully predicted 2 of 5 unbound-unbound targets to medium accuracy. Improvement over the previous method was achieved with computational mutagenesis to select decoys that match the energetics of experimentally determined hot spots. In the case of Target 21, Orc1/Sir1, this resulted in a successful docking prediction where RosettaDock alone or with simple site constraints failed. Experimental information also helped limit the interacting region of TolB/Pal, producing a successful prediction of Target 26. In addition, we docked multiple loop conformations for Target 20, and we developed a novel flexible docking algorithm to simultaneously optimize backbone conformation and rigid-body orientation to generate a wide diversity of conformations for Target 24. Continued challenges included docking of homology targets that differ substantially from their template (sequence identity <50%) and accounting for large conformational changes upon binding. Despite a larger number of unbound-unbound and homology model binding targets, Rounds 6-12 reinforced that RosettaDock is a powerful algorithm for predicting bound complex structures, especially when combined with experimental data.     

The role of proline and glycine in determining the backbone flexibility of a channel-forming peptide

Alamethicin is a helical 20-amino acid voltage-gated channel-forming peptide, which is known to exhibit segmental flexibility in solution along its backbone near alpha-methylalanine (MeA)-10 and Gly-11. In an alpha-helical configuration, MeA at position 10 would normally hydrogen-bond with position 14, but the presence of proline at this position prevents the formation of this interhelical hydrogen bond. To determine whether the presence of proline at position 14 contributes to the flexibility of this helix, two analogs of alamethicin were synthesized, one with proline 14 replaced by alanine and another with both proline 14 and glycine 11 replaced by alanine. The C-termini of these peptides were derivatized with a proxyl nitroxide, and paramagnetic enhancements produced by the nitroxide on the Calpha protons were used to estimate r-6 weighted distances between the nitroxide and the backbone protons. When compared to native alamethicin, the analog lacking proline 14 exhibited similar C-terminal to Calpha proton distances, indicating that substitution of proline alone does not alter the flexibility of this helix; however, the subsequent removal of glycine 11 resulted in a significant increase in the averaged distances between the C- and N-termini. Thus, the G-X-X-P motif found in alamethicin appears to be largely responsible for mediating high-amplitude bending motions that have been observed in the central helical domain of alamethicin in methanol. To determine whether these substitutions alter the channel behavior of alamethicin, the macroscopic and single-channel currents produced by these analogs were compared. Although the substitution of the G-X-X-P motif produces channels with altered characteristics, this motif is not essential to achieve voltage-dependent gating or alamethicin-like behavior.     

Modeling backbone flexibility to achieve sequence diversity: the design of novel alpha-helical ligands for Bcl-xL

Computational protein design can be used to select sequences that are compatible with a fixed-backbone template. This strategy has been used in numerous instances to engineer novel proteins. However, the fixed-backbone assumption severely restricts the sequence space that is accessible via design. For challenging problems, such as the design of functional proteins, this may not be acceptable. Here, we present a method for introducing backbone flexibility into protein design calculations and apply it to the design of diverse helical BH3 ligands that bind to the anti-apoptotic protein Bcl-xL, a member of the Bcl-2 protein family. We demonstrate how normal mode analysis can be used to sample different BH3 backbones, and show that this leads to a larger and more diverse set of low-energy solutions than can be achieved using a native high-resolution Bcl-xL complex crystal structure as a template. We tested several of the designed solutions experimentally and found that this approach worked well when normal mode calculations were used to deform a native BH3 helix structure, but less well when they were used to deform an idealized helix. A subsequent round of design and testing identified a likely source of the problem as inadequate sampling of the helix pitch. In all, we tested 17 designed BH3 peptide sequences, including several point mutants. Of these, eight bound well to Bcl-xL and four others showed weak but detectable binding. The successful designs showed a diversity of sequences that would have been difficult or impossible to achieve using only a fixed backbone. Thus, introducing backbone flexibility via normal mode analysis effectively broadened the set of sequences identified by computational design, and provided insight into positions important for binding Bcl-xL.     

Interdependence of backbone flexibility,residue conservation,and enzyme function: a case study on beta1,4-galactosyltransferase-I

Beta1,4-galactosyltransferase-I (beta4Gal-T1) catalyzes the transfer of a galactose from UDP-galactose to N-acetylglucosamine. A recent crystal structure determination of the substrate-bound enzyme reveals a large conformational change, which creates binding sites for the oligosaccharide and alpha-lactalbumin, when compared to the ligand-free structure. The conformational changes take place in a 21-residue-long loop (I345-H365) and in a smaller loop containing a tryptophan residue (W314) flanked by glycines (Y311-G316; Trp loop). A series of molecular dynamics simulations carried out with an implicit solvent model and with explicit water successfully identify flexibility in the two loops and in another interacting loop. These observations are confirmed by limited proteolysis experiments that reveal an intrinsic flexibility of the long loop. The multiple simulation runs starting with the substrate-free structure show that the long loop moves toward its conformation in the ligand-bound structure; however, it gets stabilized in an intermediate position. The Trp loop moves in the opposite direction to that of the long loop, making contacts with residues in the long loop. Remarkably, when the Trp loop is restrained in its starting conformation, no large conformational change takes place in the long loop, indicating residue communication of flexibility. Sequence and structural analysis of the beta4Gal-T1 family with 37 known sequences reveals that in contrast to the unconserved long loop, which undergoes a much larger conformational change, the Trp loop including the glycines is highly conserved. These observations lead us to propose a new functional mechanism that may be conserved by evolution to perform a variety of functions.     

Pyrrolidine PNA-DNA chimeric oligonucleotides with extended backbone

Cis-D-2-hydroxy-4-thymin-1-yl-pyrrolidine propionic acid unit is used to make PNA-DNA dimer block that is incorporated in DNA sequences at selected positions. Since the amide linkage is shorter than phosphodiester linkage, insertion of an extra atom in the backbone with amide linkage seems to be better accommodated for internucleotide distance-complementarity.     

New synthetic catecholate-type siderophores with triamine backbone

New analogues of triscatecholate siderophores based on linear or tripodal triamines with or without spacer groups or lipophilic and hydrophilic substituents were synthesized. The catecholate moieties were prepared in OH-forms, as acetylated compounds or masked as 8-methoxycarbonyloxy-2,4-dioxo-1,3-benzoxazine derivatives. Some of the new compounds were active as siderophores tested by growth promotion assays using various Gram-negative bacteria and mycobacteria under iron limitation and by CAS-assay. Structure-activity-correlations have been studied.     

NMR studies of the backbone flexibility and structure of human growth hormone: a comparison of high and low pH conformations

(15)N NMR relaxation parameters and amide (1)H/(2)H-exchange rates have been used to characterize the structural flexibility of human growth hormone (rhGH) at neutral and acidic pH. Our results show that the rigidity of the molecule is strongly affected by the solution conditions. At pH 7.0 the backbone dynamics parameters of rhGH are uniform along the polypeptide chain and their values are similar to those of other folded proteins. In contrast, at pH 2.7 the overall backbone flexibility increases substantially compared to neutral pH and the average order parameter approaches the lower limit expected for a folded protein. However, a significant variation of the backbone dynamics through the molecule indicates that under acidic conditions the mobility of the residues becomes more dependent on their location within the secondary structure units. In particular, the order parameters of certain loop regions decrease dramatically and become comparable to those found in unfolded proteins. Furthermore, the HN-exchange rates at low pH reveal that the residues most protected from exchange are clustered at one end of the helical bundle, forming a stable nucleus. We suggest that this nucleus maintains the overall fold of the protein under destabilizing conditions. We therefore conclude that the acid state of rhGH consists of a structurally conserved, but dynamically more flexible helical core surrounded by an aura of highly mobile, unstructured loops. However, in spite of its prominent flexibility the acid state of rhGH cannot be considered a "molten globule" state because of its high stability. It appears from our work that under certain conditions, a protein can tolerate a considerable increase in flexibility of its backbone, along with an increased penetration of water into its core, while still maintaining a stable folded conformation.     

Statistical theory for protein combinatorial libraries. Packing interactions, backbone flexibility, and the sequence variability of a main-chain structure

Combinatorial experiments provide new ways to probe the determinants of protein folding and to identify novel folding amino acid sequences. These types of experiments, however, are complicated both by enormous conformational complexity and by large numbers of possible sequences. Therefore, a quantitative computational theory would be helpful in designing and interpreting these types of experiment. Here, we present and apply a statistically based, computational approach for identifying the properties of sequences compatible with a given main-chain structure. Protein side-chain conformations are included in an atom-based fashion. Calculations are performed for a variety of similar backbone structures to identify sequence properties that are robust with respect to minor changes in main-chain structure. Rather than specific sequences, the method yields the likelihood of each of the amino acids at preselected positions in a given protein structure. The theory may be used to quantify the characteristics of sequence space for a chosen structure without explicitly tabulating sequences. To account for hydrophobic effects, we introduce an environmental energy that it is consistent with other simple hydrophobicity scales and show that it is effective for side-chain modeling. We apply the method to calculate the identity probabilities of selected positions of the immunoglobulin light chain-binding domain of protein L, for which many variant folding sequences are available. The calculations compare favorably with the experimentally observed identity probabilities.     

31P NMR spectra of oligodeoxyribonucleotide duplex lac operator-repressor headpiece complexes: importance of phosphate ester backbone flexibility in protein-DNA recognition.

The 31P NMR spectra of various 14-base-pair lac operators bound to both wild-type and mutant lac repressor headpiece proteins were analyzed to provide information on the backbone conformation in the complexes. The 31P NMR spectrum of a wild-type symmetrical operator, d(TGTGAGCGCTCACA)2, bound to the N-terminal 56-residue headpiece fragment of a Y7I mutant repressor was nearly identical to the spectrum of the same operator bound to the wild-type repressor headpiece. In contrast, the 31P NMR spectrum of the mutant operator, d(TATAGAGCGCTCATA)2, wild-type headpiece complex was significantly perturbed relative to the wild-type repressor-operator complex. The 31P chemical shifts of the phosphates of a second mutant operator, d(TGTGTGCGCACACA)2, showed small but specific changes upon complexation with either the wild-type or mutant headpiece. The 31P chemical shifts of the phosphates of a third mutant operator, d(TCTGAGCGCTCAGA)2, showed no perturbations upon addition of the wild-type headpiece. The 31P NMR results provide further evidence for predominant recognition of the 5'-strand of the 5'-TGTGA/3'-ACACT binding site in a 2:1 protein to headpiece complex. It is proposed that specific, strong-binding operator-protein complexes retain the inherent phosphate ester conformational flexibility of the operator itself, whereas the phosphate esters are conformationally restricted in the weak-binding operator-protein complexes. This retention of backbone torsional freedom in strong complexes is entropically favorable and provides a new (and speculative) mechanism for protein discrimination of different operator binding sites. It demonstrates the potential importance of phosphate geometry and flexibility on protein recognition and binding.