Aline Thomas

Analysis of domain motions in large proteins

We present a new approach for determining dynamical domains in large proteins, either based on a comparison of different experimental structures, or on a simplified normal mode calculation for a single conformation. In a first step, a deformation measure is evaluated for all residues in the protein; a high deformation indicates highly flexible interdomain regions. The sufficiently rigid parts of the protein are then classified into rigid domains and low‐deformation interdomain regions on the basis of their global motion. We demonstrate the techniques on three proteins: citrate synthase, which has been the subject of earlier domain analyses, HIV‐1 reverse transcriptase, which has a rather complex domain structure, and aspartate transcarbamylase as an example of a very large protein. These examples show that the comparison of conformations and the normal mode analysis lead to essentially the same domain identification, except for cases where the experimental conformations differ by the presence of a large ligand, such as a DNA strand. Normal mode analysis has the advantage of requiring only one experimental structure and of providing a more detailed picture of domain movements, e.g. the splitting of domains into subdomains at higher frequencies. Proteins 1999;34:369–382.

The dynamo library for molecular simulations using hybrid quantum mechanical and molecular mechanical potentials

The Dynamo module library has been developed for the simulation of molecular systems using hybrid quantum mechanical (QM) and molecular mechanical (MM) potentials. Dynamo is not a program package but is a library of Fortran 90 modules that can be employed by those interested in writing their own programs for performing molecular simulations. The library supports a range of different types of molecular calculation including geometry optimizations, reaction‐path determinations and molecular dynamics and Monte Carlo simulations. This article outlines the general structure and capabilities of the library and describes in detail Dynamos semiempirical QM/MM hybrid potential. Results are presented to indicate three particular aspects of this implementation—the handling of long‐range nonbonding interactions, the nature of the boundary between the quantum mechanical and molecular mechanical atoms and how to perform path‐integral hybrid‐potential molecular dynamics simulations.

Tertiary and quaternary conformational changes in aspartate transcarbamylase: a normal mode study

Aspartate transcarbamylase (ATCase) initiates the pyrimidine biosynthetic pathway in Escherichia coli. Binding of aspartate to this allosteric enzyme induces a cooperative transition between the tensed (T) and relaxed (R) states of the enzyme which involves large quaternary and tertiary rearrangements. The mechanisms of the transmission of the regulatory signal to the active site (60 Å away) and that of the cooperative transition are not known in detail, although a large number of single, double, and triple site‐specific mutants and chimeric forms of ATCase have been obtained and kinetically characterized. A previous analysis of the very low‐frequency normal modes of both the T and R state structures of ATCase identified some of the large‐amplitude motions mediating the intertrimer elongation and rotation that occur during the cooperative transition (Thomas et al., J. Mol. Biol. 257:1070–1087, 1996; Thomas et al., J. Mol. Biol. 261:490–506, 1996). As a complement to that study, the deformation of the quaternary and tertiary structure of ATCase by normal modes below 5 cm−1 is investigated in this article. The ability of the modes to reproduce the domain motions occurring during the transition is analyzed, with special attention to the interdomain closure in the catalytic chain, which has been shown to be critical for homotropic cooperativity. The calculations show a coupling between the quaternary motions and more localized motions involving specific residues. The particular dynamic behavior of these residues is examined in the light of biochemical results to obtain insights into their role in the transmission of the allosteric signal. Proteins 1999;34:96–112.

Langerin-Heparin Interaction: Two Binding Sites for Small and Large Ligands as revealed by a combination of NMR Spectroscopy and Cross-Linking Mapping Experiments

Langerin is a C-type lectin present on Langerhans cells that mediates capture of pathogens in a carbohydrate-dependent manner, leading to subsequent internalization and elimination in the cellular organelles called Birbeck granules. This mechanism mediated by langerin was shown to constitute a natural barrier for HIV-1 particle transmission. Besides interacting specifically with high mannose and fucosylated neutral carbohydrate structures, langerin has the ability to bind sulfated carbohydrate ligands as 6-sulfated galactosides in the Ca(2+)-dependent binding site. Very recently langerin was demonstrated to interact with sulfated glycosaminoglycans (GAGs), in a Ca(2+)-independent way, resulting in the proposal of a new binding site for GAGs. On the basis of those results, we have conducted a structural study of the interactions of small heparin (HEP)-like oligosaccharides with langerin in solution. Heparin bead cross-linking experiments, an approach specifically designed to identify HEP/heparan sulfate binding sites in proteins were first carried out and experimentally validated the previously proposed model for the interaction of langerin extracellular domain with 6 kDa HEP. High-resolution NMR studies of a set of eight synthetic HEP-like trisaccharides harboring different sulfation patterns demonstrated that all of them bound to langerin in a Ca(2+)-dependent way. The binding epitopes were determined by saturation transfer difference NMR and the bound conformations by transferred NOESY experiments. These experimental data were combined with docking and molecular dynamics and resulted in the proposal of a binding mode characterized by the coordination of calcium by the two equatorial hydroxyl groups, OH3 and OH4, at the non-reducing end. The binding also includes the carboxylate group at the adjacent iduronate residue. This epitope is shared by all eight ligands, explaining the absence of any impact on binding from differences in their substitution patterns. Finally, in contrast to the small trisaccharides, we demonstrated that a longer HEP-like hexasaccharide, bearing an additional O-sulfate group at the non-reducing end, which precludes binding to the Ca(2+) site, interacts with langerin in the previously identified Ca(2+)-independent binding site.

Insights into the Mechanism by Which Interferon-γ Basic Amino Acid Clusters Mediate Protein Binding to Heparan Sulfate

The extensive functional repertoire of heparin and heparan sulfate, which relies on their ability to interact with a large number of proteins, has recently emerged. To understand the forces that drive such interactions the binding of heparin to interferon-γ (IFNγ), used as a model system, was investigated. NMR-based titration experiments demonstrated the involvement of two adjacent cationic domains (D1: KTGKRKR and D2: RGRR), both of which are present within the carboxy-terminal sequence of the cytokine. Kinetic analysis showed that these two domains contribute differently to the interaction: D1 is required to form a complex and constitutes the actual binding site, whereas D2, although unable to associate with heparin by itself, increased the association rate of the binding. These data are consistent with the view that D2, through nonspecific electrostatic forces, places the two molecules in favorable orientations for productive binding within the encounter complex. This mechanism was supported by electrostatic potential analysis and thermodynamic investigations. They showed that D1 association to heparin is driven by both favorable enthalpic and entropic contributions, as expected for a binding sequence, but that D2 gives rise to entropic penalty, which opposes binding in a thermodynamic sense. The binding mechanism described herein, by which the D2 domain kinetically drives the interaction, has important functional consequences and gives a structural framework to better understand how specific are the interactions between proteins and heparin.

Deciphering the Glycan Preference of Bacterial Lectins by Glycan Array and Molecular Docking with Validation by Microcalorimetry and Crystallography

Recent advances in glycobiology revealed the essential role of lectins for deciphering the glycocode by specific recognition of carbohydrates. Integrated multiscale approaches are needed for characterizing lectin specificity: combining on one hand high-throughput analysis by glycan array experiments and systematic molecular docking of oligosaccharide libraries and on the other hand detailed analysis of the lectin/oligosaccharide interaction by x-ray crystallography, microcalorimetry and free energy calculations. The lectins LecB from Pseudomonas aeruginosa and BambL from Burkholderia ambifaria are part of the virulence factors used by the pathogenic bacteria to invade the targeted host. These two lectins are not related but both recognize fucosylated oligosaccharides such as the histo-blood group oligosaccharides of the ABH(O) and Lewis epitopes. The specificities were characterized using semi-quantitative data from glycan array and analyzed by molecular docking with the Glide software. Reliable prediction of protein/oligosaccharide structures could be obtained as validated by existing crystal structures of complexes. Additionally, the crystal structure of BambL/Lewis x was determined at 1.6 Å resolution, which confirms that Lewis x has to adopt a high-energy conformation so as to bind to this lectin. Free energies of binding were calculated using a procedure combining the Glide docking protocol followed by free energy rescoring with the Prime/Molecular Mechanics Generalized Born Surface Area (MM-GBSA) method. The calculated data were in reasonable agreement with experimental free energies of binding obtained by titration microcalorimetry. The established predictive protocol is proposed to rationalize large sets of data such as glycan arrays and to help in lead discovery projects based on such high throughput technology.

Reduction of Lectin Valency Drastically Changes Glycolipid Dynamics in Membranes, But not Surface Avidity.

Multivalency is proposed to play a role in the strong avidity of lectins for glycosylated cell surfaces and also in their ability to affect membrane dynamics by clustering glycosphingolipids. Lectins with modified valency were designed from the β-propeller fold of Ralstonia solanacearum lectin (RSL) that presents six fucose binding sites. After identification of key amino acids by molecular dynamics calculations, two mutants with reduced valency were produced. Isothermal titration calorimetry confirmed the loss of three high affinity binding sites for both mutants. Crystal structures indicated that residual low affinity binding occurred in W76A but not in R17A. The trivalent R17A mutant presented unchanged avidity toward fucosylated surfaces, when compared to hexavalent RSL. However, R17A is not able anymore to induce formation of membrane invaginations on giant unilamellar vesicules, indicating the crucial role of number of binding sites for clustering of glycolipids. In the human lung epithelial cell line H1299, wt-RSL is internalized within seconds whereas the kinetics of R17A uptake is largely delayed. Neolectins with tailored valency are promising tools to study membrane dynamics.

Probing the acceptor active site organization of the human recombinant β1,4-galactosyltransferase 7 and design of xyloside-based inhibitors

Background: Glycosyltransferase inhibitors have important applications in therapeutics and as chemical biology tools. Results: The human β1,4-galactosyltransferase 7 enzyme active site was mapped by modeling, mutagenesis, and in vitro/in cellulo assays, and novel inhibitors were synthesized. Conclusion: An efficient inhibitor of β1,4-galactosyltransferase 7 and glycosaminoglycan synthesis was obtained. Significance: This inhibitory molecule can be exploited to investigate glycosaminoglycan biology and modulate glycosaminoglycan synthesis in therapeutics. Among glycosaminoglycan (GAG) biosynthetic enzymes, the human β1,4-galactosyltransferase 7 (hβ4GalT7) is characterized by its unique capacity to take over xyloside derivatives linked to a hydrophobic aglycone as substrates and/or inhibitors. This glycosyltransferase is thus a prime target for the development of regulators of GAG synthesis in therapeutics. Here, we report the structure-guided design of hβ4GalT7 inhibitors. By combining molecular modeling, in vitro mutagenesis, and kinetic measurements, and in cellulo analysis of GAG anabolism and decorin glycosylation, we mapped the organization of the acceptor binding pocket, in complex with 4-methylumbelliferone-xylopyranoside as prototype substrate. We show that its organization is governed, on one side, by three tyrosine residues, Tyr194, Tyr196, and Tyr199, which create a hydrophobic environment and provide stacking interactions with both xylopyranoside and aglycone rings. On the opposite side, a hydrogen-bond network is established between the charged amino acids Asp228, Asp229, and Arg226, and the hydroxyl groups of xylose. We identified two key structural features, i.e. the strategic position of Tyr194 forming stacking interactions with the aglycone, and the hydrogen bond between the His195 nitrogen backbone and the carbonyl group of the coumarinyl molecule to develop a tight binder of hβ4GalT7. This led to the synthesis of 4-deoxy-4-fluoroxylose linked to 4-methylumbelliferone that inhibited hβ4GalT7 activity in vitro with a Ki 10 times lower than the Km value and efficiently impaired GAG synthesis in a cell assay. This study provides a valuable probe for the investigation of GAG biology and opens avenues toward the development of bioactive compounds to correct GAG synthesis disorders implicated in different types of malignancies.

A comparative QM/MM simulation study of the reaction mechanisms of human and Plasmodium falciparum HG(X)PRTases.

QM/MM hybrid potential free-energy simulations are performed to compare the reaction mechanisms of human hypoxanthine guanine phosphoribosyl transferase (HGPRTase) and the corresponding enyzme from Plasmodium falciparum (Pf), hypoxanthine guanine xanthine phosphoribosyl transferase (HGXPRTase). These enzymes share 44% of sequence identity but display very different affinities for xanthine. The calculations show that in both enzymes phosphoribosyl transfer proceeds via a dissociative mechanism from an anionic form of the substrate. Nevertheless, there are significant differences in the geometries of critical structures along the reaction paths which it may be possible to exploit for the design of specific inhibitors against the Pf enzyme.