Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites.

Zn2+ is one of the most versatile biologically available metal ions, but accurate modeling of Zn2+ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn2+ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn2+ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-ß-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn2+ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn2+ site geometries that could guide the further development of biomolecular Zn2+ models.

Structural and Functional Insight into the Mechanism of the Fe-S Cluster-Dependent Dehydratase from Paralcaligenes ureilyticus.

Enzyme-catalyzed reaction cascades play an increasingly important role for the sustainable manufacture of diverse chemicals from renewable feedstocks. For instance, dehydratases from the ilvD/EDD superfamily have been embedded into a cascade to convert glucose via pyruvate to isobutanol, a platform chemical for the production of aviation fuels and other valuable materials. These dehydratases depend on the presence of both a Fe-S cluster and a divalent metal ion for their function. However, they also represent the rate-limiting step in the cascade. Here, catalytic parameters and the crystal structure of the dehydratase from Paralcaligenes ureilyticus (PuDHT, both in presence of Mg2+ and Mn2+ ) were investigated. Rate measurements demonstrate that the presence of stoichiometric concentrations Mn2+ promotes higher activity than Mg2+ , but at high concentrations the former inhibits the activity of PuDHT. Molecular dynamics simulations identify the position of a second binding site for the divalent metal ion. Only binding of Mn2+ (not Mg2+ ) to this site affects the ligand environment of the catalytically essential divalent metal binding site, thus providing insight into an inhibitory mechanism of Mn2+ at higher concentrations. Furthermore, in silico docking identified residues that play a role in determining substrate binding and selectivity. The combined data inform engineering approaches to design an optimal dehydratase for the cascade.

Structure-Guided Modulation of the Catalytic Properties of [2Fe-2S]-Dependent Dehydratases.

The FeS cluster-dependent dihydroxyacid dehydratases (DHADs) and sugar acid-specific dehydratases (DHTs) from the ilvD/EDD superfamily are key enzymes in the bioproduction of a wide variety of chemicals. We analyzed [2Fe-2S]-dependent dehydratases in silico and in vitro, deduced functionally relevant sequence, structure, and activity relationships within the ilvD/EDD superfamily, and we propose a new classification based on their evolutionary relationships and substrate profiles. In silico simulations and analyses identified several key positions for specificity, which were experimentally investigated with site-directed and saturation mutagenesis. We thus increased the promiscuity of DHAD from Fontimonas thermophila (FtDHAD), showing >10-fold improved activity toward D-gluconate, and shifted the substrate preference of DHT from Paralcaligenes ureilyticus (PuDHT) toward shorter sugar acids (recording a six-fold improved activity toward the non-natural substrate D-glycerate). The successful elucidation of the role of important active site residues of the ilvD/EDD superfamily will further guide developments of this important biocatalyst for industrial applications.

DynaBiS: A hierarchical sampling algorithm to identify flexible binding sites for large ligands and peptides.

Knowing the ligand or peptide binding site in proteins is highly important to guide drug discovery, but experimental elucidation of the binding site is difficult. Therefore, various computational approaches have been developed to identify potential binding sites in protein structures. However, protein and ligand flexibility are often neglected in these methods due to efficiency considerations despite the recognition that protein-ligand interactions can be strongly affected by mutual structural adaptations. This is particularly true if the binding site is unknown, as the screening will typically be performed based on an unbound protein structure. Herein we present DynaBiS, a hierarchical sampling algorithm to identify flexible binding sites for a target ligand with explicit consideration of protein and ligand flexibility, inspired by our previously presented flexible docking algorithm DynaDock. DynaBiS applies soft-core potentials between the ligand and the protein, thereby allowing a certain protein-ligand overlap resulting in efficient sampling of conformational adaptation effects. We evaluated DynaBiS and other commonly used binding site identification algorithms against a diverse evaluation set consisting of 26 proteins featuring peptide as well as small ligand binding sites. We show that DynaBiS outperforms the other evaluated methods for the identification of protein binding sites for large and highly flexible ligands such as peptides, both with a holo or apo structure used as input.

Thiazoline-Specific Amidohydrolase PurAH Is the Gatekeeper of Bottromycin Biosynthesis.

The ribosomally synthesized and post-translationally modified peptide (RiPP) bottromycin A2 possesses potent antimicrobial activity. Its biosynthesis involves the enzymatic formation of a macroamidine, a process previously suggested to require the concerted efforts of a YcaO enzyme (PurCD) and an amidohydrolase (PurAH) in vivo. In vitro, PurCD alone is sufficient to catalyze formation of the macroamidine, but the process is reversible. We set out to probe the role of PurAH in macroamidine formation in vitro. We demonstrate that PurAH is highly selective for macroamidine-containing precursor peptides and cleaves C-terminal of a thiazoline, thus removing the follower peptide. After follower cleavage, macroamidine formation is irreversible, indicating PurAH as the gatekeeper of bottromycin biosynthesis. The structure of PurAH suggests residues involved in catalysis, which were probed through mutagenesis.

NewProt - a protein engineering portal.

The NewProt protein engineering portal is a one-stop-shop for in silico protein engineering. It gives access to a large number of servers that compute a wide variety of protein structure characteristics supporting work on the modification of proteins through the introduction of (multiple) point mutations. The results can be inspected through multiple visualizers. The HOPE software is included to indicate mutations with possible undesired side effects. The Hotspot Wizard software is embedded for the design of mutations that modify a proteins' activity, specificity, or stability. The NewProt portal is freely accessible at http://newprot.cmbi.umcn.nl/ and http://newprot.fluidops.net/.