Docking Stage Research Articles

Replica Exchange Simulations of the Thermodynamics of Aβ Fibril Growth

Replica exchange molecular dynamics and an all-atom implicit solvent model are used to probe the thermodynamics of deposition of Alzheimer's Aβ monomers on preformed amyloid fibrils. Consistent with the experiments, two deposition stages have been identified. The docking stage occurs over a wide temperature range, starting with the formation of the first peptide-fibril interactions at 500 K. Docking is completed when a peptide fully adsorbs on the fibril edge at the temperature of 380 K. The docking transition appears to be continuous, and occurs without free energy barriers or intermediates. During docking, incoming Aβ monomer adopts a disordered structure on the fibril edge. The locking stage occurs at the temperature of ≈360 K and is characterized by the rugged free energy landscape. Locking takes place when incoming Aβ peptide forms a parallel β-sheet structure on the fibril edge. Because the β-sheets formed by locked Aβ peptides are typically off-registry, the structure of the locked phase differs from the structure of the fibril interior. The study also reports that binding affinities of two distinct fibril edges with respect to incoming Aβ peptides are different. The peptides bound to the concave edge have significantly lower free energy compared to those bound on the convex edge. Comparison with the available experimental data is discussed.

Principles of flexible protein–protein docking

Treating flexibility in molecular docking is a major challenge in cell biology research. Here we describe the background and the principles of existing flexible protein-protein docking methods, focusing on the algorithms and their rational. We describe how protein flexibility is treated in different stages of the docking process: in the preprocessing stage, rigid and flexible parts are identified and their possible conformations are modeled. This preprocessing provides information for the subsequent docking and refinement stages. In the docking stage, an ensemble of pre-generated conformations or the identified rigid domains may be docked separately. In the refinement stage, small-scale movements of the backbone and side-chains are modeled and the binding orientation is improved by rigid-body adjustments. For clarity of presentation, we divide the different methods into categories. This should allow the reader to focus on the most suitable method for a particular docking problem.

A Flexible Approach to Induced Fit Docking

We present Fleksy, a new approach to consider both ligand and receptor flexibility in small molecule docking. Pivotal to our method is the use of a receptor ensemble to describe protein flexibility. To construct these ensembles, we use a backbone-dependent rotamer library and implement the concept of interaction sampling. The latter allows the evaluation of different orientations of ambivalent interaction partners. The docking stage consists of an ensemble-based soft-docking experiment using FlexX-Ensemble, followed by an effective flexible receptor-ligand complex optimization using Yasara. Fleksy produces a set of receptor-ligand complexes ranked using a consensus scoring function combining docking scores and force field energies. Averaged over three cross-docking datasets, containing 35 different receptor-ligand complexes in total, Fleksy reproduces the observed binding mode within 2.0 A for 78% of the complexes. This compares favorably to the rigid receptor FlexX program, which on average reaches a success rate of 44% for these datasets.

Assays of vacuole fusion resolve the stages of docking, lipid mixing, and content mixing

Membrane fusion entails organelle docking and subsequent mixing of membrane bilayers and luminal compartments. We now present an in vitro assay of fusion, using yeast vacuoles bearing domains of either Fos or Jun fused to complementary halves of beta-lactamase. Upon fusion, these proteins associate to yield beta-lactamase activity. This assay complements the standard fusion assay (activation of pro-Pho8p in protease-deficient vacuoles by proteases from pho8Delta vacuoles). Both the beta-lactamase and pro-Pho8p activation assays of fusion show the same long kinetic delay between SNARE pairing and luminal compartment mixing. Lipid-mixing occurs rapidly after SNARE pairing but well before aqueous compartment mixing. These results support a model in which SNARE pairing leads to rapid hemifusion, followed by slow further lipid rearrangement and aqueous compartment mixing.

PyDock: Electrostatics and desolvation for effective scoring of rigid‐body protein–protein docking

The accurate scoring of rigid-body docking orientations represents one of the major difficulties in protein-protein docking prediction. Other challenges are the development of faster and more efficient sampling methods and the introduction of receptor and ligand flexibility during simulations. Overall, good discrimination of near-native docking poses from the very early stages of rigid-body protein docking is essential step before applying more costly interface refinement to the correct docking solutions. Here we explore a simple approach to scoring of rigid-body docking poses, which has been implemented in a program called pyDock. The scheme is based on Coulombic electrostatics with distance dependent dielectric constant, and implicit desolvation energy with atomic solvation parameters previously adjusted for rigid-body protein-protein docking. This scoring function is not highly dependent on specific geometry of the docking poses and therefore can be used in rigid-body docking sets generated by a variety of methods. We have tested the procedure in a large benchmark set of 80 unbound docking cases. The method is able to detect a near-native solution from 12,000 docking poses and place it within the 100 lowest-energy docking solutions in 56% of the cases, in a completely unrestricted manner and without any other additional information. More specifically, a near-native solution will lie within the top 20 solutions in 37% of the cases. The simplicity of the approach allows for a better understanding of the physical principles behind protein-protein association, and provides a fast tool for the evaluation of large sets of rigid-body docking poses in search of the near-native orientation.

The Calcium-binding C-terminal Domain of Escherichia coli α-Hemolysin Is a Major Determinant in the Surface-active Properties of the Protein

alpha-Hemolysin (HlyA) from Escherichia coli is a protein toxin (1024 amino acids) that targets eukaryotic cell membranes, causing loss of the permeability barrier. HlyA consists of two main regions, an N-terminal domain rich in amphipathic helices, and a C-terminal Ca(2+)-binding domain containing a Gly- and Asp-rich nonapeptide repeated in tandem 11-17 times. The latter is called the RTX domain and gives its name to the RTX protein family. It had been commonly assumed that membrane interaction occurred mainly if not exclusively through the amphipathic helix domain. However, we have cloned and expressed the C-terminal region of HlyA, containing the RTX domain plus a few stabilizing sequences, and found that it is a potent surface-active molecule. The isolated domain binds Ca(2+) with about the same affinity (apparent K(0.5) approximately 150 microM) as the parent protein HlyA, and Ca(2+) binding induces in turn a more compact folding with an increased proportion of beta-sheet structure. Both with and without Ca(2+) the C-terminal region of HlyA can interact with lipid monolayers spread at an air-water interface. However, the C-terminal domain by itself is devoid of membrane lytic properties. The present results can be interpreted in the light of our previous studies that involved in receptor binding a peptide in the C-terminal region of HlyA. We had also shown experimentally the distinction between reversible membrane adsorption and irreversible lytic insertion of the toxin. In this context, the present data allow us to propose that both major domains of HlyA are directly involved in membrane-toxin interaction, the nonapeptide repeat, calcium-binding RTX domain being responsible for the early stages of HlyA docking to the target membrane.

Docking Ligands into Flexible and Solvated Macromolecules. 1. Development and Validation of FITTED 1.0

We report the development and validation of a novel suite of programs, FITTED 1.0, for the docking of flexible ligands into flexible proteins. This docking tool is unique in that it can deal with both the flexibility of macromolecules (side chains and main chains) and the presence of bridging water molecules while treating protein/ligand complexes as realistically dynamic systems. This software relies on a genetic algorithm to account for the flexibility of the two molecules as well as the location of bridging water molecules. In addition, FITTED 1.0 features a novel application of a switching function to retain or displace key water molecules from the protein-ligand complexes. Two independent modules, ProCESS and SMART, were developed to set up the proteins and the ligands prior to the docking stage. Validation of the accuracy of the software was achieved via the application of FITTED 1.0 to the docking of inhibitors of HIV-1 protease, thymidine kinase, trypsin, factor Xa, and MMP to their respective proteins.

Dissecting Multiple Steps of GLUT4 Trafficking and Identifying the Sites of Insulin Action

Insulin-stimulated GLUT4 translocation is central to glucose homeostasis. Functional assays to distinguish individual steps in the GLUT4 translocation process are lacking, thus limiting progress toward elucidation of the underlying molecular mechanism. Here we have developed a robust method, which relies on dynamic tracking of single GLUT4 storage vesicles (GSVs) in real time, for dissecting and systematically analyzing the docking, priming, and fusion steps of GSVs with the cell surface in vivo. Using this method, we have shown that the preparation of GSVs for fusion competence after docking at the surface is a key step regulated by insulin, whereas the docking step is regulated by PI3K and its downstream effector, the Rab GAP AS160. These data show that Akt-dependent phosphorylation of AS160 is not the major regulated step in GLUT4 trafficking, implicating alternative Akt substrates or alternative signaling pathways downstream of GSV docking at the cell surface as the major regulatory node.

Ion Regulation of Homotypic Vacuole Fusion in Saccharomyces cerevisiae

Biological membrane fusion employs divalent cations as protein cofactors or as signaling ligands. For example, Mg2+ is a cofactor for the N-ethylmaleimide-sensitive factor (NSF) ATPase, and the Ca2+ signal from neuronal membrane depolarization is required for synaptotagmin activation. Divalent cations also regulate liposome fusion, but the role of such ion interactions with lipid bilayers in Rab- and soluble NSF attachment protein receptor (SNARE)-dependent biological membrane fusion is less clear. Yeast vacuole fusion requires Mg2+ for Sec18p ATPase activity, and vacuole docking triggers an efflux of luminal Ca2+. We now report distinct reaction conditions where divalent or monovalent ions interchangeably regulate Rab- and SNARE-dependent vacuole fusion. In reactions with 5 mm Mg2+, other free divalent ions are not needed. Reactions containing low Mg2+ concentrations are strongly inhibited by the rapid Ca2+ chelator BAPTA. However, addition of the soluble SNARE Vam7p relieves BAPTA inhibition as effectively as Ca2+ or Mg2+, suggesting that Ca2+ does not perform a unique signaling function. When the need for Mg2+, ATP, and Sec18p for fusion is bypassed through the addition of Vam7p, vacuole fusion does not require any appreciable free divalent cations and can even be stimulated by their chelators. The similarity of these findings to those with liposomes, and the higher ion specificity of the regulation of proteins, suggests a working model in which ion interactions with bilayer lipids permit Rab- and SNARE-dependent membrane fusion.

Identification of Novel Parasitic Cysteine Protease Inhibitors Using Virtual Screening. 1. The ChemBridge Database

Trypanosomiasis, leishmaniasis, and malaria are major parasitic diseases in developing countries. The existing chemotherapy of these diseases suffers from lack of safe and effective drugs and/or the presence of widespread drug resistance. Cysteine proteases are exciting novel targets for antiparasitic drug design. Virtual screening was performed in an attempt to identify novel druglike nonpeptide inhibitors of parasitic cysteine proteases. The ChemBridge database consisting of approximately 241 000 compounds was screened against homology models of falcipain-2 and falcipain-3 in three consecutive stages of docking. A total of 24 diverse inhibitors were identified from an initial group of 84, of which 12 compounds appeared to be dual inhibitors of falcipain-2 and falcipain-3. Four compounds showed inhibition of both the malarial cysteine proteases as well as Leishmania donovani cysteine protease.

The Cytoplasmic Domain of Vamp4 and Vamp5 Is Responsible for Their Correct Subcellular Targeting: THE N-TERMINAL EXTENSION OF VAMP4 CONTAINS A DOMINANT AUTONOMOUS TARGETING SIGNAL FOR THE TRANS-GOLGI NETWORK

SNAREs represent a superfamily of proteins responsible for the last stage of docking and subsequent fusion in diverse intracellular membrane transport events. The Vamp subfamily of SNAREs contains 7 members (Vamp1, Vamp2, Vamp3/cellubrevin, Vamp4, Vamp5, Vamp7/Ti-Vamp, and Vamp8/endobrevin) that are distributed in various post-Golgi structures. Vamp4 and Vamp5 are distributed predominantly in the trans-Golgi network (TGN) and the plasma membrane, respectively. When C-terminally tagged with enhanced green fluorescent protein, the majority of Vamp4 and Vamp5 is correctly targeted to the TGN and plasma membrane, respectively. Swapping the N-terminal cytoplasmic region and the C-terminal membrane anchor domain between Vamp4 and Vamp5 demonstrates that the N-terminal cytoplasmic region of these two SNAREs contains the correct subcellular targeting information. As compared with Vamp5, Vamp4 contains an N-terminal extension of 51 residues. Appending this 51-residue N-terminal extension onto the N terminus of Vamp5 results in targeting of the chimeric protein to the TGN, suggesting that this N-terminal extension of Vamp4 contains a dominant and autonomous targeting signal for the TGN. Analysis of deletion mutants of this N-terminal region suggests that this TGN-targeting signal is encompassed within a smaller region consisting of a di-Leu motif followed by two acidic clusters. The essential role of the di-Leu motif and the second acidic cluster was then established by site-directed mutagenesis.

Rho1p and Cdc42p act after Ypt7p to regulate vacuole docking.

Rho GTPases, which control polarized cell growth through cytoskeletal reorganization, have recently been implicated in the control of endo- and exocytosis. We now report that both Rho1p and Cdc42p have a direct role in mediating the docking stage of homotypic vacuole fusion. Vacuoles prepared from strains with temperature-sensitive alleles of either Rho1p or Cdc42p are thermolabile for fusion. RhoGDI (Rdi1p), which extracts Rho1p and Cdc42p from the vacuole membrane, blocks vacuole fusion. The Rho GTPases can not fulfill their function as long as priming and Ypt7p-dependent tethering are inhibited. However, reactions that are reversibly blocked after docking by the calcium chelator BAPTA have passed the point of sensitivity to Rdi1p. Extraction and removal of Ypt7p, Rho1p and Cdc42p from docked vacuoles (by Gdi1p, Gyp7p and Rdi1p) does not impede subsequent membrane fusion, which is still sensitive to GTPgammaS. Thus, multiple GTPases act in a defined sequence to regulate the docking steps of vacuole fusion.

Cdc42p functions at the docking stage of yeast vacuole membrane fusion.

Membrane fusion reactions have been considered to be primarily regulated by Rab GTPases. In the model system of homotypic vacuole fusion in the yeast Saccharomyces cerevisiae, we show that Cdc42p, a member of the Rho family of GTPases, has a direct role in membrane fusion. Genetic evidence suggested a relationship between Cdc42p and Vtc1p/Nrf1p, a central part of the vacuolar membrane fusion machinery. Vacuoles from cdc42 temperature-sensitive mutants are deficient for fusion at the restrictive temperature. Specific amino acid changes on the Cdc42p protein surface in these mutants define the putative interaction domain that is crucial for its function in membrane fusion. Affinity-purified antibodies to this domain inhibited the in vitro fusion reaction. Using these antibodies in kinetic analyses and assays for subreactions of the priming, docking and post-docking phase of the reaction, we show that Cdc42p action follows Ypt7p-dependent tethering, but precedes the formation of trans-SNARE complexes. Thus, our data define an effector binding domain of Cdc42p by which it regulates the docking reaction of vacuole fusion.

Dsl1p, an essential protein required for membrane traffic at the endoplasmic reticulum/Golgi interface in yeast.

To identify novel factors required for ER to Golgi transport in yeast we performed a screen for genes that, when mutated, confer a dependence on a dominant mutant form of the ER to Golgi vesicle docking factor Sly1p, termed Sly1-20p. DSL1, a novel gene isolated in the screen, encodes an essential protein with a predicted molecular mass of 88 kDa. DSL1 is required for transport between the ER and the Golgi because strains bearing mutant alleles of this gene accumulate the pre-Golgi form of transported proteins at the restrictive temperature. Two strains bearing temperature-sensitive alleles of DSL1 display distinct phenotypes as observed by electron microscopy at the restrictive temperature; although both strains accumulate ER membrane, one also accumulates vesicles. Interestingly, the inviability of strains bearing several mutant alleles of DSL1 can be suppressed by expression of either Erv14p (a protein required for the movement of a specific protein from the ER to the Golgi), Sec21p (the gamma-subunit of the COPI coat protein complex coatomer), or Sly1-20p. Because the strongest suppressor is SEC21, we proposed that Dsl1p functions primarily in retrograde Golgi to ER traffic, although it is possible that Dsl1p functions in anterograde traffic as well, perhaps at the docking stage, as suggested by the suppression by SLY1-20.

A new role for a SNARE protein as a regulator of the Ypt7/Rab-dependent stage of docking.

The homotypic fusion of yeast vacuoles occurs in an ordered cascade of priming, docking, and fusion. The linkage between these steps has so far remained unclear. We now report that Vam7p (the vacuolar SNAP-23/25 homolog) signals from the cis-SNARE complex to Ypt7p (the vacuolar Rab/Ypt) to initiate the docking process. After Vam7p has been released from the cis-SNARE complex by Sec18p-mediated priming, it is still required for Ypt7p-dependent docking and it needs Ypt7p to remain on the vacuole. Thus, after priming, Vam7p is released from the vacuole altogether if Ypt7p has been extracted by Gdi1p or inactivated by antibody but is not released if docking is blocked simply by vacuole dilution; it is therefore Ypt7p function, and not docking per se, that retains Vam7p. In accord with this finding, cells deleted for the gene encoding Ypt7 have normal amounts of Vam7p but have little Vam7p on their isolated vacuoles. Interaction of Vam7p and Ypt7p is further indicated by two-hybrid analysis [Uetz, P., Giot, L., Cagney, G., Mansfield, T. A., Judson, R. S., Knight, J. R., Lockshon, D., Narayan, V., Srinivasan, M., Pochart, P., et al. (2000) Nature (London) 403, 623-627] and by the effect of Vam7p on the association of the Rab/Ypt-effector HOPS complex (homotypic fusion and vacuole protein sorting; Vam2p and Vam6p plus four vacuole protein sorting class C proteins) with Ypt7p. Vam7p provides a functional link between the priming step, which releases it from the cis-SNARE complex, and docking.

Conflict resolution in AGV systems

Currently, conflict-free routing in AGV systems is established by means of one of the following three approaches: (i) the problem elimination through the adoption of a segmented path flow or tandem queue configuration; (ii) the identification of imminent collisions through forward sensing and their aversion through vehicle backtracking and/or rerouting; or (iii) the imposition of zone control and extensive route pre-planning, typically based on deterministic timing of the vehicle traveling and docking stages. Among these three approaches, the segmented path flow-based approach presents the highest robustness to the system stochas-ticities/randomness, but at the cost of restricted vehicle routings and the need for complicated handling operations. This paper proposes an alternative conflict resolution strategy that will ensure robust AGV conflict resolution, while maintaining the operational flexibility provided by free vehicle travel on arbitrarily structured guidepath networks. Specifically, the approach advocated in this paper also employs zone control, but it determines vehicle routes incrementally, one zone at a time. Routing decisions are the result of a sequence of safety and performance considerations, with the former being primarily based on structural/logical rather than timing aspects of the system behavior. The resulting control problem is characterized as the AGV structural control. After defining the notion of AGV structural control, the paper proceeds to the formal characterization and analysis of the problem, and to the development of a structural control policy appropriate for the class of AGV resource allocation systems. The paper concludes with some discussion on the accommodation of emerging AGV operational features in the proposed modeling and analysis framework, and the integration of AGV structural control with the broader control of material-flow among the shop-floor workstations.

The Docking of Primed Vacuoles Can Be Reversibly Arrested by Excess Sec17p (α-SNAP)

Homotypic vacuole fusion occurs in ordered stages of priming, docking, and fusion. Priming, which prepares vacuoles for productive association, requires Sec17p (the yeast homolog of alpha-SNAP), Sec18p (the yeast NSF, an ATP-driven chaperone), and ATP. Sec17p is initially an integral part of the cis-SNARE complex together with vacuolar SNARE proteins and Sec18p (NSF). Previous studies have shown that Sec17p is rapidly released from the vacuole membrane during priming as the cis-SNARE complex is disassembled, but the order and causal relationship of these subreactions has not been known. We now report that the addition of excess recombinant his(6)-Sec17p to primed vacuoles can block subsequent docking. This inhibition is reversible by Sec18p, but the reaction cannot proceed to the tethering and trans-SNARE pairing steps of docking while the Sec17p block is in place. Once docking has occurred, excess Sec17p does not inhibit membrane fusion per se. Incubation of cells with thermosensitive Sec17-1p at nonpermissive temperature causes SNARE complex disassembly. These data suggest that Sec17p can stabilize vacuolar cis-SNARE complexes and that the release of Sec17p by Sec18p and ATP allows disassembly of this complex and activates its components for docking.

The docking stage of yeast vacuole fusion requires the transfer of proteins from a cis-SNARE complex to a Rab/Ypt protein.

The homotypic fusion of yeast vacuoles requires Sec18p (NSF)-driven priming to allow vacuole docking, but the mechanism that links priming and docking is unknown. We find that a large multisubunit protein called the Vam2/6p complex is bound to cis-paired SNAP receptors (SNAREs) on isolated vacuoles. This association of the Vam2/6p complex with the cis-SNARE complex is disrupted during priming. The Vam2/6p complex then binds to Ypt7p, a guanosine triphosphate binding protein of the Rab family, to initiate productive contact between vacuoles. Thus, cis-SNARE complexes can contain Rab/Ypt effectors, and these effectors can be mobilized by NSF/Sec18p-driven priming, allowing their direct association with a Rab/Ypt protein to activate docking.

Conflict resolution in AGV systems

Currently, conflict-free routing in AGV systems is established by means of one of the following three approaches: (i) the problem elimination through the adoption of a segmented path flow or tandem queue configuration; (ii) the identification of imminent collisions through forward sensing and their aversion through vehicle backtracking and/or rerouting; or (iii) the imposition of zone control and extensive route pre-planning, typically based on deterministic timing of the vehicle traveling and docking stages. Among these three approaches, the segmented path flow-based approach presents the highest robustness to the system stochasticities/randomness, but at the cost of restricted vehicle routings and the need for complicated handling operations. This paper proposes an alternative conflict resolution strategy that will ensure robust AGV conflict resolution, while maintaining the operational flexibility provided by free vehicle travel on arbitrarily structured guidepath networks. Specifically, the approach advocated in this paper also employs zone control, but it determines vehicle routes incrementally, one zone at a time. Routing decisions are the result of a sequence of safety and performance considerations, with the former being primarily based on structural/logical rather than timing aspects of the system behavior. The resulting control problem is characterized as the AGV structural control. After defining the notion of AGV structural control, the paper proceeds to the formal characterization and analysis of the problem, and to the development of a structural control policy appropriate for the class of AGV resource allocation systems. The paper concludes with some discussion on the accommodation of emerging AGV operational features in the proposed modeling and analysis framework, and the integration of AGV structural control with the broader control of material-flow among the shop-floor workstations.

Two-Stage Method for Protein−Ligand Docking

A two-stage method for the computational prediction of the structure of protein-ligand complexes is proposed. Given an experimentally determined structure of the protein, in the first stage a large number of plausible ligand conformations is generated using the fast docking algorithm FlexX. In the second stage these conformations are minimized and reranked using a method based on a classical force field. The two-stage method is tested for 10 different protein-ligand complexes. For 9 of them experimentally determined structures are known. It turns out that the two-stage method strongly improves the predictive power as compared to that of the fast docking stage alone. The tenth case is a bona fide prediction of a complex of thrombin with a new inhibitor for which no experimentally determined structure is available so far.