Molecular motions in drug design: the coming age of the metadynamics method

Metadynamics is emerging as a useful free energy method in physics, chemistry and biology. Recently, it has been applied also to investigate ligand binding to biomolecules of pharmacological interest. Here, after introducing the basic idea of the method, we review applications to challenging targets for pharmaceutical intervention. We show that this methodology, especially when combined with a variety of other computational approaches such as molecular docking and/or molecular dynamics simulation, may be useful to predict structure and energetics of ligand/target complexes even when the targets lack a deep binding cavity, such as DNA and proteins undergoing fibrillation in neurodegenerative diseases. Furthermore, the method allows investigating the routes of molecular recognition and the associated binding energy profiles, providing a molecular interpretation to experimental data.     

Structure-energy relations in methylcobalamin with and without bound axial base

The properties of the Co-C bond in methylcobalamin (MeCbl) are analyzed by means of first-principles molecular dynamics. The optimized structure is in very good agreement with experiments, reproducing the bent-up deformation of the corrin ring as well as the metal-ligand bond distances. The analysis of the binding energies, bond orders, and vibrational stretching frequencies shows that the axial base slightly weakens the Co-C bond (by 4%), while the alkyl ligand substantially reinforces the Co-axial base bond (by 90%). These findings support several experiments and provide insight into the conversion between the base-on and base-off forms of the MeCbl cofactor.     

Structural Basis of Chitin Oligosaccharide Deacetylation

Cell signaling and other biological activities of chitooligosaccharides (COSs) seem to be dependent not only on the degree of polymerization, but markedly on the specific de‐N‐acetylation pattern. Chitin de‐N‐acetylases (CDAs) catalyze the hydrolysis of the acetamido group in GlcNAc residues of chitin, chitosan, and COS. A major challenge is to understand how CDAs specifically define the distribution of GlcNAc and GlcNH₂ moieties in the oligomeric chain. We report the crystal structure of the Vibrio cholerae CDA in four relevant states of its catalytic cycle. The two enzyme complexes with chitobiose and chitotriose represent the first 3D structures of a CDA with its natural substrates in a productive mode for catalysis, thereby unraveling an induced‐fit mechanism with a significant conformational change of a loop closing the active site. We propose that the deacetylation pattern exhibited by different CDAs is governed by critical loops that shape and differentially block accessible subsites in the binding cleft of CE4 enzymes.     

A Native Ternary Complex Trapped in a Crystal Reveals the Catalytic Mechanism of a Retaining Glycosyltransferase

Glycosyltransferases (GTs) comprise a prominent family of enzymes that play critical roles in a variety of cellular processes, including cell signaling, cell development, and host–pathogen interactions. Glycosyl transfer can proceed with either inversion or retention of the anomeric configuration with respect to the reaction substrates and products. The elucidation of the catalytic mechanism of retaining GTs remains a major challenge. A native ternary complex of a GT in a productive mode for catalysis is reported, that of the retaining glucosyl‐3‐phosphoglycerate synthase GpgS from M. tuberculosis in the presence of the sugar donor UDP‐Glc, the acceptor substrate phosphoglycerate, and the divalent cation cofactor. Through a combination of structural, chemical, enzymatic, molecular dynamics, and quantum‐mechanics/molecular‐mechanics (QM/MM) calculations, the catalytic mechanism was unraveled, thereby providing a strong experimental support for a front–side substrate‐assisted S_Ni‐type reaction.     

The conformational free energy landscape of beta-D-glucopyranose. Implications for substrate preactivation in beta-glucoside hydrolases

Using ab initio metadynamics we have computed the conformational free energy landscape of beta-D-glucopyranose as a function of the puckering coordinates. We show that the correspondence between the free energy and the Stoddard's pseudorotational itinerary for the system is rather poor. The number of free energy minima (9) is smaller than the number of ideal structures (13). Moreover, only six minima correspond to a canonical conformation. The structural features, the electronic properties, and the relative stability of the predicted conformers permit the rationalization of the occurrence of distorted sugar conformations in all the available X-ray structures of beta-glucoside hydrolase Michaelis complexes. We show that these enzymes recognize the most stable distorted conformers of the isolated substrate and at the same time the ones better prepared for catalysis in terms of bond elongation/shrinking and charge distribution. This suggests that the factors governing the distortions present in these complexes are largely dictated by the intrinsic properties of a single glucose unit.     

Multidimensional view of amyloid fibril nucleation in atomistic detail

Starting from a disordered aggregate, we have simulated the formation of ordered amyloid-like beta structures in a system formed by 18 polyvaline chains in explicit solvent, employing molecular dynamics accelerated by bias-exchange metadynamics. We exploited 8 different collective variables to compute the free energy of hundreds of putative aggregate structures, with variable content of parallel and antiparallel β-sheets and different packing among the sheets. This allowed characterizing in detail a possible nucleation pathway for the formation of amyloid fibrils: first the system forms a relatively large ordered nucleus of antiparallel β-sheets, and then a few parallel sheets start appearing. The relevant nucleation process culminates at this point: when a sufficient number of parallel sheets is formed, the free energy starts to decrease toward a new minimum in which this structure is predominant. The complex nucleation pathway we found cannot be described within classical nucleation theory, namely employing a unique simple reaction coordinate like the total content of β-sheets.     

The conformational free-energy landscape of β-D-mannopyranose: evidence for a (1)S(5) → B(2,5) → (O)S(2) catalytic itinerary in β-mannosidases

The mechanism of glycosidic bond cleavage by glycosidases involves substrate ring distortions in the Michaelis complex that favor catalysis. Retaining β-mannosidases bind the substrate in a (1)S(5) conformation, and recent experiments have proposed an unusual substrate conformational pathway ((1)S(5) → B(2,5) → (O)S(2)) for the hydrolysis reaction. By means of Car-Parrinello metadynamics simulations, we have obtained the conformational free-energy surface (FES) of a β-d-mannopyranose molecule associated with the ideal Stoddart conformational diagram. We have found that (1)S(5) is among the most stable conformers and simultaneously is the most preactivated conformation in terms of elongation/shortening of the C1-O1/C1-O5 bonds, C1-O1 orientation, and charge development at the anomeric carbon. Analysis of the computed FES gives support to the proposed (1)S(5) → B(2,5) → (O)S(2) catalytic itinerary, showing that the degree of preactivation of the substrate in glycoside hydrolases (GHs) is related to the properties of an isolated sugar ring. We introduce a simple preactivation index integrating several structural, electronic, and energetic properties that can be used to predict the conformation of the substrate in the Michaelis complex of any GH.     

Catalytic itinerary in 1,3-1,4-β-glucanase unraveled by QM/MM metadynamics. Charge is not yet fully developed at the oxocarbenium ion-like transition state

Retaining glycoside hydrolases (GHs), key enzymes in the metabolism of polysaccharides and glycoconjugates and common biocatalysts used in chemoenzymatic oligosaccharide synthesis, operate via a double-displacement mechanism with the formation of a glycosyl-enzyme intermediate. However, the degree of oxocarbenium ion character of the reaction transition state and the precise conformational itinerary of the substrate during the reaction, pivotal in the design of efficient inhibitors, remain elusive for many GHs. By means of QM/MM metadynamics, we unravel the catalytic itinerary of 1,3-1,4-β-glucanase, one of the most active GHs, belonging to family 16. We show that, in the Michaelis complex, the enzyme environment restricts the conformational motion of the substrate to stabilize a (1,4)B/(1)S(3) conformation of the saccharide ring at the -1 subsite, confirming that this distortion preactivates the substrate for catalysis. The metadynamics simulation of the enzymatic reaction captures the complete conformational itinerary of the substrate during the glycosylation reaction ((1,4)B/(1)S(3) -(4)E/(4)H(3) - (4)C(1)) and shows that the transition state is not the point of maximum charge development at the anomeric carbon. The overall catalytic mechanism is of dissociative type, and proton transfer to the glycosidic oxygen is a late event, clarifying previous kinetic studies of this enzyme.