High-resolution structure prediction of a circular permutation loop

Methods for rapid and reliable design and structure prediction of linker loops would facilitate a variety of protein engineering applications. Circular permutation, in which the existing termini of a protein are linked by the polypeptide chain and new termini are created, is one such application that has been employed for decreasing proteolytic susceptibility and other functional purposes. The length and sequence of the linker can impact the expression level, solubility, structure and function of the permuted variants. Hence it is desirable to achieve atomic‐level accuracy in linker design. Here, we describe the use of RosettaRemodel for design and structure prediction of circular permutation linkers on a model protein. A crystal structure of one of the permuted variants confirmed the accuracy of the computational prediction, where the all‐atom rmsd of the linker region was 0.89 Å between the model and the crystal structure. This result suggests that RosettaRemodel may be generally useful for the design and structure prediction of protein loop regions for circular permutations or other structure‐function manipulations.     

MollDE: a homology modeling framework you can click with

SUMMARY: Molecular Integrated Development Environment (MolIDE) is an integrated application designed to provide homology modeling tools and protocols under a uniform, user-friendly graphical interface. Its main purpose is to combine the most frequent modeling steps in a semi-automatic, interactive way, guiding the user from the target protein sequence to the final three-dimensional protein structure. The typical basic homology modeling process is composed of building sequence profiles of the target sequence family, secondary structure prediction, sequence alignment with PDB structures, assisted alignment editing, side-chain prediction and loop building. All of these steps are available through a graphical user interface. MolIDE's user-friendly and streamlined interactive modeling protocol allows the user to focus on the important modeling questions, hiding from the user the raw data generation and conversion steps. MolIDE was designed from the ground up as an open-source, cross-platform, extensible framework. This allows developers to integrate additional third-party programs to MolIDE. AVAILABILITY: http://dunbrack.fccc.edu/molide/molide.php CONTACT: rl_dunbrack@fccc.edu.     

Structure of yeast triosephosphate isomerase at 1.9-A resolution

The structure of yeast triosephosphate isomerase (TIM) has been solved at 3.0-A resolution and refined at 1.9-A resolution to an R factor of 21.0%. The final model consists of all non-hydrogen atoms in the polypeptide chain and 119 water molecules, a number of which are found in the interior of the protein. The structure of the active site clearly indicates that the carboxylate of the catalytic base, Glu 165, is involved in a hydrogen-bonding interaction with the hydroxyl of Ser 96. In addition, the interactions of the other active site residues, Lys 12 and His 95, are also discussed. For the first time in any TIM structure, the "flexible loop" has well-defined density; the conformation of the loop in this structure is stabilized by a crystal contact. Analysis of the subunit interface of this dimeric enzyme hints at the source of the specificity of one subunit for another and allows us to estimate an association constant of 10(14)-10(16) M-1 for the two monomers. The analysis also suggests that the interface may be a particularly good target for drug design. The conserved positions (20%) among sequences from 13 sources ranging on the evolutionary scale from Escherichia coli to humans reveal the intense pressure to maintain the active site structure.     

A computational method for the design of nested proteins by loop‐directed domain insertion

The computational design of novel nested proteins—in which the primary structure of one protein domain (insert) is flanked by the primary structure segments of another (parent)—would enable the generation of multifunctional proteins. Here we present a new algorithm, called Loop‐Directed Domain Insertion (LooDo), implemented within the Rosetta software suite, for the purpose of designing nested protein domain combinations connected by flexible linker regions. Conformational space for the insert domain is sampled using large libraries of linker fragments for linker‐to‐parent domain superimposition followed by insert‐to‐linker superimposition. The relative positioning of the two domains (treated as rigid bodies) is sampled efficiently by a grid‐based, mutual placement compatibility search. The conformations of the loop residues, and the identities of loop as well as interface residues, are simultaneously optimized using a generalized kinematic loop closure algorithm and Rosetta EnzymeDesign, respectively, to minimize interface energy. The algorithm was found to consistently sample near‐native conformations and interface sequences for a benchmark set of structurally similar but functionally divergent domain‐inserted enzymes from the α/β hydrolase superfamily, and discriminates well between native and nonnative conformations and sequences, although loop conformations tended to deviate from the native conformations. Furthermore, in cross‐domain placement tests, native insert‐parent domain combinations were ranked as the best‐scoring structures compared to nonnative domain combinations. This algorithm should be broadly applicable to the design of multi‐domain protein complexes with any combination of inserted or tandem domain connections.     

Atomic-Accuracy Prediction of Protein Loop Structures through an RNA-Inspired Ansatz

Consistently predicting biopolymer structure at atomic resolution from sequence alone remains a difficult problem, even for small sub-segments of large proteins. Such loop prediction challenges, which arise frequently in comparative modeling and protein design, can become intractable as loop lengths exceed 10 residues and if surrounding side-chain conformations are erased. Current approaches, such as the protein local optimization protocol or kinematic inversion closure (KIC) Monte Carlo, involve stages that coarse-grain proteins, simplifying modeling but precluding a systematic search of all-atom configurations. This article introduces an alternative modeling strategy based on a ‘stepwise ansatz’, recently developed for RNA modeling, which posits that any realistic all-atom molecular conformation can be built up by residue-by-residue stepwise enumeration. When harnessed to a dynamic-programming-like recursion in the Rosetta framework, the resulting stepwise assembly (SWA) protocol enables enumerative sampling of a 12 residue loop at a significant but achievable cost of thousands of CPU-hours. In a previously established benchmark, SWA recovers crystallographic conformations with sub-Angstrom accuracy for 19 of 20 loops, compared to 14 of 20 by KIC modeling with a comparable expenditure of computational power. Furthermore, SWA gives high accuracy results on an additional set of 15 loops highlighted in the biological literature for their irregularity or unusual length. Successes include cis-Pro touch turns, loops that pass through tunnels of other side-chains, and loops of lengths up to 24 residues. Remaining problem cases are traced to inaccuracies in the Rosetta all-atom energy function. In five additional blind tests, SWA achieves sub-Angstrom accuracy models, including the first such success in a protein/RNA binding interface, the YbxF/kink-turn interaction in the fourth ‘RNA-puzzle’ competition. These results establish all-atom enumeration as an unusually systematic approach to ab initio protein structure modeling that can leverage high performance computing and physically realistic energy functions to more consistently achieve atomic accuracy.     

Modeling protein loops using a phi i + 1, psi i dimer database.

We present an automated method for modeling backbones of protein loops. The method samples a database of phi i + 1 and psi i angles constructed from a nonredundant version of the Protein Data Bank (PDB). The dihedral angles phi i + 1 and psi i completely define the backbone conformation of a dimer when standard bond lengths, bond angles, and a trans planar peptide configuration are used. For the 400 possible dimers resulting from 20 natural amino acids, a list of allowed phi i + 1, psi i pairs for each dimer is created by pooling all such pairs from the loop segments of each protein in the nonredundant version of the PDB. Starting from the N-terminus of the loop sequence, conformations are generated by assigning randomly selected pairs of phi i + 1, psi i for each dimer from the respective pool using standard bond lengths, bond angles, and a trans peptide configuration. We use this database to simulate protein loops of lengths varying from 5 to 11 amino acids in five proteins of known three-dimensional structures. Typically, 10,000-50,000 models are simulated for each protein loop and are evaluated for stereochemical consistency. Depending on the length and sequence of a given loop, 50-80% of the models generated have no stereochemical strain in the backbone atoms. We demonstrate that, when simulated loops are extended to include flanking residues from homologous segments, only very few loops from an ensemble of sterically allowed conformations orient the flanking segments consistent with the protein topology. The presence of near-native backbone conformations for loops from five different proteins suggests the completeness of the dimeric database for use in modeling loops of homologous proteins. Here, we take advantage of this observation to design a method that filters near-native loop conformations from an ensemble of sterically allowed conformations. We demonstrate that our method eliminates the need for a loop-closure algorithm and hence allows for the use of topological constraints of the homologous proteins or disulfide constraints to filter near-native loop conformations.     

A Unified Conformational Selection and Induced Fit Approach to Protein-Peptide Docking

Protein-peptide interactions are vital for the cell. They mediate, inhibit or serve as structural components in nearly 40% of all macromolecular interactions, and are often associated with diseases, making them interesting leads for protein drug design. In recent years, large-scale technologies have enabled exhaustive studies on the peptide recognition preferences for a number of peptide-binding domain families. Yet, the paucity of data regarding their molecular binding mechanisms together with their inherent flexibility makes the structural prediction of protein-peptide interactions very challenging. This leaves flexible docking as one of the few amenable computational techniques to model these complexes. We present here an ensemble, flexible protein-peptide docking protocol that combines conformational selection and induced fit mechanisms. Starting from an ensemble of three peptide conformations (extended, a-helix, polyproline-II), flexible docking with HADDOCK generates 79.4% of high quality models for bound/unbound and 69.4% for unbound/unbound docking when tested against the largest protein-peptide complexes benchmark dataset available to date. Conformational selection at the rigid-body docking stage successfully recovers the most relevant conformation for a given protein-peptide complex and the subsequent flexible refinement further improves the interface by up to 4.5 Å interface RMSD. Cluster-based scoring of the models results in a selection of near-native solutions in the top three for ∼75% of the successfully predicted cases. This unified conformational selection and induced fit approach to protein-peptide docking should open the route to the modeling of challenging systems such as disorder-order transitions taking place upon binding, significantly expanding the applicability limit of biomolecular interaction modeling by docking.     

Multiple copy sampling in protein loop modeling: computational efficiency and sensitivity to dihedral angle perturbations.

Multiple copy sampling and the bond scaling-relaxation technique are combined to generate 3-dimensional conformations of protein loop segments. The computational efficiency and sensitivity to initial loop copy dispersion are analyzed. The multicopy loop modeling method requires approximately 20-50% of the computational time required by the single-copy method for the various protein segments tested. An analytical formula is proposed to estimate the computational gain prior to carrying out a multicopy simulation. When 7-residue loops within flexible proteins are modeled, each multicopy simulation can sample a set of loop conformations with initial dispersions up to +/- 15 degrees for backbone and +/- 30 degrees for side-chain rotatable dihedral angles. The dispersions are larger for shorter and smaller for longer and/or surface loops. The degree of convergence of loop copies during a simulation can be used to complement commonly used target functions (such as potential energy) for distinguishing between native and misfolded conformations. Furthermore, this convergence also reflects the conformational flexibility of the modeled protein segment. Application to simultaneously building all 6 hypervariable loops of an antibody is discussed.     

Structural analysis of lipid complexes of GM2-activator protein

The GM2-activator protein (GM2-AP) is a small lysosomal lipid transfer protein essential for the hydrolytic conversion of ganglioside GM2 to GM3 by beta-hexosaminidase A. The crystal structure of human apo-GM2-AP is known to consist of a novel beta-cup fold with a spacious hydrophobic interior. Here, we present two new structures of GM2-AP with bound lipids, showing two different lipid-binding modes within the apolar pocket. The 1.9A structure with GM2 bound shows the position of the ceramide tail and significant conformational differences among the three molecular copies in the asymmetric unit. The tetrasaccharide head group is not visible and is presumed to be disordered. However, its general position could be established through modeling. The structure of a low-pH crystal, determined at 2.5A resolution, has a significantly enlarged hydrophobic channel that merges with the apolar pocket. Electron density inside the pocket and channel suggests the presence of a trapped phospholipid molecule. Structure alignments among the four crystallographically unique monomers provide information on the potential role for lipid binding of flexible chain segments at the rim of the cavity opening. Two discrete orientations of the S130-T133 loop define an open and a closed configuration of the hydrophobic channel that merges with the apolar pocket. We propose: (i) that the low-pH structure represents an active membrane-binding conformation; (ii) that the mobile S130-T133 loop serves as a gate for passage of ligand into the apolar pocket; and (iii) that this loop and the adjacent apolar V59-W63 loop form a surface patch with two exposed tryptophan residues that could interface with lipid bilayers.     

Crystal structure of reovirus attachment protein sigma1 reveals evolutionary relationship to adenovirus fiber.

Reovirus attaches to cellular receptors with the sigma1 protein, a fiber-like molecule protruding from the 12 vertices of the icosahedral virion. The crystal structure of a receptor-binding fragment of sigma1 reveals an elongated trimer with two domains: a compact head with a new beta-barrel fold and a fibrous tail containing a triple beta-spiral. Numerous structural and functional similarities between reovirus sigma1 and the adenovirus fiber suggest an evolutionary link in the receptor-binding strategies of these two viruses. A prominent loop in the sigma1 head contains a cluster of residues that are conserved among reovirus serotypes and are likely to form a binding site for junction adhesion molecule, an integral tight junction protein that serves as a reovirus receptor. The fibrous tail is mainly responsible for sigma1 trimer formation, and it contains a highly flexible region that allows for significant movement between the base of the tail and the head. The architecture of the trimer interface and the observed flexibility indicate that sigma1 is a metastable structure poised to undergo conformational changes upon viral attachment and cell entry.     

ModLoop: automated modeling of loops in protein structures

SUMMARY: ModLoop is a web server for automated modeling of loops in protein structures. The input is the atomic coordinates of the protein structure in the Protein Data Bank format, and the specification of the starting and ending residues of one or more segments to be modeled, containing no more than 20 residues in total. The output is the coordinates of the non-hydrogen atoms in the modeled segments. A user provides the input to the server via a simple web interface, and receives the output by e-mail. The server relies on the loop modeling routine in MODELLER that predicts the loop conformations by satisfaction of spatial restraints, without relying on a database of known protein structures. For a rapid response, ModLoop runs on a cluster of Linux PC computers. AVAILABILITY: The server is freely accessible to academic users at http://salilab.org/modloop     

Accounting for loop flexibility during protein-protein docking

Although reliable docking can now be achieved for systems that do not undergo important induced conformational change upon association, the presence of flexible surface loops, which must adapt to the steric and electrostatic properties of a partner, generally presents a major obstacle. We report here the first docking method that allows large loop movements during a systematic exploration of the possible arrangements of the two partners in terms of position and rotation. Our strategy consists in taking into account an ensemble of possible loop conformations by a multi-copy representation within a reduced protein model. The docking process starts from regularly distributed positions and orientations of the ligand around the whole receptor. Each starting configuration is submitted to energy minimization during which the best-fitting loop conformation is selected based on the mean-field theory. Trials were carried out on proteins with significant differences in the main-chain conformation of the binding loop between isolated form and complexed form, which were docked to their partner considered in their bound form. The method is able to predict complexes very close to the crystal complex both in terms of relative position of the two partners and of the geometry of the flexible loop. We also show that introducing loop flexibility on the isolated protein form during systematic docking largely improves the predictions of relative position of the partners in comparison with rigid-body docking.     

RosettaScripts: a scripting language interface to the Rosetta macromolecular modeling suite

Macromolecular modeling and design are increasingly useful in basic research, biotechnology, and teaching. However, the absence of a user-friendly modeling framework that provides access to a wide range of modeling capabilities is hampering the wider adoption of computational methods by non-experts. RosettaScripts is an XML-like language for specifying modeling tasks in the Rosetta framework. RosettaScripts provides access to protocol-level functionalities, such as rigid-body docking and sequence redesign, and allows fast testing and deployment of complex protocols without need for modifying or recompiling the underlying C++ code. We illustrate these capabilities with RosettaScripts protocols for the stabilization of proteins, the generation of computationally constrained libraries for experimental selection of higher-affinity binding proteins, loop remodeling, small-molecule ligand docking, design of ligand-binding proteins, and specificity redesign in DNA-binding proteins.     

Anchored design of protein-protein interfaces

Background

Few existing protein-protein interface design methods allow for extensive backbone rearrangements during the design process. There is also a dichotomy between redesign methods, which take advantage of the native interface, and de novo methods, which produce novel binders.

Methodology

Here, we propose a new method for designing novel protein reagents that combines advantages of redesign and de novo methods and allows for extensive backbone motion. This method requires a bound structure of a target and one of its natural binding partners. A key interaction in this interface, the anchor, is computationally grafted out of the partner and into a surface loop on the design scaffold. The design scaffold''s surface is then redesigned with backbone flexibility to create a new binding partner for the target. Careful choice of a scaffold will bring experimentally desirable characteristics into the new complex. The use of an anchor both expedites the design process and ensures that binding proceeds against a known location on the target. The use of surface loops on the scaffold allows for flexible-backbone redesign to properly search conformational space.

Conclusions and Significance

This protocol was implemented within the Rosetta3 software suite. To demonstrate and evaluate this protocol, we have developed a benchmarking set of structures from the PDB with loop-mediated interfaces. This protocol can recover the correct loop-mediated interface in 15 out of 16 tested structures, using only a single residue as an anchor.     

Efficient algorithms to explore conformation spaces of flexible protein loops

Several applications in biology - e.g., incorporation of protein flexibility in ligand docking algorithms, interpretation of fuzzy X-ray crystallographic data, and homology modeling - require computing the internal parameters of a flexible fragment (usually, a loop) of a protein in order to connect its termini to the rest of the protein without causing any steric clash. One must often sample many such conformations in order to explore and adequately represent the conformational range of the studied loop. While sampling must be fast, it is made difficult by the fact that two conflicting constraints - kinematic closure and clash avoidance - must be satisfied concurrently. This paper describes two efficient and complementary sampling algorithms to explore the space of closed clash-free conformations of a flexible protein loop. The "seed sampling" algorithm samples broadly from this space, while the "deformation sampling" algorithm uses seed conformations as starting points to explore the conformation space around them at a finer grain. Computational results are presented for various loops ranging from 5 to 25 residues. More specific results also show that the combination of the sampling algorithms with a functional site prediction software (FEATURE) makes it possible to compute and recognize calcium-binding loop conformations. The sampling algorithms are implemented in a toolkit (LoopTK), which is available at https://simtk.org/home/looptk.     

The FEATURE framework for protein function annotation: modeling new functions,improving performance,and extending to novel applications

Structural genomics efforts contribute new protein structures that often lack significant sequence and fold similarity to known proteins. Traditional sequence and structure-based methods may not be sufficient to annotate the molecular functions of these structures. Techniques that combine structural and functional modeling can be valuable for functional annotation. FEATURE is a flexible framework for modeling and recognition of functional sites in macromolecular structures. Here, we present an overview of the main components of the FEATURE framework, and describe the recent developments in its use. These include automating training sets selection to increase functional coverage, coupling FEATURE to structural diversity generating methods such as molecular dynamics simulations and loop modeling methods to improve performance, and using FEATURE in large-scale modeling and structure determination efforts.     

Engineered extrahelical base destabilization enhances sequence discrimination of DNA methyltransferase M.HhaI

Improved sequence specificity of the DNA cytosine methyltransferase HhaI was achieved by disrupting interactions at a hydrophobic interface between the active site of the enzyme and a highly conserved flexible loop. Transient fluorescence experiments show that mutations disrupting this interface destabilize the positioning of the extrahelical, "flipped" cytosine base within the active site. The ternary crystal structure of the F124A M.HhaI bound to cognate DNA and the cofactor analogue S-adenosyl-l-homocysteine shows an increase in cavity volume between the flexible loop and the core of the enzyme. This cavity disrupts the interface between the loop and the active site, thereby destabilizing the extrahelical target base. The favored partitioning of the base-flipped enzyme-DNA complex back to the base-stacked intermediate results in the mutant enzyme discriminating better than the wild-type enzyme against non-cognate sites. Building upon the concepts of kinetic proofreading and our understanding of M.HhaI, we describe how a 16-fold specificity enhancement achieved with a double mutation at the loop/active site interface is acquired through destabilization of intermediates prior to methyltransfer rather than disruption of direct interactions between the enzyme and the substrate for M.HhaI.     

Elucidation of the ATP7B N-domain Mg2+-ATP coordination site and its allosteric regulation

The diagnostic of orphan genetic disease is often a puzzling task as less attention is paid to the elucidation of the pathophysiology of these rare disorders at the molecular level. We present here a multidisciplinary approach using molecular modeling tools and surface plasmonic resonance to study the function of the ATP7B protein, which is impaired in the Wilson disease. Experimentally validated in silico models allow the elucidation in the Nucleotide binding domain (N-domain) of the Mg(2+)-ATP coordination site and answer to the controversial role of the Mg(2+) ion in the nucleotide binding process. The analysis of protein motions revealed a substantial effect on a long flexible loop branched to the N-domain protein core. We demonstrated the capacity of the loop to disrupt the interaction between Mg(2+)-ATP complex and the N-domain and propose a role for this loop in the allosteric regulation of the nucleotide binding process.     

Zinc ion coordination as a modulating factor of the ZnuA histidine-rich loop flexibility: A molecular modeling and fluorescence spectroscopy study

ZnuA is the soluble component of the high-affinity ZnuABC zinc transporter belonging to the ATP-binding cassette-type periplasmic Zn-binding proteins. The zinc transporter ZnuABC is composed by three proteins: ZnuB, the membrane permease, ZnuC, the ATPase component and ZnuA, the soluble periplasmic metal-binding protein which captures Zn and delivers it to ZnuB.The ZnuA protein contains a charged flexible loop, rich in histidines and acidic residues, showing significant species-specific differences. Various studies have established that this loop contributes to the formation of a secondary zinc binding site, which has been proposed to be important in the acquisition of periplasmic Zn for its delivery to ZnuB or for regulation of zinc uptake. Due to its high mobility the structure of the histidine-rich loop has never been solved by X-ray diffraction studies. In this paper, through a combined use of molecular modeling, mutagenesis and fluorescence spectroscopy, we confirm the presence of two zinc binding sites characterized by different affinities for the metal ion and show that the flexibility of the loop is modulated by the binding of the zinc ions to the protein. The data obtained by fluorescence spectroscopy have then be used to validate a 3D model including the unsolved histidine-rich loop.