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.     

Inherent protein structural flexibility at the RNA-binding interface of L30e

The Saccharomyces cerevisiae ribosomal protein L30 autoregulates its own expression by binding to a purine-rich internal loop in its pre-mRNA and mRNA. NMR studies of L30 and its RNA complex showed that both the internal loop of the RNA as well as a region of the protein become substantially more ordered upon binding. A crystal structure of a maltose binding protein (MBP)-L30 fusion protein with two copies in the asymmetric unit has been determined. The flexible RNA-binding region in the L30 copies has two distinct conformations, one resembles the RNA bound form solved by NMR and the other is unique. Structure prediction algorithms also had difficulty accurately predicting this region, which is consistent with conformational flexibility seen in the NMR and X-ray crystallography studies. Inherent conformational flexibility may be a hallmark of regions involved in intermolecular interactions.     

Protein flexibility in ligand docking and virtual screening to protein kinases

The main complicating factor in structure-based drug design is receptor rearrangement upon ligand binding (induced fit). It is the induced fit that complicates cross-docking of ligands from different ligand-receptor complexes. Previous studies have shown the necessity to include protein flexibility in ligand docking and virtual screening. Very few docking methods have been developed to predict the induced fit reliably and, at the same time, to improve on discriminating between binders and non-binders in the virtual screening process.We present an algorithm called the ICM-flexible receptor docking algorithm (IFREDA) to account for protein flexibility in virtual screening. By docking flexible ligands to a flexible receptor, IFREDA generates a discrete set of receptor conformations, which are then used to perform flexible ligand-rigid receptor docking and scoring. This is followed by a merging and shrinking step, where the results of the multiple virtual screenings are condensed to improve the enrichment factor. In the IFREDA approach, both side-chain rearrangements and essential backbone movements are taken into consideration, thus sampling adequately the conformational space of the receptor, even in cases of large loop movements.As a preliminary step, to show the importance of incorporating protein flexibility in ligand docking and virtual screening, and to validate the merging and shrinking procedure, we compiled an extensive small-scale virtual screening benchmark of 33 crystal structures of four different protein kinases sub-families (cAPK, CDK-2, P38 and LCK), where we obtained an enrichment factor fold-increase of 1.85±0.65 using two or three multiple experimental conformations. IFREDA was used in eight protein kinase complexes and was able to find the correct ligand conformation and discriminate the correct conformations from the “misdocked” conformations solely on the basis of energy calculation. Five of the generated structures were used in the small-scale virtual screening stage and, by merging and shrinking the results with those of the original structure, we show an enrichment factor fold increase of 1.89±0.60, comparable to that obtained using multiple experimental conformations.Our cross-docking tests on the protein kinase benchmark underscore the necessity of incorporating protein flexibility in both ligand docking and virtual screening. The methodology presented here will be extremely useful in cases where few or no experimental structures of complexes are available, while some binders are known.     

Reaction pattern of Photosystem II: oxidative water cleavage and protein flexibility

This short communication addresses three topics of photosynthetic water cleavage in Photosystem II (PS II): (a) effect of protonation in the acidic range on the extent of the ‘fast’ ns kinetics of P680^+· reduction by Y_Z, (b) mechanism of O–O bond formation and (c) role of protein flexibility in the functional integrity of PS II. Based on measurements of light-induced absorption changes and quasielastic neutron scattering in combination with mechanistic considerations, evidence is presented for the protein acting as a functionally active constituent of the water cleavage machinery, in particular, for directed local proton transfer. A specific flexibility emerging above a threshold of about 230 K is an indispensable prerequisite for oxygen evolution and plastoquinol formation.     

Future directions in folding: The multi-state nature of protein structure

All possible protein folding intermediates exist in equilibrium with the native protein at native as well as non-native conditions, with occupation determined by their free energy level. The study of these forms can illuminate the fundamental principles of protein structure and folding. Hydrogen exchange methods can be used to detect and characterize these partially unfolded forms at native conditions and as a function of mild denaturant and temperature. This information illuminates the requirements that govern the ability of kinetic and equilibrium methods to study folding intermediates.     

Understanding protein-ligand interactions: the price of protein flexibility

In order to design selective, high-affinity ligands to a target protein, it is advantageous to understand the structural determinants for protein-ligand complex formation at the atomic level. In a model system, we have successively mapped the factor Xa binding site onto trypsin, showing that certain mutations influence both protein structure and inhibitor specificity. Our previous studies have shown that introduction of the 172SSFI175 sequence of factor Xa into rat or bovine trypsin results in the destabilisation of the intermediate helix with burial of Phe174 (the down conformation). Surface exposure of the latter residue (the up conformation) is critical for the correct formation of the aromatic box found in factor Xa-ligand complexes. In the present study, we investigate the influence of aromatic residues in position 174. Replacement with the bulky tryptophan (SSWI) shows reduced affinity for benzamidine-based inhibitors (1) and (4), whereas removal of the side-chain (alanine, SSAI) or exchange with a hydrophilic residue (arginine, SSRI) leads to a significant loss in affinity for all inhibitors studied. The variants could be crystallised in the presence of different inhibitors in multiple crystal forms. Structural characterisation of the variants revealed three different conformations of the intermediate helix and 175 loop in SSAI (down, up and super-up), as well as a complete disorder of this region in one crystal form of SSRI, suggesting that the compromised affinity of these variants is related to conformational flexibility. The influence of Glu217, peripheral to the ligand-binding site in factor Xa, was investigated. Introduction of Glu217 into trypsin variants containing the SSFI sequence exhibited enhanced affinity for the factor Xa ligands (2) and (3). The crystal structures of these variants also exhibited the down and super-up conformations, the latter of which could be converted to up upon soaking and binding of inhibitor (2). The improved affinity of the Glu217-containing variants appears to be due to a shift towards the up conformation. Thus, the reduction in affinity caused by conformational variability of the protein target can be partially or wholly offset by compensatory binding to the up conformation. The insights provided by these studies will be helpful in improving our understanding of ligand binding for the drug design process.     

Disulfide crosslinks to probe the structure and flexibility of a designed four-helix bundle protein.

The introduction of disulfide crosslinks is a generally useful method by which to identify regions of a protein that are close together in space. Here we describe the use of disulfide crosslinks to investigate the structure and flexibility of a family of designed 4-helix bundle proteins. The results of these analyses lend support to our working model of the proteins' structure and suggest that the proteins have limited main-chain flexibility.     

Relative flexibility of DNA and RNA: a molecular dynamics study

State of the art molecular dynamics simulations are used to study the structure, dynamics, molecular interaction properties and flexibility of DNA and RNA duplexes in aqueous solution. Special attention is paid to the deformability of both types of structures, revisiting concepts on the relative flexibility of DNA and RNA duplexes. Our simulations strongly suggest that the concepts of flexibility, rigidity and deformability are much more complex than usually believed, and that it is not always true that DNA is more flexible than RNA.     

Conformational flexibility in calcitonin: The dynamic properties of human and salmon calcitonin in solution

We have studied the dynamic properties of human (h) and salmon (s) calcitonin (CT) in solution. For both hormones, distance geometry in torsion-angle space has been used to generate three-dimensional structures consistent with NMR data obtained in sodium dodecyl sulfate micelles. For sCT and hCT we used, respectively, 356 and 275 interproton distances together with hydrogen-bonds as restraints. To better characterize their flexibility and dynamic properties two fully unrestrained 1100-ps molecular dynamics (MD) simulations in methanol were performed on the lowest-energy structures of both hormones. Statistical analyses of average geometric parameters and of their fluctuations performed in the last 1000 ps of the MD run show typical helical values for residues 9–19 of sCT during the whole trajectory. For hCT a shorter helix was observed involving residues 13–21, with a constant helical region in the range 13–19. Angular order parameters S() and S() indicate that hCT exhibits a higher flexibility, distributed along the whole chain, including the helix, while the only flexible amino acid residues in sCT connect three well-defined domains. Finally, our study shows that simulated annealing in torsion-angle space can efficiently be extended to NMR-based three-dimensional structure calculations of helical polypeptides. Furthermore, provided that a sufficient number of NMR restraints describes the system, the method allows the detection of equilibria in solution. This identification occurs through the generation of 'spurious' high-energy structures, which, for right-handed -helices, are likely to be represented by left-handed -helices.     

Incorporation of side chain flexibility into protein binding pockets using MTflex

The plasticity of active sites plays a significant role in drug recognition and binding, but the accurate incorporation of ‘receptor flexibility’ remains a significant computational challenge. Many approaches have been put forward to address receptor flexibility in docking studies by generating relevant ensembles on the energy surface. Herein, we describe the Movable Type with flexibility (MT_flex) method that generates ensembles on the more relevant free energy surface in a computationally tractable manner. This novel approach enumerates conformational states based on side chain flexibility and then estimates their relative free energies using the MT methodology. The resultant conformational states can then be used in subsequent docking and scoring exercises. In particular, we demonstrate that using the MT_flex ensembles improves MT’s ability to predict binding free energies over docking only to the crystal structure.     

Frataxin from Psychromonas ingrahamii as a model to study stability modulation within the CyaY protein family

Adaptation of life to low temperatures influences both protein stability and flexibility. Thus, proteins from psychrophilic organisms are excellent models to study relations between these properties. Here we focused on frataxin from Psychromonas ingrahamii (pFXN), an extreme psychrophilic sea ice bacterium that can grow at temperatures as low as − 12 °C. This α/β protein is highly conserved and plays a key role in iron homeostasis as an iron chaperone. In contrast to other frataxin homologs, chemical and temperature unfolding experiments showed that the thermodynamic stability of pFXN is strongly modulated by pHs: ranging from 5.5 ± 0.9 (pH 6.0) to 0.9 ± 0.3 kcal mol^− 1 (pH 8.0). This protein was crystallized and its X-ray structure solved at 1.45 Å. Comparison of B-factor profiles between Escherichia coli and P. ingrahamii frataxin variants (51% of identity) suggests that, although both proteins share the same structural features, their flexibility distribution is different. Molecular dynamics simulations showed that protonation of His44 or His67 in pFXN lowers the mobility of regions encompassing residues 20–30 and the C-terminal end, probably through favorable electrostatic interactions with residues Asp27, Glu42 and Glu99. Since the C-terminal end of the protein is critical for the stabilization of the frataxin fold, the predictions presented may be reporting on the microscopic origin of the decrease in global stability produced near neutral pH in the psychrophilic variant. We propose that suboptimal electrostatic interactions may have been an evolutionary strategy for the adaptation of frataxin flexibility and function to cold environments.     

Comparative analysis of essential collective dynamics and NMR-derived flexibility profiles in evolutionarily diverse prion proteins

Collective motions on ns-µs time scales are known to have a major impact on protein folding, stability, binding and enzymatic efficiency. It is also believed that these motions may have an important role in the early stages of prion protein misfolding and prion disease. In an effort to accurately characterize these motions and their potential influence on the misfolding and prion disease transmissibility we have conducted a combined analysis of molecular dynamic simulations and NMR-derived flexibility measurements over a diverse range of prion proteins. Using a recently developed numerical formalism, we have analyzed the essential collective dynamics (ECD) for prion proteins from eight different species including human, cow, elk, cat, hamster, chicken, turtle and frog. We also compared the numerical results with flexibility profiles generated by the random coil index (RCI) from NMR chemical shifts. Prion protein backbone flexibility derived from experimental NMR data and from theoretical computations show strong agreement with each other, demonstrating that it is possible to predict the observed RCI profiles employing the numerical ECD formalism. Interestingly, flexibility differences in the loop between second b strand (S2) and the second a helix (HB) appear to distinguish prion proteins from species that are susceptible to prion disease and those that are resistant. Our results show that the different levels of flexibility in the S2-HB loop in various species are predictable via the ECD method, indicating that ECD may be used to identify disease resistant variants of prion proteins, as well as the influence of prion proteins mutations on disease susceptibility or misfolding propensity.     

The dynamic properties of the Hepatitis C Virus E2 envelope protein unraveled by molecular dynamics

Hepatitis C Virus (HCV) is one of the most persistent human viruses. Although effective therapeutic approaches have been recently discovered, their use is limited by the elevated costs. Therefore, the development of alternative/complementary strategies is an urgent need. The E2 glycoprotein, the most immunogenic HCV protein, and its variants represent natural candidates to achieve this goal. Here we report an extensive molecular dynamics (MD) analysis of the intrinsic properties of E2. Our data provide interesting clues on the global and local intrinsic dynamic features of the protein. Present MD data clearly indicate that E2 combines a flexible structure with a network of covalent bonds. Moreover, the analysis of the two most important antigenic regions of the protein provides some interesting insights into their intrinsic structural and dynamic properties. Our data indicate that a fluctuating β-hairpin represents a populated state by the region E2^412?423. Interestingly, the analysis of the epitope E2^427?446 conformation, that undergoes a remarkable rearrangement in the simulation, has significant similarities with the structure that the E2^430?442 fragment adopts in complex with a neutralizing antibody. Present data also suggest that the strict conservation of Gly436 in E2 protein of different HCV genotypes is likely dictated by structural restraints. Moreover, the analysis of the E2^412?423 flexibility provides insights into the mechanisms that some antibodies adopt to anchor Trp437 that is fully buried in E2. Finally, the present investigation suggests that MD simulations should systematically complement crystallographic studies on flexible proteins that are studied in combination with antibodies.     

Conformational flexibility of N-glycans in solution studied by REMD simulations

Protein–glycan recognition regulates a wide range of biological and pathogenic processes. Conformational diversity of glycans in solution is apparently incompatible with specific binding to their receptor proteins. One possibility is that among the different conformational states of a glycan, only one conformer is utilized for specific binding to a protein. However, the labile nature of glycans makes characterizing their conformational states a challenging issue. All-atom molecular dynamics (MD) simulations provide the atomic details of glycan structures in solution, but fairly extensive sampling is required for simulating the transitions between rotameric states. This difficulty limits application of conventional MD simulations to small fragments like di- and tri-saccharides. Replica-exchange molecular dynamics (REMD) simulation, with extensive sampling of structures in solution, provides a valuable way to identify a family of glycan conformers. This article reviews recent REMD simulations of glycans carried out by us or other research groups and provides new insights into the conformational equilibria of N-glycans and their alteration by chemical modification. We also emphasize the importance of statistical averaging over the multiple conformers of glycans for comparing simulation results with experimental observables. The results support the concept of “conformer selection” in protein–glycan recognition.     

Two‐stage dynamic DNA quality check by xeroderma pigmentosum group C protein

Xeroderma pigmentosum group C (XPC) protein initiates the DNA excision repair of helix‐distorting base lesions. To understand how this versatile subunit searches for aberrant sites within the vast background of normal genomic DNA, the real‐time redistribution of fluorescent fusion constructs was monitored after high‐resolution DNA damage induction. Bidirectional truncation analyses disclosed a surprisingly short recognition hotspot, comprising ～15% of human XPC, that includes two β‐hairpin domains with a preference for non‐hydrogen‐bonded bases in double‐stranded DNA. However, to detect damaged sites in living cells, these DNA‐attractive domains depend on the partially DNA‐repulsive action of an adjacent β‐turn extension that promotes the mobility of XPC molecules searching for lesions. The key function of this dynamic interaction surface is shown by a site‐directed charge inversion, which results in increased affinity for native DNA, retarded nuclear mobility and diminished repair efficiency. These studies reveal a two‐stage discrimination process, whereby XPC protein first deploys a dynamic sensor interface to rapidly interrogate the double helix, thus forming a transient recognition intermediate before the final installation of a more static repair‐initiating complex.     

Evolutionary computer programming of protein folding and structure predictions

In order to understand the mechanism of protein folding and to assist the rational de-novo design of fast-folding, non-aggregating and stable artificial enzymes it is very helpful to be able to simulate protein folding reactions and to predict the structures of proteins and other biomacromolecules. Here, we use a method of computer programming called "evolutionary computer programming" in which a program evolves depending on the evolutionary pressure exerted on the program. In the case of the presented application of this method on a computer program for folding simulations, the evolutionary pressure exerted was towards faster finding deep minima in the energy landscape of protein folding. Already after 20 evolution steps, the evolved program was able to find deep minima in the energy landscape more than 10 times faster than the original program prior to the evolution process.     

A single residue exchange between two HLA-B27 alleles triggers increased peptide flexibility

For more than 30 years, human leukocyte antigen B27 (HLA-B27) has been known to be closely related to the autoimmune disease ankylosing spondylitis, yet little is known about the molecular mechanisms of pathogenesis. Crystal structures of two closely related, but differently disease-associated, subtypes (B*2705 and B*2709) also did not resolve this situation as they revealed the bound nonapeptide in essentially identical conformations. As the peptide is part of putative binding epitopes for the T cell receptor, we performed molecular dynamics simulations to gain deeper insight into the dynamic behaviour of HLA-B27 molecules. We find increased flexibility of the peptide in the binding groove of subtype B*2709 due to weaker interactions in the F pocket. Possible implications of this flexibility for T cell recognition and signalling are discussed.Abbreviations ₂m ₂-microglobulin - AS ankylosing spondylitis - CDR complementarity determining region - HC heavy chain - HLA human leukocyte antigen - MD molecular dynamics - MHC major histocompatibility complex - pMHC peptide-loaded MHC - RMSD root mean square deviation - RMSF root mean square fluctuation - TCR T cell receptor An erratum to this article can be found at