Xiaolin Cheng

Characterizing Loop Dynamics and Ligand Recognition in Human- and Avian-Type Influenza Neuraminidases via Generalized Born Molecular Dynamics and End-Point Free Energy Calculations

The comparative dynamics and inhibitor binding free energies of group-1 and group-2 pathogenic influenza A subtype neuraminidase (NA) enzymes are of fundamental biological interest and relevant to structure-based drug design studies for antiviral compounds. In this work, we present seven generalized Born molecular dynamics simulations of avian (N1)- and human (N9)-type NAs in order to probe the comparative flexibility of the two subtypes, both with and without the inhibitor oseltamivir bound. The enhanced sampling obtained through the implicit solvent treatment suggests several provocative insights into the dynamics of the two subtypes, including that the group-2 enzymes may exhibit similar motion in the 430-binding site regions but different 150-loop motion. End-point free energy calculations elucidate the contributions to inhibitor binding free energies and suggest that entropic considerations cannot be neglected when comparing across the subtypes. We anticipate the findings presented here will have broad implications for the development of novel antiviral compounds against both seasonal and pandemic influenza strains.

Intramembrane Proton Binding Site Linked to Activation of Bacterial Pentameric Ion Channel

Background: Bacterial pentameric ion channels provide model systems to delineate structures and mechanisms of eukaryotic Cys-loop receptor channels. Results: In the bacterial channel GLIC, histidine protonation promotes hydrogen bonding between transmembrane (TM) helices. Conclusion: Protonation of intramembrane binding site promotes ion channel opening in GLIC. Significance: Interhelical hydrogen bonds promote concerted tilting of TM helices associated with the open channel conformation. The TM site may contribute to phototaxis in bacteria and allosteric potentiation in eukaryotes. Prokaryotic orthologs of eukaryotic Cys-loop receptor channels recently emerged as structural and mechanistic surrogates to investigate this superfamily of intercellular signaling proteins. Here, we examine proton activation of the prokaryotic ortholog GLIC using patch clamp electrophysiology, mutagenesis, and molecular dynamics (MD) simulations. Whole-cell current recordings from human embryonic kidney (HEK) 293 cells expressing GLIC show half-maximal activation at pH 6, close to the pKa of histidine, implicating the three native His residues in proton sensing linked to activation. The mutation H235F abolishes proton activation, H277Y is without effect, and all nine mutations of His-127 prevent expression on the cell surface. In the GLIC crystal structure, His-235 on transmembrane (TM) α-helix 2, hydrogen bonds to the main chain carbonyl oxygen of Ile-259 on TM α-helix 3. MD simulations show that when His-235 is protonated, the hydrogen bond persists, and the channel remains in the open conformation, whereas when His-235 is deprotonated, the hydrogen bond dissociates, and the channel closes. Mutations of the proximal Tyr-263, which also links TM α-helices 2 and 3 via a hydrogen bond, alter proton sensitivity over a 1.5 pH unit range. MD simulations show that mutations of Tyr-263 alter the hydrogen bonding capacity of His-235. The overall findings show that His-235 in the TM region of GLIC is a novel proton binding site linked to channel activation.

The molecular structure of a phosphatidylserine bilayer determined by scattering and molecular dynamics simulations.

Phosphatidylserine (PS) lipids play essential roles in biological processes, including enzyme activation and apoptosis. We report on the molecular structure and atomic scale interactions of a fluid bilayer composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS). A scattering density profile model, aided by molecular dynamics (MD) simulations, was developed to jointly refine different contrast small-angle neutron and X-ray scattering data, which yielded a lipid area of 62.7 Å(2) at 25 °C. MD simulations with POPS lipid area constrained at different values were also performed using all-atom and aliphatic united-atom models. The optimal simulated bilayer was obtained using a model-free comparison approach. Examination of the simulated bilayer, which agrees best with the experimental scattering data, reveals a preferential interaction between Na(+) ions and the terminal serine and phosphate moieties. Long-range inter-lipid interactions were identified, primarily between the positively charged ammonium, and the negatively charged carboxylic and phosphate oxygens. The area compressibility modulus KA of the POPS bilayer was derived by quantifying lipid area as a function of surface tension from area-constrained MD simulations. It was found that POPS bilayers possess a much larger KA than that of neutral phosphatidylcholine lipid bilayers. We propose that the unique molecular features of POPS bilayers may play an important role in certain physiological functions.

An Adaptive Fast Multipole Boundary Element Method for Poisson-Boltzmann Electrostatics

The numerical solution of the Poisson−Boltzmann (PB) equation is a useful but a computationally demanding tool for studying electrostatic solvation effects in chemical and biomolecular systems. Recently, we have described a boundary integral equation-based PB solver accelerated by a new version of the fast multipole method (FMM). The overall algorithm shows an order N complexity in both the computational cost and memory usage. Here, we present an updated version of the solver by using an adaptive FMM for accelerating the convolution type matrix-vector multiplications. The adaptive algorithm, when compared to our previous nonadaptive one, not only significantly improves the performance of the overall memory usage but also remarkably speeds the calculation because of an improved load balancing between the local- and far-field calculations. We have also implemented a node-patch discretization scheme that leads to a reduction of unknowns by a factor of 2 relative to the constant element method without sacrificing accuracy. As a result of these improvements, the new solver makes the PB calculation truly feasible for large-scale biomolecular systems such as a 30S ribosome molecule even on a typical 2008 desktop computer.

AFMPB: An adaptive fast multipole Poisson Boltzmann solver for calculating electrostatics in biomolecular systems

A Fortran program package is introduced for rapid evaluation of the electrostatic potentials and forces in biomolecular systems modeled by the linearized Poisson-Boltzmann equation. The numerical solver utilizes a well-conditioned boundary integral equation (BIE) formulation, a node-patch discretization scheme, a Krylov subspace iterative solver package with reverse communication protocols, and an adaptive new version of fast multipole method in which the exponential expansions are used to diagonalize the multipole to local translations. The program and its full description, as well as several closely related libraries and utility tools are available at http://lsec.cc.ac.cn/lubz/afmpb.html and a mirror site at http://mccammon.ucsd.edu/. This paper is a brief summary of the program: the algorithms, the implementation and the usage.

Mechanical Properties of Nanoscopic Lipid Domains

The lipid raft hypothesis presents insights into how the cell membrane organizes proteins and lipids to accomplish its many vital functions. Yet basic questions remain about the physical mechanisms that lead to the formation, stability, and size of lipid rafts. As a result, much interest has been generated in the study of systems that contain similar lateral heterogeneities, or domains. In the current work we present an experimental approach that is capable of isolating the bending moduli of lipid domains. This is accomplished using neutron scattering and its unique sensitivity to the isotopes of hydrogen. Combining contrast matching approaches with inelastic neutron scattering, we isolate the bending modulus of ∼13 nm diameter domains residing in 60 nm unilamellar vesicles, whose lipid composition mimics the mammalian plasma membrane outer leaflet. Importantly, the bending modulus of the nanoscopic domains differs from the modulus of the continuous phase surrounding them. From additional structural measurements and all-atom simulations, we also determine that nanoscopic domains are in-register across the bilayer leaflets. Taken together, these results inform a number of theoretical models of domain/raft formation and highlight the fact that mismatches in bending modulus must be accounted for when explaining the emergence of lateral heterogeneities in lipid systems and biological membranes.

A Solvent-Free Coarse Grain Model for Crystalline and Amorphous Cellulose Fibrils.

Understanding biomass structure and dynamics on a range of time and length scales is important for the development of cellulosic biofuels. Here, to enable length and time scale extension, we develop a coarse grain (CG) model for molecular dynamics (MD) simulations of cellulose. For this purpose, we use distribution functions from fully atomistic MD simulations as target observables. A single bead per monomer level coarse graining is found to be sufficient to successfully reproduce structural features of crystalline cellulose. Without the use of constraints the CG crystalline fibril is found to remain stable over the maximum simulation length explored in this study (>1 μs). We also extend the CG representation to model fully amorphous cellulose fibrils. This is done by using an atomistic MD simulation of fully solvated individual cellulose chains as a target for developing the corresponding fully amorphous CG force field. Fibril structures with different degrees of crystallinity are obtained using force fields derived using a parameter coupling the crystalline and amorphous potentials. The method provides an accurate and constraint-free approach to derive CG models for cellulose with a wide range of crystallinity, suitable for incorporation into large-scale models of lignocellulosic biomass.

Molecular-Dynamics Simulations of ELIC—a Prokaryotic Homologue of the Nicotinic Acetylcholine Receptor

The ligand-gated ion channel from Erwinia chrysanthemi (ELIC) is a prokaryotic homolog of the eukaryotic nicotinic acetylcholine receptor (nAChR) that responds to the binding of neurotransmitter acetylcholine and mediates fast signal transmission. ELIC is similar to the nAChR in its primary sequence and overall subunit organization, but despite their structural similarity, it is not clear whether these two ligand-gated ion channels operate in a similar manner. Further, it is not known to what extent mechanistic insights gleaned from the ELIC structure translate to eukaryotic counterparts such as the nAChR. Here we use molecular-dynamics simulations to probe the conformational dynamics and hydration of the transmembrane pore of ELIC. The results are compared with those from our previous simulation of the human alpha7 nAChR. Overall, ELIC displays increased stability compared to the nAChR, whereas the two proteins exhibit remarkable similarity in their global motion and flexibility patterns. The majority of the increased stability of ELIC does not stem from the deficiency of the models used in the simulations, and but rather seems to have a structural basis. Slightly altered dynamical correlation features are also observed among several loops within the membrane region. In sharp contrast to the nAChR, ELIC is completely dehydrated from the pore center to the extracellular end throughout the simulation. Finally, the simulation of an ELIC mutant substantiates the important role of F246 on the stability, hydration and possibly function of the ELIC channel.

Simulation of a cellulose fiber in ionic liquid suggests a synergistic approach to dissolution

Ionic liquids dissolve cellulose in a more efficient and environmentally acceptable way than conventional methods in aqueous solution. An understanding of how ionic liquids act on cellulose is essential for improving pretreatment conditions and thus detailed knowledge of the interactions between the cations, anions and cellulose is necessary. Here, to explore ionic liquid effects, we perform all-atom molecular dynamics simulations of a cellulose microfibril in 1-butyl-3-methylimidazolium chloride and analyze site–site interactions and cation orientations at the solute–solvent interface. The results indicate that Cl− anions predominantly interact with cellulose surface hydroxyl groups but with differences between chains of neighboring cellulose layers, referred to as center and origin chains; Cl− binds to C3-hydroxyls on the origin chains but to C2- and C6-hydroxyls on the center chains, thus resulting in a distinct pattern along glucan chains of the hydrophilic fiber surfaces. In particular, Cl− binding disrupts intrachain O3H–O5 hydrogen bonds on the origin chains but not those on the center chains. In contrast, Bmim+ cations stack preferentially on the hydrophobic cellulose surface, governed by non-polar interactions with cellulose. Complementary to the polar interactions between Cl− and cellulose, the stacking interaction between solvent cation rings and cellulose pyranose rings can compensate the interaction between stacked cellulose layers, thus stabilizing detached cellulose chains. Moreover, a frequently occurring intercalation of Bmim+ on the hydrophilic surface is observed, which by separating cellulose layers can also potentially facilitate the initiation of fiber disintegration. The results provide a molecular description why ionic liquids are ideal cellulose solvents, the concerted action of anions and cations on the hydrophobic and hydrophilic surfaces being key to the efficient dissolution of the amphiphilic carbohydrate.

Lateral organization, bilayer asymmetry, and inter-leaflet coupling of biological membranes

Understanding of cell membrane organization has evolved significantly from the classic fluid mosaic model. It is now recognized that biological membranes are highly organized structures, with differences in lipid compositions between inner and outer leaflets and in lateral structures within the bilayer plane, known as lipid rafts. These organizing principles are important for protein localization and function as well as cellular signaling. However, the mechanisms and biophysical basis of lipid raft formation, structure, dynamics and function are not clearly understood. One key question, which we focus on in this review, is how lateral organization and leaflet compositional asymmetry are coupled. Detailed information elucidating this question has been sparse because of the small size and transient nature of rafts and the experimental challenges in constructing asymmetric bilayers. Resolving this mystery will require advances in both experimentation and modeling. We discuss here the preparation of model systems along with experimental and computational approaches that have been applied in efforts to address this key question in membrane biology. We seek to place recent and future advances in experimental and computational techniques in context, providing insight into in-plane and transverse organization of biological membranes.