Christopher N. Rowley

Ion selectivity in channels and transporters

A multitude of biological processes requires the participation of specific cations, such as H+, Na+, K+, Ca2+, and Mg2+. Many of these processes can take place only when proteins have the ability to discriminate between different ions with a very high fidelity. How this is possible is a fundamental

The Solvation Structure of Na(+) and K(+) in Liquid Water Determined from High Level ab Initio Molecular Dynamics Simulations.

Knowledge of the hydration structure of Na(+) and K(+) in the liquid phase has wide ranging implications in the field of biological chemistry. Despite numerous experimental and computational studies, even basic features such as the coordination number of these alkali ions in liquid water, thought to play a critical role in selectivity, continue to be the subject of intensive debates. Simulations based on accurate potential energy surfaces offer one approach to resolve these issues by providing reliable results on ion hydration. In this article, we report the results from molecular dynamics simulations of Na(+) and K(+) hydration based on a novel and rigorous strategy designed to overcome the challenges of QM/MM simulations of solvent molecules in the liquid phase. In this method, which we call Flexible Inner Region Ensemble Separator (FIRES), the ion and a fixed number of nearest water molecules form a dynamical and flexible inner region that is represented with high level ab initio quantum mechanical (QM) methods, while the water molecules from the surrounding bulk form an outer region that is represented by a polarizable molecular mechanical (MM) force field. Simulations yield rigorously correct thermodynamic averages as long as the solvent molecules in the flexible inner and outer regions are not allowed to exchange. Extensive FIRES simulations were carried out based on a QM/MM model in which the Na(+) or K(+) ion and the 12 nearest water molecules were represented by high level ab initio methods (RI-MP2/def2-TZVP and density functional theory with PBE/def2-TZVP), while the surrounding MM water molecules were represented by the polarizable SWM4-NDP potential. On the basis of these results, the ion coordination numbers are estimated to be within the range of 5.7-5.8 for Na(+) and 6.9-7.0 for K(+).

Solution conformation of C-linked antifreeze glycoprotein analogues and modulation of ice recrystallization.

Antifreeze glycoproteins (AFGPs) are a unique class of proteins that are found in many organisms inhabiting subzero environments and ensure their survival by preventing ice growth in vivo. During the last several years, our laboratory has synthesized functional C-linked AFGP analogues (3 and 5) that possess custom-tailored antifreeze activity suitable for medical, commercial, and industrial applications. These compounds are potent inhibitors of ice recrystallization and do not exhibit thermal hysteresis. The current study explores how changes in the length of the amide-containing side chain between the carbohydrate moiety and the polypeptide backbone in 5 influences ice recrystallization inhibition (IRI) activity. Analogue 5 (n = 3, where n is the number of carbons in the side chain) was a potent inhibitor of ice recrystallization, while 4, 6, and 7 (n = 4, 2, and 1, respectively) exhibited no IRI activity. The solution conformation of the polypeptide backbone in C-linked AFGP analogues 4-7 was examined using circular dichroism (CD) spectroscopy. The results suggested that all of the analogues exhibit a random coil conformation in solution and that the dramatic increase in IRI activity observed with 5 is not due to a change in long-range solution conformation. Variable-temperature (1)H NMR studies on truncated analogues 26-28 failed to elucidate the presence of persistent intramolecular bonds between the amide in the side chain and the peptide backbone. Molecular dynamics simulations performed on these analogues also failed to show persistent intramolecular hydrogen bonds. However, the simulations did indicate that the side chain of IRI-active analogue 26 (n = 3) adopts a unique short-range solution conformation in which it is folded back onto the peptide backbone, orienting the more hydrophilic face of the carbohydrate moiety away from the bulk solvent. In contrast, the solution conformation of IRI-inactive analogues 25, 27, and 28 had fully extended side chains, with the carbohydrate moiety being exposed to bulk solvent. These results illustrate how subtle changes in conformation and carbohydrate orientation dramatically influence IRI activity in C-linked AFGP analogues.

Why Can Hydrogen Sulfide Permeate Cell Membranes

The high membrane permeability of H2S was studied using polarizable molecular dynamics simulations of a DPPC lipid bilayer. The solubility-diffusion model predicts permeability coefficients of H2S and H2O that are in good agreement with experiment. The computed diffusion coefficient profile shows H2S to diffuse at a lower rate than H2O, but the barrier for H2S permeation on the Gibbs energy profile is negligible. The hydrophobicity of H2S allows it to partition into the paraffinic interior of the membrane readily.

Amidolithium and Amidoaluminum Catalyzed Synthesis of Substituted Guanidines : An Interplay of DFT Modeling and Experiment

The synthesis of substituted guanidines is of significant interest for their use as versatile ligands and for the synthesis of bioactive molecules. Lithium amides supported by tetramethylethylenediamine have recently been shown to catalyze the guanylation of amines with carbodiimide. In this report, density functional theory (DFT) calculations are used to provide insight into the mechanism of this transformation. The mechanism identified through our calculations is a carbodiimide insertion into the lithium-amide bond to form a lithium guanidinate, followed by a proton transfer from the amine. The proton transfer transition state requires the dissociation of one of the chelating nitrogen centers of the lithium guanidinate, proton abstraction from the amine, and bond formation between the lithium center and the amine nitrogen. On the basis of this mechanism, further calculations predicted that aluminum amides would also function as active catalysts for the guanylation of amines. We confirm this experimentally and report the development of aluminum amides as a new main group catalyst for the guanylation of a range of electron-poor amines with carbodiimide.

Simulation-Based Approaches for Determining Membrane Permeability of Small Compounds.

Predicting the rate of nonfacilitated permeation of solutes across lipid bilayers is important to drug design, toxicology, and signaling. These rates can be estimated using molecular dynamics simulations combined with the inhomogeneous solubility-diffusion model, which requires calculation of the potential of mean force and position-dependent diffusivity of the solute along the transmembrane axis. In this paper, we assess the efficiency and accuracy of several methods for the calculation of the permeability of a model DMPC bilayer to urea, benzoic acid, and codeine. We compare umbrella sampling, replica exchange umbrella sampling, adaptive biasing force, and multiple-walker adaptive biasing force for the calculation of the transmembrane PMF. No definitive advantage for any of these methods in their ability to predict the membrane permeability coefficient Pm was found, provided that a sufficiently long equilibration is performed. For diffusivities, a Bayesian inference method was compared to a generalized Langevin method, both being sensitive to chosen parameters and the slow relaxation of membrane defects. Agreement within 1.5 log units of the computed Pm with experiment is found for all permeants and methods. Remaining discrepancies can likely be attributed to limitations of the force field as well as slowly relaxing collective movements within the lipid environment. Numerical calculations based on model profiles show that Pm can be reliably estimated from only a few data points, leading to recommendations for calculating Pm from simulations.

Molecular simulation of nonfacilitated membrane permeation

This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Range-Separated DFT Functionals are Necessary to Model Thio-Michael Additions.

The textbook mechanism for the addition of a thiol to an olefin is the Michael-type addition, which involves a nucleophilic attack of a thiolate anion on an alkene to form a carbanion intermediate. Previous computational models of these reactions have proposed alternative mechanisms, as no minimum corresponding to the carbanion intermediate was present on the potential energy surface. We show that many popular pure and hybrid DFT functionals, such as PBE and B3LYP, erroneously predict that the carbanion is not an intermediate, favoring a noncovalent charge-transfer complex stabilized spuriously by delocalization error. Range-separated DFT functionals correct this problem and predict stable carbanion structures and energies. In particular, calculations using the ωB97X-D functional are in close agreement with CCSD(T) data for the structures and energies of a series of thio-carbanions. Range-separated functionals will make it possible to model the reaction mechanisms of Michael-type additions that occur in biochemistry, such as the covalent modification of a cysteine side chain by drugs containing an electrophilic double bond.

Analysis of the Critical Step in Catalytic Carbodiimide Transformation : Proton Transfer from Amines, Phosphines, and Alkynes to Guanidinates, Phosphaguanidinates, and Propiolamidinates with Li and Al Catalysts

While lithium amides supported by tetramethylethylenediamine (TMEDA) are efficient catalysts in the synthesis of substituted guanidines via the guanylation of an amine with carbodiimide, as well as the guanylation of phosphines and conversion of alkynes into propiolamidines, aluminum amides are only efficient catalysts for the guanylation of amides. Density functional theory (DFT) calculations were used to explain this difference in activity. The origin of this behavior is apparent in the critical step where a proton is transferred from the substrate to a metal guanidinate. The activation energies of these steps are modest for amines, phosphines, and alkynes when a lithium catalyst was used, but are prohibitively high for the analogous reactions with phosphines and alkynes for aluminum amide catalysts. Energy decomposition analysis (EDA) indicates that these high activations energies are due to the high energetic cost of the detachment of a chelating guanidinate nitrogen from the aluminum in the proton transfer transition state. Amines are able to adopt an ideal geometry for facile proton transfer to the aluminum guanidinate and concomitant Al-N bond formation, while phosphines and alkynes are not.

The CHARMM–TURBOMOLE interface for efficient and accurate QM/MM molecular dynamics, free energies, and excited state properties

The quantum mechanical (QM)/molecular mechanical (MM) interface between Chemistry at HARvard Molecular Mechanics (CHARMM) and TURBOMOLE is described. CHARMM provides an extensive set of simulation algorithms, like molecular dynamics (MD) and free energy perturbation, and support for mature nonpolarizable and Drude polarizable force fields. TURBOMOLE provides fast QM calculations using density functional theory or wave function methods and excited state properties. CHARMM–TURBOMOLE is well‐suited for extended QM/MM MD simulations using first principles methods with large (triple‐ζ) basis sets. We demonstrate these capabilities with a QM/MM simulation of Mg2+(aq), where the MM outer sphere water molecules are represented using the SWM4‐NDP Drude polarizable force field and the ion and inner coordination sphere are represented using QM PBE, PBE0, and MP2 methods. The relative solvation free energies of Mg2+ and Zn2+ were calculated using thermodynamic integration. We also demonstrate the features for excited state properties. We calculate the time‐averaged solution absorption spectrum of indole, the emission spectrum of the indole 1La excited state, and the electronic circular dichroism spectrum of an oxacepham.