Oliver S. Smart

HOLE: A program for the analysis of the pore dimensions of ion channel structural models

A method (HOLE) that allows the analysis of the dimensions of the pore running through a structural model of an ion channel is presented. The algorithm uses a Monte Carlo simulated annealing procedure to find the best route for a sphere with variable radius to squeeze through the channel. Results can be displayed in a graphical fashion or visualized with most common molecular graphical packages. Advances include a method to analyze the anisotropy within a pore. The method can also be used to predict the conductance of channels using a simple empirically corrected ohmic model. As an example the program is applied to the cholera toxin B-subunit pentamer. The compatibility of the crystal structure and conductance data is established.

The pore dimensions of gramicidin A.

The ion channel forming peptide gramicidin A adopts a number of distinct conformations in different environments. We have developed a new method to analyze and display the pore dimensions of ion channels. The procedure is applied to two x-ray crystal structures of gramicidin that adopt distinct antiparallel double helical dimer conformations and a nuclear magnetic resonance (NMR) structure for the beta6.3 NH2-terminal to NH2-terminal dimer. The results are discussed with reference to ion conductance properties and dependence of pore dimensions on the environment.

A novel method for structure-based prediction of ion channel conductance properties.

A rapid and easy-to-use method of predicting the conductance of an ion channel from its three-dimensional structure is presented. The method combines the pore dimensions of the channel as measured in the HOLE program with an Ohmic model of conductance. An empirically based correction factor is then applied. The method yielded good results for six experimental channel structures (none of which were included in the training set) with predictions accurate to within an average factor of 1.62 to the true values. The predictive r2 was equal to 0.90, which is indicative of a good predictive ability. The procedure is used to validate model structures of alamethicin and phospholamban. Two genuine predictions for the conductance of channels with known structure but without reported conductances are given. A modification of the procedure that calculates the expected results for the effect of the addition of nonelectrolyte polymers on conductance is set out. Results for a cholera toxin B-subunit crystal structure agree well with the measured values. The difficulty in interpreting such studies is discussed, with the conclusion that measurements on channels of known structure are required.

A new method to calculate reaction paths for conformation transitions of large molecules

Abstract Path energy minimization (PEM), a novel method for the generation of a reaction path linking two known conformers of a molecule, is presented. The technique is based on optimizing a function which closely approximates the peak potential energy of a quasi-continuous path between the fixed end points. A transition involving the change in the pucker angle of α-D-xylulofuranose is used as a test case. The method is shown to, be capable of identifying transition state structures and energy barries. The utility of the method is demonstrated by an application to substantial conformational transition of the ion-channel forming polypeptide gramicidin A.

Channels formed by the transmembrane helix of phospholamban: a simulation study

Phospholamban is a small membrane protein which can form cation selective ion channels in lipid bilayers. Each subunit contains a single, largely hydrophobic transmembrane helix. The helices are thought to assemble as a pentameric and approximately parallel bundle surrounding a central pore. A model of this assembly (PDB code IPSL) has been used as the starting point for molecular dynamics (MD) simulations of a system consisting of the pentameric helix bundle, plus 217 water molecules located within and at either mouth of the pore. Interhelix distance restraints were employed to maintain the integrity of the helix bundle during a 500 ps MD simulation. Water molecules within the pore exhibited reduced diffusional and rotational mobility. Interactions between the alpha-helix dipoles and the water dipoles, the latter aligned anti-parallel to the former, contribute to the stability of the system. Analysis of the potential energy of interaction of a K+ ion as it was moved through the pore suggested that unfavourable interactions of the cation with the aligned helix dipoles at the N-terminal mouth were overcome by favourable ion-water interactions. Comparable analysis for a Cl ion revealed that the ion-(pore + water) interactions were unfavourable along the whole of the pore, increasingly so from the N- to the C-terminal mouth. Overall, the interaction energy profiles were consistent with a pore selective for cations over anions. Pore radius profiles were used to predict a channel conductance of 50 to 70 ps in 0.2 M KCl, which compares well with an experimental value of 100 ps.

On the Simulation of Conformational Transitions: Smoothing Path Energy Minimization Results

Abstract Path energy minimization is a method for calculating a quasi-continuous reaction path linking two known conformers of a molecule. Although the method locates the transition state configuration lower energy regions of the path are not optimized. This work examines four methods of following routes downhill from the transition state. Results are compared for the test system of a pucker angle change for the sugar α-D-xylulofuranose. Restrained energy minimization is found to be the superior method. The steepest descents path obtained by this method is shown to be distinct from the route found by adiabatic mapping.

New methods for the analysis of the protein-solvent interface

Abstract The protein-solvent interface is complex and may include solvent channels and cavities as well as the normal surface water molecules. We describe several algorithms for investigating the intra- and inter-molecular interactions of proteins in general but with the aim of developing methods to accurately and definitively characterise the interactions of water and other small ligands with proteins. Specifically, we present the methods which underlie three programs (AQUARIUS2, HOLE and PRO_ACT) which can be used to to look at different aspects of these interactions.