Features of Large Hinge-Bending Conformational Transitions. Prediction of Closed Structure from Open State

We performed a detailed analysis of conformational transition pathways for a set of 10 proteins, which undergo large hinge-bending-type motions with 4–12 Å RMSD (root mean-square distance) between open and closed crystal structures. Anisotropic network model-Monte Carlo (ANM-MC) algorithm generates a targeted pathway between two conformations, where the collective modes from the ANM are used for deformation at each iteration and the conformational energy of the deformed structure is minimized via an MC algorithm. The target structure was approached successfully with an RMSD of 0.9–4.1 Å when a relatively low cutoff radius of 10 Å was used in ANM. Even though one predominant mode (first or second) directed the open-to-closed conformational transition, changes in the dominant mode character were observed for most cases along the transition. By imposing radius of gyration constraint during mode selection, it was possible to predict the closed structure for eight out of 10 proteins (with initial 4.1–7.1 Å and final 1.7–2.9 Å RMSD to target). Deforming along a single mode leads to most successful predictions. Based on the previously reported free energy surface of adenylate kinase, deformations along the first mode produced an energetically favorable path, which was interestingly facilitated by a change in mode shape (resembling second and third modes) at key points. Pathway intermediates are provided in our database of conformational transitions (http://safir.prc.boun.edu.tr/anmmc/method/1).     

Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes

Conformational transitions between open/closed or free/bound states in proteins possess functional importance. We propose a technique in which the collective modes obtained from an anisotropic network model (ANM) are used in conjunction with a Monte Carlo (MC) simulation approach, to investigate conformational transition pathways and pathway intermediates. The ANM-MC technique is applied to adenylate kinase (AK) and hemoglobin. The iterative method, in which normal modes are continuously updated during the simulation, proves successful in accomplishing the transition between open-closed conformations of AK and tense-relaxed forms of hemoglobin (C^α− root mean square deviations between two end structures of 7.13 Å and 3.55 Å, respectively). Target conformations are reached by root mean-square deviations of 2.27 Å and 1.90 Å for AK and hemoglobin, respectively. The intermediate conformations overlap with crystal structures from the AK family within a 3.0-Å root mean-square deviation. In the case of hemoglobin, the transition of tense-to-relaxed passes through the relaxed state. In both cases, the lowest-frequency modes are effective during transitions. The targeted Monte Carlo approach is used without the application of collective modes. Both the ANM-MC and targeted Monte Carlo techniques can explore sequences of events in transition pathways with an efficient yet realistic conformational search.     

Ligand Docking to Intermediate and Close-To-Bound Conformers Generated by an Elastic Network Model Based Algorithm for Highly Flexible Proteins

Incorporating receptor flexibility in small ligand-protein docking still poses a challenge for proteins undergoing large conformational changes. In the absence of bound structures, sampling conformers that are accessible by apo state may facilitate docking and drug design studies. For this aim, we developed an unbiased conformational search algorithm, by integrating global modes from elastic network model, clustering and energy minimization with implicit solvation. Our dataset consists of five diverse proteins with apo to complex RMSDs 4.7–15 Å. Applying this iterative algorithm on apo structures, conformers close to the bound-state (RMSD 1.4–3.8 Å), as well as the intermediate states were generated. Dockings to a sequence of conformers consisting of a closed structure and its “parents” up to the apo were performed to compare binding poses on different states of the receptor. For two periplasmic binding proteins and biotin carboxylase that exhibit hinge-type closure of two dynamics domains, the best pose was obtained for the conformer closest to the bound structure (ligand RMSDs 1.5–2 Å). In contrast, the best pose for adenylate kinase corresponded to an intermediate state with partially closed LID domain and open NMP domain, in line with recent studies (ligand RMSD 2.9 Å). The docking of a helical peptide to calmodulin was the most challenging case due to the complexity of its 15 Å transition, for which a two-stage procedure was necessary. The technique was first applied on the extended calmodulin to generate intermediate conformers; then peptide docking and a second generation stage on the complex were performed, which in turn yielded a final peptide RMSD of 2.9 Å. Our algorithm is effective in producing conformational states based on the apo state. This study underlines the importance of such intermediate states for ligand docking to proteins undergoing large transitions.     

Homology-Modelling Protein-Ligand Interactions: Allowing for Ligand-Induced Conformational Change

Current homology-modelling methods do not consider small molecules in their automated processes. Therefore, the development of a reliable tool for protein-ligand homology modelling is an important next step in generating plausible models for molecular interactions. Two automated protein-ligand homology-modelling strategies, requiring no expert knowledge from the user, are investigated here. Both employ the “induced fit” concept with flexibility in side chains and ligand. The most successful strategy superimposes the new ligand over the original ligand before homology modelling, allowing the new ligand to be taken into consideration during protein modelling (rather than after), facilitating conformational change in the local backbone if necessary. We show that this approach results in successful modelling of the ligand and key binding-site residues of angiotensin-converting enzyme 2 (ACE2) from its homologue ACE, which is not possible via conventional homology modelling or by homology modelling followed by docking. Several other difficult target complexes are also successfully modelled, reproducing native protein-ligand contacts with significantly different biological substrates and different binding-site conformations. These include the modelling of Cdk5 (cyclin-dependent kinase 5) from Cdk2, thymidine phosphorylase from a bacterial homologue, and dihydrofolate reductase from a recombinant variant with a markedly different inhibitor. In terms of average modelling quality across 82 targets, the ligand RMSD with respect to the experimental structure is 1.4 Å (and 2.0 Å for the protein binding site) for “easy” cases and 2.9 Å for the ligand (and 2.7 Å for the protein binding site) in “hard” cases. This demonstrates the importance of selecting an optimal template. Ligand-modelling accuracy is strongly dependent on target-template ligand structural similarity, rather than target-template sequence identity. However, protein-modelling accuracy is dependent on both. Our automated protein-ligand homology-modelling strategy generates a higher degree of accuracy than homology modelling followed by docking, generating an average ligand RMSD that is 1-2 Å better than docking with homology models.     

Conformational Changes and Ligand Recognition of Escherichia colid-Xylose Binding Protein Revealed

ATP binding cassette transport systems account for most import of necessary nutrients in bacteria. The periplasmic binding component (or an equivalent membrane-anchored protein) is critical to recognizing cognate ligand and directing it to the appropriate membrane permease. Here we report the X-ray structures of d-xylose binding protein from Escherichia coli in ligand-free open form, ligand-bound open form, and ligand-bound closed form at 2.15 Å, 2.2 Å, and 2.2 Å resolutions, respectively. The ligand-bound open form is the first such structure to be reported at high resolution; the combination of the three different forms from the same protein furthermore gives unprecedented details concerning the conformational changes involved in binding protein function. As is typical of the structural family, the protein has two similar globular domains, which are connected by a three-stranded hinge region. The open liganded structure shows that xylose binds first to the C-terminal domain, with only very small conformational changes resulting. After a 34° closing motion, additional interactions are formed with the N-terminal domain; changes in this domain are larger and serve to make the structure more ordered near the ligand. An analysis of the interactions suggests why xylose is the preferred ligand. Furthermore, a comparison with the most closely related proteins in the structural family shows that the conformational changes are distinct in each type of binding protein, which may have implications for how the individual proteins act in concert with their respective membrane permeases.     

Structural basis of poly(3-hydroxybutyrate) hydrolysis by PhaZ7 depolymerase from Paucimonas lemoignei

The crystal structure of poly(3-hydroxybutyrate) (PHB) depolymerase PhaZ7 purified from Paucimonas lemoignei was determined at 1.90 Å resolution. The structure consists of a single domain with an α/β hydrolase fold in its core. The active site is analogous to that of serine esterases/lipases and is characterized by the presence of a catalytic triad comprising Ser136, Asp242, and His306. Comparison with other structures in the Protein Data Bank showed a high level of similarity with the Bacillus subtilis lipase LipA (RMSD, 1.55 Å). Structural comparison with Penicillium funiculosum PHB depolymerase, the only PHB depolymerase whose structure is already known, revealed significant differences, resulting in an RMSD of 2.80-3.58 Å. The two enzymes appear to utilize different types of solvent-exposed residues for biopolymer binding, with aliphatic and hydroxyl residues used in P. funiculosum PHB depolymerase and aromatic residues in PhaZ7. Moreover, the active site of P. funiculosum PHB depolymerase is accessible to the substrate in contrast to the active site of PhaZ7, which is buried. Hence, considerable conformational changes are required in PhaZ7 for the creation of a channel leading to the active site. Taken together, the structural data suggest that PhaZ7 and P. funiculosum PHB depolymerase have adopted different strategies for effective substrate binding in response to their diverse substrate specificity and the lack of a substrate-binding domain.     

Models with energy penalty on interresidue rotation address insufficiencies of conventional elastic network models

In this study, I present a new elastic network model, to our knowledge, that addresses insufficiencies of two conventional models—the Gaussian network model (GNM) and the anisotropic network model (ANM). It has been shown previously that the GNM is not rotation-invariant due to its energy, which penalizes rigid-body rotation (external rotation). As a result, GNM models are found contaminated with rigid-body rotation, especially in the most collective ones. A new model (EPIRM) is proposed to remove such external component in modes. The extracted internal motions result from a potential that penalizes interresidue stretching and rotation in a protein. The new model is shown to pertinently describe crystallographic temperature factors (B-factors) and protein open↔closed transitions. Also, the capability of separating internal and external motions in GNM slow modes permits reexamining important mechanochemical properties in enzyme active sites. The results suggest that catalytic residues stay closer to rigid-body rotation axes than their immediate backbone neighbors. I show that the cumulative density of states for EPIRM and ANM follow different power laws as functions of low-mode frequencies. When using a cutoff distance of 7.5 Å, The cumulative density of states of EPIRM scales faster than that of all-atom normal mode analysis and slower than that of simple lattices.     

Conformational Transitions upon Ligand Binding: Holo-Structure Prediction from Apo Conformations

Biological function of proteins is frequently associated with the formation of complexes with small-molecule ligands. Experimental structure determination of such complexes at atomic resolution, however, can be time-consuming and costly. Computational methods for structure prediction of protein/ligand complexes, particularly docking, are as yet restricted by their limited consideration of receptor flexibility, rendering them not applicable for predicting protein/ligand complexes if large conformational changes of the receptor upon ligand binding are involved. Accurate receptor models in the ligand-bound state (holo structures), however, are a prerequisite for successful structure-based drug design. Hence, if only an unbound (apo) structure is available distinct from the ligand-bound conformation, structure-based drug design is severely limited. We present a method to predict the structure of protein/ligand complexes based solely on the apo structure, the ligand and the radius of gyration of the holo structure. The method is applied to ten cases in which proteins undergo structural rearrangements of up to 7.1 Å backbone RMSD upon ligand binding. In all cases, receptor models within 1.6 Å backbone RMSD to the target were predicted and close-to-native ligand binding poses were obtained for 8 of 10 cases in the top-ranked complex models. A protocol is presented that is expected to enable structure modeling of protein/ligand complexes and structure-based drug design for cases where crystal structures of ligand-bound conformations are not available.     

Predicting Ligand Binding Affinity with Alchemical Free Energy Methods in a Polar Model Binding Site

We present a combined experimental and modeling study of organic ligand molecules binding to a slightly polar engineered cavity site in T4 lysozyme (L99A/M102Q). For modeling, we computed alchemical absolute binding free energies. These were blind tests performed prospectively on 13 diverse, previously untested candidate ligand molecules. We predicted that eight compounds would bind to the cavity and five would not; 11 of 13 predictions were correct at this level. The RMS error to the measurable absolute binding energies was 1.8 kcal/mol. In addition, we computed “relative” binding free energies for six phenol derivatives starting from two known ligands: phenol and catechol. The average RMS error in the relative free energy prediction was 2.5 kcal/mol (phenol) and 1.1 kcal/mol (catechol). To understand these results at atomic resolution, we obtained x-ray co-complex structures for nine of the diverse ligands and for all six phenol analogs. The average RMSD of the predicted pose to the experiment was 2.0 Å (diverse set), 1.8 Å (phenol-derived predictions), and 1.2 Å (catechol-derived predictions). We found that predicting accurate affinities and rank-orderings required near-native starting orientations of the ligand in the binding site. Unanticipated binding modes, multiple ligand binding, and protein conformational change all proved challenging for the free energy methods. We believe that these results can help guide future improvements in physics-based absolute binding free energy methods.     

Arginine kinase: joint crystallographic and NMR RDC analyses link substrate-associated motions to intrinsic flexibility

The phosphagen kinase family, including creatine and arginine kinases (AKs), catalyzes the reversible transfer of a “high-energy” phosphate between ATP and a phosphoguanidino substrate. They have become a model for the study of both substrate-induced conformational change and intrinsic protein dynamics. Prior crystallographic studies indicated large substrate-induced domain rotations, but differences among a recent set of AK structures were interpreted as a plastic deformation. Here, the structure of Limulus substrate-free AK is refined against high-resolution crystallographic data and compared quantitatively with NMR chemical shifts and residual dipolar couplings (RDCs). This demonstrates the feasibility of this type of RDC analysis of proteins that are large by NMR standards (42 kDa) and illuminates the solution structure, free from crystal-packing constraints. Detailed comparison of the 1.7 Å resolution substrate-free crystal structure against the 1.7 Å transition-state analog complex shows large substrate-induced domain motions that can be broken down into movements of smaller quasi-rigid bodies. The solution-state structure of substrate-free AK is most consistent with an equilibrium of substrate-free and substrate-bound structures, with the substrate-free form dominating, but with varying displacements of the quasi-rigid groups. Rigid-group rotations evident from the crystal structures are about axes previously associated with intrinsic millisecond dynamics using NMR relaxation dispersion. Thus, “substrate-induced” motions are along modes that are intrinsically flexible in the substrate-free enzyme and likely involve some degree of conformational selection.     

Structural formation of huntingtin exon 1 aggregates probed by small-angle neutron scattering

In several neurodegenerative disorders, including Huntington's disease, aspects concerning the earliest of protein structures that form along the aggregation pathway have increasingly gained attention because these particular species are likely to be neurotoxic. We used time-resolved small-angle neutron scattering to probe in solution these transient structures formed by peptides having the N-terminal sequence context of mutant huntingtin exon 1. We obtained snapshots of the formed aggregates as the kinetic reaction ensued to yield quantitative information on their size and mass. At the early stage, small precursor species with an initial radius of gyration of 16.1 ± 5.9 Å and average mass of a dimer to trimer were monitored. Structural growth was treated as two modes with a transition from three-dimensional early aggregate formation to two-dimensional fibril growth and association. Our small-angle neutron scattering results on the internal structure of the mature fibrils demonstrate loose packing with ∼1 peptide per 4.75 Å β-sheet repeat distance, which is shown to be quantitatively consistent with a β-helix model. This research provides what we believe to be new insights into the structures forming along the pathway of huntingtin exon 1 aggregation and should assist in determining the role that precursors play in neuronal toxicity.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

A Free-Energy Approach for All-Atom Protein Simulation

All-atom free-energy methods offer a promising alternative to kinetic molecular mechanics simulations of protein folding and association. Here we report an accurate, transferable all-atom biophysical force field (PFF02) that stabilizes the native conformation of a wide range of proteins as the global optimum of the free-energy landscape. For 32 proteins of the ROSETTA decoy set and six proteins that we have previously folded with PFF01, we find near-native conformations with an average backbone RMSD of 2.14 Å to the native conformation and an average Z-score of −3.46 to the corresponding decoy set. We used nonequilibrium sampling techniques starting from completely extended conformations to exhaustively sample the energy surface of three nonhomologous hairpin-peptides, a three-stranded β-sheet, the all-helical 40 amino-acid HIV accessory protein, and a zinc-finger ββα motif, and find near-native conformations for the minimal energy for each protein. Using a massively parallel evolutionary algorithm, we also obtain a near-native low-energy conformation for the 54 amino-acid engrailed homeodomain. Our force field thus stabilized near-native conformations for a total of 20 proteins of all structure classes with an average RMSD of only 3.06 Å to their respective experimental conformations.     

Computational protein design and large-scale assessment by I-TASSER structure assembly simulations

Protein design aims at designing new protein molecules of desired structure and functionality. One of the major obstacles to large-scale protein design are the extensive time and manpower requirements for experimental validation of designed sequences. Recent advances in protein structure prediction have provided potentials for an automated assessment of the designed sequences via folding simulations. We present a new protocol for protein design and validation. The sequence space is initially searched by Monte Carlo sampling guided by a public atomic potential, with candidate sequences selected by the clustering of sequence decoys. The designed sequences are then assessed by I-TASSER folding simulations, which generate full-length atomic structural models by the iterative assembly of threading fragments. The protocol is tested on 52 nonhomologous single-domain proteins, with an average sequence identity of 24% between the designed sequences and the native sequences. Despite this low sequence identity, three-dimensional models predicted for the first designed sequence have an RMSD of < 2 Å to the target structure in 62% of cases. This percentage increases to 77% if we consider the three-dimensional models from the top 10 designed sequences. Such a striking consistency between the target structure and the structural prediction from nonhomologous sequences, despite the fact that the design and folding algorithms adopt completely different force fields, indicates that the design algorithm captures the features essential to the global fold of the target. On average, the designed sequences have a free energy that is 0.39 kcal/(mol residue) lower than in the native sequences, potentially affording a greater stability to synthesized target folds.     

The PDB is a covering set of small protein structures

Structure comparisons of all representative proteins have been done. Employing the relative root mean square deviation (RMSD) from native enables the assessment of the statistical significance of structure alignments of different lengths in terms of a Z-score. Two conclusions emerge: first, proteins with their native fold can be distinguished by their Z-score. Second and somewhat surprising, all small proteins up to 100 residues in length have significant structure alignments to other proteins in a different secondary structure and fold class; i.e. 24.0% of them have 60% coverage by a template protein with a RMSD below 3.5 Å and 6.0% have 70% coverage. If the restriction that we align proteins only having different secondary structure types is removed, then in a representative benchmark set of proteins of 200 residues or smaller, 93% can be aligned to a single template structure (with average sequence identity of 9.8%), with a RMSD less than 4 Å, and 79% average coverage. In this sense, the current Protein Data Bank (PDB) is almost a covering set of small protein structures. The length of the aligned region (relative to the whole protein length) does not differ among the top hit proteins, indicating that protein structure space is highly dense. For larger proteins, non-related proteins can cover a significant portion of the structure. Moreover, these top hit proteins are aligned to different parts of the target protein, so that almost the entire molecule can be covered when combined. The number of proteins required to cover a target protein is very small, e.g. the top ten hit proteins can give 90% coverage below a RMSD of 3.5 Å for proteins up to 320 residues long. These results give a new view of the nature of protein structure space, and its implications for protein structure prediction are discussed.     

Robust probabilistic superposition and comparison of protein structures

Background  

Protein structure comparison is a central issue in structural bioinformatics. The standard dissimilarity measure for protein structures is the root mean square deviation (RMSD) of representative atom positions such as α-carbons. To evaluate the RMSD the structures under comparison must be superimposed optimally so as to minimize the RMSD. How to evaluate optimal fits becomes a matter of debate, if the structures contain regions which differ largely - a situation encountered in NMR ensembles and proteins undergoing large-scale conformational transitions.     

Structural characterization of the active and inactive states of Src kinase in solution by small-angle X-ray scattering

Src kinase plays an important role in several signaling and regulation mechanisms in vivo. Enzymatic activity is tightly regulated through the phosphorylation and dephosphorylation of tyrosine 527, which is placed at the C-terminal tail. Here, we have addressed domain rearrangements involved in the regulation mechanism of Src kinase in solution using small-angle X-ray scattering. In the phosphorylated wild-type form of Src kinase corresponding to the inactive state of the protein, a single conformation compatible with a closed crystallographic structure was found in solution. In the Y527F point mutant representing the active state, analysis of scattering data reveals an equilibrium between two differently populated conformations differing in the radius of gyration by 5 Å. The major species (85% of the total population) presents a closed conformation indistinguishable from the crystallographic structure of the inactive state. The minor species (15% of the total population) is an open conformation similar to the crystallographic structure in the active state. The latter structure has the SH3, SH2, and SH2-catalytic domain linker assembled as a pseudo-two-domain protein. The regulation model emerging from this study, including at least three different conformational states, allows the tight regulation of the enzyme without compromising fast response in the presence of natural targets.     

Determination of Ensemble-Average Pairwise Root Mean-Square Deviation from Experimental B-Factors

Root mean-square deviation (RMSD) after roto-translational least-squares fitting is a measure of global structural similarity of macromolecules used commonly. On the other hand, experimental x-ray B-factors are used frequently to study local structural heterogeneity and dynamics in macromolecules by providing direct information about root mean-square fluctuations (RMSF) that can also be calculated from molecular dynamics simulations. We provide a mathematical derivation showing that, given a set of conservative assumptions, a root mean-square ensemble-average of an all-against-all distribution of pairwise RMSD for a single molecular species, <RMSD²>^1/2, is directly related to average B-factors (<B>) and <RMSF²>^1/2. We show this relationship and explore its limits of validity on a heterogeneous ensemble of structures taken from molecular dynamics simulations of villin headpiece generated using distributed-computing techniques and the Folding@Home cluster. Our results provide a basis for quantifying global structural diversity of macromolecules in crystals directly from x-ray experiments, and we show this on a large set of structures taken from the Protein Data Bank. In particular, we show that the ensemble-average pairwise backbone RMSD for a microscopic ensemble underlying a typical protein x-ray structure is ∼1.1 Å, under the assumption that the principal contribution to experimental B-factors is conformational variability.     

Effective handling of induced-fit motion in flexible docking

For structure-based drug design, where various ligand structures need to be docked to a target protein structure, a docking method that can handle conformational flexibility of not only the ligand, but also the protein, is indispensable. We have developed a simple and effective approach for dealing with the local induced-fit motion of the target protein, and implemented it in our docking tool, ADAM. Our approach efficiently combines the following two strategies: a vdW-offset grid in which the protein cavity is enlarged uniformly, and structure optimization allowing the motion of ligand and protein atoms. To examine the effectiveness of our approach, we performed docking validation studies, including redocking in 18 test cases and foreign-docking, in which various ligands from foreign crystal structures of complexes are docked into a target protein structure, in 22 cases (on five target proteins). With the original ADAM, the correct docking modes (RMSD < 2.0 A) were not present among the top 20 models in one case of redocking and four cases of foreign-docking. When the handling of induced-fit motion was implemented, the correct solutions were acquired in all 40 test cases. In foreign-docking on thymidine kinase, the correct docking modes were obtained as the top-ranked solutions for all 10 test ligands by our combinatorial approach, and this appears to be the best result ever reported with any docking tool. The results of docking validation have thus confirmed the effectiveness of our approach, which can provide reliable docking models even in the case of foreign-docking, where conformational change of the target protein cannot be ignored. We expect that this approach will contribute substantially to actual drug design, including virtual screening.