Sophie Barbe

Insights into lid movements of Burkholderia cepacia lipase inferred from molecular dynamics simulations

The interfacial activation of many lipases at water/lipid interface is mediated by large conformational changes of a so‐called lid subdomain that covers up the enzyme active site. Here we investigated using molecular dynamic simulations in different explicit solvent environments (water, octane and water/octane interface) the molecular mechanism by which the lid motion of Burkholderia cepacia lipase might operate. Although B. cepacia lipase has so far only been crystallized in open conformation, this study reveals for the first time the major conformational rearrangements that the enzyme undergoes under the influence of the solvent, which either exposes or shields the active site from the substrate. In aqueous media, the lid switches from an open to a closed conformation while the reverse motion occurs in organic environment. In particular, the role of a subdomain facing the lid on B. cepacia lipase conformational rearrangements was investigated using position‐restrained MD simulations. Our conclusions indicate that the sole mobility of α9 helix side‐chains of B. cepacia lipase is required for the full completion of the lid conformational change which is essentially driven by α5 helix movement. The role of selected α5 hydrophobic residues on the lid movement was further examined. In silico mutations of two residues, V138 and F142, were shown to drastically modify the conformational behavior of B. cepacia lipase. Overall, our results provide valuable insight into the role played by the surrounding environment on the lid conformational rearrangement and the activation of B. cepacia lipase. Proteins 2009.

Control of Lipase Enantioselectivity by Engineering the Substrate Binding Site and Access Channel

Lipase from Burkholderia cepacia (BCL) has proven to be a very useful biocatalyst for the resolution of 2‐substituted racemic acid derivatives, which are important chiral building blocks. Our previous work showed that enantioselectivity of the wild‐type BCL could be improved by chemical engineering of the substrates molecular structure. From this earlier study, three amino acids (L17, V266 and L287) were proposed as targets for mutagenesis aimed at tailoring enzyme enantioselectivity. In the present work, a small library of 57 BCL single mutants targeted on these three residues was constructed and screened for enantioselectivity towards (R,S)‐2‐chloro ethyl 2‐bromophenylacetate. This led to the fast isolation of three single mutants with a remarkable tenfold enhanced or reversed enantioselectivity. Analysis of substrate docking and access trajectories in the active site was then performed. From this analysis, the construction of 13 double mutants was proposed. Among them, an outstanding improved mutant of BCL was isolated that showed an E value of 178 and a 15‐fold enhanced specific activity compared to the parental enzyme; thus, this study demonstrates the efficiency of the semirational engineering strategy.

Applying pairwise combinations of amino Acid mutations for sorting out highly efficient glucosylation tools for chemo-enzymatic synthesis of bacterial oligosaccharides.

Iterative saturation mutagenesis and combinatorial active site saturation focused on vicinal amino acids were used to alter the acceptor specificity of amylosucrase from Neisseria polysaccharea , a sucrose-utilizing α-transglucosidase, and sort out improved variants. From the screening of three semirational sublibraries accounting in total for 20,000 variants, we report here the isolation of three double mutants of N. polysaccharea amylosucrase displaying a spectacular specificity enhancement toward both sucrose, the donor substrate, and the allyl 2-acetamido-2-deoxy-α-D-glucopyranoside acceptor as compared to the wild-type enzyme. Such levels of activity improvement have never been reported before for this class of carbohydrate-active enzymes. X-ray structure of the best performing enzymes supported by molecular dynamics simulations showed local rigidity of the -1 subsite as well as flexibility of loops involved in active site topology, which both account for the enhanced catalytic performances of the mutants. The study well illustrates the importance of taking into account the local conformation of catalytic residues as well as protein dynamics during the catalytic process, when designing enzyme libraries.

Crystal Structures of Bacterial Peptidoglycan Amidase Ampd and an Unprecedented Activation Mechanism.

AmpD is a cytoplasmic peptidoglycan (PG) amidase involved in bacterial cell-wall recycling and in induction of β-lactamase, a key enzyme of β-lactam antibiotic resistance. AmpD belongs to the amidase_2 family that includes zinc-dependent amidases and the peptidoglycan-recognition proteins (PGRPs), highly conserved pattern-recognition molecules of the immune system. Crystal structures of Citrobacter freundii AmpD were solved in this study for the apoenzyme, for the holoenzyme at two different pH values, and for the complex with the reaction products, providing insights into the PG recognition and the catalytic process. These structures are significantly different compared with the previously reported NMR structure for the same protein. The NMR structure does not possess an accessible active site and shows the protein in what is proposed herein as an inactive “closed” conformation. The transition of the protein from this inactive conformation to the active “open” conformation, as seen in the x-ray structures, was studied by targeted molecular dynamics simulations, which revealed large conformational rearrangements (as much as 17 Å) in four specific regions representing one-third of the entire protein. It is proposed that the large conformational change that would take the inactive NMR structure to the active x-ray structure represents an unprecedented mechanism for activation of AmpD. Analysis is presented to argue that this activation mechanism might be representative of a regulatory process for other intracellular members of the bacterial amidase_2 family of enzymes.

A new framework for computational protein design through cost function network optimization

MOTIVATION The main challenge for structure-based computational protein design (CPD) remains the combinatorial nature of the search space. Even in its simplest fixed-backbone formulation, CPD encompasses a computationally difficult NP-hard problem that prevents the exact exploration of complex systems defining large sequence-conformation spaces. RESULTS We present here a CPD framework, based on cost function network (CFN) solving, a recent exact combinatorial optimization technique, to efficiently handle highly complex combinatorial spaces encountered in various protein design problems. We show that the CFN-based approach is able to solve optimality a variety of complex designs that could often not be solved using a usual CPD-dedicated tool or state-of-the-art exact operations research tools. Beyond the identification of the optimal solution, the global minimum-energy conformation, the CFN-based method is also able to quickly enumerate large ensembles of suboptimal solutions of interest to rationally build experimental enzyme mutant libraries. AVAILABILITY The combined pipeline used to generate energetic models (based on a patched version of the open source solver Osprey 2.0), the conversion to CFN models (based on Perl scripts) and CFN solving (based on the open source solver toulbar2) are all available at http://genoweb.toulouse.inra.fr/~tschiex/CPD

Computational protein design as an optimization problem

Proteins are chains of simple molecules called amino acids. The three-dimensional shape of a protein and its amino acid composition define its biological function. Over millions of years, living organisms have evolved a large catalog of proteins. By exploring the space of possible amino acid sequences, protein engineering aims at similarly designing tailored proteins with specific desirable properties. In Computational Protein Design (CPD), the challenge of identifying a protein that performs a given task is defined as the combinatorial optimization of a complex energy function over amino acid sequences. In this paper, we introduce the CPD problem and some of the main approaches that have been used by structural biologists to solve it, with an emphasis on the exact method embodied in the dead-end elimination/A? algorithm (DEE/A?). The CPD problem is a specific form of binary Cost Function Network (CFN, aka Weighted CSP). We show how DEE algorithms can be incorporated and suitably modified to be maintained during search, at reasonable computational cost. We then evaluate the efficiency of CFN algorithms as implemented in our solver toulbar2, on a set of real CPD instances built in collaboration with structural biologists. The CPD problem can be easily reduced to 0/1 Linear Programming, 0/1 Quadratic Programming, 0/1 Quadratic Optimization, Weighted Partial MaxSAT and Graphical Model optimization problems. We compare toulbar2 with these different approaches using a variety of solvers. We observe tremendous differences in the difficulty that each approach has on these instances. Overall, the CFN approach shows the best efficiency on these problems, improving by several orders of magnitude against the exact DEE/A? approach. The introduction of dead-end elimination before or during search allows to further improve these results.

Structural Investigation of the Thermostability and Product Specificity of Amylosucrase from the Bacterium Deinococcus geothermalis

Background: Amylosucrases (AS) hold great potential for glycodiversification. Results: The first three-dimensional structure of AS from Deinococcus geothermalis solved here revealed an unusual dimer organization. Structures of complex of AS with turanose were also determined. Conclusion: Dimerization may contribute to thermostability. Turanose versus trehalulose formation is controlled by residues from subsite +1. Significance: This study improves the comprehension of AS properties and provides new insight for AS design. Amylosucrases are sucrose-utilizing α-transglucosidases that naturally catalyze the synthesis of α-glucans, linked exclusively through α1,4-linkages. Side products and in particular sucrose isomers such as turanose and trehalulose are also produced by these enzymes. Here, we report the first structural and biophysical characterization of the most thermostable amylosucrase identified so far, the amylosucrase from Deinoccocus geothermalis (DgAS). The three-dimensional structure revealed a homodimeric quaternary organization, never reported before for other amylosucrases. A sequence signature of dimerization was identified from the analysis of the dimer interface and sequence alignments. By rigidifying the DgAS structure, the quaternary organization is likely to participate in the enhanced thermal stability of the protein. Amylosucrase specificity with respect to sucrose isomer formation (turanose or trehalulose) was also investigated. We report the first structures of the amylosucrases from Deinococcus geothermalis and Neisseria polysaccharea in complex with turanose. In the amylosucrase from N. polysaccharea (NpAS), key residues were found to force the fructosyl moiety to bind in an open state with the O3′ ideally positioned to explain the preferential formation of turanose by NpAS. Such residues are either not present or not similarly placed in DgAS. As a consequence, DgAS binds the furanoid tautomers of fructose through a weak network of interactions to enable turanose formation. Such topology at subsite +1 is likely favoring other possible fructose binding modes in agreement with the higher amount of trehalulose formed by DgAS. Our findings help to understand the inter-relationships between amylosucrase structure, flexibility, function, and stability and provide new insight for amylosucrase design.

CAZyme discovery and design for sweet dreams

Development of synthetic routes to complex carbohydrates and glyco-conjugates is often hampered by the lack of enzymes with requisite properties or specificities. Indeed, assembly or degradation of carbohydrates requires carbohydrate-active enzymes (CAZymes) able to act on a vast range of glycosidic monomers, oligomers or polymers in a regio-specific or stereo-specific manner in order to produce the desired structure. Sequence-based analyses allow finding the most original enzymes. Novel screening methods have emerged that enable a more efficient exploitation of the CAZyme diversity found in the microbial world or generated by protein engineering. Computational biology methods also play a prominent role in the success of CAZyme design. Such progress allows circumventing current limitations of carbohydrate synthesis and opens new opportunities related to the synthetic biology field.

Enantioselectivity of Candida rugosa lipases (Lip1, Lip3, and Lip4) towards 2-bromo phenylacetic acid octyl esters controlled by a single amino acid

Enantiomer discrimination by enzymes is a very accurate mechanism, which often involves few amino acids located at the active site. Lipase isoforms from Candida rugosa are very good enzymatic models to study this phenomenon as they display high sequence homology (>80%) and their enantioselectivity is often pointed out. In the present work, we investigated three lipases from C. rugosa (Lip1, Lip3, and Lip4, respectively) towards the resolution of 2‐bromo‐arylacetic acid esters, an important class of chemical intermediates in the pharmaceutical industry. All exhibited a high enantioselectivity, with Lip4 preferring the R‐enantiomer (E‐value = 15), while Lip1 and Lip3 showed an S‐enantioselectivity >200. A combination of sequence and structure analysis of the three C. rugosa lipases suggested that position 296 could play a role in S‐ or R‐enantiomer preference of C. rugosa lipases. This led to the construction by site‐directed mutagenesis of Lip1 and Lip4 variants in which position 296 was, respectively, exchanged by a Gly, Ala, Leu, or Phe amino acid. Screening of these variants for their enantioselectivity toward 2‐bromo phenyl acetic acid octyl esters revealed that steric hindrance of the amino acid residue introduced at position 296 controls both the enantiopreference and the enantioselectivity value of the lipase: bulkier is the amino acid at position 296, larger is the selectivity towards the S‐enantiomer. To investigate further these observations at an atomic level, we carried out a preliminary modeling study of the tetrahedral intermediates formed by Lip1 and Lip4 with the (R, S)‐2‐bromo‐phenylacetic acid octyl ester enantiomers that provides some insight regarding the determinants responsible for lipase enantiodiscrimination. Biotechnol. Bioeng. 2011; 108:1749–1756.

Guaranteed Discrete Energy Optimization on Large Protein Design Problems

In Computational Protein Design (CPD), assuming a rigid backbone and amino-acid rotamer library, the problem of finding a sequence with an optimal conformation is NP-hard. In this paper, using Dunbracks rotamer library and Talaris2014 decomposable energy function, we use an exact deterministic method combining branch and bound, arc consistency, and tree-decomposition to provenly identify the global minimum energy sequence-conformation on full-redesign problems, defining search spaces of size up to 10(234). This is achieved on a single core of a standard computing server, requiring a maximum of 66GB RAM. A variant of the algorithm is able to exhaustively enumerate all sequence-conformations within an energy threshold of the optimum. These proven optimal solutions are then used to evaluate the frequencies and amplitudes, in energy and sequence, at which an existing CPD-dedicated simulated annealing implementation may miss the optimum on these full redesign problems. The probability of finding an optimum drops close to 0 very quickly. In the worst case, despite 1,000 repeats, the annealing algorithm remained more than 1 Rosetta unit away from the optimum, leading to design sequences that could differ from the optimal sequence by more than 30% of their amino acids.