Src Kinase Conformational Activation: Thermodynamics, Pathways, and Mechanisms

Tyrosine kinases of the Src-family are large allosteric enzymes that play a key role in cellular signaling. Conversion of the kinase from an inactive to an active state is accompanied by substantial structural changes. Here, we construct a coarse-grained model of the catalytic domain incorporating experimental structures for the two stable states, and simulate the dynamics of conformational transitions in kinase activation. We explore the transition energy landscapes by constructing a structural network among clusters of conformations from the simulations. From the structural network, two major ensembles of pathways for the activation are identified. In the first transition pathway, we find a coordinated switching mechanism of interactions among the αC helix, the activation-loop, and the β strands in the N-lobe of the catalytic domain. In a second pathway, the conformational change is coupled to a partial unfolding of the N-lobe region of the catalytic domain. We also characterize the switching mechanism for the αC helix and the activation-loop in detail. Finally, we test the performance of a Markov model and its ability to account for the structural kinetics in the context of Src conformational changes. Taken together, these results provide a broad framework for understanding the main features of the conformational transition taking place upon Src activation

Brownian dynamics simulations of ions channels: A general treatment of electrostatic reaction fields for molecular pores of arbitrary geometry

Copyright 2001 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in the Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1390507.A general method has been developed to include the electrostatic reaction field in Brownian dynamics (BD) simulations of ions diffusing through complex molecular channels of arbitrary geometry. Assuming that the solvent is represented as a featureless continuum dielectric medium, a multipolar basis-set expansion is developed to express the reaction field Green’s function. A reaction field matrix, which provides the coupling between generalized multipoles, is calculated only once and stored before the BD simulations. The electrostatic energy and forces are calculated at each time step by updating the generalized multipole moments. The method is closely related to the generalized solvent boundary potential [Im et al., J. Chem. Phys. 114, 2924 (2001)] which was recently developed to include the influence of distant atoms on a small region part of a large macromolecular system in molecular dynamics simulations. It is shown that the basis-set expansion is accurate and computationally inexpensive for three simple models such as a spherical ionic system, an impermeable membrane system, and a cylindrical pore system as well as a realistic system such as OmpF porin with all atomic details. The influence of the static field and the reaction field on the ion distribution and conductance in the OmpF channel is studied and discussed

A Grand Canonical Monte Carlo–Brownian Dynamics Algorithm for Simulating Ion Channels

This is the published version. Copyright 2000 by Elsevier B. V.A computational algorithm based on Grand Canonical Monte Carlo (GCMC) and Brownian Dynamics (BD) is described to simulate the movement of ions in membrane channels. The proposed algorithm, GCMC/BD, allows the simulation of ion channels with a realistic implementation of boundary conditions of concentration and transmembrane potential. The method is consistent with a statistical mechanical formulation of the equilibrium properties of ion channels (Roux, B. 1999; Biophys. J. 77:139–153). The GCMC/BD algorithm is illustrated with simulations of simple test systems and of the OmpF porin of Escherichia coli. The approach provides a framework for simulating ion permeation in the context of detailed microscopic models

Statistical Parsing of Spanish and Data Driven Lemmatization

International audienceAlthough parsing performances have greatly improved in the last years, grammar inference from treebanks for morphologically rich lan- guages, especially from small treebanks, is still a challenging task. In this paper we in- vestigate how state-of-the-art parsing perfor- mances can be achieved on Spanish, a lan- guage with a rich verbal morphology, with a non-lexicalized parser trained on a treebank containing only around 2,800 trees. We rely on accurate part-of-speech tagging and data- driven lemmatization in order to cope with lexical data sparseness. Providing state-of- the-art results on Spanish, our methodology is applicable to other languages

Ion Permeation through the α-Hemolysin Channel: Theoretical Studies Based on Brownian Dynamics and Poisson-Nernst-Plank Electrodiffusion Theory

This is the published version. Copyright 2004 by Elsevier.Identification of the molecular interaction governing ion conduction through biological pores is one of the most important goals of modern electrophysiology. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) and three-dimensional Poisson-Nernst-Plank (3d-PNP) electrodiffusion algorithms offer powerful and general approaches to study of ion permeation through wide molecular pores. A detailed analysis of ion flows through the staphylococcal α-hemolysin channel based on series of simulations at different concentrations and transmembrane potentials is presented. The position-dependent diffusion coefficient is approximated on the basis of a hydrodynamic model. The channel conductance calculated by GCMC/BD is ∼10% higher than (electrophysiologically measured) experimental values, whereas results from 3d-PNP are always 30–50% larger. Both methods are able to capture all important electrostatic interactions in equilibrium conditions. The asymmetric conductance upon the polarity of the transmembrane potential observed experimentally is reproduced by GCMC/BD and 3d-PNP. The separation of geometrical and energetic influence of the channel on ion conduction reveals that such asymmetries arise from the permanent charge distribution inside the pore. The major determinant of the asymmetry is unbalanced charge in the triad of polar residues D127, D128, and K131. The GCMC/BD or 3d-PNP calculations reproduce also experimental reversal potentials and permeability rations in asymmetric ionic solutions. The weak anionic selectivity of the channel results from the presence of the salt bridge between E111 and K147 in the constriction zone. The calculations also reproduce the experimentally derived dependence of the reversible potential to the direction of the salt gradient. The origin of such effect arises from the asymmetrical distribution of energetic barriers along the channel axis, which modulates the preferential ion passage in different directions

Architecture and assembly of the Gram-positive cell wall

The bacterial cell wall is a mesh polymer of peptidoglycan – linear glycan strands cross-linked by flexible peptides – that determines cell shape and provides physical protection. While the glycan strands in thin ‘Gram-negative’ peptidoglycan are known to run circumferentially around the cell, the architecture of the thicker ‘Gram-positive’ form remains unclear. Using electron cryotomography, here we show that Bacillus subtilis peptidoglycan is a uniformly dense layer with a textured surface. We further show it rips circumferentially, curls and thickens at free edges, and extends longitudinally when denatured. Molecular dynamics simulations show that only atomic models based on the circumferential topology recapitulate the observed curling and thickening, in support of an ‘inside-to-outside’ assembly process. We conclude that instead of being perpendicular to the cell surface or wrapped in coiled cables (two alternative models), the glycan strands in Gram-positive cell walls run circumferentially around the cell just as they do in Gram-negative cells. Together with providing insights into the architecture of the ultimate determinant of cell shape, this study is important because Gram-positive peptidoglycan is an antibiotic target crucial to the viability of several important rod-shaped pathogens including Bacillus anthracis, Listeria monocytogenes, and Clostridium difficile

Long-pore Electrostatics in Inward-rectifier Potassium Channels

Inward-rectifier potassium (Kir) channels differ from the canonical K+ channel structure in that they possess a long extended pore (∼85 Å) for ion conduction that reaches deeply into the cytoplasm. This unique structural feature is presumably involved in regulating functional properties specific to Kir channels, such as conductance, rectification block, and ligand-dependent gating. To elucidate the underpinnings of these functional roles, we examine the electrostatics of an ion along this extended pore. Homology models are constructed based on the open-state model of KirBac1.1 for four mammalian Kir channels: Kir1.1/ROMK, Kir2.1/IRK, Kir3.1/GIRK, and Kir6.2/KATP. By solving the Poisson-Boltzmann equation, the electrostatic free energy of a K+ ion is determined along each pore, revealing that mammalian Kir channels provide a favorable environment for cations and suggesting the existence of high-density regions in the cytoplasmic domain and cavity. The contribution from the reaction field (the self-energy arising from the dielectric polarization induced by the ion's charge in the complex geometry of the pore) is unfavorable inside the long pore. However, this is well compensated by the electrostatic interaction with the static field arising from the protein charges and shielded by the dielectric surrounding. Decomposition of the static field provides a list of residues that display remarkable correspondence with existing mutagenesis data identifying amino acids that affect conduction and rectification. Many of these residues demonstrate interactions with the ion over long distances, up to 40 Å, suggesting that mutations potentially affect ion or blocker energetics over the entire pore. These results provide a foundation for understanding ion interactions in Kir channels and extend to the study of ion permeation, block, and gating in long, cation-specific pores