Quantum mechanical/molecular mechanical approaches to transition state structure: mechanism of sialidase action

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Bacterial Quorum-Sensing Molecules as Promising Natural Inhibitors of Candida albicans Virulence Dimorphism: An In Silico and In Vitro Study

Farnesol, a self-secreted quorum-sensing molecule (QSM) of Candida albicans, has been known to limit yeast-to-hyphal transition by blocking the RAS1-cAMP-PKA pathway. In a similar fashion, certain bacterial QSMs have also been reported to be successful in attenuating C. albicans biofilm and hyphal formation at relatively high cell density. This prompted us to investigate the antihyphal efficacy of certain bacterial QSMs through virtual docking against seminal drug targets, viz., CYCc and RAS1, that have been reported to be the hallmark players in C. albicans dimorphic virulence cascade. Against this backdrop, 64 QSMs belonging to five different bacterial QS signaling systems were subjected to initial virtual screening with farnesol as reference. Data of the virtual screening unveiled QSMs belonging to diketopiperazines (DKPs), i.e., 3-benzyl-6-isobutylidene-2,5-piperazinedione (QSSM 1157) and cyclo(l-Pro-l-Leu) (QSSM 1112), as potential inhibitors of CYCc and RAS1 with binding energies of -8.2 and -7.3 kcal mol-1, respectively. Further, the molecular dynamics simulations (for 50 ns) of CYCc-QSSM 1157 and RAS1-QSSM 1112 complexes revealed the mean ligand root mean square deviation (RMSD) values of 0.35 and 0.27 Å, respectively, which endorsed the rigid nature, less fluctuation in binding stiffness, and conformation of binding complexes. Furthermore, the identified two QSMs were found to be good in solubility, absorption, and permeation and less toxic in nature, as revealed by pharmacokinetics and toxicity analyses. In addition, the in vitro antihyphal assays using liquid and solid media, germ-tube experiment, and microscopic analysis strongly validated DKP-QSSM 1112 as a promising inhibitor of hyphal transition. Taken together, the present study unequivocally proves that DKPs can be used as potent inhibitors of C. albicans virulence dimorphism.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

QM/MM through the 1990s: The First Twenty Years of Method Development and Applications

The 2013 Nobel Prize in Chemistry was awarded to the authors of the first two publications utilizing the concept of combined quantum mechanical and molecular mechanical (QM/MM) methods. In celebrating this great event in computational chemistry, we review the early development of combined QM/MM techniques and the associated events that took place through the mid-1990s. We also offer some prospects for the future development of quantum mechanical techniques for macromolecular systems.

Kinetic Isotope Effects from QM/MM Subset Hessians: "Cut-Off" Analysis for SN2 Methyl Transfer in Solution

Isotopic partition-function ratios and kinetic isotope effects for reaction of S-adenosylmethionine with catecholate in water are evaluated using a subset of 324 atoms within its surrounding aqueous environment at the AM1/TIP3P level. Two alternative methods for treating motion in the six librational degrees of freedom of the subset atoms relative to their environment are compared. A series of successively smaller subset Hessians are generated by cumulative deletion of rows and columns from the initial 972 × 972 Hessian. We find that it is better to treat these librations as vibrations than as translations and rotations and that there is no need to invoke the Teller-Redlich product rule. The validity of "cut-off" procedures for computation of isotope effects with truncated atomic subsets is assessed: to ensure errors in ln(KIE) < 1% (or 2% for the quantum-corrected KIE) for all isotopic substitutions considered, it is necessary to use a less-restrictive procedure than is suggested by the familiar two-bond cutoff rule.

Structure and Mechanism of Action of an Inverting Mutant Sialidase

Mutagenesis of the conserved tyrosine (Y370) of the Micromonospora viridifaciens sialidase to small amino acids changes the mechanism of catalysis from retention of anomeric configuration to inversion [Watson, J. N., et al. (2003) Biochemistry 42, 12682-12690]. For the Y370G mutant enzyme-catalyzed hydrolysis of a series of aryl sialosides and 3'-sialyllactose, the derived Brønsted parameters (beta(lg)) on k(cat) and k(cat)/K(m) are -0.63 +/- 0.05 and -0.80 +/- 0.08, respectively. Thus, for the Y370G enzyme, glycosidic C-O bond cleavage is rate-determining. Analysis of the activity of the Y370G mutant and wild-type enzymes against a substrate [3,4-dihydro-2H-pyrano[3,2-c]pyridinium alpha-d-N-acetylneuraminide (DHP-alphaNeu5Ac)] whose hydrolysis cannot be accelerated by acid catalysis is consistent with these reactions proceeding via S(N)1 and S(N)2 mechanisms, respectively. The overall structure of the Y370G mutant sialidase active site is very similar to the previously reported wild-type structure [Gaskell, A., et al. (1995) Structure 3, 1197-1205], although removal of the tyrosine residue creates two significant changes to the active site. First, the anomeric oxygen atom of the hydrolysis product (beta-N-acetylneuraminic acid) and four water molecules bind in the large cavity created by the Y370G mutation. Second, the side chain of Asn310 moves to make a strong hydrogen bond to one of the bound water molecules.

Computational enzymology

Numerical simulations of enzyme reaction mechanisms are beginning to provide quantitative as well as qualitative insights. Methods based on hybrid quantum mechanical/molecular mechanical technique permit the natural inclusion of protein solvation effects. Coupled with modern experimental techniques, the numerical simulations are providing details at the atomic level about how enzyme structure influences its function.