The role of bulk protein in local models of ion-binding to proteins: comparative study of KcsA, its semisynthetic analog with a locked-in binding site, and valinomycin

Mini-grand canonical ensemble: Chemical potential in the solvation shell

Quantifying the statistics of occupancy of solvent molecules in the vicinity of solutes is central to our understanding of solvation phenomena. Number fluctuations in small solvation shells around solutes cannot be described within the macroscopic grand canonical framework using a single chemical potential that represents the solvent bath. In this communication, we hypothesize that molecular-sized observation volumes such as solvation shells are best described by coupling the solvation shell with a mixture of particle baths each with its own chemical potential. We confirm our hypotheses by studying the enhanced fluctuations in the occupancy statistics of hard sphere solvent particles around a distinguished hard sphere solute particle. Connections with established theories of solvation are also discussed.

Detecting temperature fluctuations at equilibrium

Equilibrium and dynamical properties of small systems are captured accurately when their temperature is allowed to vary.

K(+) and Na(+) conduction in selective and nonselective ion channels via molecular dynamics simulations

Generations of scientists have been captivated by ion channels and how they control the workings of the cell by admitting ions from one side of the cell membrane to the other. Elucidating the molecular determinants of ion conduction and selectivity are two of the most fundamental issues in the field of biophysics. Combined with ongoing progress in structural studies, modeling and simulation have been an integral part of the development of the field. As of this writing, the relentless growth in computational power, the development of new algorithms to tackle the so-called rare events, improved force-field parameters, and the concomitant increasing availability of membrane protein structures, allow simulations to contribute even further, providing more-complete models of ion conduction and selectivity in ion channels. In this report, we give an overview of the recent progress made by simulation studies on the understanding of ion permeation in selective and nonselective ion channels.

How does overcoordination create ion selectivity?

Some biological molecules can distinguish between ions of similar nature, which may be achieved by enforcing specific coordination numbers on ions in the binding site. It is suggested that when this number is favourable for one ion type, but too large for another, this creates ion selectivity through the proposed mechanism of 'overcoordination'. Much debate has occurred about the role overcoordination plays, and suggestions made as to how molecules can enforce particular coordination numbers, but there has not been an examination of the microscopic underpinning of ion selectivity by overcoordination. Here we use molecular-dynamics to systematically investigate how the number of ligands affects the ion-ligand and ligand-ligand interaction energies, and thus the thermodynamic ion selectivity, of a combination of model systems: three ions (Li(+)/Na(+)/K(+)) with three different ligands (water/formaldehyde/formamide). We find that the ligand-ligand repulsion controls the changes in geometry of each system with changing ligand number. Ion selectivity by overcoordination is achieved as smaller ions exhibit anomalous geometrical changes with the addition of extra ligands, whilst larger ions do not.

Separating the role of protein restraints and local metal-site interaction chemistry in the thermodynamics of a zinc finger protein

We express the effective Hamiltonian of an ion-binding site in a protein as a combination of the Hamiltonian of the ion-bound site in vacuum and the restraints of the protein on the site. The protein restraints are described by the quadratic elastic network model. The Hamiltonian of the ion-bound site in vacuum is approximated as a generalized Hessian around the minimum energy configuration. The resultant of the two quadratic Hamiltonians is cast into a pure quadratic form. In the canonical ensemble, the quadratic nature of the resultant Hamiltonian allows us to express analytically the excess free energy, enthalpy, and entropy of ion binding to the protein. The analytical expressions allow us to separate the roles of the dynamic restraints imposed by the protein on the binding site and the temperature-independent chemical effects in metal-ligand coordination. For the consensus zinc-finger peptide, relative to the aqueous phase, the calculated free energy of exchanging Zn(2+) with Fe(2+), Co(2+), Ni(2+), and Cd(2+) are in agreement with experiments. The predicted excess enthalpy of ion exchange between Zn(2+) and Co(2+) also agrees with the available experimental estimate. The free energy of applying the protein restraints reveals that relative to Zn(2+), the Co(2+), and Cd(2+)-site clusters are more destabilized by the protein restraints. This leads to an experimentally testable hypothesis that a tetrahedral metal binding site with minimal protein restraints will be less selective for Zn(2+) over Co(2+) and Cd(2+) compared to a zinc finger peptide. No appreciable change is expected for Fe(2+) and Ni(2+). The framework presented here may prove useful in protein engineering to tune metal selectivity.