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.