PubMed Health⌕ Search

Biomedical subjects

Carol K Hall

Publications and source records attributed to Carol K Hall.

At least 19 recordsLinked to original sources
PubMed2007

Effects of chain length on the aggregation of model polyglutamine peptides: molecular dynamics simulations.

Aggregation in the brain of polyglutamine-containing proteins is either a cause or an associated symptom of nine hereditary neurodegenerative disorders including Huntington's disease. The molecular level mechanisms by which these proteins aggregate are still unclear. In an effort to shed light on this important phenomenon, we are investigating the aggregation of model polyglutamine peptides using molecular-level computer simulation with a simplified model of polyglutamine that we have developed. This model accounts for the most important types of intra- and inter-molecular interactions-hydrogen bonding and hydrophobic interactions-while allowing the folding process to be simulated in a reasonable time frame. The model is used to examine the folding of isolated polyglutamine peptides 16, 32, and 48 residues long and the folding and aggregation of systems of 24 model polyglutamine peptides 16, 24, 32, 36, 40, and 48 residues long. Although the isolated polyglutamine peptides did form some alpha and beta backbone-backbone hydrogen bonds they did not have as many of these bonds as they would have if they had folded into a complete alpha helix or beta sheet. In one of the simulations on the isolated polyglutamine peptide 48 residues long, we observed a structure that resembles a beta helix. In the multi-chain simulations we observed amorphous aggregates at low temperatures, ordered aggregates with significant beta sheet character at intermediate temperatures, and random coils at high temperatures. We have found that the temperature at which the model peptides undergo the transition from amorphous aggregates to ordered aggregates and the temperature at which the model peptides undergo the transition from ordered aggregates to random coils increase with increasing chain length. Our finding that the stability of the ordered aggregates increases as the peptide chain length increases may help to explain the experimentally observed relation between polyglutamine tract length and aggregation in vitro and disease progression in vivo. We have also observed in our simulations that the optimal temperature for the formation of beta sheets increases with chain length up to 36 glutamine residues but not beyond. Equivalently, at fixed temperature we find a transition from a region dominated by random coils at chain lengths less than 36 to a region dominated by relatively ordered beta sheet structures at chain lengths greater than 36. Our finding of this critical chain length of 36 glutamine residues is interesting because a critical chain length of 37 glutamine residues has been observed experimentally.

Computer Simulation↗
PubMed2006

Computer simulation study of molecular recognition in model DNA microarrays.

DNA microarrays have been widely adopted by the scientific community for a variety of applications. To improve the performance of microarrays there is a need for a fundamental understanding of the interplay between the various factors that affect microarray sensitivity and specificity. We use lattice Monte Carlo simulations to study the thermodynamics and kinetics of hybridization of single-stranded target genes in solution with complementary probe DNA molecules immobilized on a microarray surface. The target molecules in our system contain 48 segments and the probes tethered on a hard surface contain 8-24 segments. The segments on the probe and target are distinct and each segment represents a sequence of nucleotides ( approximately 11 nucleotides). Each probe segment interacts exclusively with its unique complementary target segment with a single hybridization energy; all other interactions are zero. We examine how the probe length, temperature, or hybridization energy, and the stretch along the target that the probe segments complement, affect the extent of hybridization. For systems containing single probe and single target molecules, we observe that as the probe length increases, the probability of binding all probe segments to the target decreases, implying that the specificity decreases. We observe that probes 12-16 segments ( approximately 132-176 nucleotides) long gave the highest specificity and sensitivity. This agrees with the experimental results obtained by another research group, who found an optimal probe length of 150 nucleotides. As the hybridization energy increases, the longer probes are able to bind all their segments to the target, thus improving their specificity. The hybridization kinetics reveals that the segments at the ends of the probe are most likely to start the hybridization. The segments toward the center of the probe remain bound to the target for a longer time than the segments at the ends of the probe.

Computer Simulation↗
PubMed2006

Spontaneous formation of annular structures observed in molecular dynamics simulations of polyglutamine peptides.

Annular structures have been observed experimentally in aggregates of polyglutamine-containing proteins and other proteins associated with diseases of the brain. Here we report the observation of annular structures in molecular-level simulations of large systems of model polyglutamine peptides. A system of 24 polyglutamine chains 16 residues long at a concentration of 5 mM spontaneously formed large beta sheets which curved to form tube-like annular structures that resemble beta barrels. This work was accomplished by extending the PRIME model to polyglutamine. PRIME is an off-lattice, unbiased, intermediate-resolution protein model based on an amino acid representation of between three and seven united atoms depending on the residue being modeled. Our results are interesting not only because of the recent discovery of tubular protofibrils in experiments on aggregation of mutant huntingtin fragments containing expanded polyglutamine tracts but also because Perutz predicted that polyglutamine forms water filled nanotubes.

Brain Diseases↗
PubMed2006

Side-chain interactions determine amyloid formation by model polyglutamine peptides in molecular dynamics simulations.

The pathological manifestation of nine hereditary neurodegenerative diseases is the presence within the brain of aggregates of disease-specific proteins that contain polyglutamine tracts longer than a critical length. To improve our understanding of the processes by which polyglutamine-containing proteins misfold and aggregate, we have conducted molecular dynamics simulations of the aggregation of model polyglutamine peptides. This work was accomplished by extending the PRIME model to polyglutamine. PRIME is an off-lattice, unbiased, intermediate-resolution protein model based on an amino acid representation of between three and seven united atoms, depending on the residue being modeled. The effects of hydrophobicity on the system are studied by varying the strength of the hydrophobic interaction from 12.5% to 5% of the hydrogen-bonding interaction strength. In our simulations, we observe the spontaneous formation of aggregates and annular structures that are made up of beta-sheets starting from random configurations of random coils. This result was interesting because tubular protofibrils were recently found in experiments on polyglutamine aggregation and because of Perutz's prediction that polyglutamine would form water-filled nanotubes.

Amino Acid Sequence↗
PubMed2006

Spontaneous fibril formation by polyalanines; discontinuous molecular dynamics simulations.

Fibrillary protein aggregates rich in beta-sheet structure have been implicated in the pathology of several neurodegenerative diseases. In this work, we investigate the formation of fibrils by performing discontinuous molecular dynamics simulations on systems containing 12 to 96 model Ac-KA(14)K-NH(2) peptides using our newly developed off-lattice, implicit-solvent, intermediate-resolution model, PRIME. We find that, at a low concentration, random-coil peptides assemble into alpha-helices at low temperatures. At intermediate concentrations, random-coil peptides assemble into alpha-helices at low temperatures and large beta-sheet structures at high temperatures. At high concentrations, the system forms beta-sheets over a wide range of temperatures. These assemble into fibrils above a critical temperature which decreases with concentration and exceeds the isolated peptide's folding temperature. At very high temperatures and all concentrations, the system is in a random-coil state. All of these results are in good qualitative agreement with those by Blondelle and co-workers on Ac-KA(14)K-NH(2) peptides. The fibrils observed in our simulations mimic the structural characteristics observed in experiments in terms of the number of sheets formed, the values of the intra- and intersheet separations, and the parallel peptide arrangement within each beta-sheet. Finally, we find that when the strength of the hydrophobic interaction between nonpolar side chains is high compared to the strength of hydrogen bonding, amorphous aggregates, rather than fibrillar aggregates, are formed.

Amyloid↗
PubMed2006

Computational approaches to fibril structure and formation.

Assembly of normally soluble proteins into amyloid fibrils is a cause or associated symptom of numerous human disorders. Although some progress toward understanding the molecular-level details of fibril structure has been made through in vitro experiments, the insoluble nature of fibrils make them difficult to study experimentally. We describe two computational approaches used to investigate fibril formation and structure: intermediate-resolution discontinuous molecular dynamics simulations and atomistic molecular dynamics simulations. Each method has its strengths and weaknesses, but taken together the two approaches provide a useful molecular-level picture of fibril structure and formation.

Amyloid↗
PubMed2005

Parametric studies of interaction strengths in polymer/CO2 systems: discontinuous molecular dynamics simulations.

Discontinuous molecular dynamics simulations are performed on homopolymer/solvent and surfactant/solvent systems. The homopolymer and surfactant molecules are modeled as freely jointed square-well chains. Solvent molecules are modeled as both hard spheres and square-well spheres. We explore how the various interaction parameters affect the types of phase behavior and micellization observed in the homopolymer/solvent and surfactant/solvent systems. Increasing the packing fraction of homopolymers in both hard-sphere solvents and square-well solvents increases the solvent's ability to dissolve homopolymers only when the segment-solvent interaction strength exceeds a critical value. Although only upper critical solution temperature (UCST) behavior is observed for homopolymers in hard-sphere solvents, both UCST and lower critical solution temperature (LCST) behavior are observed for homopolymers in square-well solvents, depending upon the interaction strengths and chain length. This indicates that it is necessary to account for the solvent-solvent attraction to model LCST behavior in supercritical CO2. Our simulation results on surfactants in hard-sphere solvents show that it is necessary to account for the interactions experienced by both the head and tail blocks in order to capture the essential features of surfactant/supercritical CO2 systems.

Journal Article↗
PubMed2005

The dynamics of single chains within a model polymer melt.

Discontinuous molecular dynamics simulations are performed on a system containing 32 hard chains of length 192 at a volume fraction of phi = 0.45 to explore the idea that localized entanglements have a significant effect on the dynamics of the individual chains within an entangled polymer melt. Anomalous behavior can still be observed when studying the dynamics of the individual chains, although increased time averaging causes the anomalous relaxation-memory-release behavior that was observed previously in the system to smooth out. First, the individual chain mean squared displacements and apparent diffusion coefficients are calculated, and a wide distribution of diffusive behavior is found. Although the apparent diffusion coefficient curve averaged over all chains displays the predicted long-time diffusive behavior, the curves for the individual chains differ both qualitatively and quantitatively. They display superdiffusive, diffusive, and subdiffusive behavior, with the largest percentage of chains exhibiting superdiffusive behavior and the smallest percentage exhibiting the predicted diffusive behavior. Next, the individual chain end-to-end vector autocorrelation functions and relaxation times are determined, and a wide distribution of stress relaxation behavior is found. The times when the end-to-end vector autocorrelation functions relax completely span almost an order of magnitude in reduced time. For some chains, the end-to-end vector autocorrelation function relaxes smoothly toward zero similar to the system average; however, for other chains the relaxation is slowed greatly, indicating the presence of additional entanglements. Almost half of the chains exhibit the anomalous behavior in the end-to-end vector autocorrelation function. Finally, the dynamic properties are displayed for a single chain exhibiting anomalous relaxation-memory-release behavior, supporting the idea that the relaxation-memory-release behavior is a single-chain property.

Journal Article↗
PubMed2005

Designing pattern-recognition surfaces for selective adsorption of copolymer sequences using lattice monte carlo simulation.

We describe a simulation method to design surfaces for recognizing specific monomer sequences in copolymers. We fix the monomer sequence statistics of the AB copolymers on a surface containing two types of sites and allow the simulation to iterate towards an optimal surface pattern that can recognize and selectively adsorb the sequence in the copolymer. During the simulation the surface pattern is designed by switching identities of two randomly picked sites. For copolymers with less blocky sequences the designed surfaces recognize the correct sequence well when the segment-surface interactions dominate over the intersegment interactions. For copolymers with more blocky sequences recognition is good when the segment-surface interactions are only slightly stronger than the intersegment interactions.

Adsorption↗
PubMed2004

Kinetics of fibril formation by polyalanine peptides.

Ordered beta-sheet complexes, termed amyloid fibrils, are the underlying structural components of the intra- and extracellular fibrillar protein deposits that are associated with a variety of human diseases, including Alzheimer's, Parkinson's, and the prion diseases. In this work, we investigated the kinetics of fibril formation using our newly developed off-lattice intermediate resolution model, PRIME. The model is simple enough to allow the treatment of large multichain systems while maintaining a fairly realistic description of protein dynamics without built-in bias toward any conformation when used in conjunction with constant temperature discontinuous molecular dynamics, a fast alternative to conventional molecular dynamics. Simulations were performed on systems containing 48-96 model Ac-KA14K-NH2 peptides. We found that fibril formation for polyalanines incorporate features that are characteristic of three models, the templated assembly, nucleated polymerization, and nucleated conformational conversion models, but that none of them gave a completely satisfactory description of the simulation kinetics. Fibril formation was nucleation-dependent, occurring after a lag time that decreased with increasing peptide concentration and increased with increasing temperature. Fibril formation appeared to be a conformational conversion process in which small amorphous aggregates --> beta-sheets --> ordered nucleus --> subsequent rapid growth of a small stable fibril or protofilament. Fibril growth in our simulations involved both beta-sheet elongation, in which the fibril grew by adding individual peptides to the end of each beta-sheet, and lateral addition, in which the fibril grew by adding already formed beta-sheets to its side. The initial rate of fibril formation increased with increasing concentration and decreased with increasing temperature.

Amyloid↗
PubMed2004

Molecular dynamics simulations of spontaneous fibril formation by random-coil peptides.

Assembly of normally soluble proteins into amyloid fibrils is a cause or associated symptom of numerous human disorders, including Alzheimer's and the prion diseases. We report molecular-level simulation of spontaneous fibril formation. Systems containing 12-96 model polyalanine peptides form fibrils at temperatures greater than a critical temperature that decreases with peptide concentration and exceeds the peptide's folding temperature, consistent with experimental findings. Formation of small amorphous aggregates precedes ordered nucleus formation and subsequent rapid fibril growth through addition of beta-sheets laterally and monomeric peptides at fibril ends. The fibril's structure is similar to that observed experimentally.

Amyloid↗
PubMed2004

Phase diagrams describing fibrillization by polyalanine peptides.

Amyloid fibrils are the structural components underlying the intra- and extracellular protein deposits that are associated with a variety of human diseases, including Alzheimer's, Parkinson's, and the prion diseases. In this work, we examine the thermodynamics of fibril formation using our newly-developed off-lattice intermediate-resolution protein model, PRIME. The model is simple enough to allow the treatment of large multichain systems while maintaining a fairly realistic description of protein dynamics when used in conjunction with constant-temperature discontinuous molecular dynamics, a fast alternative to conventional molecular dynamics. We conduct equilibrium simulations on systems containing 96 Ac-KA14K-NH2 peptides over a wide range of temperatures and peptide concentrations using the replica-exchange method. Based on measured values of the heat capacity, radius of gyration, and percentage of peptides that form the various structures, a phase diagram in the temperature-concentration plane is constructed delineating the regions where each structure is stable. There are four distinct single-phase regions: alpha-helices, fibrils, nonfibrillar beta-sheets, and random coils; and four two-phase regions: random coils/nonfibrillar beta-sheets, random coils/fibrils, fibrils/nonfibrillar beta-sheets, and alpha-helices/nonfibrillar beta-sheets. The alpha-helical region is at low temperature and low concentration. The nonfibrillar beta-sheet region is at intermediate temperatures and low concentrations and expands to higher temperatures as concentration is increased. The fibril region occurs at intermediate temperatures and intermediate concentrations and expands to lower as the peptide concentration is increased. The random-coil region is at high temperatures and all concentrations; this region shifts to higher temperatures as the concentration is increased.

Amyloid↗
PubMed2004

Phase behavior in model homopolymer/CO2 and surfactant/CO2 systems: discontinuous molecular dynamics simulations.

Discontinuous molecular dynamics simulations are performed on surfactant (HmTn)/solvent systems modeled as a mixture of single-sphere solvent molecules and freely jointed surfactant chains composed of m slightly solvent-philic head spheres (H) and n solvent-philic tail spheres (T), all of the same size. We use a square-well potential to account for the head-head, head-solvent, tail-tail, and tail-solvent interactions and a hard-sphere potential for the head-tail and solvent-solvent interactions. We first simulate homopolymer/supercritical CO2 (scCO2) systems to establish the appropriate interaction parameters for a surfactant/scCO2 system. Next, we simulate surfactant/scCO2 systems and explore the effect of the surfactant volume fraction, packing fraction, and temperature on the phase behavior. The transition from the two-phase region to the one-phase region is located by monitoring the contrast structure factor of the equilibrated surfactant/scCO2 system, and the micelle to unimer transition is located by monitoring the aggregate size distribution of the equilibrated surfactant/scCO2 system. We find a two-phase region, a micelle phase, and a unimer phase with increasing packing fraction at fixed temperature or with increasing temperature at fixed packing fraction. The phase diagram for the surfactant/scCO2 system in the surfactant volume fraction-packing fraction plane and the density dependence of the critical micelle concentration are in qualitative agreement with experimental observations. The phase behavior of a surfactant/scCO2 system can be directly related to the solubilities of the corresponding homopolymers that serve as the head and tail blocks for the surfactant. The influence of surfactant structure (head and tail lengths) on the phase transitions is explored.

Journal Article↗
PubMed2004

Lattice Monte Carlo simulations of phase separation and micellization in supercritical CO2/surfactant systems: effect of CO2 density.

Lattice Monte Carlo simulations are used to study the effect of nonionic surfactant concentration and CO2 density on the micellization and phase equilibria of supercritical CO2/surfactant systems. The interaction parameter for carbon dioxide is obtained by matching the critical temperature of the model fluid with the experimental critical temperature. Various properties such as the critical micelle concentration and the size, shape, and structure ofmicelles are calculated, and the phase diagram in the surfactant concentration-CO2 density space is constructed. On increasing the CO2 density, we find an increase in the critical micelle concentration and a decrease in the micellar size; this is consistent with existing experimental results. The variation of the micellar shape and structure with CO2 density shows that the micelles are spherical and that the extension of the micellar core increases with increasing micellar size, while the extension of the micellar corona increases with increasing CO2 density. The predicted phase diagram is in qualitative agreement with experimental phase diagrams for nonionic surfactants in carbon dioxide.

Journal Article↗
PubMed2004

Assembly and kinetic folding pathways of a tetrameric beta-sheet complex: molecular dynamics simulations on simplified off-lattice protein models.

We have performed discontinuous molecular dynamic simulations of the assembly and folding kinetics of a tetrameric beta-sheet complex that contains four identical four-stranded antiparallel beta-sheet peptides. The potential used in the simulation is a hybrid Go-type potential characterized by the bias gap parameter g, an artificial measure of a model protein's preference for its native state, and the intermolecular contact parameter eta, which measures the ratio of intermolecular to intramolecular native attractions. The formation of the beta-sheet complex and its equilibrium properties strongly depend on the size of the intermolecular contact parameter eta. The ordered beta-sheet complex in the folded state and nonaligned beta-sheets or tangled chains in the misfolded state are distinguished by measuring the squared radius of gyration Rg2 and the fraction of native contacts Q. The folding yield for the folded state is high at intermediate values of eta, but is low at both small and large values of eta. The folded state at small eta is liquid-like, but is solid-like at both intermediate and large eta. The misfolded state at small eta contains nonaligned beta-sheets and tangled chains with poor secondary structure at large eta. Various folding pathways via dimeric and trimeric intermediates are observed, depending on eta. Comparison with experimental results on protein aggregation indicates that intermediate eta values are most appropriate for modeling fibril formation and small eta values are most appropriate for modeling the formation of amorphous aggregates.

Binding Sites↗
PubMed2004

Thermodynamics and stability of a beta-sheet complex: molecular dynamics simulations on simplified off-lattice protein models.

We have performed discontinuous molecular dynamics simulations of the thermodynamics and stability of a tetrameric beta-sheet complex that contains four identical four-stranded antiparallel beta-sheet peptides. The potential used in the simulation is a hybrid Go-type potential characterized by the bias gap parameter g, an artificial measure of the preference of a model protein for its native state, and the intermolecular contact parameter eta, which measures the ratio of intermolecular to intramolecular native attractions. Despite the simplicity of the model, a complex set of thermodynamic transitions for the beta-sheet complex is revealed that shows there are three distinct oligomer (partially ordered, ordered, and highly ordered beta-sheet complex) states and four noninteracting monomers phases. The thermodynamic properties of the three oligomer states strongly depend on both the size of the intermolecular contact parameter eta and the temperature. The partially ordered beta-sheet complex is made up of four ordered globules and is observed at intermediate to large eta at high temperatures. The ordered beta-sheet complex contains four native beta-sheets and is located at small to intermediate eta at low temperatures in the phase diagram. The highly ordered beta-sheet complex has fully-stiff beta-sheet strands, the same as the global energy minimum structure, and is observed for all eta at low temperatures.

Computer Simulation↗
PubMed2004

Solvent effects on the conformational transition of a model polyalanine peptide.

We have investigated the folding of polyalanine by combining discontinuous molecular dynamics simulation with our newly developed off-lattice intermediate-resolution protein model. The thermodynamics of a system containing a single Ac-KA(14)K-NH(2) molecule has been explored by using the replica exchange simulation method to map out the conformational transitions as a function of temperature. We have also explored the influence of solvent type on the folding process by varying the relative strength of the side-chain's hydrophobic interactions and backbone hydrogen bonding interactions. The peptide in our simulations tends to mimic real polyalanine in that it can exist in three distinct structural states: alpha-helix, beta-structures (including beta-hairpin and beta-sheet-like structures), and random coil, depending upon the solvent conditions. At low values of the hydrophobic interaction strength between nonpolar side-chains, the polyalanine peptide undergoes a relatively sharp transition between an alpha-helical conformation at low temperatures and a random-coil conformation at high temperatures. As the hydrophobic interaction strength increases, this transition shifts to higher temperatures. Increasing the hydrophobic interaction strength even further induces a second transition to a beta-hairpin, resulting in an alpha-helical conformation at low temperatures, a beta-hairpin at intermediate temperatures, and a random coil at high temperatures. At very high values of the hydrophobic interaction strength, polyalanines become beta-hairpins and beta-sheet-like structures at low temperatures and random coils at high temperatures. This study of the folding of a single polyalanine-based peptide sets the stage for a study of polyalanine aggregation in a forthcoming paper.

Hydrogen Bonding↗