Tunnels and channels facilitate the transport of small molecules, ions and water solvent in a large variety of proteins. Characteristics of individual transport pathways, including their geometry, physico-chemical properties and dynamics are instrumental for understanding of structure-function relationships of these proteins, for the design of new inhibitors and construction of improved biocatalysts. CAVER is a software tool widely used for the identification and characterization of transport pathways in static macromolecular structures. Herein we present a new version of CAVER enabling automatic analysis of tunnels and channels in large ensembles of protein conformations. CAVER 3.0 implements new algorithms for the calculation and clustering of pathways. A trajectory from a molecular dynamics simulation serves as the typical input, while detailed characteristics and summary statistics of the time evolution of individual pathways are provided in the outputs. To illustrate the capabilities of CAVER 3.0, the tool was applied for the analysis of molecular dynamics simulation of the microbial enzyme haloalkane dehalogenase DhaA. CAVER 3.0 safely identified and reliably estimated the importance of all previously published DhaA tunnels, including the tunnels closed in DhaA crystal structures. Obtained results clearly demonstrate that analysis of molecular dynamics simulation is essential for the estimation of pathway characteristics and elucidation of the structural basis of the tunnel gating. CAVER 3.0 paves the way for the study of important biochemical phenomena in the area of molecular transport, molecular recognition and enzymatic catalysis. The software is freely available as a multiplatform command-line application at http://www.caver.cz.

Virtual screening is becoming an important tool for drug discovery. However, the application of virtual screening has been limited by the lack of accurate scoring functions. Here, we present a novel scoring function, MedusaScore, for evaluating protein-ligand binding. MedusaScore is based on models of physical interactions that include van der Waals, solvation and hydrogen bonding energies. To ensure the best transferability of the scoring function, we do not use any protein-ligand experimental data for parameter training. We then test the MedusaScore for docking decoy recognition and binding affinity prediction and find superior performance compared to other widely used scoring functions. Statistical analysis indicates that one source of inaccuracy of MedusaScore may arise from the unaccounted entropic loss upon ligand binding, which suggests avenues of approach for further MedusaScore improvement.

Allylation of aromatic aldehydes 1a-m with allyl- and crotyl-trichlorosilanes 2- 4, catalyzed by the chiral N-oxide QUINOX (9), has been found to exhibit a significant dependence on the electronics of the aldehyde, with p-(trifluoromethyl)benzaldehyde 1g and its p-methoxy counterpart 1h affording the corresponding homoallylic alcohols 6g, h in 96 and 16% ee, respectively, at -40 degrees C. The kinetic and computational data indicate that the reaction is likely to proceed via an associative pathway involving neutral, octahedral silicon complex 22 with only one molecule of the catalyst involved in the rate- and selectivity-determining step. The crotylation with (E) and (Z)-crotyltrichlorosilanes 3 and 4 is highly diastereoselective, suggesting the chairlike transition state 5, which is supported by computational data. High-level quantum chemical calculations further suggest that attractive aromatic interactions between the catalyst 9 and the aldehyde 1 contribute to the enantiodifferentiation and that the dramatic drop in enantioselectivity, observed with the electron-rich aldehyde 1h, originates from narrowing the energy gap between the (R)- and (S)-reaction channels in the associative mechanism (22). Overall, a good agreement between the theoretically predicted enantioselectivities for 1a and 1h and the experimental data allowed to understand the specific aspects of the reaction mechanism.

Background: Tunnel properties affect ligand passage in enzymes with buried active sites. Results: A tunnel mutation from leucine to tryptophan changes the mechanism of bromide ion release from haloalkane dehalogenase LinB. Conclusion:Interactions of the bromide ion with the tryptophan increase free energy barrier for its passage, causing the reaction mechanism change. Significance: The results provide guidelines for enzyme engineering.

No abstract

Proline-tryptophan complexes derived from experimental structures are investigated by quantum chemical procedures known to properly describe the London dispersion energy. We study two geometrical arrangements: the "L-shaped", stabilized by an H-bond, and the "stacked-like", where the two residues are in parallel orientation without any H-bond. Interestingly, the interaction energies in both cases are comparable and very large ( approximately 7 kcal mol(-1)). The strength of stabilization in the stacked arrangement is rather surprising considering the fact that only one partner has an aromatic character. The interaction energy decomposition using the SAPT method further demonstrates the very important role of dispersion energy in such arrangement. To elucidate the structural features responsible for this unexpectedly large stabilization we examined the role of the nitrogen heteroatom and the importance of the cyclicity of the proline residue. We show that the electrostatic interaction due to the presence of the dipole, caused by the nitrogen heteroatom, contributes largely to the strength of the interaction. Nevertheless, the cyclic arrangement of proline, which allows for the largest amount of dispersive contact with the aromatic partner, also has a notable-effect. Geometry optimizations carried out for the "stacked-like" complexes show that the arrangements derived from protein structure are close to their gas phase optimum geometry, suggesting that the environment has only a minor effect on the geometry of the interaction. We conclude that the strength of proline non-covalent interactions, combined with this residue's rigidity, might be the explanation for its prominent role in protein stabilization and recognition processes.

Recombinant ligands derived from small protein scaffolds show promise as robust research and diagnostic reagents and next generation protein therapeutics. Here, we derived high-affinity binders of human interferon gamma (hIFNγ) from the three helix bundle scaffold of the albumin-binding domain (ABD) of protein G from Streptococcus G148. Computational interaction energy mapping, solvent accessibility assessment, and in silico alanine scanning identified 11 residues from the albumin-binding surface of ABD as suitable for randomization. A corresponding combinatorial ABD scaffold library was synthesized and screened for hIFNγ binders using in vitro ribosome display selection, to yield recombinant ligands that exhibited K(d) values for hIFNγ from 0.2 to 10 nM. Molecular modeling, computational docking onto hIFNγ, and in vitro competition for hIFNγ binding revealed that four of the best ABD-derived ligands shared a common binding surface on hIFNγ, which differed from the site of human IFNγ receptor 1 binding. Thus, these hIFNγ ligands provide a proof of concept for design of novel recombinant binding proteins derived from the ABD scaffold.

Rutinosidases (α‐l‐rhamnosyl‐β‐d‐glucosidases) catalyze the cleavage of the glycosidic bond between the aglycone and the disaccharide rutinose (α‐l‐rhamnopyranosyl‐(1→6)‐β‐d‐glucopyranose) of specific flavonoid glycosides such as rutin (quercetin 3‐O‐rutinoside). Microbial rutinosidases are part of the rutin catabolic pathway, enabling the microorganism to utilize rutin and related plant phenolic glycosides. Here, we report the first three‐dimensional structure of a rutinosidase determined at 1.27‐Å resolution. The rutinosidase from Aspergillus niger K2 (AnRut), a member of glycoside hydrolase family GH‐5, subfamily 23, was heterologously produced in Pichia pastoris. The X‐ray structure of AnRut is represented by a distorted (β/α)8 barrel fold with its closest structural homologue being an exo‐β‐(1,3)‐glucanase from Candida albicans (CaExg). The catalytic site is located in a deep pocket with a striking structural similarity to CaExg. However, the entrance to the active site of AnRut has been found to be different from that of CaExg – a mostly unstructured section of ~ 40 residues present in CaExg is missing in AnRut, whereas an additional loop of 13 amino acids partially covers the active site of AnRut. NMR analysis of reaction products provided clear evidence for a retaining reaction mechanism of AnRut. Unexpectedly, quercetin 3‐O‐glucoside was found to be a better substrate than rutin, and thus, AnRut cannot be considered a typical diglycosidase. Mutational analysis of conserved active site residues in combination with in silico modeling allowed identification of essential interactions for enzyme activity and helped to reveal further details of substrate binding. The protein sequence of AnRut has been revised.DatabasesThe nucleotide sequence of the rutinosidase‐encoding gene is available in the GenBank database under the accession number MN393234. Structural data are available in the PDB database under the accession number 6I1A.Enzymeα‐l‐Rhamnosyl‐β‐d‐glucosidase (EC 3.2.1.168).