Computing infrared spectra of proteins using the exciton model.

The ability to compute from first principles the infrared spectrum of a protein in solution phase representing a biological system would provide a useful connection to atomistic models of protein structure and dynamics. Indeed, such calculations are a vital complement to 2DIR experimental measurements, allowing the observed signals to be interpreted in terms of detailed structural and dynamical information. In this article, we have studied nine structurally and spectroscopically well-characterized proteins, representing a range of structural types. We have simulated the equilibrium conformational dynamics in an explicit point charge water model. Using the resulting trajectories based on MD simulations, we have computed the one and two dimensional infrared spectra in the Amide I region, using an exciton approach, in which a local mode basis of carbonyl stretches is considered. The role of solvent in shifting the Amide I band (by 30 to 50 cm-1 ) is clearly evident. Similarly, the conformational dynamics contribute to the broadening of peaks in the spectrum. The inhomogeneous broadening in both the 1D and 2D spectra reflects the significant conformational diversity observed in the simulations. Through the computed 2D cross-peak spectra, we show how different pulse schemes can provide additional information on the coupled vibrations. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Hit-to-lead and lead optimization binding free energy calculations for G protein-coupled receptors.

We apply the hit-to-lead ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and lead-optimization TIES (thermodynamic integration with enhanced sampling) methods to compute the binding free energies of a series of ligands at the A1 and A2A adenosine receptors, members of a subclass of the GPCR (G protein-coupled receptor) superfamily. Our predicted binding free energies, calculated using ESMACS, show a good correlation with previously reported experimental values of the ligands studied. Relative binding free energies, calculated using TIES, accurately predict experimentally determined values within a mean absolute error of approximately 1 kcal mol-1. Our methodology may be applied widely within the GPCR superfamily and to other small molecule-receptor protein systems.

Ensemble-Based Steered Molecular Dynamics Predicts Relative Residence Time of A2A Receptor Binders.

Drug-target residence time, the length of time for which a small molecule stays bound to its receptor target, has increasingly become a key property for optimization in drug discovery programs. However, its in silico prediction has proven difficult. Here we describe a method, using atomistic ensemble-based steered molecular dynamics (SMD), to observe the dissociation of ligands from their target G protein-coupled receptor in a time scale suitable for drug discovery. These dissociation simulations accurately, precisely, and reproducibly identify ligand-residue interactions and quantify the change in ligand energy values for both protein and water. The method has been applied to 17 ligands of the A2A adenosine receptor, all with published experimental kinetic binding data. The residues that interact with the ligand as it dissociates are known experimentally to have an effect on binding affinities and residence times. There is a good correlation ( R2 = 0.79) between the computationally calculated change in water-ligand interaction energy and experimentally determined residence time. Our results indicate that ensemble-based SMD is a rapid, novel, and accurate semi-empirical method for the determination of drug-target relative residence time.

Examining the role of protein structural dynamics in drug resistance in Mycobacterium tuberculosis.

Antimicrobial resistance represents a growing global health problem. The emergence of novel resistance mechanisms necessitates the development of alternative approaches to investigate the molecular fundamentals of resistance, leading ultimately to new strategies for counteracting them. To gain deeper insight into antibiotic-target interactions, the binding of the frontline anti-tuberculosis drug isoniazid (INH) to a target enzyme, InhA, from Mycobacterium tuberculosis was studied using ultrafast two-dimensional infrared (2D-IR) spectroscopy and molecular simulations. Comparing wild-type InhA with a series of single point mutations, it was found that binding of the INH-NAD inhibitor to susceptible forms of the enzyme increased the vibrational coupling between residues located in the Rossmann fold co-factor binding site of InhA and suppressed dynamic fluctuations of the enzyme structure. The effect correlated with biochemical assay data, being reduced in the INH-resistant S94A mutant and absent in the biochemically-inactive P193A control. Molecular dynamics simulations and calculations of inter-residue couplings indicate that the changes in coupling and dynamics are not localised to the co-factor binding site, but permeate much of the protein. We thus propose that the resistant S94A mutation circumvents subtle changes in global structural dynamics caused by INH upon binding to the wild-type enzyme that may impact upon the formation of important protein-protein complexes in the fatty acid synthase pathway of M. tuberculosis.

Simulation of Two-Dimensional Infrared Spectroscopy of Peptides Using Localized Normal Modes.

Nonlinear two-dimensional infrared spectroscopy (2DIR) is most commonly simulated within the framework of the exciton method. The key parameters for these calculations include the frequency of the oscillators within their molecular environments and coupling constants that describe the strength of coupling between the oscillators. It is shown that these quantities can be obtained directly from harmonic frequency calculations by exploiting a procedure that localizes the normal modes. This approach is demonstrated using the amide I modes of polypeptides. For linear and cyclic diamides and hexapeptide Z-Aib-L-Leu-(Aib)2-Gly-Aib-OtBu, the computed parameters are compared with those from existing schemes, and the resulting 2DIR spectra are consistent with experimental observations. The incorporation of conformational averaging of structures from molecular dynamics simulations is discussed, and a hybrid scheme wherein the Hamiltonian matrix from the quantum chemical local-mode approach is combined with fluctuations from empirical schemes is shown to be consistent with experiment. The work demonstrates that localized vibrational modes can provide a foundation for the calculation of 2DIR spectra that does not rely on extensive parametrization and can be applied to a wide range of systems. For systems that are too large for quantum chemical harmonic frequency calculations, the local-mode approach provides a convenient platform for the development of site frequency and coupling maps.