Structural Origin of Metal Specificity in Isatin Hydrolase from Labrenzia aggregata Investigated by Computer Simulations.

We performed quantum-chemical calculations, ab initio molecular dynamics, hybrid quantum mechanics/molecular mechanics (QM/MM) and enhanced sampling metadynamics simulations to investigate the origin of metal specificity in isatin hydrolase from Labrenzia aggregata. The peculiar octahedral binding geometry of the Mn2+ ion in the Michaelis complex includes both the isatin substrate and the catalytic water within the first coordination shell of the cation. Our calculations show that the same arrangement of the ligands cannot be efficiently achieved in the presence of other small divalent metal cations such as Zn2+ or Cu2+ . On the contrary, bulkier alkaline-earth cations such as Mg2+ , which allow octahedral coordination, are not able to activate the catalytic water into the stronger OH- nucleophile required to attack the stable N-aryl-amide moiety of isatin.

Mechanical unfolding pathway of a model ß-peptide foldamer.

Foldamers constructed from oligomers of ß-peptides form stable secondary helix structures already for small chain lengths, which makes them ideal candidates for the investigation of the (un)folding of polypeptides. Here, the results of molecular simulations of the mechanical unfolding of a ß-heptapeptide in methanol solvent revealing the detailed unfolding pathway are reported. The unfolding process is shown to proceed via a stable intermediate even for such a small system. This result is arrived at performing non-equilibrium force ramp simulations employing different pulling velocities and also using standard calculations of the potential of mean force, i.e., the free energy as a function of the helix elongation. It is thus demonstrated that even with the rather large pulling velocities employed in the force ramp simulations relevant information about the equilibrium kinetics can be obtained. The smallness of the system allows a detailed analysis of the unfolding pathway, which is characterized by an opening of the terminal loops followed by the unfolding of the center. This sequence is in accord with the configurational preferences of the system that also are responsible for the stability of the 314-helix. From an analysis of the distributions of rupture forces and the force spectra, the kinetic rates for both transitions were determined and common models were used to extract geometric quantities describing the free energy landscape of the system.

A fundamental catalytic difference between zinc and manganese dependent enzymes revealed in a bacterial isatin hydrolase.

The catalytic mechanism of the cyclic amidohydrolase isatin hydrolase depends on a catalytically active manganese in the substrate-binding pocket. The Mn2+ ion is bound by a motif also present in other metal dependent hydrolases like the bacterial kynurenine formamidase. The crystal structures of the isatin hydrolases from Labrenzia aggregata and Ralstonia solanacearum combined with activity assays allow for the identification of key determinants specific for the reaction mechanism. Active site residues central to the hydrolytic mechanism include a novel catalytic triad Asp-His-His supported by structural comparison and hybrid quantum mechanics/classical mechanics simulations. A hydrolytic mechanism for a Mn2+ dependent amidohydrolases that disfavour Zn2+ as the primary catalytically active site metal proposed here is supported by these likely cases of convergent evolution. The work illustrates a fundamental difference in the substrate-binding mode between Mn2+ dependent isatin hydrolase like enzymes in comparison with the vast number of Zn2+ dependent enzymes.

Determining Factors for the Unfolding Pathway of Peptides, Peptoids, and Peptidic Foldamers.

We present a study of the mechanical unfolding pathway of five different oligomers (α-peptide, ß-peptide, δ-aromatic-peptides, α/γ-peptides, and ß-peptoids), adopting stable helix conformations. Using force-probe molecular dynamics, we identify the determining structural factors for the unfolding pathways and reveal the interplay between the hydrogen bond strength and the backbone rigidity in the stabilization of their helix conformations. On the basis of their behavior, we classify the oligomers in four groups and deduce a set of rules for the prediction of the unfolding pathways of small foldamers.

Comparative Study of the Mechanical Unfolding Pathways of α- and ß-Peptides.

Using molecular simulations, we analyze the unfolding pathways of various peptides. We compare the mechanical unfolding of a ß-alanine's octamer (ß-HAla8) and an α-alanine's decamer (α-Ala10). Using force-probe molecular-dynamics simulations, to induce unfolding, we show that the 3(14)-helix formed by ß-HAla8 is mechanically more stable than the α-helix formed by α-Ala10, although both structures are stabilized by six hydrogen bonds. Additionally, computations of the potential of mean force validate this result and show that also the thermal stability of the 3(14)-helix is higher. It is demonstrated that ß-HAla8 unfolds in a two-step fashion with a stable intermediate. This is contrasted with the known single-step scenario of the unfolding of α-Ala10. Furthermore, we present a study of the chain-length dependence of the mechanical unfolding pathway of the 3(14)-helix. The calculation of the dynamic strength for oligomers with chain lengths ranging from 6 to 18 monomers shows that the unfolding pathway of helices with an integer and noninteger number of turns has m + 1 and m energy barriers, respectively, with m being the number of complete turns. The additional barrier for helices with an integer number of turns is shown to be related to the breaking of the N-terminus' hydrogen bond.