Biomedical subjects

Alexandre M J J Bonvin

Publications and source records attributed to Alexandre M J J Bonvin.

At least 19 recordsLinked to original sources

PubMed2006

Activity-structure correlations in divergent lectin evolution: fine specificity of chicken galectin CG-14 and computational analysis of flexible ligand docking for CG-14 and the closely related CG-16.

Gene duplication and sequence divergence are driving forces toward establishing protein families. To examine how sequence changes affect carbohydrate specificity, the two closely related proto-type chicken galectins CG-14 and CG-16 were selected as models. Binding properties were analyzed using a highly sensitive solid-phase assay. We tested 56 free saccharides and 34 well-defined glycoproteins. The two galectins share preference for the II (Galbeta1-4GlcNAc) versus I (Galbeta1-3GlcNAc) version of beta-galactosides. A pronounced difference is found owing to the reactivity of CG-14 with histo-blood group ABH active oligosaccharides and A/B active glycoproteins. These experimental results prompted to determine activity-structure correlations by modeling. Computational analysis included consideration of the flexibility of binding partners and the presence of water molecules. It provided a comparative description of complete carbohydrate recognition domains, which had so far not been characterized in animal galectins. The structural models assigned II, I selectivity to a region downstream of the central Trp moiety. Docking revealed that the tetrasaccharides can be accommodated in their free-state low-energy conformations. CG-14's preference for A versus B epitopes could be attributed to a contact between His124 and the N-acetyl group of GalNAc. Regarding intergalectin comparison, the Ala53/Cys51 exchange affects the interaction potential of His54/His52. Close inspection of simulated dynamic interplay revealed reorientation of His124 at the site of the His124/Glu123 substitution, with potential impact on ligand dissociation. In summary, this study identifies activity differences and provides information on their relation to structural divergence, epitomizing the value of this combined approach beyond galectins.

Amino Acid Sequence↗

PubMed2006

Solvated docking: introducing water into the modelling of biomolecular complexes.

MOTIVATION: Interfacial water, which plays an important role in mediating biomolecular interactions, has been neglected in the modelling of biomolecular complexes. METHODS: We present a solvated docking approach that explicitly accounts for the presence of water in protein-protein complexes. Our solvated docking protocol is based on the concept of the first encounter complex in which a water layer is present in-between the molecules. It mimics the pathway from this initial complex towards the final assembly in which most waters have been expelled from the interface. Docking is performed from solvated biomolecules and waters are removed in a biased Monte Carlo procedure based on water-mediated contact propensities obtained from an analysis of high-resolution crystal structures. RESULTS: We demonstrate the feasibility of this approach for protein-protein complexes representing both 'wet' and 'dry' interfaces. Solvated docking leads to improvements both in quality and scoring. Water molecules are recovered that closely match the ones in the crystal structures. AVAILABILTY: Solvated docking will be made available in the future release of HADDOCK version 2.0 (http://www.nmr.chem.uu.nl/haddock).

Binding Sites↗

PubMed2006

Binding site structure of one LRP-RAP complex: implications for a common ligand-receptor binding motif.

The low-density lipoprotein receptor-related protein (LRP) interacts with more than 30 ligands of different sizes and structures that can all be replaced by the receptor-associated protein (RAP). The double module of complement type repeats, CR56, of LRP binds many ligands including all three domains of RAP and alpha2-macroglobulin, which promotes the catabolism of the Abeta-peptide implicated in Alzheimer's disease. To understand the receptor-ligand cross-talk, the NMR structure of CR56 has been solved and ligand binding experiments with RAP domain 1 (RAPd1) have been performed. From chemical shift perturbations of both binding partners upon complex formation, a HADDOCK model of the complex between CR56 and RAPd1 has been obtained. The binding residues are similar to a common binding motif suggested from alpha2-macroglobulin binding studies and provide evidence for an understanding of their mutual cross-competition pattern. The present structural results convey a simultaneous description of both binding partners of an LRP-ligand complex and open a route to a broader understanding of the binding specificity of the LRP receptor, which may involve a general four-residue receptor-ligand recognition motif common to all LRP ligands. The present result may be beneficial in the design of antagonists of ligand binding to the LDL receptor family, and especially of drugs for treatment of Alzheimer's disease.

Amino Acid Motifs↗

PubMed2006

Comparative NMR study on the impact of point mutations on protein stability of Pseudomonas mendocina lipase.

In this work we compare the dynamics and conformational stability of Pseudomonas mendocina lipase enzyme and its F180P/S205G mutant that shows higher activity and stability for use in washing powders. Our NMR analyses indicate virtually identical structures but reveal remarkable differences in local dynamics, with striking correspondence between experimental data (i.e., (15)N relaxation and H/D exchange rates) and data from Molecular Dynamics simulations. While overall the cores of both proteins are very rigid on the pico- to nanosecond timescale and are largely protected from H/D exchange, the two point mutations stabilize helices alpha1, alpha4, and alpha5 and locally destabilize the H-bond network of the beta-sheet (beta7-beta9). In particular, it emerges that helix alpha5, undergoing some fast destabilizing motions (on the pico- to nanosecond timescale) in wild-type lipase, is substantially rigidified by the mutation of Phe180 for a proline at its N terminus. This observation could be explained by the release of some penalizing strain, as proline does not require any "N-capping" hydrogen bond acceptor in the i+3 position. The combined experimental and simulated data thus indicate that reduced molecular flexibility of the F180P/S205G mutant lipase underlies its increased stability, and thus reveals a correlation between microscopic dynamics and macroscopic thermodynamic properties. This could contribute to the observed altered enzyme activity, as may be inferred from recent studies linking enzyme kinetics to their local molecular dynamics.

Amino Acid Sequence↗

PubMed2006

Information-driven protein-DNA docking using HADDOCK: it is a matter of flexibility.

Intrinsic flexibility of DNA has hampered the development of efficient protein-DNA docking methods. In this study we extend HADDOCK (High Ambiguity Driven DOCKing) [C. Dominguez, R. Boelens and A. M. J. J. Bonvin (2003) J. Am. Chem. Soc. 125, 1731-1737] to explicitly deal with DNA flexibility. HADDOCK uses non-structural experimental data to drive the docking during a rigid-body energy minimization, and semi-flexible and water refinement stages. The latter allow for flexibility of all DNA nucleotides and the residues of the protein at the predicted interface. We evaluated our approach on the monomeric repressor-DNA complexes formed by bacteriophage 434 Cro, the Escherichia coli Lac headpiece and bacteriophage P22 Arc. Starting from unbound proteins and canonical B-DNA we correctly predict the correct spatial disposition of the complexes and the specific conformation of the DNA in the published complexes. This information is subsequently used to generate a library of pre-bent and twisted DNA structures that served as input for a second docking round. The resulting top ranking solutions exhibit high similarity to the published complexes in terms of root mean square deviations, intermolecular contacts and DNA conformation. Our two-stage docking method is thus able to successfully predict protein-DNA complexes from unbound constituents using non-structural experimental data to drive the docking.

Computational Biology↗

PubMed2006

Direct use of unassigned resonances in NMR structure calculations with proxy residues.

We present a method that significantly enhances the robustness of (automated) NMR structure determination by allowing the NOE data corresponding to unassigned NMR resonances to be used directly in the calculations. The unassigned resonances are represented by additional atoms or groups of atoms that have no interaction with the regular protein atoms except through distance restraints. These so-called "proxy" residues can be used to generate NOE-based distance restraints in a similar fashion as for the assigned part of the protein. If sufficient NOE information is available, the restraints are expected to place the proxies at positions close to the correct atoms for the unassigned resonance, which can facilitate subsequent assignment. Convergence can be further improved by supplying additional information about the possible identities of the unassigned resonances. We have implemented this approach in the widely used automated assignment and structure calculation protocols ARIA and CANDID. We find that it significantly increases the robustness of structure calculations with regard to missing assignments and yields structures of higher quality. Our approach is still able to find correctly folded structures with up to 30% randomly missing resonance assignments, and even when only backbone and beta resonances are present! This should be of significant value to NMR-based structural proteomics initiatives.

Journal Article↗

PubMed2006

Intramolecular surface contacts contain information about protein-protein interface regions.

MOTIVATION: Some amino acids clearly show preferences over others in protein-protein interfaces. These preferences, or so-called interface propensities can be used for a priori interface prediction. We investigated whether the prediction accuracy could be improved by considering not single but pairs of residues in an interface. Here we present the first systematic analysis of intramolecular surface contacts in interface prediction. RESULTS: We show that preferences do exist for contacts within and around an interface region within one molecule: specific pairs of amino acids are more often occurring than others. Using intramolecular contact propensities in a blind test, higher average scores were assigned to interface residues than to non-interface residues. This effect persisted as small but significant when the contact propensities were corrected to eliminate the influence of single amino acid interface propensity. This indicates that intramolecular contact propensities may replace interface propensities in protein-protein interface prediction. AVAILABILITY: The source code is available on request from the authors.

Algorithms↗

PubMed2006

WHISCY: what information does surface conservation yield? Application to data-driven docking.

Protein-protein interactions play a key role in biological processes. Identifying the interacting residues is a first step toward understanding these interactions at a structural level. In this study, the interface prediction program WHISCY is presented. It combines surface conservation and structural information to predict protein-protein interfaces. The accuracy of the predictions is more than three times higher than a random prediction. These predictions have been combined with another interface prediction program, ProMate [Neuvirth et al. J Mol Biol 2004;338:181-199], resulting in an even more accurate predictor. The usefulness of the predictions was tested using the data-driven docking program HADDOCK [Dominguez et al. J Am Chem Soc 2003;125:1731-1737] in an unbound docking experiment, with the goal of generating as many near-native structures as possible. Unrefined rigid body docking solutions within 10 A ligand RMSD from the true structure were generated for 22 out of 25 docked complexes. For 18 complexes, more than 100 of the 8000 generated models were correct. Our results demonstrates the potential of using interface predictions to drive protein-protein docking.

Binding Sites↗

PubMed2006

Flexible protein-protein docking.

Predicting the structure of protein-protein complexes using docking approaches is a difficult problem whose major challenges include identifying correct solutions, and properly dealing with molecular flexibility and conformational changes. Flexibility can be addressed at several levels: implicitly, by smoothing the protein surfaces or allowing some degree of interpenetration (soft docking) or by performing multiple docking runs from various conformations (cross or ensemble docking); or explicitly, by allowing sidechain and/or backbone flexibility. Although significant improvements have been achieved in the modeling of sidechains, methods for the explicit inclusion of backbone flexibility in docking are still being developed. A few novel approaches have emerged involving collective degrees of motion, multicopy representations and multibody docking, which should allow larger conformational changes to be modeled.

Binding Sites↗

PubMed2006

Kinetics and thermodynamics of type VIII beta-turn formation: a CD, NMR, and microsecond explicit molecular dynamics study of the GDNP tetrapeptide.

We report an experimental and theoretical study on type VIII beta-turn using a designed peptide of sequence GDNP. CD and NMR studies reveal that this peptide exists in equilibrium between type VIII beta-turn and extended conformations. Extensive MD simulations give a description of the free energy landscape of the peptide in which we retrieve the same two main conformations suggested by the experiments. The free energy difference between the two conformational states is very small and the transition between them occurs within a few kT at 300 K on a nanosecond timescale. The equilibrium is mainly driven by entropic contribution, which favors extended conformations over beta-turns. This confirms other theoretical studies showing that beta-turns are marginally stable in water solution because of the larger entropy of the extended state unless some stabilizing interactions exist. Our observations may be extended to any type of beta-turn and have important consequences for protein folding. A comparison of our MD and CD results also suggests a possible type VIII beta-turn CD signature indicated by one main band at 200 nm, close to that of random coil, and a fairly large shoulder at 220 nm. Last, our results clearly show that the XXXP motif can only fold into a type VIII beta-turn, which is consistent with its fairly strong propensity for this type of turn. This important finding may help for peptide design and is in line with recent studies on bioactive elastin peptides.

Circular Dichroism↗

PubMed2006

The alpha-to-beta conformational transition of Alzheimer's Abeta-(1-42) peptide in aqueous media is reversible: a step by step conformational analysis suggests the location of beta conformation seeding.

Current views of the role of beta-amyloid (Abeta) peptide fibrils range from regarding them as the cause of Alzheimer's pathology to having a protective function. In the last few years, it has also been suggested that soluble oligomers might be the most important toxic species. In all cases, the study of the conformational properties of Abeta peptides in soluble form constitutes a basic approach to the design of molecules with "antiamyloid" activity. We have experimentally investigated the conformational path that can lead the Abeta-(1-42) peptide from the native state, which is represented by an alpha helix embedded in the membrane, to the final state in the amyloid fibrils, which is characterized by beta-sheet structures. The conformational steps were monitored by using CD and NMR spectroscopy in media of varying polarities. This was achieved by changing the composition of water and hexafluoroisopropanol (HFIP). In the presence of HFIP, beta conformations can be observed in solutions that have very high water content (up to 99 % water; v/v). These can be turned back to alpha helices simply by adding the appropriate amount of HFIP. The transition of Abeta-(1-42) from alpha to beta conformations occurs when the amount of water is higher than 80 % (v/v). The NMR structure solved in HFIP/H2O with high water content showed that, on going from very apolar to polar environments, the long N-terminal helix is essentially retained, whereas the shorter C-terminal helix is lost. The complete conformational path was investigated in detail with the aid of molecular-dynamics simulations in explicit solvent, which led to the localization of residues that might seed beta conformations. The structures obtained might help to find regions that are more affected by environmental conditions in vivo. This could in turn aid the design of molecules able to inhibit fibril deposition or revert oligomerization processes.

Algorithms↗

PubMed2006

The solution structure of the AppA BLUF domain: insight into the mechanism of light-induced signaling.

The transcriptional antirepressor AppA from the photosynthetic bacterium Rhodobacter sphaeroides senses both the light climate and the intracellular redox state. Under aerobic conditions in the dark, AppA binds to and thereby blocks the function of PpsR, a transcriptional repressor. Absorption of a blue photon dissociates AppA from PpsR and allows the latter to repress photosynthesis gene expression. The N terminus of AppA contains sequence homology to flavin-containing photoreceptors that belong to the BLUF family. Structural and chemical aspects of signal transduction mediated by AppA are still largely unknown. Here we present NMR studies of the N-terminal flavin-binding BLUF domain of AppA. Its solution structure adopts an alpha/beta-sandwich fold with a beta alpha beta beta alpha beta beta topology, which represents a new flavin-binding fold. Considerable disorder is observed for residues near the chromophore due to conformational exchange. This disorder is observed both in the dark and in the light-induced signaling state of AppA. Furthermore, we detect light-induced structural changes in a patch of surface residues that provide a structural link between light absorption and signal-transduction events.

Bacterial Proteins↗

PubMed2006

Combining NMR relaxation with chemical shift perturbation data to drive protein-protein docking.

The modeling of biomolecular complexes by computational docking using the known structures of their constituents is developing rapidly to become a powerful tool in structural biology. It is especially useful in combination with even limited experimental information describing the interface. Here we demonstrate for the first time the use of diffusion anisotropy in combination with chemical shift perturbation data to drive protein-protein docking. For validation purposes we make use of simulated diffusion anisotropy data. Inclusion of this information, which can be derived from NMR relaxation rates and reports on the orientation of the components of a complex with respect to the rotational diffusion tensor, substantially improves the docking results.

Anisotropy↗

PubMed2005

Model for RNA binding and the catalytic site of the RNase Kid of the bacterial parD toxin-antitoxin system.

The toxin Kid and antitoxin Kis are encoded by the parD operon of Escherichia coli plasmid R1. Kid and its chromosomal homologues MazF and ChpBK have been shown to inhibit protein synthesis in cell extracts and to act as ribosome-independent endoribonucleases in vitro. Kid cleaves RNA preferentially at the 5' side of the A residue in the nucleotide sequence 5'-UA(A/C)-3' of single-stranded regions. Here, we show that RNA cleavage by Kid yields two fragments with a 2':3'-cyclic phosphate group and a free 5'-OH group, respectively. The cleavage mechanism is similar to that of RNases A and T1, involving the uracil 2'-OH group. Via NMR titration studies with an uncleavable RNA mimic, we demonstrate that residues of both monomers of the Kid dimer together form a concatenated RNA-binding surface. Docking calculations based on the NMR chemical shifts, the cleavage mechanism and previously reported mutagenesis data provide a detailed picture of the position of the AUACA fragment within the binding pocket. We propose that residues D75, R73 and H17 form the active site of the Kid toxin, where D75 and R73 are the catalytic base and acid, respectively. The RNA sequence specificity is defined by residues T46, S47, A55, F57, T69, V71 and R73. Our data show the importance of these residues for Kid function, and the implications of our results for related toxins, such as MazF, CcdB and RelE, are discussed.

Bacterial Proteins↗

PubMed2005

Various strategies of using residual dipolar couplings in NMR-driven protein docking: application to Lys48-linked di-ubiquitin and validation against 15N-relaxation data.

When classical, Nuclear Overhauser Effect (NOE)-based approaches fail, it is possible, given high-resolution structures of the free molecules, to model the structure of a complex in solution based solely on chemical shift perturbation (CSP) data in combination with orientational restraints from residual dipolar couplings (RDCs) when available. RDCs can be incorporated into the docking following various strategies: as direct restraints and/or as intermolecular intervector projection angle restraints (Meiler et al., J Biomol NMR 2000;16:245-252). The advantage of the latter for docking is that they directly define the relative orientation of the molecules. A combined protocol in which RDCs are first introduced as intervector projection angle restraints and at a later stage as direct restraints is shown here to give the best performance. This approach, implemented in our information-driven docking approach HADDOCK (Dominguez et al., J Am Chem Soc 2003;125:1731-1737), is used to determine the solution structure of the Lys48-linked di-ubiquitin, for which chemical shift mapping, RDCs, and (15)N-relaxation data have been previously obtained (Varadan et al., J Mol Biol 2002;324:637-647). The resulting structures, derived from CSP and RDC data, are cross-validated using (15)N-relaxation data. The solution structure differs from the crystal structure by a 20 degrees rotation of the two ubiquitin units relative to each other.

Algorithms↗

PubMed2005

RECOORD: a recalculated coordinate database of 500+ proteins from the PDB using restraints from the BioMagResBank.

State-of-the-art methods based on CNS and CYANA were used to recalculate the nuclear magnetic resonance (NMR) solution structures of 500+ proteins for which coordinates and NMR restraints are available from the Protein Data Bank. Curated restraints were obtained from the BioMagResBank FRED database. Although the original NMR structures were determined by various methods, they all were recalculated by CNS and CYANA and refined subsequently by restrained molecular dynamics (CNS) in a hydrated environment. We present an extensive analysis of the results, in terms of various quality indicators generated by PROCHECK and WHAT_CHECK. On average, the quality indicators for packing and Ramachandran appearance moved one standard deviation closer to the mean of the reference database. The structural quality of the recalculated structures is discussed in relation to various parameters, including number of restraints per residue, NOE completeness and positional root mean square deviation (RMSD). Correlations between pairs of these quality indicators were generally low; for example, there is a weak correlation between the number of restraints per residue and the Ramachandran appearance according to WHAT_CHECK (r = 0.31). The set of recalculated coordinates constitutes a unified database of protein structures in which potential user- and software-dependent biases have been kept as small as possible. The database can be used by the structural biology community for further development of calculation protocols, validation tools, structure-based statistical approaches and modeling. The RECOORD database of recalculated structures is publicly available from http://www.ebi.ac.uk/msd/recoord.

Databases, Protein↗

PubMed2005

Characterization and structural analyses of nonspecific lipid transfer protein 1 from mung bean.

Plant nonspecific lipid transfer proteins (nsLTPs) are thermal stable proteins that are capable of transferring lipid molecules between bilayers in vitro. This family of proteins, abundant in plants, is proposed to be involved in defense, pollination, and germination; the in vivo biological function remains, however, elusive. Here we report the purification and sequencing of an nsLTP1 from mung bean sprouts. We have also determined the solution structure of this nsLTP1, which represents the first 3D structure of the dicotyledonous nsLTP1 family. The global fold of mung bean nsLTP1 is similar to those of the monocotyledonous nsLTP1 structures and consists of four alpha-helices stabilized by four disulfide bonds. There are, however, some notable differences in the C-terminal tails and internal hydrophobic cavities. Circular dichroism and fluorescence spectroscopy were used to compare the thermodynamics and lipid transfer properties of mung bean nsLTP1 with those of rice nsLTP1. Docking of a lipid molecule into the solution structure of mung bean nsLTP1 reveals similar binding cavities and hydrophobic interactions as in rice nsLTP1, consistent with their comparable lipid transfer properties measured experimentally.

Amino Acid Sequence↗

PubMed2005

The orientations of cytochrome c in the highly dynamic complex with cytochrome b5 visualized by NMR and docking using HADDOCK.

The interaction of bovine microsomal ferricytochrome b5 with yeast iso-1-ferri and ferrocytochrome c has been investigated using heteronuclear NMR techniques. Chemical-shift perturbations for 1H and 15N nuclei of both cytochromes, arising from the interactions with the unlabeled partner proteins, were used for mapping the interacting surfaces on both proteins. The similarity of the binding shifts observed for oxidized and reduced cytochrome c indicates that the complex formation is not influenced by the oxidation state of the cytochrome c. Protein-protein docking simulations have been performed for the binary cytochrome b5-cytochrome c and ternary (cytochrome b5)-(cytochrome c)2 complexes using a novel HADDOCK approach. The docking procedure, which makes use of the experimental data to drive the docking, identified a range of orientations assumed by the proteins in the complex. It is demonstrated that cytochrome c uses a confined surface patch for interaction with a much more extensive surface area of cytochrome b5. Taken together, the experimental data suggest the presence of a dynamic ensemble of conformations assumed by the proteins in the complex.

Animals↗