Tatyana Mamonova | PubMed Health

PubMed2006

Structure and energetics of channel-forming protein-polysaccharide complexes inferred via computational statistical thermodynamics.

The ion channel protein alpha-hemolysin (alphaHL) forms supramolecular complexes with the polysaccharide beta-cyclodextrin (betaCD). This system has potential uses in nanoscale device engineering. It has been found recently that betaCD formed longer- or shorter-lived complexes with some engineered alphaHL mutants then with a wild type protein (Gu et al. J. Gen. Physiol. 2001, 118, 481-493). However, how changes in the protein sequence affect complex lifetime was not completely understood in part due to the lack of knowledge of structures of these metastable complexes. In this paper, we present an extensive molecular modeling study of the betaCD-alphaHL and selected mutant complexes to gain insights into the betaCD-alphaHL interaction mechanisms and to predict possible structures and energetics of the complexes. Thermodynamic integration (TI) and umbrella sampling (US) techniques (with the weighted histogram analysis method (WHAM)) were used to calculate the relative binding affinities of the complexes formed with the wild type alphaHL and the M113N, M113E, M113A, and M113V mutants. Our results are in excellent agreement with experiment. While betaCD-M113N and betaCD-M113A complexes were stable in the configuration of the wild type complex, the equilibrium configuration of the betaCD-M113V and betaCD-M113E complexes was significantly different. In these cases, TI alone was insufficient to accurately calculate the corresponding free energy differences. By utilizing a TI/US combination in a novel manner, we were able to accurately calculate free energy changes in these flexible systems. The betaCD-M113A and betaCD-M113E complexes, which exhibited shorter lifetimes than other complexes in an experiment, in simulations exhibited greater flexibility and higher water solvation of the betaCD adapter. MD simulations of the betaCD-M113N complex with betaCD in a downward orientation were also performed.

Bacterial Proteins↗

PubMed2005

Protein flexibility using constraints from molecular dynamics simulations.

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.

Algorithms↗