Biomedical subjects

Donald Hamelberg

Publications and source records attributed to Donald Hamelberg.

17 recordsLinked to original sources

PubMed2006

On the application of accelerated molecular dynamics to liquid water simulations.

Our group recently proposed a robust bias potential function that can be used in an efficient all-atom accelerated molecular dynamics (MD) approach to simulate the transition of high energy barriers without any advance knowledge of the potential-energy landscape. The main idea is to modify the potential-energy surface by adding a bias, or boost, potential in regions close to the local minima, such that all transitions rates are increased. By applying the accelerated MD simulation method to liquid water, we observed that this new simulation technique accelerates the molecular motion without losing its microscopic structure and equilibrium properties. Our results showed that the application of a small boost energy on the potential-energy surface significantly reduces the statistical inefficiency of the simulation while keeping all the other calculated properties unchanged. On the other hand, although aggressive acceleration of the dynamics simulation increases the self-diffusion coefficient of water molecules greatly and dramatically reduces the correlation time of the simulation, configurations representative of the true structure of liquid water are poorly sampled. Our results also showed the strength and robustness of this simulation technique, which confirm this approach as a very useful and promising tool to extend the time scale of the all-atom simulations of biological system with explicit solvent models. However, we should keep in mind that there is a compromise between the strength of the boost applied in the simulation and the reproduction of the ensemble average properties.

Algorithms↗

PubMed2006

Insight into the role of hydration on protein dynamics.

The potential energy surface of a protein is rough. This intrinsic energetic roughness affects diffusion, and hence the kinetics. The dynamics of a system undergoing Brownian motion on this surface in an implicit continuum solvent simulation can be tuned via the frictional drag or collision frequency to be comparable to that of experiments or explicit solvent simulations. We show that the kinetic rate constant for a local rotational isomerization in stochastic simulations with continuum solvent and a collision frequency of 2 ps(-1) is about 10(4) times faster than that in explicit water and experiments. A further increase in the collision frequency to 60 ps(-1) slows down the dynamics, but does not fully compensate for the lack of explicit water. We also show that the addition of explicit water does not only slow down the dynamics by increasing the frictional drag, but also increases the local energetic roughness of the energy landscape by as much as 1.0 kcal/mol.

Computer Simulation↗

PubMed2006

Elasticity of peptide omega bonds.

We calculated the changes of the free energy profile of the peptidyl-prolyl torsional angle of the dipeptide valine-proline under pulling forces by simulations. Using a dynamic model built on the equilibrium properties of this system and previously studied dynamic properties of cis-trans isomerization of other dipeptides, we calculated the dynamic viscoelasticity of this degree of freedom. The results show significant differences between how thermal and mechanical forces alter the equilibrium and the dynamics of the isomerization transition. The former does not change the barrier heights but changes the prefactor of the kinetics owing to temperature effects, while the latter changes minima and thus the population. The force that is required to "excite" this degree of freedom is small. Compared to other systems, we found that this degree of freedom becomes already quite rigid at several hertz, which is a much lower value due to the high barrier of the cis-trans isomerization. We also found that the tensile elastic modulus of densely packed omega bonds is at the order of GPa, which is comparable to that of polymer materials. These results give mechanical properties of polyproline elasticity of a local nature and provide guidance for future experimental designs.

Binding Sites↗

PubMed2005

Mapping the tropomyosin isoform 5 binding site on human erythrocyte tropomodulin: further insights into E-Tmod/TM5 interaction.

Actin protofilaments in the erythrocyte membrane skeleton are uniformly approximately 37nm. This length may be in part attributed to a "molecular ruler" made of erythrocyte tropomodulin (E-Tmod) and tropomyosin (TM) isoforms 5 or 5b. We previously mapped the E-Tmod binding site to TM5 N-terminal heptad repeat residues "a" (I(7), I(14)), "d" (V(10)) and "f" (R(12)). We now map the TM5 binding site to E-Tmod residues at L(116), E(117) and/or E(118) by identifying among 35 deletion clones and a series of point mutations that no longer bind to human TM5 and rat TM5b. Upstream residues 71-104 contain an actin binding site. The N-terminal "KRK ring" may participate in balancing electrostatic force with hydrophobic interaction in dimerization of TM and its binding to E-Tmod.

Amino Acid Sequence↗

PubMed2005

Fast peptidyl cis-trans isomerization within the flexible Gly-rich flaps of HIV-1 protease.

The catalytic aspartyl protease of the HIV-1 virus is a homodimer with two flaps that control access to the active site and are known to be flexible. However, knowledge of the atomistic mechanism of the flexibility is lacking. We show that the Gly-Gly omega-bond in the glycine-rich flap tips undergoes fast cis-trans isomerization on the microsecond to millisecond time scale rather than in the usual seconds. Further study reveals that the unexpectedly fast isomerization is a direct consequence of the beta-hairpin loop structure of the flap tips, which appears to be counterintuitive. After loop formation of a linear peptide containing the Gly-Gly motif, the rate of isomerization is shown to increase by many orders of magnitude.

Glycine↗

PubMed2005

Relating kinetic rates and local energetic roughness by accelerated molecular-dynamics simulations.

We show that our accelerated molecular-dynamics (MD) approach can extend the time scale in all-atom MD simulations of biopolymers. We also show that this technique allows for the kinetic rate information to be recaptured. In deducing the kinetic rates, the relationship between the local energetic roughness of the potential-energy landscape and the effective diffusion coefficient is established. These are demonstrated on a very slow but important biomolecular process: the dynamics of cis-trans-isomerization of Ser-Pro motifs. We do not only recapture the slow kinetic rates, which is difficult in traditional MD, but also obtain the underlying roughness of the energy landscape of proteins at atomistic resolution.

Algorithms↗

PubMed2005

Phosphorylation effects on cis/trans isomerization and the backbone conformation of serine-proline motifs: accelerated molecular dynamics analysis.

The presence of serine/threonine-proline motifs in proteins provides a conformational switching mechanism of the backbone through the cis/trans isomerization of the peptidyl-prolyl (omega) bond. The reversible phosphorylation of the serine/threonine modulates this switching in regulatory proteins to alter signaling and transcription. However, the mechanism is not well understood. This is partly because cis/trans isomerization is a very slow process and, hence, difficult to study. We have used our accelerated molecular dynamics method to study the cis/trans proline isomerization, preferred backbone conformation of a serine-proline motif, and the effects of phosphorylation of the serine residue. We demonstrate that, unlike normal molecular dynamics, the accelerated molecular dynamics allows for the system to escape very easily from the trans isomer to cis isomer, and vice versa. Moreover, for both the unphosphorylated and phosphorylated peptides, the statistical thermodynamic properties are recaptured, and the results are consistent with experimental values. Isomerization of the proline omega bond is shown to be asymmetric and strongly dependent on the psi backbone angle before and after phosphorylation. The rates of escape decrease after phosphorylation. Also, the alpha-helical backbone conformation is more favored after phosphorylation. This accelerated molecular dynamics approach provides a general approach for enhancing the conformational transitions of molecular systems without having prior knowledge of the location of the minima and barriers on the potential-energy landscape.

Amino Acid Motifs↗

PubMed2005

The folding energy landscape and phosphorylation: modeling the conformational switch of the NFAT regulatory domain.

An energy landscape approach predicts the conformational changes of the configurations of the regulatory domain of the protein nuclear factor of activated T cells (NFAT) caused by phosphorylation of specific multiple sites. Structurally local effects and secondary structural changes are modeled using all-atom Brownian dynamics to investigate the changes of the backbone torsional distributions upon phosphorylation. For tertiary and global changes, we employ a coarse-grained model to sample ensembles of conformations both with and without phosphorylation. At the secondary structure level, phosphorylation moderately increases the helical propensity and gives a more rigid local backbone conformation. The tertiary effects of phosphorylation caused by the extensive charge modification are more pronounced and collectively change the conformation of the regulatory domain of NFAT from a flexible globular ensemble to a rather rigid helical bundle, blocking access to the nuclear localization sequence. These studies give computational support to one scenario conjectured from experiments.

Amino Acid Sequence↗

PubMed2004

Structural selectivity of aromatic diamidines.

Competition dialysis was used to study the interactions of 13 substituted aromatic diamidine compounds with 13 nucleic acid structures and sequences. The results show a striking selectivity of these compounds for the triplex structure poly dA:(poly dT)2, a novel aspect of their interaction with nucleic acids not previously described. The triplex selectivity of selected compounds was confirmed by thermal denaturation studies. Triplex selectivity was found to be modulated by the location of amidine substituents on the core phenyl-furan-phenyl ring scaffold. Molecular models were constructed to rationalize the triplex selectivity of DB359, the most selective compound in the series. Its triplex selectivity was found to arise from optimal ring stacking on base triplets, along with proper positioning of its amidine substituents to occupy the minor and the major-minor grooves of the triplex. New insights into the molecular recognition of nucleic acid structures emerged from these studies, adding to the list of available design principles for selectively targeting DNA and RNA.

Benzamidines↗

PubMed2004

Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method.

Localized water molecules in the binding pockets of proteins play an important role in noncovalent association of proteins and small drug compounds. At times, the dominant contribution to the binding free energy comes from the release of localized water molecules in the binding pockets of biomolecules. Therefore, to quantify the energetic importance of these water molecules for drug design purposes, we have used the double-decoupling approach to calculate the standard free energy of tying up a water molecule in the binding pockets of two protein complexes. The double-decoupling approach is based on the underlying principle of statistical thermodynamics. We have calculated the standard free energies of tying up the water molecule in the binding pockets of these complexes to be favorable. These water molecules stabilize the protein-drug complexes by interacting with the ligands and binding pockets. Our results offer ideas that could be used in optimizing protein-drug interactions, by designing ligands that are capable of targeting localized water molecules in protein binding sites. The resulting free energy of ligand binding could benefit from the potential free energy gain accompanying the release of these water molecules. Furthermore, we have examined the theoretical background of the double-decoupling method and its connection to the molecular dynamics thermodynamic integration techniques.

Binding Sites↗

PubMed2004

Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules.

Many interesting dynamic properties of biological molecules cannot be simulated directly using molecular dynamics because of nanosecond time scale limitations. These systems are trapped in potential energy minima with high free energy barriers for large numbers of computational steps. The dynamic evolution of many molecular systems occurs through a series of rare events as the system moves from one potential energy basin to another. Therefore, we have proposed a robust bias potential function that can be used in an efficient accelerated molecular dynamics approach to simulate the transition of high energy barriers without any advance knowledge of the location of either the potential energy wells or saddle points. In this method, the potential energy landscape is altered by adding a bias potential to the true potential such that the escape rates from potential wells are enhanced, which accelerates and extends the time scale in molecular dynamics simulations. Our definition of the bias potential echoes the underlying shape of the potential energy landscape on the modified surface, thus allowing for the potential energy minima to be well defined, and hence properly sampled during the simulation. We have shown that our approach, which can be extended to biomolecules, samples the conformational space more efficiently than normal molecular dynamics simulations, and converges to the correct canonical distribution.

Binding Sites↗

PubMed2004

DNA sequence dependent monomer-dimer binding modulation of asymmetric benzimidazole derivatives.

A number of studies indicate that DNA sequences such as AATT and TTAA have significantly different physical and interaction properties. To probe these interaction differences in detail and determine the influence of charge, we have synthesized three bisbenzimidazole derivatives, a diamidine, DB185, and monoamidines, DB183 and DB210, that are related to the well-known minor groove agent, Hoechst 33258. Footprinting studies with several natural and designed DNA fragments indicate that the synthetic compounds bind at AT sequences in the minor groove and interact more weakly at sites with TpA steps relative to sites without such steps. Circular dichroism spectroscopy also indicates that the compounds bind in the DNA minor groove. Surprisingly, Tm studies as a function of ratio indicate that the monoamidines bind to TTAA sequences as dimers, whereas the diamidine binds as a monomer. Biosensor-surface plasmon resonance (SPR) studies allowed us to quantitate the interaction differences in more detail. SPR results clearly show that the monoamidine compounds bind to the TTAA sequence in a cooperative 2:1 complex but bind as monomers to AATT. The dication binds to both sequences in monomer complexes but the binding to AATT is significantly stronger than binding to TTAA. Molecular dynamics simulations indicate that the AATT sequence has a narrow time-average minor groove width that is a very good receptor site for the bisbenzimidazole compounds. The groove in TTAA sequences is wider and the width must be reduced to form a favorable monomer complex. The monocations thus form cooperative dimers that stack in an antiparallel orientation and closely fit the structure of the TTAA minor groove. The amidine groups in the dimer are oriented in the 5' direction of the strand to which they are closest. Charge repulsion in the dication apparently keeps it from forming the dimer. It instead reduces the TTAA groove width, in an induced fit process, sufficiently to form a minor groove complex. The dimer-binding mode of DB183 and DB210 is a new DNA recognition motif and offers novel design concepts for selective targeting of DNA sequences with a wider minor groove, including those with TpA steps.

Benzimidazoles↗

PubMed2004

Characterization of a novel DNA minor-groove complex.

Many dicationic amidine compounds bind in the DNA minor groove and have excellent biological activity against a range of infectious diseases. Para-substituted aromatic diamidines such as furamidine, which is currently being tested against trypanosomiasis in humans, and berenil, which is used in animals, are typical examples of this class. Recently, a meta-substituted diamidine, CGP 40215A, has been found to have excellent antitrypanosomal activity. The compound has a linear, conjugated linking group that can be protonated under physiological conditions when the compound interacts with DNA. Structural and molecular dynamics analysis of the DNA complex indicated an unusual AT-specific complex that involved water-mediated H-bonds between one amidine of the compound and DNA bases at the floor of the minor groove. To investigate this unique system in more detail DNase I footprinting, surface plasmon resonance biosensor techniques, linear dichroism, circular dichroism, ultraviolet-visible spectroscopy, and additional molecular dynamics simulations have been conducted. Spectrophotometric titrations of CGP 40215A binding to poly(dAT)(2) have characteristics of DNA-binding-induced spectral changes as well as effects due to binding-induced protonation of the compound linker. Both footprinting and surface plasmon resonance results show that this compound has a high affinity for AT-rich sequences of DNA but very weak binding to GC sequences. The dissociation kinetics of the CGP 40215A-DNA complex are much slower than with similar diamidines such as berenil. The linear dichroism results support a minor-groove complex for the compound in AT DNA sequences. Molecular dynamics studies complement the structural analysis and provide a clear picture of the importance of water in mediating the dynamic interactions between the ligand and the DNA bases in the minor groove.

Base Sequence↗

PubMed2004

An empirical study of the universal chemical key algorithm for assigning unique keys to chemical compounds.

In this paper, we introduce an algorithm that assigns an essentially unique key called the Universal Chemical Key (UCK) to molecular structures. The molecular structures are represented as labeled graphs whose nodes abstract atoms and whose edges abstract bonds. The algorithm was tested on 236,917 compounds obtained from the National Cancer Institute (NCI) database of chemical compounds. On this database, the UCK algorithm assigned unique keys for chemicals with distinct molecular structures.

Algorithms↗

PubMed2003

Experimental studies of the Universal Chemical Key (UCK) algorithm on the NCI database of chemical compounds.

We have developed an algorithm called the Universal Chemical Key (UCK) algorithm that constructs a unique key for a molecular structure. The molecular structures are represented as undirected labeled graphs with the atoms representing the vertices of the graph and the bonds representing the edges. The algorithm was tested on 236,917 compounds obtained from the National Cancer Institute (NCI) database of chemical compounds. In this paper we present the algorithm,some examples and the experimental results on the NCI database. On the NCI database, the UCK algorithm provided distinct unique keys for chemicals with different molecular structures.

Algorithms↗

PubMed2002

Strong binding in the DNA minor groove by an aromatic diamidine with a shape that does not match the curvature of the groove.

A combination of biophysical techniques has been used to characterize the interaction of an antitrypanosomal agent, CGP 40215A, with DNA. The results from a broad array of methods (DNase I footprinting, surface plasmon resonance, X-ray crystallography, and molecular dynamics) indicate that this compound binds to the minor groove of AT DNA sequences. Despite its unusual linear shape that is not complementary to that of the DNA groove, a high binding affinity was observed in comparison with other similar but more curved diamidine compounds. The amidine groups at both ends of the ligand and the -NH groups on the linker are involved in extensive and dynamic H-bonds to the DNA bases. Complementary and consistent results were obtained from both the X-ray and molecular dynamics studies; both of these methods reveal direct and water-mediated H-bonds between the ligand and the DNA.

DNA↗

PubMed2002

Effect of a neutralized phosphate backbone on the minor groove of B-DNA: molecular dynamics simulation studies.

Alternative models have been presented to provide explanations for the sequence-dependent variation of the DNA minor groove width. In a structural model groove narrowing in A-tracts results from direct, short-range interactions among DNA bases. In an electrostatic model, the narrow minor groove of A-tracts is proposed to respond to sequence-dependent localization of water and cations. Molecular dynamics simulations on partially methylphosphonate substituted helical chains of d(TATAGGCCTATA) and d(CGCGAATTCGCG) duplexes have been carried out to help evaluate the effects of neutralizing DNA phosphate groups on the minor groove width. The results show that the time-average minor groove width of the GGCC duplex becomes significantly more narrow on neutralizing the phosphate backbone with methylphosphonates. The minor groove of the AATT sequence is normally narrow and the methylphosphonate substitutions have a smaller but measurable affect on this sequence. These results and models provide a system that can be tested by experiment and they support the hypothesis that the electrostatic environment around the minor groove affects the groove width in a sequence-dependent dynamic and time-average manner.

Base Sequence↗