SEARCH · PubMed Health

Results for “trajectory analysis”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 559 records · Page 31Linked to original sources

PubMed2002

Homology modeling and molecular dynamics study of NAD-dependent glycerol-3-phosphate dehydrogenase from Trypanosoma brucei rhodesiense, a potential target enzyme for anti-sleeping sickness drug development.

Sleeping sickness and Chagas disease are among the most severe diseases in Africa as well as Latin America. These two diseases are caused by Trypanosoma spp. Recently, an enzyme of a glycolytic pathway, NAD-dependent glycerol-3-phosphate dehydrogenase, of Leishmania mexicana was crystallized and its structure determined by x-ray crystallography. This structure has offered an excellent template for modeling of the homologous enzymes from another Trypanosoma species. Here, a homology model of the T. brucei enzyme based on the x-ray structure of LmGPDH has been generated. This model was used as the starting point for molecular dynamics simulation in a water box. The analysis of the molecular dynamics trajectory indicates that the functionally important motifs have both a very stable secondary structure and tertiary arrangement.

Amino Acid Sequence↗

PubMed1996

Dynamics of proteins in different solvent systems: analysis of essential motion in lipases.

We have investigated the effect of different solvents on the dynamics of Rhizomucor miehei lipase. Molecular dynamics simulations were performed in water, methyl hexanoate, and cyclohexane. Analysis of the 400-ps trajectories showed that the solvent has a pronounced effect on the geometrical properties of the protein. The radius of gyration and total accessibility surface decrease in organic solvents, whereas the number of hydrogen bonds increases. The essential motions of the protein in different solvents can be described in a low-dimensional "essential subspace," and the dynamic behavior in this subspace correlates with the polarity of the solvent. Methyl hexanoate, which is a substrate for R. miehei lipase, significantly increases the fluctuations in the active-site loop. During the simulation, a methyl hexanoate entered the active-site groove. This observation provides insight into the possible docking mechanism of the substrate.

Fungal Proteins↗

PubMed2003

Conformations of terminal sialyloligosaccharide fragments--a molecular dynamics study.

Molecular dynamics simulations have been performed to understand the conformational features of the terminal sialyloligosaccharide fragments NeuNAc alpha(2-3)Gal, NeuNAc alpha(2-6)Gal, NeuNAc alpha(2-8)NeuNAc and NeuNAc alpha(2-9)NeuNAc. The conformational regions A(i), B(i) and C(i) were identified in the Ramachandran plot. Analysis of the 1000 ps trajectories collected through simulation (2000 ps in the case of NeuNAc alpha (2-9)NeuNAc) revealed that these molecules have conformational propensity in region B(i). The occurrence of these molecules in the common conformational space leads to a structural similarity between them. This structural similarity may be an essential requirement for the neuraminidase activity towards sialyloligosaccharides. The local change in the conformation of the active site residues of neuraminidases may contribute for the specificity differences between different linkages of sialyloligosaccharides. A highly conserved water-mediated hydrogen bond observed in these structures between the sugar residues, acts as an additional stabilizing force.

Carbohydrate Conformation↗

PubMed2002

Deuterium spin probes of side-chain dynamics in proteins. 2. Spectral density mapping and identification of nanosecond time-scale side-chain motions.

In the previous paper in this issue we have demonstrated that it is possible to measure the five different relaxation rates of a deuteron in (13)CH(2)D methyl groups of (13)C-labeled, fractionally deuterated proteins. The extensive set of data acquired in these experiments provides an opportunity to investigate side-chain dynamics in proteins at a level of detail that heretofore was not possible. The data, acquired on the B1 domain of peptostreptococcal protein L, include 16 (9) relaxation measurements at 4 (2) different magnetic field strengths, 25 degrees C (5 degrees C). These data are shown to be self-consistent and are analyzed using a spectral density mapping procedure which allows extraction of values of the spectral density function at a number of frequencies with no assumptions about the underlying dynamics. Dynamics data from 31 of 35 methyls in the protein for which data could be obtained were well-fitted using the two-parameter Lipari-Szabo model (Lipari, G.; Szabo, A. J. Am. Chem. Soc. 1982, 104, 4546). The data from the remaining 4 methyls can be fitted using a three-parameter version of the Lipari-Szabo model that takes into account, in a simple manner, additional nanosecond time-scale local dynamics. This interpretation is supported by analysis of a molecular dynamics trajectory where spectral density profiles calculated for side-chain methyl sites reflect the influence of slower (nanosecond) time-scale motions involving jumps between rotameric wells. A discussion of the minimum number of relaxation measurements that are necessary to extract the full complement of dynamics information is presented along with an interpretation of the extracted dynamics parameters.

Bacterial Proteins↗

PubMed2005

Theoretical study of the truncated hemoglobin HbN: exploring the molecular basis of the NO detoxification mechanism.

Mycobacterium tuberculosis is the causative agent of human tuberculosis. The nitric oxide reaction with oxy-truncated hemoglobin N (trHbN) has been proposed to be responsible for the resistance mechanism by which this microorganism can evade the toxic effects of NO. In this work, we explore the molecular basis of the NO detoxification mechanism using a combination of classical and hybrid quantum-classical (QM-MM) simulation techniques. We have investigated the structural flexibility of the protein, the ligand affinity properties, and the nitric oxide reaction with coordinated O2. The analysis of the classical MD trajectory allowed us to identify Phe62 as the gate of the main channel for ligand diffusion to the active site. Moreover, the opening of the channel stems from the interplay between collective backbone motions and local rearrangements in the side chains of the residues that form the bottleneck of the tunnel. Even though the protein environment is not found to make a significant contribution to the heme moiety catalyzed reaction, the binding site influences the physiological function of the enzyme at three different levels. First, by isolating the intermediates formed in the reaction, it prevents nondesired reactions from proceeding. Second, it modulates the ligand (O2, NO) affinity of the protein, which can be ascribed to both distal and proximal effects. Finally, the stabilization of the Tyr33-Gln58 pair upon O2 binding might alter the essential dynamics of the protein, leading in turn to a mechanism for ligand-induced regulation.

Heme↗

PubMed2005

Single molecule spectroscopy as a probe for dye-polymer interactions.

Experimental (Single Molecule Spectroscopy) and theoretical (quantum-chemical calculations and Monte Carlo and molecular dynamics simulations) techniques are combined to investigate the behavior and dynamics of a polymer-dye molecule system. It is shown that the dye molecule of interest (1,1'-dioctadecyl-3,3,3',3'-tetramethylindo-dicarbocyanine) adopts two classes of conformations, namely planar and nonplanar ones, when embedded in a poly(styrene) matrix. From an in-depth analysis of the fluorescence lifetime trajectories, the planar conformers can be further classified according to the way their alkyl side chains interact with the surrounding poly(styrene) chains.

Journal Article↗

PubMed1996

Investigation of the proton-assisted pathway to formation of the catalytically active, ferryl species of P450s by molecular dynamics studies of P450eryF.

The recently determined crystal structure of cytochrome P450eryF (6-deoxyerythronolide B hydroxylase; CYP107A1) in its ferric heme substrate-bound form has been used to address one of the most fundamental unresolved aspects of the mechanism of oxidation common to this ubiquitous family of metabolizing heme proteins, the pathway from the twice reduced dioxygen species to the putative catalytically active ferryl oxygen species. Both of these species are too transient to have been characterized experimentally, and the transformation from one to the other has been only partially characterized. The observed requirement of two protons and the formation of water in this transformation suggests a proton-assisted dioxygen bond cleavage as a plausible pathway. However, this pathway is difficult to establish by experiment alone, and the source of the protons in the largely hydrophobic binding pocket of the P450s remains unclear. In this work we have performed molecular dynamics simulations of the twice reduced dioxygen substrate-bound form of this isozyme in order to (i) determine the plausibility of the proposed pathway to compound I formation, a proton-assisted cleavage of the dioxygen bond, and (ii) investigate the possible source of these protons. The analysis of the molecular dynamics trajectories of this species does indeed provide further evidence for this pathway and points to a source of protons. Specifically, two dynamically stable hydrogen bonds to the distal oxygen atom of the dioxygen ligand, one by the substrate and the other by a bound water, are found, consistent with the proposed proton-assisted cleavage of the bond and formation of water. In addition, an extensive dynamically stable hydrogen bond network is formed that connects the distal oxygen to Glu 360, a well-conserved residue in a channel accessible to solvent that could be the ultimate source of protons. The simulations were done for both a protonated and unprotonated Glu and led to a proposed mechanism of proton transfer by it to the distal oxygen atom. In order to validate the procedures used for the simulation of this transient twice-reduced species, we have used these same procedures to perform molecular dynamics simulations of two other forms of P450eryF, the ferric and ferryl substrate-bound species, and compared the results with experiment. The results for the ferric substrate-bound species were assessed by comparisons to the experimentally determined X-ray structure and fluctuations, and good agreement was found. The simulations performed for the ferryl substrate-bound species led to the correct prediction of the observed regio- and stereospecific hydroxylation of its natural substrate, 6-deoxyerythronolide B (6-DEB) at the 6S position. The results of these two additional studies lend credibility to the important mechanistic inferences from the simulations of the transient twice reduced dioxygen species: further evidence for a proton-assisted pathway from it to the catalytically active ferryl species and a possible source of the protons.

Bacterial Proteins↗

PubMed2005

n-Heptane under pressure: structure and dynamics from molecular simulations.

Atomistic molecular dynamics simulations have been performed in the isothermal-isobaric ensemble to explore the phase behavior of n-heptane. Motivated by recent high-pressure spectroscopic experiments on n-heptane, the present work aims at understanding the liquid-solid and the alluded to solid-solid transitions upon increasing pressure. Starting from the stabilized solid phase at 300 K and 10 kbar, we have investigated the range of these two transitions by a gradual decrease and increase of pressure, respectively. Although the solid-liquid transition has clear signatures such as the formation of gauche defects along the molecular backbone, the present model does not show any sign of a first-order solid-solid transition at high pressures. However, interesting changes in the environment around methyl groups and in their dynamics are observed. These have been substantiated by calculations of the vibrational density of states obtained from a normal-mode analysis and from the simulation trajectory.

Journal Article↗

PubMed2005

Subsecond luminescence intensity fluctuations of single CdSe quantum dots.

Photoluminescence (PL) intermittency characteristics are examined for single quantum dots (QDs) in a CdSe QD sample synthesized at a slow rate at 75 degrees C. Although the PL quantum efficiency was relatively low ( approximately 0.25), we noticed that the PL intensity of single CdSe QDs fluctuated on a subsecond time scale with short-lived "on" and "off" states. The subsecond PL intensity fluctuations of CdSe QDs are different from "on" and "off" PL blinking generally observed for QDs fluctuating on a millisecond to minute time scale. We characterized single QDs by identifying polarized excitations, topographic imaging using atomic force microscopy (AFM), and transmission electron microscopy (TEM). From analysis of the PL intensity trajectories from >100 single CdSe QDs, the average intermittency time was 213 ms. From the PL quantum efficiency, slow growth of QDs, intensity trajectory analyses, and previous reports relating surface trap states and PL properties of QDs, we attribute the subsecond PL intensity fluctuations of single CdSe QDs and short-lived "on" and "off" states to a high-density distribution of homogeneous surface trap states.

Cadmium Compounds↗

PubMed2004

Kinetic control of dimer structure formation in amyloid fibrillogenesis.

Amyloid fibril formation involves nonfibrillar oligomeric intermediates, which are important as possible cytotoxic species in neurodegenerative diseases. However, their transient nature and polydispersity have made it difficult to identify their formation mechanism or structure. We have investigated the dimerization process, the first step in aggregate formation, by multiple molecular dynamics simulations of five beta-sheet-forming peptides. Contrary to the regular beta-sheet structure of the amyloid fibril, the dimers exhibit all possible combinations of beta-sheets, with an overall preference for antiparallel arrangements. Through statistical analysis of 1,000 dimerization trajectories, each 1 ns in length, we have demonstrated that the observed distribution of dimer configurations is kinetically determined; hydrophobic interactions orient the peptides so as to minimize the solvent accessible surface area, and the dimer structures become trapped in energetically unfavorable conformations. Once the hydrophobic contacts are present, the backbone hydrogen bonds form rapidly by a zipper-like mechanism. The initial nonequilibrium structures formed are stable during the 1-ns simulation time for all five peptides at room temperature. In contrast, at higher temperatures, where rapid equilibration among different configurations occurs, the distribution follows the global energies. The relaxation time of dimers at room temperature was estimated to be longer than the time for diffusional encounters with other oligomers at typical concentrations. These results suggest that kinetic trapping could play a role in the structural evolution of early aggregates in amyloid fibrillogenesis.

Amyloid beta-Peptides↗

PubMed2002

Analysis of a library of macaque nuclear mitochondrial sequences confirms macaque origin of divergent sequences from old oral polio vaccine samples.

Nuclear mtDNA sequences (numts) are a widespread family of paralogs evolving as pseudogenes in chromosomal DNA [Zhang, D. E. & Hewitt, G. M. (1996) TREE 11, 247-251 and Bensasson, D., Zhang, D., Hartl, D. L. & Hewitt, G. M. (2001) TREE 16, 314-321]. When trying to identify the species origin of an unknown DNA sample by way of an mtDNA locus, PCR may amplify both mtDNA and numts. Indeed, occasionally numts dominate confounding attempts at species identification [Bensasson, D., Zhang, D. X. & Hewitt, G. M. (2000) Mol. Biol. Evol. 17, 406-415; Wallace, D. C., et al. (1997) Proc. Natl. Acad. Sci. USA 94, 14900-14905]. Rhesus and cynomolgus macaque mtDNA haplotypes were identified in a study of oral polio vaccine samples dating from the late 1950s [Blancou, P., et al. (2001) Nature (London) 410, 1045-1046]. They were accompanied by a number of putative numts. To confirm that these putative numts were of macaque origin, a library of numts corresponding to a small segment of 12S rDNA locus has been made by using DNA from a Chinese rhesus macaque. A broad distribution was found with up to 30% sequence variation. Phylogenetic analysis showed that the evolutionary trajectories of numts and bona fide mtDNA haplotypes do not overlap with the signal exception of the host species; mtDNA fragments are continually crossing over into the germ line. In the case of divergent mtDNA sequences from old oral polio vaccine samples [Blancou, P., et al. (2001) Nature (London) 410, 1045-1046], all were closely related to numts in the Chinese macaque library.

Animals↗

PubMed2003

Revealing competitive Forster-type resonance energy-transfer pathways in single bichromophoric molecules.

We demonstrate measurements of the efficiency of competing Förster-type energy-transfer pathways in single bichromophoric systems by monitoring simultaneously the fluorescence intensity, fluorescence lifetime, and the number of independent emitters with time. Peryleneimide end-capped fluorene trimers, hexamers, and polymers with interchromophore distances of 3.4, 5.9, and on average 42 nm, respectively, served as bichromophoric systems. Because of different energy-transfer efficiencies, variations in the interchromophore distance enable the switching between homo-energy transfer (energy hopping), singlet-singlet annihilation, and singlet-triplet annihilation. The data suggest that similar energy-transfer pathways have to be considered in the analysis of single-molecule trajectories of donor/acceptor pairs as well as in natural and synthetic multichromophoric systems such as light-harvesting antennas, oligomeric fluorescent proteins, and dendrimers. Here we report selectively visualization of different energy-transfer pathways taking place between identical fluorophores in individual bichromophoric molecules.

Fluorescence Resonance Energy Transfer↗

PubMed2005

Stability and change of neuroticism in aging.

Data from the Longitudinal Aging Study Amsterdam were used to study the relationship between neuroticism and aging. At baseline, cross-sectional analyses of data from 2,117 respondents (aged 55-85 years, M = 70) showed no significant age differences. The magnitude of the 3- and 6-year stability coefficients was high, and 12% of the elderly participants showed a clinically relevant mean level change. Longitudinal multilevel analyses showed a small but statistical significant change with aging, but the mean change was not considered clinically relevant. A U-formed course was found, showing a slight decrease until respondents reached the age of 70. Adjusting the model for physical health-related variables slightly increased the stability. An additional interaction analysis showed that the individual trajectory of neuroticism was not affected by the physical health status. In conclusion, neuroticism remains rather stable in middle and older adulthood, with some apparent increase in late life.

Aged↗

PubMed2002

Closed loops: persistence of the protein chain returns.

It has recently been discovered that globular proteins are universally built from standard loop-n-lock units of about 30 amino acid residues. The hypothesis has been put forward on the loop stage in the protein evolution when the units were autonomous. Later they joined together making longer chains. One would expect that the early individual loop-n-lock elements might still be detected in modern protein sequences as remnants of the hypothetical 30-residue sequence prototypes. Among several strong sequence motifs, extracted from protein sequences of 23 complete bacterial proteomes, one 32-residue prototype was studied here in detail. Numerous sequence segments related to the prototype are identified in the crystal structures of proteins of a PDB_SELECT database. Analysis of the respective chain trajectories for the cases with different degrees of sequence conservation confirms that the majority of the segments correspond to the closed loops. In the evolutionary diversification of the prototypes the secondary structure yields first, while the sequence is still moderately conserved. The last feature to go is the chain return property. Apparently, the opening of the loops would severely destabilize the protein fold, which explains their conservation.

Amino Acid Motifs↗

PubMed1996

Molecular modelling of Staphylococcal delta-toxin ion channels by restrained molecular dynamics.

Delta-Toxin is a 26-residue channel-forming peptide from Staphylococcus aureus which forms an amphipathic alpha-helix in a membrane environment. Channel formation in planar bilayers suggests that an average of six delta-toxin helices self-assemble to form transbilayer pores. Molecular models for channels formed by delta-toxin and by a synthetic analogue have been generated using a simulated annealing protocol applied via restrained molecular dynamics. These models are analysed in terms of the predicted geometric and energetic properties of the transbilayer pores. Pore radius calculations of the models demonstrate that rings of channel-lining residues contribute a series of constrictions along the pore. Electrostatic properties of the pores are determined both by pore-lining charged side chains and by the aligned helix dipoles of the parallel helix bundle. Molecular dynamics simulations (100 ps) of delta-toxin models containing intra-pore water were performed. Analysis of the resultant dynamics trajectories further supports the proposal that alternative conformations of pore-constricting side chains may be responsible for the observed conductance heterogeneity of delta-toxin ion channels.

Amino Acid Sequence↗

PubMed2002

Using Lyapunov exponents to predict the onset of chaos in nonlinear oscillators.

An analytic technique for predicting the emergence of chaotic instability in nonlinear nonautonomous dissipative oscillators is proposed. The method is based on the Lyapunov-type stability analysis of an arbitrary phase trajectory and the standard procedure of calculating the Lyapunov characteristic exponents. The concept of temporally local Lyapunov exponents is then utilized for specifying the area in the phase space where any trajectory is asymptotically stable, and, therefore, the existence of chaotic attractors is impossible. The procedure of linear coordinate transform optimizing the linear part of the vector field is developed for the purpose of maximizing the stability area in the vicinity of a stable fixed point. By considering the inverse conditions of asymptotic stability, this approach allows formulating a necessary condition for chaotic motion in a broad class of nonlinear oscillatory systems, including many cases of practical interest. The examples of externally excited one- and two-well Duffing oscillators and a planar pendulum demonstrate efficiency of the proposed method, as well as a good agreement of the theoretical predictions with the results of numerical experiments. The comparison of the proposed method with Melnikov's criterion shows a potential advantage of using the former one at high values of dissipation parameter and/or multifrequency type of excitation in dynamical systems.

Journal Article↗

PubMed1991

Conformation and inhibitory properties of peptides based on the tissue kallikrein-aprotinin complex.

A series of peptides encompassing the primary binding segment (residues 12-19) of aprotinin has been synthesized and tested for their ability to inhibit porcine pancreatic kallikrein. A minimum sequence of five amino acids spanning residues 12-16 of aprotinin is necessary for inhibition of porcine pancreatic kallikrein. An octapeptide homologous with the binding segment of aprotinin has a Ki-value of 1.2 x 10(-4) M. The solution structure of the octapeptide was studied by one- and two-dimensional NMR methods for comparison with the known structure of the segment of aprotinin that contacts tissue kallikrein. NMR experiments suggest that the peptide is either a random coil or that it samples several conformations on the NMR time scale. Analysis of the molecular dynamics trajectory of the octapeptide also suggests that the peptide is highly flexible. Thus, inhibition by the octapeptide occurs because of its homology with residues 12-19 of aprotinin. Moreover, the absence of a stable solution conformation similar to that of the binding segment of aprotinin is consistent with the 150,000-fold increase in Ki of the octapeptide compared to intact aprotinin.

Amino Acid Sequence↗

PubMed1997

Investigating locomotion of dairy cows by use of high speed cinematography.

The longterm influence of management systems on the locomotion of 17 dairy cows was investigated by high speed cinematography (100 frames/s) and kinematic analysis. Angular patterns and hoof trajectories of the left fore- and hindlimbs are presented and statistics made of occurring minimum and maximum angles. At the recording, 3 cows had been kept in tie-stalls (TI) and 6 cows in cubicles (CI) for a consecutive time of about 2.5 years while 8 cows had been kept on grass for about 3 months. Four of the grazing cows had earlier been kept in cubicles (CG) and 4 in tie-stalls (TG) during earlier off grazing seasons together with TI and CI cows. The CI cows had a smaller maximum angle of the elbow joint compared to TI, TG and CG cows. The hock joint angle of the CI cows was less flexed during the stance phase than in TI and CG cows while the minimum angle during the swing phase was greater in the TI and CI cows compared to TG and CG cows. Pastured cows (TG and CG) had a less pronounced flexion of the fetlock joint angle during the stance compared to cows kept indoors (TI and CI). The results suggest that slatted floor and lack of exercise during summer grazing may affect locomotion. This is indicated by restrictions in the movements of the elbow and hock joints and in less fetlock joint flexion at full support.

Animals↗