Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “dynamics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5Linked to original sources

[Comparison of contrast resolution between dynamic MRI and dynamic CT in liver scanning].

Dynamic studies employing contrast agents have been used in the detection of hepatocellular carcinoma (HCC). Some studies have indicated that the detectability of small nodules in liver is significantly better with MRI dynamic study than with CT dynamic study. However, the reasons for this are unclear. In an attempt to clarify this difference, we compared contrast resolution using two contrast phantoms made for MRI and CT. Both phantoms were constructed using the same ratios of dilution of the contrast materials Gd-DTPA and iodine. The signal intensity of different dilution ratios of contrast medium were measured in order to compare the signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and contrast. The results showed that the contrast of phantoms in MRI and CT dynamic studies was almost the same except for certain MRI scan techniques. In terms of SNR and CNR, MRI dynamic study was superior to CT dynamic study. In addition, T(1) relaxation of the contrast agent should be considered for suitable scanning, because the contrast of MRI was affected by scan parameters.

Contrast Media↗

Contrasting dynamics of Bartonella spp. in cyclic field vole populations: the impact of vector and host dynamics.

Many zoonotic disease agents are transmitted between hosts by arthropod vectors, including fleas, but few empirical studies of host-vector-microparasite dynamics have investigated the relative importance of hosts and vectors. This study investigates the dynamics of 4 closely related Bartonella species and their flea vectors in cyclic populations of field voles (Microtus agrestis) over 3 years. The probability of flea infestation was positively related to field vole density 12 months previously in autumn, but negatively related to more recent host densities, suggesting a dilution effect. The 4 Bartonella species exhibited contrasting dynamics. Only B. grahamii, showed a distinct seasonal pattern. Infection probability increased with field vole density for B. doshiae, B. taylorii and BGA (a previously unidentified species) and with density of coexisting wood mice for B. doshiae and B. grahamii. However, only the infection probability of BGA in spring was related to flea prevalence. B. doshiae and BGA were most common in older animals, but the other 2 were most common in non-reproductive hosts. Generally, host density rather than vector abundance appears most important for the dynamics of flea-transmitted Bartonella spp., possibly reflecting the importance of flea exchange between hosts. However, even closely related species showed quite different dynamics, emphasising that other factors such as population age structure can impact on zoonotic risk.

Animals↗

Exploring the dynamic information content of a protein NMR structure: comparison of a molecular dynamics simulation with the NMR and X-ray structures of Escherichia coli ribonuclease HI.

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7-ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high-resolution (1.48 A) crystal structure than with either the lower-resolution (2.05 A) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side-chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn-pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean-square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non-rotameric/syn-pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X-ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed.

Cluster Analysis↗

Comparison of the dynamics of bovine and human angiogenin: a molecular dynamics study.

Molecular dynamics simulations have been carried out for 1 ns on human and bovine angiogenin systems in an effort to compare and contrast their dynamics. An analysis of their dynamics is done by examining the rms deviations, following hydrogen-bonding interactions and looking at the role of water in and around the protein. The C-terminus of bovine angiogenin moves appreciably during dynamics suggesting a better structure for ligand binding. However, we do not find any evidence of a conformation where the glutamate residue that obstructs the active site takes on a different conformation. We observe a differential hydrogen-bonding pattern in the active site regions of bovine and human angiogenins, which could have a bearing on the different catalytic activities of the proteins. We also propose that the differential binding of the monoclonal antibody toward the two proteins might be due sequential and not conformational differences. Water molecules might play an important functional role in both proteins given their subtle functional differences. A simple computation on the molecular dynamics data has been carried out to identify locations in and around the protein that are invariably occupied by water. The locations of nearly half the waters we have identified from the simulation as being invariant in bovine angiogenin occupy similar locations in the bovine angiogenin crystal structure. The positions of the waters identified in human angiogenin differ considerably from that of bovine angiogenin.

Animals↗

Molecular dynamics simulations of papilloma virus E2 DNA sequences: dynamical models for oligonucleotide structures in solution.

The specificity of papilloma virus E2 protein-DNA binding depends critically upon the sequence of a region of the DNA not in direct contact with the protein, and represents one of the simplest known examples of indirect readout. A detailed characterization of this system in solution is important to the further investigation hypothesis of a structural code for DNA recognition by regulatory proteins. In the crystalline state, the E2 DNA oligonucleotide sequence, d(ACCGAATTCGGT), exhibits three different structural forms. We report herein studies of the structure of E2 DNA in solution based on a series of molecular dynamics (MD) simulations including counterions and water, utilizing both the canonical and various crystallographic structures as initial points of departure. All MDs converged on a single dynamical structure of d(ACCGAATTCGGT) in solution. The predicted structure is in close accord with two of the three crystal structures, and indicates that a significant kink in the double helix at the central ApT step in the other crystal molecule may be a packing effect. The dynamical fine structure was analyzed on the basis of helicoidal parameters. The calculated curvature in the sequence was found to originate primarily from YPR steps in the regions flanking the central AATT tract. In order to study the role of structural adaptation of the DNA in the binding process, a subsequent simulation on the 16-mer cognate sequence d(CAACCGAATTCGGTTG) was initiated from the crystallographic coordinates of the bound DNA in the crystal structure of the protein DNA complex. MD simulations starting with the protein-bound form relaxed rapidly back to the dynamical structure predicted from the previous simulations on the uncomplexed DNA. The MD results show that the bound form E2 DNA is a dynamically unstable structure in the absence of protein, and arises as a consequence of both structural changes intrinsic to the sequence and induced by the interaction with protein.

Animals↗

Dynamics of microtubules visualized by darkfield microscopy: treadmilling and dynamic instability.

Individual microtubules undergoing treadmilling in vitro were visualized by darkfield light microscopy, and the relationship between treadmilling and dynamic instability was studied as a function of microtubule-associated proteins (MAPs). In order to demonstrate treadmilling directly by real-time observation, we constructed three-block microtubules, the center-block of which was decorated with Tetrahymena dynein. The decorated block can easily be distinguished from undecorated blocks in the darkfield microscope because the decorated one appears much thicker. At steady-state conditions, the length of an undecorated block at one end increased and that at another end decreased, while the decorated center-block did not change in its length. The results from these direct observations show that calf brain 3X-microtubules exhibit a treadmilling flux of 0.9 micron/h. Using a similar microscopy technique, we previously demonstrated that phosphocellulose PC-microtubules existed in either the growing or the shortening phase and alternated quite frequently at steady-state conditions (dynamic instability). How does treadmilling relate to dynamic instability? An image recording of individual 3X-microtubules containing MAPs revealed that the microtubules undergo treadmilling and do not exhibit any dynamic instability. This evidence shows that MAPs suppress the dynamic instability of microtubules. That is, treadmilling can take place in the steady state only after microtubules have been stabilized by MAPs.

Animals↗

Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: calculation of structural and dynamic properties.

Long dynamics simulations were carried out on the B1 immunoglobulin-binding domain of streptococcal protein G (ProtG) and bovine pancreatic trypsin inhibitor (BPTI) using atomistic descriptions of the proteins and a continuum representation of solvent effects. To mimic frictional and random collision effects, Langevin dynamics (LD) were used. The main goal of the calculations was to explore the stability of tens-of-nanosecond trajectories as generated by this molecular mechanics approximation and to analyze in detail structural and dynamical properties. Conformational fluctuations, order parameters, cross correlation matrices, residue solvent accessibilities, pKa values of titratable groups, and hydrogen-bonding (HB) patterns were calculated from all of the trajectories and compared with available experimental data. The simulations comprised over 40 ns per trajectory for ProtG and over 30 ns per trajectory for BPTI. For comparison, explicit water molecular dynamics simulations (EW/MD) of 3 ns and 4 ns, respectively, were also carried out. Two continuum simulations were performed on each protein using the CHARMM program, one with the all-atom PAR22 representation of the protein force field (here referred to as PAR22/LD simulations) and the other with the modifications introduced by the recently developed CMAP potential (CMAP/LD simulations). The explicit solvent simulations were performed with PAR22 only. Solvent effects are described by a continuum model based on screened Coulomb potentials (SCP) reported earlier, i.e., the SCP-based implicit solvent model (SCP-ISM). For ProtG, both the PAR22/LD and the CMAP/LD 40-ns trajectories were stable, yielding C(alpha) root mean square deviations (RMSD) of about 1.0 and 0.8 A respectively along the entire simulation time, compared to 0.8 A for the EW/MD simulation. For BPTI, only the CMAP/LD trajectory was stable for the entire 30-ns simulation, with a C(alpha) RMSD of approximately 1.4 A, while the PAR22/LD trajectory became unstable early in the simulation, reaching a C(alpha) RMSD of about 2.7 A and remaining at this value until the end of the simulation; the C(alpha) RMSD of the EW/MD simulation was about 1.5 A. The source of the instabilities of the BPTI trajectories in the PAR22/LD simulations was explored by an analysis of the backbone torsion angles. To further validate the findings from this analysis of BPTI, a 35-ns SCP-ISM simulation of Ubiquitin (Ubq) was carried out. For this protein, the CMAP/LD simulation was stable for the entire simulation time (C(alpha) RMSD of approximately 1.0 A), while the PAR22/LD trajectory showed a trend similar to that in BPTI, reaching a C(alpha) RMSD of approximately 1.5 A at 7 ns. All the calculated properties were found to be in agreement with the corresponding experimental values, although local deviations were also observed. HB patterns were also well reproduced by all the continuum solvent simulations with the exception of solvent-exposed side chain-side chain (sc-sc) HB in ProtG, where several of the HB interactions observed in the crystal structure and in the EW/MD simulation were lost. The overall analysis reported in this work suggests that the combination of an atomistic representation of a protein with a CMAP/CHARMM force field and a continuum representation of solvent effects such as the SCP-ISM provides a good description of structural and dynamic properties obtained from long computer simulations. Although the SCP-ISM simulations (CMAP/LD) reported here were shown to be stable and the properties well reproduced, further refinement is needed to attain a level of accuracy suitable for more challenging biological applications, particularly the study of protein-protein interactions.

Alanine↗

Conformational and helicoidal analysis of the molecular dynamics of proteins: "curves," dials and windows for a 50 psec dynamic trajectory of BPTI.

A new procedure for the graphic analysis of molecular dynamics (MD) simulations on proteins is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the conformational and helicoidal parameters for each structure entering the analysis via the method "Curves," developed for proteins by Sklenar, Etchebest, and Lavery (Proteins: Structure, Function Genet. 6:46-60, 1989) followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels ("dials") for each torsional parameter and "windows" on the range values assumed by the linear and angular helicoidal parameters, and is present in a form isomorphous with the primary structure per se. The complete time evolution of dynamic structure can then be depicted in a set of four composite figures. Dynamic aspects of secondary and tertiary structure are also provided. The procedure is illustrated with an analysis of a 50 psec in vacuo simulation on the 58 residue protein, bovine pancreatic trypsin inhibitor (BPTI), in the vicinity of the local minimum on the energy surface corresponding to a high resolution crystal structure. The time evolution of 272 conformational and 788 helicoidal parameters for BPTI is analyzed. A number of interesting features can be discerned in the analysis, including the dynamic range of conformational and helicoidal motions, the dynamic extent of 2 degrees structure motifs, and the calculated fluctuations in the helix axis. This approach is expected to be useful for a critical analysis of the effects of various assumptions about force field parameters, truncation of potentials, solvation, and electrostatic effects, and can thus contribute to the development of more reliable simulation protocols for proteins. Extensions of the analysis to present differential changes in conformational and helicoidal parameters is expected to be valuable in MD studies of protein complexes with substrates, inhibitors, and effectors and in determining the nature of structural changes in protein-protein interactions.

Biological Evolution↗

Long timestep dynamics of peptides by the dynamics driver approach.

Previous experience with the Langevin/implicit-Euler scheme for dynamics ("LI") on model systems (butane, water) has shown that LI is numerically stable for timesteps in the 5-20 fs range but quenches high-frequency modes. To explore applications to polypeptides, we apply LI to model systems (several dipeptides, a tetrapeptide, and a 13-residue oligoalanine) and also develop a new dynamics driver approach ("DA"). The DA scheme, based on LI, addresses the important issue of proper sampling, which is unlikely to be solved by small-timestep integration methods or implicit methods with intrinsic damping at room temperature, such as LI. Equilibrium averages, time-dependent molecular properties, and sampling trends at room temperature are reported for both LI and DA dynamics simulations, which are then compared to those generated by a standard explicit discretization of the Langevin equation with a 1 fs timestep. We find that LI's quenching effects are severe on both the fast and slow (due to vibrational coupling) frequency modes of all-atom polypeptides and lead to more restricted dynamics at moderate timesteps (40 fs). The DA approach empirically counteracts these damping effects by adding random atomic perturbations to the coordinates at each step (before the minimization of a dynamics function). By restricting the energetic fluctuations and controlling the kinetic energy, we are able with a 60 fs timestep to generate continuous trajectories that sample more of the relevant conformational space and also reproduce reasonably Boltzmann statistics. Although the timescale for transition may be accelerated by the DA approach, the transitional information obtained for the alanine dipeptide and the tetrapeptide is consistent with that obtained by several other theoretical approaches that focus specifically on the determination of pathways. While the trajectory for oligoalanine by the explicit scheme over the nanosecond timeframe remains in the vicinity of the full alpha R-helix starting structure, and a high-temperature (600 degrees K) MD trajectory departs slowly from the alpha helical structure, the DA-generated trajectory for the same CPU time exhibits unfolding and refolding and reveals a range of conformations with an intermediate helix content. Significantly, this range of states is more consistent with spectroscopic experiments on small peptides, as well as the cooperative two-state model for helix-coil transition. The good, near-Boltzmann statistics reported for the smaller systems above, in combination with the interesting oligoalanine results, suggest that DA is a promising tool for efficiently exploring conformational spaces of biomolecules and exploring folding/unfolding processes of polypeptides.

Alanine↗

Structure and dynamics of ferrocytochrome c553 from Desulfovibrio vulgaris studied by NMR spectroscopy and restrained molecular dynamics.

The solution structure of Desulfovibrio vulgaris Hildenborough (DvH) ferrocytochrome c553 has been determined by nuclear magnetic resonance spectroscopy and combined simulated annealing/high temperature restrained molecular dynamics calculations. This three-stage protocol consists of an initial determination of overall fold from randomised co-ordinates, followed by a 20 picosecond exploratory stage, during which the non-bonded terms are simplified to facilitate as broad a sampling of conformational space as possible, and a 26 picosecond refinement stage, using the full AMBER force field. This latter stage systematically improved the energetic and convergence characteristics of the ensemble, while still satisfying the experimental restraints. Forty structures have been obtained from a total of 875 distance constraints for this protein of 79 amino acid residues. The root-mean-square deviation over all residues with respect to the mean is 0.70(+/- 0.12)A for the backbone (N, C alpha and C') atoms. Two conformations of the turn motif at the solvent/heme cleft interface have been identified, both fulfilling the experimental data and having equally viable energetic characteristics. The stability of the ensemble and the dynamic characteristics have been further investigated by subjecting ten of the structures to constraint-free molecular dynamics calculations (130 picoseconds) in vacuo. The structures were found to be stable to within 1.5 A of the initial backbone conformation. Comparison with the dynamic behaviour of the restrained molecular dynamics calculations has been used to identify regions of inherent flexibility in the molecule.

Cytochrome c Group↗

Essential dynamics from NMR clusters: dynamic properties of the Myb DNA-binding domain and a hinge-bending enhancing variant.

Application of the "essential dynamics" method to the NMR cluster of structures for the R2R3 DNA-binding domain of the mouse c-Myb transcriptional activator is described. Using this method, large concerted fluctuations of atoms are extracted showing a hinge-bending motion between the two (R2 and R3) Myb repeats on the basis of NMR data alone. Molecular dynamics simulation of the same protein allowed quantitative comparison of the large concerted motions calculated from experimental and theoretical data, showing a significant degree of similarity. Detailed inspection of the motions reveals a conserved proline that plays a key role in determining hinge flexibility. The proline-to-alanine mutation at this position, which has previously been characterized biochemically, was subjected to molecular dynamics and subsequent essential dynamics analysis. The hinge-bending motion between the two repeats was found to be enhanced for the mutant. The approach described should have general applications, predicting the effect of mutations on protein dynamic properties of other proteins.

Binding Sites↗

Safety and efficacy of dynamic graciloplasty for fecal incontinence: report of a prospective, multicenter trial. Dynamic Graciloplasty Therapy Study Group.

PURPOSE: Dynamic graciloplasty has been used for intractable fecal incontinence, and good results have been reported. The aim of this study was to assess prospectively the safety and efficacy of dynamic graciloplasty for intractable fecal incontinence in a prospective, multicenter trial. METHODS: A total of 123 adults were treated with dynamic graciloplasty at 20 institutions. Continence was assessed preoperatively and postoperatively by use of 14-day diaries. RESULTS: There was one treatment-related death. One hundred eighty-nine adverse events occurred in 91 patients (74 percent). Forty-nine patients (40 percent) required one or more operations to treat complications. One hundred seventy (90 percent) events were resolved. Sixty-three percent of patients without pre-existing stomas recorded a 50 percent or greater decrease in incontinent events 12 months after dynamic graciloplasty, and an additional 11 percent experienced lesser degrees of improvement. Twenty-six percent were not improved, worsened, or exited. In patients with pre-existing stomas, 33 percent achieved successful outcomes at 12 months. This number increased to 60 percent at 18 months. Seventy-eight percent of patients had increased enema retention time, and mean anal canal pressures improved significantly at 12 months. Significant changes in quality of life were also observed. CONCLUSIONS: Objective improvement can be demonstrated in the majority of patients with end-stage fecal incontinence treated with dynamic graciloplasty. Reduction in incontinence episodes can be correlated with improved quality of life. Adverse events are frequently encountered, but most resolve with treatment.

Adolescent↗

Dynamization timing and its effect on bone healing when using the Orthofix Dynamic Axial Fixator.

We report the effect of dynamization timing on healing rates when treating tibial fractures with the Dynamic Axial Factor (Orthofix, registered trademark). Sixty-nine patients with 71 fresh tibial fractures were grouped retrospectively for age, Gustilo and Anderson wound grading, site of tibial fracture and presence of an intact fibula or stable fibular fracture. Pairs were then generated chronologically to compare healing rates in those patients dynamized before 4 weeks with those dynamized after this time. This resulted in 22 pairs of patients; the other fractures could not be paired. Their healing times were compared using a paired t-test and were found to be significantly quicker in the group dynamized within the first 4 weeks (P < 0.05) and we therefore recommend that this should be routine practice.

Adolescent↗

Molecular dynamics simulations of the holo and apo forms of retinol binding protein. Structural and dynamical changes induced by retinol removal.

The effects of removing retinol from the X-ray structure of holo-retinol binding protein are studied using the molecular dynamics technique. Structural and dynamical properties emerging from an 80 ps simulation of the apo form, for which no crystallographic structure is available, are compared with the results of a 70 ps trajectory of the holo-protein. Dynamical stationarity is attained after roughly 30 ps, and the resulting average structure is proposed as a reasonable model of the apo-protein. Conformational changes are observed for the loops at the beta-barrel entrance during the non-equilibrium part of the apo-trajectory. Tryptophan labelling experiments and retinoid reconstitution experiments point towards this part of the molecule as being involved in prealbumin binding. Structural changes in this region may therefore explain the differences in prealbumin affinity between the apo and holo forms. Furthermore, a change in the position of the alpha-helix, corresponding to a pivot around its C terminus, is observed for the apo-protein. The resulting conformation of the alpha-helix is found to be similar to that in apo-beta-lactoglobulin, which also can bind retinol and for which a crystal structure exists. The results from the holo simulation are compared to the crystallographic data and show good agreement. The dynamics of the secondary and tertiary structural elements are analysed and compared for the two forms. The beta-barrel is found to be extremely cooperative in its atomic motions in both simulations, and the top and bottom beta-sheets perform collective fluctuations with respect to each other in the low-frequency limit of the simulations. The dynamics of the alpha-helical region presents clear differences between the two forms; while the holo-protein has a well-defined spectrum for the longitudinal stretching mode, the apo form displays a fairly large bending of the alpha-helix at several points of the trajectory.

Computer Simulation↗

Bimolecular reaction simulation using Weighted Ensemble Brownian dynamics and the University of Houston Brownian Dynamics program.

We discuss here the implementation of the Weighted Ensemble Brownian (WEB) dynamics algorithm of Huber and Kim in the University of Houston Brownian Dynamics (UHBD) suite of programs and its application to bimolecular association problems. WEB dynamics is a biased Brownian dynamics (BD) algorithm that is more efficient than the standard Northrup-Allison-McCammon (NAM) method in cases where reaction events are infrequent because of intervening free energy barriers. Test cases reported here include the Smoluchowski rate for association of spheres, the association of the enzyme copper-zinc superoxide dismutase with superoxide anion, and the binding of the superpotent sweetener N-(p-cyanophenyl)-N'-(diphenylmethyl)-guanidinium acetic acid to a monoclonal antibody fragment, NC6.8. Our results show that the WEB dynamics algorithm is a superior simulation method for enzyme-substrate reaction encounters with large free energy barriers.

Acetates↗

Dynamics-function correlation in Cu, Zn superoxide dismutase: a spectroscopic and molecular dynamics simulation study.

A single mutation (Val29-->Gly) at the subunit interface of a Cu, Zn superoxide dismutase dimer leads to a twofold increase in the second order catalytic rate, when compared to the native enzyme, without causing any modification of the structure or the electric field distribution. To check the role of dynamic processes in this catalytic enhancement, the flexibility of the dimeric protein at the subunit interface region has been probed by the phosphorescence and fluorescence properties of the unique tryptophan residue. Multiple spectroscopic data indicate that Trp83 experiences a very similar, and relatively hydrophobic, environment in both wild-type and mutant protein, whereas its mobility is distinctly more restrained in the latter. Molecular dynamics simulation confirms this result, and provides, at the molecular level, details of the dynamic change felt by tryptophan. Moreover, the simulation shows that the loops surrounding the active site are more flexible in the mutant than in the native enzyme, making the copper more accessible to the incoming substrate, and being thus responsible for the catalytic rate enhancement. Evidence for increased, dynamic copper accessibility also comes from faster copper removal in the mutant by a metal chelator. These results indicate that differences in dynamic, rather than structural, features of the two enzymes are responsible for the observed functional change.

Amino Acid Sequence↗

Molecular dynamics simulation of dimeric and monomeric forms of human prion protein: insight into dynamics and properties.

A central theme in prion protein research is the detection of the process that underlies the conformational transition from the normal cellular prion form (PrP(C)) to its pathogenic isoform (PrP(Sc)). Although the three-dimensional structures of monomeric and dimeric human prion protein (HuPrP) have been revealed by NMR spectroscopy and x-ray crystallography, the process underlying the conformational change from PrP(C) to PrP(Sc) and the dynamics and functions of PrP(C) remain unknown. The dimeric form is thought to play an important role in the conformational transition. In this study, we performed molecular dynamics (MD) simulations on monomeric and dimeric HuPrP at 300 K and 500 K for 10 ns to investigate the differences in the properties of the monomer and the dimer from the perspective of dynamic and structural behaviors. Simulations were also undertaken with Asp178Asn and acidic pH, which is known as a disease-associated factor. Our results indicate that the dynamics of the dimer and monomer were similar (e.g., denaturation of helices and elongation of the beta-sheet). However, additional secondary structure elements formed in the dimer might result in showing the differences in dynamics and properties between the monomer and dimer (e.g., the greater retention of dimeric than monomeric tertiary structure).

Binding Sites↗

Structural and dynamic properties of the homodimeric hemoglobin from Scapharca inaequivalvis Thr-72-->Ile mutant: molecular dynamics simulation, low temperature visible absorption spectroscopy, and resonance Raman spectroscopy studies.

Molecular dynamics simulations, low temperature visible absorption spectroscopy, and resonance Raman spectroscopy have been performed on a mutant of the Scapharca inaequivalvis homodimeric hemoglobin, where residue threonine 72, at the subunit interface, has been substituted by isoleucine. Molecular dynamics simulation indicates that in the Thr-72-->Ile mutant several residues that have been shown to play a role in ligand binding fluctuate around orientations and distances similar to those observed in the x-ray structure of the CO derivative of the native hemoglobin, although the overall structure remains in the T state. Visible absorption spectroscopy data indicate that in the deoxy form the Soret band is less asymmetric in the mutant than in the native protein, suggesting a more planar heme structure; moreover, these data suggest a similar heme-solvent interaction in both the liganded and unliganded states of the mutant protein, at variance with that observed in the native protein. The "conformation sensitive" band III of the deoxy mutant protein is shifted to lower energy by >100 cm-1 with respect to the native one, about one-half of that observed in the low temperature photoproducts of both proteins, indicating a less polar or more hydrophobic heme environment. Resonance Raman spectroscopy data show a slight shift of the iron-proximal histidine stretching mode of the deoxy mutant toward lower frequency with respect to the native protein, which can be interpreted in terms of either a change in packing of the phenyl ring of Phe-97, as also observed from the simulation, or a loss of water in the heme pocket. In line with this latter interpretation, the number of water molecules that dynamically enters the intersubunit interface, as calculated by the molecular dynamics simulation, is lower in the mutant than in the native protein. The 10-ns photoproduct for the carbonmonoxy mutant derivative has a higher iron-proximal histidine stretching frequency than does the native protein. This suggests a subnanosecond relaxation that is slowed in the mutant, consistent with a stabilization of the R structure. Taken together, the molecular dynamics and the spectroscopic data indicate that the higher oxygen affinity displayed by the Thr-72-->Ile mutant is mainly due to a local perturbation in the dimer interface that propagates to the heme region, perturbing the polarity of the heme environment and propionate interactions. These changes are consistent with a destabilization of the T state and a stabilization of the R state in the mutant relative to the native protein.

Animals↗