Relative solvent accessible surface area predicts protein conformational changes upon binding

Protein interactions are often accompanied by significant changes in conformation. We have analyzed the relationships between protein structures and the conformational changes they undergo upon binding. Based upon this, we introduce a simple measure, the relative solvent accessible surface area, which can be used to predict the magnitude of binding-induced conformational changes from the structures of either monomeric proteins or bound subunits. Applying this to a large set of protein complexes suggests that large conformational changes upon binding are common. In addition, we observe considerable enrichment of intrinsically disordered sequences in proteins predicted to undergo large conformational changes. Finally, we demonstrate that the relative solvent accessible surface area of monomeric proteins can be used as a simple proxy for protein flexibility. This reveals a powerful connection between the flexibility of unbound proteins and their binding-induced conformational changes, consistent with the conformational selection model of molecular recognition.     

Side-chain conformational changes upon Protein-Protein Association

Conformational changes upon protein-protein association are the key element of the binding mechanism. The study presents a systematic large-scale analysis of such conformational changes in the side chains. The results indicate that short and long side chains have different propensities for the conformational changes. Long side chains with three or more dihedral angles are often subject to large conformational transition. Shorter residues with one or two dihedral angles typically undergo local conformational changes not leading to a conformational transition. A relationship between the local readjustments and the equilibrium fluctuations of a side chain around its unbound conformation is suggested. Most of the side chains undergo larger changes in the dihedral angle most distant from the backbone. The frequencies of the core-to-surface interface transitions of six nonpolar residues and Tyr are larger than the frequencies of the opposite surface-to-core transitions. The binding increases both polar and nonpolar interface areas. However, the increase of the nonpolar area is larger for all considered classes of protein complexes, suggesting that the protein association perturbs the unbound interfaces to increase the hydrophobic contribution to the binding free energy. To test modeling approaches to side-chain flexibility in protein docking, conformational changes in the X-ray set were compared with those in the docking decoy sets. The results lead to a better understanding of the conformational changes in proteins and suggest directions for efficient conformational sampling in docking protocols.     

Conformational changes in DNA‐binding proteins: Relationships with precomplex features and contributions to specificity and stability

Both Proteins and DNA undergo conformational changes in order to form functional complexes and also to facilitate interactions with other molecules. These changes have direct implications for the stability and specificity of the complex, as well as the cooperativity of interactions between multiple entities. In this work, we have extensively analyzed conformational changes in DNA‐binding proteins by superimposing DNA‐bound and unbound pairs of protein structures in a curated database of 90 proteins. We manually examined each of these pairs, unified the authors' annotations, and summarized our observations by classifying conformational changes into six structural categories. We explored a relationship between conformational changes and functional classes, binding motifs, target specificity, biophysical features of unbound proteins, and stability of the complex. In addition, we have also investigated the degree to which the intrinsic flexibility can explain conformational changes in a subset of 52 proteins with high quality coordinate data. Our results indicate that conformational changes in DNA‐binding proteins contribute significantly to both the stability of the complex and the specificity of targets recognized by them. We also conclude that most conformational changes occur in proteins interacting with specific DNA targets, even though unbound protein structures may have sufficient information to interact with DNA in a nonspecific manner. Proteins 2014; 82:841–857. © 2013 Wiley Periodicals, Inc.     

Statistical analysis of structural determinants for protein–DNA‐binding specificity

DNA‐binding proteins play critical roles in biological processes including gene expression, DNA packaging and DNA repair. They bind to DNA target sequences with different degrees of binding specificity, ranging from highly specific (HS) to nonspecific (NS). Alterations of DNA‐binding specificity, due to either genetic variation or somatic mutations, can lead to various diseases. In this study, a comparative analysis of protein–DNA complex structures was carried out to investigate the structural features that contribute to binding specificity. Protein–DNA complexes were grouped into three general classes based on degrees of binding specificity: HS, multispecific (MS), and NS. Our results show a clear trend of structural features among the three classes, including amino acid binding propensities, simple and complex hydrogen bonds, major/minor groove and base contacts, and DNA shape. We found that aspartate is enriched in HS DNA binding proteins and predominately binds to a cytosine through a single hydrogen bond or two consecutive cytosines through bidentate hydrogen bonds. Aromatic residues, histidine and tyrosine, are highly enriched in the HS and MS groups and may contribute to specific binding through different mechanisms. To further investigate the role of protein flexibility in specific protein–DNA recognition, we analyzed the conformational changes between the bound and unbound states of DNA‐binding proteins and structural variations. The results indicate that HS and MS DNA‐binding domains have larger conformational changes upon DNA‐binding and larger degree of flexibility in both bound and unbound states. Proteins 2016; 84:1147–1161. © 2016 Wiley Periodicals, Inc.     

Ubiquitin Dynamics in Complexes Reveal Molecular Recognition Mechanisms Beyond Induced Fit and Conformational Selection

Protein-protein interactions play an important role in all biological processes. However, the principles underlying these interactions are only beginning to be understood. Ubiquitin is a small signalling protein that is covalently attached to different proteins to mark them for degradation, regulate transport and other functions. As such, it interacts with and is recognised by a multitude of other proteins. We have conducted molecular dynamics simulations of ubiquitin in complex with 11 different binding partners on a microsecond timescale and compared them with ensembles of unbound ubiquitin to investigate the principles of their interaction and determine the influence of complex formation on the dynamic properties of this protein. Along the main mode of fluctuation of ubiquitin, binding in most cases reduces the conformational space available to ubiquitin to a subspace of that covered by unbound ubiquitin. This behaviour can be well explained using the model of conformational selection. For lower amplitude collective modes, a spectrum of zero to almost complete coverage of bound by unbound ensembles was observed. The significant differences between bound and unbound structures are exclusively situated at the binding interface. Overall, the findings correspond neither to a complete conformational selection nor induced fit scenario. Instead, we introduce a model of conformational restriction, extension and shift, which describes the full range of observed effects.     

From Induced Fit to Conformational Selection: A Continuum of Binding Mechanism Controlled by the Timescale of Conformational Transitions

In receptor-ligand binding, a question that generated considerable interest is whether the mechanism is induced fit or conformational selection. This question is addressed here by a solvable model, in which a receptor undergoes transitions between active and inactive forms. The inactive form is favored while unbound but the active form is favored while a ligand is loosely bound. As the active-inactive transition rates increase, the binding mechanism gradually shifts from conformational selection to induced fit. The timescale of conformational transitions thus plays a crucial role in controlling binding mechanisms.     

The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking

We investigate the extent to which the conformational fluctuations of proteins in solution reflect the conformational changes that they undergo when they form binary protein-protein complexes. To do this, we study a set of 41 proteins that form such complexes and whose three-dimensional structures are known, both bound in the complex and unbound. We carry out molecular dynamics simulations of each protein, starting from the unbound structure, and analyze the resulting conformational fluctuations in trajectories of 5 ns in length, comparing with the structure in the complex. It is found that fluctuations take some parts of the molecules into regions of conformational space close to the bound state (or give information about it), but at no point in the simulation does each protein as whole sample the complete bound state. Subsequent use of conformations from a clustered MD ensemble in rigid-body docking is nevertheless partially successful when compared to docking the unbound conformations, as long as the unbound conformations are themselves included with the MD conformations and the whole globally rescored. For one key example where sub-domain motion is present, a ribonuclease inhibitor, principal components analysis of the MD was applied and was also able to produce conformations for docking that gave enhanced results compared to the unbound. The most significant finding is that core interface residues show a tendency to be less mobile (by size of fluctuation or entropy) than the rest of the surface even when the other binding partner is absent, and conversely the peripheral interface residues are more mobile. This surprising result, consistent across up to 40 of the 41 proteins, suggests different roles for these regions in protein recognition and binding, and suggests ways that docking algorithms could be improved by treating these regions differently in the docking process.     

Robust recognition of zinc binding sites in proteins

Metals play a variety of roles in biological processes, and hence their presence in a protein structure can yield vital functional information. Because the residues that coordinate a metal often undergo conformational changes upon binding, detection of binding sites based on simple geometric criteria in proteins without bound metal is difficult. However, aspects of the physicochemical environment around a metal binding site are often conserved even when this structural rearrangement occurs. We have developed a Bayesian classifier using known zinc binding sites as positive training examples and nonmetal binding regions that nonetheless contain residues frequently observed in zinc sites as negative training examples. In order to allow variation in the exact positions of atoms, we average a variety of biochemical and biophysical properties in six concentric spherical shells around the site of interest. At a specificity of 99.8%, this method achieves 75.5% sensitivity in unbound proteins at a positive predictive value of 73.6%. We also test its accuracy on predicted protein structures obtained by homology modeling using templates with 30%-50% sequence identity to the target sequences. At a specificity of 99.8%, we correctly identify at least one zinc binding site in 65.5% of modeled proteins. Thus, in many cases, our model is accurate enough to identify metal binding sites in proteins of unknown structure for which no high sequence identity homologs of known structure exist. Both the source code and a Web interface are available to the public at http://feature.stanford.edu/metals.     

What is the relationship between the global structures of apo and holo proteins?

It is well known that ligand binding and release may induce a wide range of structural changes in a receptor protein, varying from small movements of loops or side chains in the binding pocket to large‐scale domain hinge‐bending and shear motions or even partial unfolding that facilitates the capture and release of a ligand. An interesting question is what in general are the conformational changes triggered by ligand binding? The aim of this work is analyze the magnitude of structural changes in a protein resulting from ligand binding to assess if the state of ligand binding needs to be included in template‐based protein structure prediction algorithms. To address this issue, a nonredundant dataset of 521 paired protein structures in the ligand‐free and ligand‐bound form was created and used to estimate the degree of both local and global structure similarity between the apo and holo forms. In most cases, the proteins undergo relatively small conformational rearrangements of their tertiary structure upon ligand binding/release (most root‐mean‐square‐deviations from native, RMSD, are <1 Å). However, a clear difference was observed between single‐ and multiple‐domain proteins. For the latter, RMSD changes greater than 1 Å and sometimes larger were found for almost 1/3 of the cases; these are mainly associated with large‐scale hinge‐bending movements of entire domains. The changes in the mutual orientation of individual domains in multiple‐domain proteins upon ligand binding were investigated using a mechanistic model based on mass‐weighted principal axes as well as interface buried surface calculations. Some preferences toward the anticipated mechanism of protein domain movements are predictable based on the examination of just the ligand‐free structural form. These results have applications to protein structure prediction, particularly in the context of protein domain assembly, if additional information concerning ligand binding is exploited. Proteins 2008. © 2007 Wiley‐Liss, Inc.     

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.     

Evaluating conformational changes in protein structures binding RNA

Many protein-RNA recognition events are known to exhibit conformational changes from qualitative observations of individual complexes. However, a quantitative estimation of conformational changes is required if protein-RNA docking and template-based methods for RNA binding site prediction are to be developed. This study presents the first quantitative evaluation of conformational changes that occur when proteins bind RNA. The analysis of twelve RNA-binding proteins in the bound and unbound states using error-scaled difference distance matrices is presented. The binding site residues are mapped to each structure, and the conformational changes that affect these residues are evaluated. Of the twelve proteins four exhibit greater movements in nonbinding site residues, and a further four show the greatest movements in binding site residues. The remaining four proteins display no significant conformational change. When interface residues are found to be in conformationally variable regions of the protein they are typically seen to move less than 2 A between the bound and unbound conformations. The current data indicate that conformational changes in the binding site residues of RNA binding proteins may not be as significant as previously suggested, but a larger data set is required before wider conclusions may be drawn. The implications of the observed conformational changes for protein function prediction are discussed.     

19F NMR studies of the D-galactose chemosensory receptor. 1. Sugar binding yields a global structural change.

The Escherichia coli D-galactose and D-glucose receptor is an aqueous sugar-binding protein and the first component in the distinct chemosensory and transport pathways for these sugars. Activation of the receptor occurs when the sugar binds and induces a conformational change, which in turn enables docking to specific membrane proteins. Only the structure of the activated receptor containing bound D-glucose is known. To investigate the sugar-induced structural change, we have used 19F NMR to probe 12 sites widely distributed in the receptor molecule. Five sites are tryptophan positions probed by incorporation of 5-fluorotryptophan; the resulting 19F NMR resonances were assigned by site-directed mutagenesis. The other seven sites are phenylalanine positions probed by incorporation of 3-fluorophenylalanine. Sugar binding to the substrate binding cleft was observed to trigger a global structural change detected via 19F NMR frequency shifts at 10 of the 12 labeled sites. Two of the altered sites lie in the substrate binding cleft in van der Waals contact with the bound sugar molecule. The other eight altered sites, specifically two tryptophans and six phenylalanines distributed equally between the two receptor domains, are distant from the cleft and therefore experience allosteric structural changes upon sugar binding. The results are consistent with a model in which multiple secondary structural elements, known to extend between the substrate cleft and the protein surface, undergo shifts in their average positions upon sugar binding to the cleft. Such structural coupling provides a mechanism by which sugar binding to the substrate cleft can cause structural changes at one or more docking sites on the receptor surface.     

The impact of protein flexibility on protein-protein docking

Accounting for protein flexibility in protein-protein docking algorithms is challenging, and most algorithms therefore treat proteins as rigid bodies or permit side-chain motion only. While the consequences are obvious when there are large conformational changes upon binding, the situation is less clear for the modest conformational changes that occur upon formation of most protein-protein complexes. We have therefore studied the impact of local protein flexibility on protein-protein association by means of rigid body and torsion angle dynamics simulation. The binding of barnase and barstar was chosen as a model system for this study, because the complexation of these 2 proteins is well-characterized experimentally, and the conformational changes accompanying binding are modest. On the side-chain level, we show that the orientation of particular residues at the interface (so-called hotspot residues) have a crucial influence on the way contacts are established during docking from short protein separations of approximately 5 A. However, side-chain torsion angle dynamics simulations did not result in satisfactory docking of the proteins when using the unbound protein structures. This can be explained by our observations that, on the backbone level, even small (2 A) local loop deformations affect the dynamics of contact formation upon docking. Complementary shape-based docking calculations confirm this result, which indicates that both side-chain and backbone levels of flexibility influence short-range protein-protein association and should be treated simultaneously for atomic-detail computational docking of proteins.     

Elucidating the Role of Microprocessor Protein DGCR8 in Bending RNA Structures

Although conformational dynamics of RNA molecules are potentially important in microRNA (miRNA) processing, the role of the protein binding partners in facilitating the requisite structural changes is not well understood. In previous work, we and others have demonstrated that nonduplex structural elements and the conformational flexibility they support are necessary for efficient RNA binding and cleavage by the proteins associated with the two major stages of miRNA processing. However, recent studies showed that the protein DGCR8 binds primary miRNA and duplex RNA with similar affinities. Here, we study RNA binding by a small recombinant construct of the DGCR8 protein and the RNA conformation changes that result. This construct, the DGCR8 core, contains two double-stranded RNA-binding domains (dsRBDs) and a C-terminal tail. To assess conformational changes resulting from binding, we applied small-angle x-ray scattering with contrast variation to detect conformational changes of primary-miR-16-1 in complex with the DGCR8 core. This method reports only on the RNA conformation within the complex and suggests that the protein bends the RNA upon binding. Supporting work using smFRET to study the conformation of RNA duplexes bound to the core also shows bending. Together, these studies elucidate the role of DGCR8 in interacting with RNA during the early stages of miRNA processing.     

Residual Structure Accelerates Binding of Intrinsically Disordered ACTR by Promoting Efficient Folding upon Encounter

Intrinsically disordered proteins (IDPs) often fold into stable structures upon specific binding. The roles of residual structure of unbound IDPs in coupling binding and folding have been under much debate. While many studies emphasize the importance of conformational flexibility for IDP recognition, it was recently demonstrated that stabilization the N-terminal helix of intrinsically disordered ACTR accelerated its binding to another IDP, NCBD of the CREB-binding protein. To understand how enhancing ACTR helicity accelerates binding, we derived a series of topology-based coarse-grained models that mimicked various ACTR mutants with increasing helical contents and reproduced their NCBD binding affinities. Molecular dynamics simulations were then performed to sample hundreds of reversible coupled binding and folding transitions. The results show that increasing ACTR helicity does not alter the baseline mechanism of synergistic folding, which continues to follow “extended conformational selection” with multiple stages of selection and induced folding. Importantly, these coarse-grained models, while only calibrated based on binding thermodynamics, recapitulate the observed kinetic acceleration with increasing ACTR helicity. However, the residual helices do not enhance the association kinetics via more efficient seeding of productive collisions. Instead, they allow the nonspecific collision complexes to evolve more efficiently into the final bound and folded state, which is the primary source of accelerated association kinetics. Meanwhile, reduced dissociation kinetics with increasing ACTR helicity can be directly attributed to smaller entropic cost of forming the bound state. Altogether, this study provides important mechanistic insights into how residual structure may modulate thermodynamics and kinetics of IDP interactions.     

The structural and functional role of lysine residues in the binding domain of cytochrome c in the electron transfer to cytochrome c oxidase.

The interactions of yeast iso-1 cytochrome c with bovine cytochrome c oxidase were studied using cytochrome c variants in which lysines of the binding domain were substituted by alanines. Resonance Raman spectra of the fully oxidized complexes of both proteins reveal structural changes of both the heme c and the hemes a and a3. The structural changes in cytochrome c are the same as those observed upon binding to phospholipid vesicles where the bound protein exists in two conformers, B1 and B2. Whereas the structure of B1 is the same as that of the unbound cytochrome c, the formation of B2 is associated with substantial alterations of the heme pocket. In cytochrome c oxidase, the structural changes in both hemes refer to more subtle perturbations of the immediate protein environment and may be a result of a conformational equilibrium involving two states. These changes are qualitatively different to those observed for cytochrome c oxidase upon poly-l-lysine binding. The resonance Raman spectra of the various cytochrome c/cytochrome c oxidase complexes were analyzed quantitatively. The spectroscopic studies were paralleled by steady-state kinetic measurements of the same protein combinations. The results of the spectra analysis and the kinetic studies were used to determine the stability of the complexes and the conformational equilibria B2/B1 for all cytochrome c variants. The complex stability decreases in the order: wild-type WT > J72K > K79A > K73A > K87A > J72A > K86A > K73A/K79A (where J is the natural trimethyl lysine). This order is not exhibited by the conformational equilibria. The electrostatic control of state B2 formation does not depend on individual intermolecular salt bridges, but on the charge distribution in a specific region of the front surface of cytochrome c that is defined by the lysyl residues at positions 72, 73 and 79. On the other hand, the conformational changes in cytochrome c oxidase were found to be independent of the identity of the bound cytochrome c variant. The maximum rate constants determined from steady-state kinetic measurements could be related to the conformational equilibria of the bound cytochrome c using a simple model that assumes that the conformational transitions are faster than product formation. Within this model, the data analysis leads to the conclusion that the interprotein electron transfer rate constant is around two times higher in state B2 than in B1. These results can be interpreted in terms of an increase of the driving force in state B2 as a result of the large negative shift of the reduction potential.     

The solution structure of BMPR-IA reveals a local disorder-to-order transition upon BMP-2 binding

The structure of the extracellular domain of BMP receptor IA was determined in solution by NMR spectroscopy and compared to its structure when bound to its ligand BMP-2. While most parts of the secondary structure are highly conserved between the bound and unbound forms, large conformational rearrangements can be observed in the beta4beta5 loop of BMPR-IA, which is in contact with BMP-2 and harbors the main binding determinants for the BMPR-IA-BMP-2 interaction. In its unbound form, helix alpha1 in BMPR-IA, which is in the center of the binding epitope for BMP-2, is missing. Since BMP-2 also shows conformational changes in the type I receptor epitope upon binding to BMPR-IA, both binding partners pass through an induced fit mechanism to adapt their binding interfaces to a given interaction surface. The inherent flexibility of both partners possibly explains the promiscuous ligand-receptor interaction observed in the BMP protein superfamily.