The anthropogenic toxic compound 1,2,3-trichloropropane is poorly degradable by natural enzymes. We have previously constructed the haloalkane dehalogenase DhaA31 by focused directed evolution ( Pavlova, M. et al. Nat. Chem. Biol. 2009 , 5 , 727 - 733 ), which is 32 times more active than the wild-type enzyme and is currently the most active variant known against that substrate. Recent evidence has shown that the structural basis responsible for the higher activity of DhaA31 was poorly understood. Here we have undertaken a comprehensive computational study of the main steps involved in the biocatalytic hydrolysis of 1,2,3-trichloropropane to decipher the structural basis for such enhancements. Using molecular dynamics and quantum mechanics approaches we have surveyed (i) the substrate binding, (ii) the formation of the reactive complex, (iii) the chemical step, and (iv) the release of the products. We showed that the binding of the substrate and its transport through the molecular tunnel to the active site is a relatively fast process. The cleavage of the carbon-halogen bond was previously identified as the rate-limiting step in the wild-type. Here we demonstrate that this step was enhanced in DhaA31 due to a significantly higher number of reactive configurations of the substrate and a decrease of the energy barrier to the S2 reaction. C176Y and V245F were identified as the key mutations responsible for most of those improvements. The release of the alcohol product was found to be the rate-limiting step in DhaA31 primarily due to the C176Y mutation. Mutational dissection of DhaA31 and kinetic analysis of the intermediate mutants confirmed the theoretical observations. Overall, our comprehensive computational approach has unveiled mechanistic details of the catalytic cycle which will enable a balanced design of more efficient enzymes. This approach is applicable to deepen the biochemical knowledge of a large number of other systems and may contribute to robust strategies in the development of new biocatalysts.

No abstract

No abstract

The emergence of large-scale order in self-organized systems relies on local interactions between individual components. During bacterial cell division, the tubulin-homolog FtsZ polymerizes into treadmilling filaments that further assemble into a cytoskeletal ring. Although minimal in vitro assays have shown the striking self-organization capacity of FtsZ filaments, such as dynamic chiral assemblies, how these large-scale structures emerge and relate to individual filament properties remains poorly understood. To understand this quantitatively, we combined minimal chiral active matter simulations with biochemical reconstitution experiments. Using STED and TIRF microscopy as well as high-speed AFM, we imaged the behavior of FtsZ filaments on different spatial scales. Simulations and experiments revealed that filament density and flexibility define the local and global order of the system: At intermediate densities, flexible filaments organize into chiral rings and polar bands, while an effectively nematic organization dominates for high filament densities and for mutant filaments with increased rigidity. Our predicted phase diagram captured these features quantitatively, demonstrating how filament flexibility, density and chirality cooperate with activity to give rise to a large repertoire of collective behaviors. These properties are likely important for the dynamic organization of soft chiral matter, including that of treadmilling FtsZ filaments during bacterial cell division.