Bogdan Lev

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

Relative Free Energies for Hydration of Monovalent Ions from QM and QM/MM Simulations

Methods directly evaluating the hydration structure and thermodynamics of physiologically relevant cations (Na(+), K(+), Cl(-), etc.) have wide ranging applications in the fields of inorganic, physical, and biological chemistry. All-atom simulations based on accurate potential energy surfaces appear to offer a viable option for assessing the chemistry of ion solvation. Although MD and free energy simulations of ion solvation with classical force fields have proven their usefulness, a number of challenges still remain. One of them is the difficulty of force field benchmarking and validation against structural and thermodynamic data obtained for a condensed phase. Hybrid quantum mechanical/molecular mechanical (QM/MM) models combined with sampling algorithms have the potential to provide an accurate solvation model and to incorporate the effects from the surrounding, which is often missing in gas-phase ab initio computations. Herein, we report the results from QM/MM free energy simulations of Na(+)/K(+) and Cl(-)/Br(-) hydration where we simultaneously characterized the relative thermodynamics of ion solvation and changes in the solvation structure. The Flexible Inner Region Ensemble Separator (FIRES) method was used to impose a spatial separation between QM region and the outer sphere of solvent molecules treated with the CHARMM27 force field. FEP calculations based on QM/MM simulations utilizing the CHARMM/deMon2k interface were performed with different basis set combinations for K(+)/Na(+) and Cl(-)/Br(-) perturbations to establish the dependence of the computed free energies on the basis set level. The dependence of the computed relative free energies on the size of the QM and MM regions is discussed. The current methodology offers an accurate description of structural and thermodynamic aspects of the hydration of alkali and halide ions in neat solvents and can be used to obtain thermodynamic data on ion solvation in condensed phase along with underlying structural properties of the ion-solvent system.

String method solution of the gating pathways for a pentameric ligand-gated ion channel

Significance High-resolution structures of pentameric ligand-gated ion channels have created an opportunity to discover the mechanisms of rapid synaptic transduction in the brain. This study describes the mechanisms of allosteric channel gating using string method simulations, applied to a complete atomistic ion channel, combined with a transition analysis approach to extract free-energy surfaces from swarms of trajectories. We reproduce pH-modulated activity of the channel, identify the molecular interactions associated with interdomain communication, and quantify the energetics of the gating process. These results provide general mechanistic understanding of the function of pentameric ligand-gated channels, with potential applications in the design of improved anesthetics, neuromodulatory drugs, antiparasitics, and pesticides. Pentameric ligand-gated ion channels control synaptic neurotransmission by converting chemical signals into electrical signals. Agonist binding leads to rapid signal transduction via an allosteric mechanism, where global protein conformational changes open a pore across the nerve cell membrane. We use all-atom molecular dynamics with a swarm-based string method to solve for the minimum free-energy gating pathways of the proton-activated bacterial GLIC channel. We describe stable wetted/open and dewetted/closed states, and uncover conformational changes in the agonist-binding extracellular domain, ion-conducting transmembrane domain, and gating interface that control communication between these domains. Transition analysis is used to compute free-energy surfaces that suggest allosteric pathways; stabilization with pH; and intermediates, including states that facilitate channel closing in the presence of an agonist. We describe a switching mechanism that senses proton binding by marked reorganization of subunit interface, altering the packing of β-sheets to induce changes that lead to asynchronous pore-lining M2 helix movements. These results provide molecular details of GLIC gating and insight into the allosteric mechanisms for the superfamily of pentameric ligand-gated channels.

Identification of Electric-Field-Dependent Steps in the Na+,K+-Pump Cycle

The charge-transporting activity of the Na(+),K(+)-ATPase depends on its surrounding electric field. To isolate which steps of the enzymes reaction cycle involve charge movement, we have investigated the response of the voltage-sensitive fluorescent probe RH421 to interaction of the protein with BTEA (benzyltriethylammonium), which binds from the extracellular medium to the Na(+),K(+)-ATPases transport sites in competition with Na(+) and K(+), but is not occluded within the protein. We find that only the occludable ions Na(+), K(+), Rb(+), and Cs(+) cause a drop in RH421 fluorescence. We conclude that RH421 detects intramembrane electric field strength changes arising from charge transport associated with conformational changes occluding the transported ions within the protein, not the electric fields of the bound ions themselves. This appears at first to conflict with electrophysiological studies suggesting extracellular Na(+) or K(+) binding in a high field access channel is a major electrogenic reaction of the Na(+),K(+)-ATPase. All results can be explained consistently if ion occlusion involves local deformations in the lipid membrane surrounding the protein occurring simultaneously with conformational changes necessary for ion occlusion. The most likely origin of the RH421 fluorescence response is a change in membrane dipole potential caused by membrane deformation.

QM/MM Calculations with deMon2k

The density functional code deMon2k employs a fitted density throughout (Auxiliary Density Functional Theory), which offers a great speed advantage without sacrificing necessary accuracy. Powerful Quantum Mechanical/Molecular Mechanical (QM/MM) approaches are reviewed. Following an overview of the basic features of deMon2k that make it efficient while retaining accuracy, three QM/MM implementations are compared and contrasted. In the first, deMon2k is interfaced with the CHARMM MM code (CHARMM-deMon2k); in the second MM is coded directly within the deMon2k software; and in the third the Chemistry in Ruby (Cuby) wrapper is used to drive the calculations. Cuby is also used in the context of constrained-DFT/MM calculations. Each of these implementations is described briefly; pros and cons are discussed and a few recent applications are described briefly. Applications include solvated ions and biomolecules, polyglutamine peptides important in polyQ neurodegenerative diseases, copper monooxygenases and ultra-rapid electron transfer in cryptochromes.

Halothane Solvation in Water and Organic Solvents from Molecular Simulations with New Polarizable Potential Function

The partitioning of a substrate from one phase into another is a complex process with widespread applications: from chemical technology to the pharmaceutical industry. One particularly well-known and well-studied example is 2-bromo-2-chloro-1,1,1-trifluoroethane (halothane) trafficking through the lipid bilayer. Halothane is a model volatile anesthetic known to impact functions of model lipid bilayers, altering the structure and thickness upon its partitioning from the bulk phase. A number of theoretical and experimental investigations suggest the importance of electronic polarizability, determining a preference for halothane to partition in the interfacial systems as in lipid bilayers or binary solvents. The recently published protocol for the development of polarizable force fields based on the classical Drude model has provided fresh impetus to efforts directed at understanding the molecular principles governing complex thermodynamics of the hydrophobic hydration. Here, molecular simulations were combined with free energy simulations to study solvation of halothane in polarizable water and methanol. The absolute free energy of halothane solvation in different solvents (water, methanol, and n-hexane) has been evaluated for additive and polarizable models. It was found that both additive and polarizable models provide an adequate description of the halothane solvation in high-dielectric (polar) solvents such as water, but explicit accounting for electronic polarization is imperative for a correct description of the solvation thermodynamics in nonpolar systems. To study halothane dynamics in binary mixtures, all-atom molecular dynamics (MD) simulations for halothane-methanol mixtures in a wide range of concentrations were performed alongside an analysis of structural organization, dynamics, and thermodynamic properties to dissect the molecular determinants of the halothane solvation in polar and amphiphilic liquids such as methanol. Additionally, a theoretical test of the hypothesis on the weak hydrogen bonding of halothane and methanol in the condensed phase is provided, which was presented on the basis of spectroscopic analysis of the C-H vibrations in different gas-phase complexes. The simulations performed in the condensed phase suggest that hydrophobic interactions between halothane and methanol play a dominant role in preferential solvation.

The voltage-sensitive dye RH421 detects a Na+,K+-ATPase conformational change at the membrane surface

RH421 is a voltage-sensitive fluorescent styrylpyridinium dye which has often been used to probe the kinetics of Na+,K+-ATPase partial reactions. The origin of the dyes response has up to now been unclear. Here we show that RH421 responds to phosphorylation of the Na+,K+-ATPase by inorganic phosphate with a fluorescence increase. Analysis of the kinetics of the fluorescence response indicates that the probe is not detecting phosphorylation itself but rather a shift in the proteins E1/E2 conformational equilibrium induced by preferential phosphate binding to and phosphorylation of enzyme in the E2 conformation. Molecular dynamics simulations of crystal structures in lipid bilayers indicate some change in the proteins hydrophobic thickness during the E1-E2 transition, which may influence the dye response. However, the transition is known to involve significant rearrangement of the proteins highly charged lysine-rich cytoplasmic N-terminal sequence. Using poly-l-lysine as a model of the N-terminus, we show that an analogous response of RH421 to the E1→E2P conformational change is produced by poly-l-lysine binding to the surface of the Na+,K+-ATPase-containing membrane fragments. Thus, it seems that the prime origin of the RH421 fluorescence response is a change in the interaction of the proteins N-terminus with the surrounding membrane. Quantum mechanical calculations of the dyes visible absorption spectrum give further support to this conclusion. The results obtained indicate that membrane binding and release of the N-terminus of the Na+,K+-ATPase α-subunit are intimately involved in the proteins catalytic cycle and could represent an effective site of regulation.