ProDy: protein dynamics inferred from theory and experiments

SUMMARY: We developed a Python package, ProDy, for structure-based analysis of protein dynamics. ProDy allows for quantitative characterization of structural variations in heterogeneous datasets of structures experimentally resolved for a given biomolecular system, and for comparison of these variations with the theoretically predicted equilibrium dynamics. Datasets include structural ensembles for a given family or subfamily of proteins, their mutants and sequence homologues, in the presence/absence of their substrates, ligands or inhibitors. Numerous helper functions enable comparative analysis of experimental and theoretical data, and visualization of the principal changes in conformations that are accessible in different functional states. ProDy application programming interface (API) has been designed so that users can easily extend the software and implement new methods. AVAILABILITY: ProDy is open source and freely available under GNU General Public License from http://www.csb.pitt.edu/ProDy/.     

DNA basepair step deformability inferred from molecular dynamics simulations

The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.     

Osteoarthritic gait dynamics from force-plate measurements

A method is presented for calculating the displacement of the centre of body mass in lateral, horizontal and vertical directions during pathological gait; evaluations are made from measurements obtained by a single pair of small force plates. Results are presented for pre- and postoperative gait for 42 patients who, initially suffering from osteoarthritis of the hip, have had total hip replacements. The method is shown to be more convenient, more direct and less expensive than techniques in present use.     

Ribosomal dynamics explored by cryo-electron microscopy

Cryo-electron microscopy has furnished direct evidence for conformational changes of the ribosome as it proceeds, in a cyclic manner, through different functional states. Strategies to explore the ribosome dynamics include trapping of particular functionally meaningful states by chemical, genetic, or physical means. The new atomic information obtained by X-ray crystallography should make it possible to track conformational changes observed by cryo-electron microscopy to changes of the underlying structural framework of the ribosome.     

Nitrogen dynamics in subtropical fringe and basin mangrove forests inferred from stable isotopes

Mangroves exhibit low species richness compared to other tropical forests, but great structural and functional diversity. Aiming to contribute to a better understanding of the functioning of mangrove forests, we investigated nitrogen (N) dynamics in two physiographic types of mangroves (fringe and basin forests) in southeastern Brazil. Because fringe forests are under great influence of tidal flushing we hypothesized that these forests would exhibit higher N cycling rates in sediment and higher N losses to the atmosphere compared to basin forests. We quantified net N mineralization and nitrification rates in sediment and natural abundance of N stable isotopes (δ¹⁵N) in the sediment–plant–litter system. The fringe forest exhibited higher net N mineralization rates and δ¹⁵N in the sediment–plant–litter system, but net nitrification rates were similar to those of the basin forest. The results of the present study suggest that fringe forests exhibit higher N availability and N cycling in sediment compared to basin forests.     

Cell cycle dynamics inferred from the static properties of cells in balanced growth

The duration of a morphological phase of the cell cycle is reflected in the steady state distribution of the sizes of cells in that phase. Relationships presented here provide a method for estimating the timing and variability of any cell cycle phase. It is shown that the mean size of cells initiating and finishing any phase can be estimated from (1) the frequency of cells exhibiting the distinguishing morphological or autoradiographic features of the phase; (2) the mean size of cells in the phase; and (3) their coefficient of variation. The calculations are based on a submodel of the Koch-Schaechter Growth Controlled Model which assumes that (i) the distribution of division sizes is Gaussian; (ii) there is no correlation in division sizes between successive generations; and (iii) every cell division gives rise to two daughter cells of equal size. The calculations should be useful for a wider range of models, however, because the extrapolation factors are not sensitive to the chosen model. Criteria are proposed to allow the user to check the method's applicability for any experimental case. The method also provides a more efficient test of the dependence of growth on cell size than does the Collins-Richmond method. This is because the method uses the mean and coefficient of variation of the size of the total population, in conjunction with those of the cells in a final phase of the cell cycle, to test potential growth laws. For Escherichia coli populations studied by electron microscopy, an exponential growth model provided much better agreement than did a linear growth model. The computer simulations were used to generate rules for three types of cell phases: those that end at cell division, those that start at cell division, and those totally contained within a single cell cycle. For the last type, additional criteria are proposed to establish if the phase is well enough contained for the formulae and graphs to be used. The most useful rule emerging from these computer studies is that the fraction of the cell cycle time occupied by a phase is the product of the frequency of the phase and the ratio of the mean size of cells in that phase to the mean size of all cells in the population. A further advantage of the techniques presented here is that they use the 'extant' distributions that were actually measured, and not hypothesized distributions nor the special distributions needed for Collins-Richmond method that can only be calculated from the observed distributions of dividing or newborn cells on the basis of an assumed growth law.     

Boolean dynamics of genetic regulatory networks inferred from microarray time series data

MOTIVATION: Methods available for the inference of genetic regulatory networks strive to produce a single network, usually by optimizing some quantity to fit the experimental observations. In this article we investigate the possibility that multiple networks can be inferred, all resulting in similar dynamics. This idea is motivated by theoretical work which suggests that biological networks are robust and adaptable to change, and that the overall behavior of a genetic regulatory network might be captured in terms of dynamical basins of attraction. RESULTS: We have developed and implemented a method for inferring genetic regulatory networks for time series microarray data. Our method first clusters and discretizes the gene expression data using k-means and support vector regression. We then enumerate Boolean activation-inhibition networks to match the discretized data. Finally, the dynamics of the Boolean networks are examined. We have tested our method on two immunology microarray datasets: an IL-2-stimulated T cell response dataset and a LPS-stimulated macrophage response dataset. In both cases, we discovered that many networks matched the data, and that most of these networks had similar dynamics. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Structural features of the binding site of cholera toxin inferred from fluorescence measurements

The dependence on pH of the fluorescence of cholera toxin and its A and B subunits has been studied at 25 degrees C. The fluorescence intensity of cholera toxin is highly pH-dependent. In the pH range 7-9.5 it reaches a maximum corresponding to a quantum yield of 0.076. In the pH range 4-7 a strong increase in fluorescence intensity is observed (delta Q/Qmax = 0.64). Evaluation of the pH sensitivity of the fluorescence intensity of the A and B subunits reveals that the B subunit is mainly responsible for the observed pH effect (delta Q/Qmax for B subunit = 0.64). The intensity changes are paralleled by similar although less pronounced changes in the average fluorescence excited state life-time tau (delta tau/tau max = 0.33 for cholera toxin). Fluorimetric titration of the B subunit, which is related to the indole fluorescence of the lone Trp-88, reveals that the fluorescence intensity changes in the pH range 4-7 are due to reaction of two types of ionizable quencher displaying apparent pKa values of 4.4 and 6.2, respectively. It is suggested that the increase in fluorescence intensity with a midpoint at pH 6.2 is the result of deionization of the imidazolium side-chain of one or two out of the four histidine residues present in each beta-polypeptide chain, whereas a deionized carboxyl group is responsible for the quenching with midpoint at pH 4.4. Complex formation of cholera toxin or B subunit with the monosialoganglioside GM1 or the oligosaccharide moiety of GM1 (oligo-GM1) completely prevents the quenching by both quenchers. Addition of 6 M urea also eliminates the pH effect. The quenching is not the result of the dissociation of the B subunit into its constituent monomers. Upon fluorimetric titration of cholera toxin or B subunit above pH 9, a progressive drop in both fluorescence intensity and tau occurs. This decrease could be due to energy transfer from the indole moiety of Trp-88 to ionized tyrosines or by quenching through an unprotonated epsilon-amino group of lysine. Fluorimetric titration of the A subunit indicates that the tryptophan fluorescence is only moderately altered by ionizable groups displaying a pKa in the range 4 to 9. Activation of A subunit does not affect this lack of pH sensitivity. Above pH 9, however, a much more significant drop in the fluorescence intensity of activated A subunit occurs. The structural implications of the results are discussed.     

An allosteric mechanism inferred from molecular dynamics simulations on phospholamban pentamer in lipid membranes

Phospholamban functions as a regulator of Ca(2+) concentration of cardiac muscle cells by triggering the bioactivity of sarcoplasmic reticulum Ca(2+)-ATPase. In order to understand its dynamic mechanism in the environment of bilayer surroundings, we performed long time-scale molecular dynamic simulations based on the high-resolution NMR structure of phospholamban pentamer. It was observed from the molecular dynamics trajectory analyses that the conformational transitions between the "bellflower" and "pinwheel" modes were detected for phospholamban. Particularly, the two modes became quite similar to each other after phospholamban was phosphorylated at Ser16. Based on these findings, an allosteric mechanism was proposed to elucidate the dynamic process of phospholamban interacting with Ca(2+)-ATPase.     

Insights into lid movements of Burkholderia cepacia lipase inferred from molecular dynamics simulations

The interfacial activation of many lipases at water/lipid interface is mediated by large conformational changes of a so‐called lid subdomain that covers up the enzyme active site. Here we investigated using molecular dynamic simulations in different explicit solvent environments (water, octane and water/octane interface) the molecular mechanism by which the lid motion of Burkholderia cepacia lipase might operate. Although B. cepacia lipase has so far only been crystallized in open conformation, this study reveals for the first time the major conformational rearrangements that the enzyme undergoes under the influence of the solvent, which either exposes or shields the active site from the substrate. In aqueous media, the lid switches from an open to a closed conformation while the reverse motion occurs in organic environment. In particular, the role of a subdomain facing the lid on B. cepacia lipase conformational rearrangements was investigated using position‐restrained MD simulations. Our conclusions indicate that the sole mobility of α9 helix side‐chains of B. cepacia lipase is required for the full completion of the lid conformational change which is essentially driven by α5 helix movement. The role of selected α5 hydrophobic residues on the lid movement was further examined. In silico mutations of two residues, V138 and F142, were shown to drastically modify the conformational behavior of B. cepacia lipase. Overall, our results provide valuable insight into the role played by the surrounding environment on the lid conformational rearrangement and the activation of B. cepacia lipase. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural features of luliberin (luteinising hormone-releasing factor) inferred from fluorescence measurements.

The fluorescence and excitation spectra of luliberin (luteinizing hormone-releasing factor) in 0.005 M aqueous ammonium acetate are identical in shape to those of N-acetyltryptophan amide and are related to the indole side chain of Trp3. The change of fluoresecence intensity of luliberin with pH was measured in the range of pH 4-11. The increase of pH from 4 to 7.5 is followed by about 50% increase in fluorescence intensity due to deprotonation of the imidazolium side chain of His2. The fluorimetric titration curve in this pH region reveals a pK value for His2 of 5.95. Increasing of pH from 8 to 11 results in about 40% quenching of the fluorescence due to electronic energy transfer from the excited indole of Trp3 to the phenolate side chain of Tyr5. The pK value of Tyr5, obtained independently from the fluorimetric and photometric titrations indicate that at pH 7-8 luliberin contains only one charged residue, Arg8, which is in close vicinity to both His2 and Tyr5. The side chains of His2, Tyr5 and Arg8 presumably form a combined unit which may play an active role in the hormone action. Trp3 is at a maximal distance from this unit and may thus act as an independent active unit.     

Estimating net primary productivity from grassland biomass dynamics measurements

To address the need for a high quality data set based upon field observations suitable for parameterization, calibration, and validation of terrestrial biosphere models, we have developed a comprehensive global database on net primary productivity (NPP). We have compiled field measurements of biomass and associated environmental data for multiple study sites in major grassland types worldwide. Where sufficient data were available, we compared aboveground and total NPP estimated by six computational methods (algorithms) for 31 grassland sites. As has been found previously, NPP estimates were 2–5 times higher using methods which accounted for the dynamics of dead matter, compared with what is still the most commonly applied estimate of NPP (maximum peak live biomass). It is suggested that assumptions such as the use of peak biomass as an indicator of NPP in grasslands may apply only within certain subbiomes, e.g. temperate steppe grasslands. Additional data on belowground dynamics, or other reliable estimates of belowground productivity, are required if grasslands are to be fully appreciated for their role in the global carbon cycle.     

Thermodynamic interpretation of protein dynamics from NMR relaxation measurements

Protein dynamics and thermodynamics can be characterized through measurements of relaxation rates of side chain (2)H and (13)C, and backbone (15)N nuclei using NMR spectroscopy. The rates reflect protein motions on timescales from picoseconds to milliseconds. Backbone and methyl side chain NMR relaxation measurements for several proteins are beginning to reveal the role of protein dynamics in protein stability and ligand binding.     

Evolutionary dynamics of the human endogenous retrovirus family HERV-K inferred from full-length proviral genomes

Several distinct families of endogenous retroviruses exist in the genomes of primates. Most of them are remnants of ancient germ-line infections. The human endogenous retrovirus family HERV-K represents the unique known case of endogenous retrovirus that amplified in the human genome after the divergence of human and chimpanzee lineages. There are two types of HERV-K proviral genomes differing by the presence or absence of 292 bp in the pol-env boundary. Human-specific insertions exist for both types. The analyses shown in the present work reveal that several lineages of type 1 and type 2 HERV-K proviruses remained transpositionally active after the human/chimpanzee split. The data also reflect the important role of mosaic evolution (either by recombination or gene conversion) during the evolutionary history of HERV-K. Received: 5 February 2001 / Accepted: 22 March 2001     

Ribosomal DNA variations in Erianthus, a wild sugarcane relative (Andropogoneae-Saccharinae)

Variation at the 18S+26S and 5S ribosomal DNA loci was assessed on 62 Erianthus Michx. clones, representing 11 species, and 15 clones from two Saccharum L. species used as a reference. Genus-specific markers for Erianthus Michx. sect. Ripidium Henrard (Old World species) were identified. Ribosomal DNA units in Erianthus sect. Ripidium exhibited an additional BamHI site compared to Saccharum, and 5S units showed length and restriction-site differences between Erianthus and Saccharum. These markers will be useful to follow introgression in Saccharum x Erianthus hybrids. Six ribosomal units (for 18+26S genes) were revealed in Erianthus sect. Ripidium, differing by restriction-site positions and/or length. These results provided new information on species relationships and evolution within the genus Erianthus. The Indonesian and Indian forms of E. arundinaceus (Retz.) Jeswiet gave different restriction patterns, which were similar to those of E. bengalense (Retz.) R. C. Bharadwaja and E. procerus (Roxb.) Raizade, respectively. The two 2n=20 species, E. ele-phantinus Hook.f. and E. ravennae (L.) P. Beauv., could also be differentiated at this locus. Two of the New World Erianthus species studied, E. rufipilus (Steud.) Griseb. and E. longisetosus Andersson, appeared more like Erianthus sect. Ripidium, whereas E. trinii Hack, and E. brevibardis Michx. showed patterns consistent with Miscanthus sinensis Andersson and S. spontaneum L., respectively. Finally, the comparison of rDNA restriction maps among Erianthus sect. Ripidium, Saccharum, sorghum and maize, led to unexpected conclusions concerning the relationships between the different genera and the position of Erianthus in the Saccharum complex.     

Ribosomal DNA variations in wild and cultivated rubber tree (Hevea brasiliensis).

Ribosomal DNA variations were surveyed to assess the genetic variability among Hevea brasiliensis genetic resources. One hundred and sixty-eight individuals, including 73 cultivated Wickham clones and 95 wild clones from a prospection, were analyzed. Restriction mapping of rDNA units showed that RFLP variations are the result of both length and site (EcoRI) variations in the intergenic spacer (IGS). These variations can be revealed between as well as within individuals. A total of 12 spacer length variants is scored in the whole population, as well as two different ribosomal units (refered as type I and type II), defined by the presence or absence of an EcoRI site. Particular associations between spacer length variants and unit types can be revealed, leading to complex RFLP patterns. Cultivated clones appear to be less variable than prospections but show, however, a relatively high level of variability despite their narrow genetic base. Furthermore, IGS variations allowed a structuring within wild clones to be drawn, based mainly on their geographical origin. Some interesting discrepancies with previous work on isozyme variations are discussed and show the interest of surveying different genetic markers for diversity studies.     

Evolutionary dynamics of ecological niche in three Rhinogobio fishes from the upper Yangtze River inferred from morphological traits

In the past decades, it has been debated whether ecological niche should be conserved among closely related species (phylogenetic niche conservatism, PNC) or largely divergent (traditional ecological niche theory and ecological speciation) and whether niche specialist and generalist might remain in equilibrium or niche generalist could not appear. In this study, we employed morphological traits to describe ecological niche and test whether different niche dimensions exhibit disparate evolutionary patterns. We conducted our analysis on three Rhinogobio fish species (R. typus, R. cylindricus, and R. ventralis) from the upper Yangtze River, China. Among the 32 measured morphological traits except body length, PCA extracted the first four principal components with their loading scores >1.000. To find the PNC among species, Mantel tests were conducted with the Euclidean distances calculated from the four principal components (representing different niche dimensions) against the pairwise distances calculated from mitochondrial cytochrome b sequence variations. The results showed that the second and the third niche dimension, both related to swimming ability and behavior, exhibited phylogenetic conservatism. Further comparison on niche breadth among these three species revealed that the fourth dimension of R. typus showed the greatest width, indicating that this dimension exhibited niche generalism. In conclusion, our results suggested that different niche dimensions could show different evolutionary dynamic patterns: they may exhibit PNC or not, and some dimensions may evolve generalism.     

Correlation between local structural dynamics of proteins inferred from NMR ensembles and evolutionary dynamics of homologues of known structure

Conformational changes in proteins are extremely important for their biochemical functions. Correlation between inherent conformational variations in a protein and conformational differences in its homologues of known structure is still unclear. In this study, we have used a structural alphabet called Protein Blocks (PBs). PBs are used to perform abstraction of protein 3-D structures into a 1-D strings of 16 alphabets (a–p) based on dihedral angles of overlapping pentapeptides. We have analyzed the variations in local conformations in terms of PBs represented in the ensembles of 801 protein structures determined using NMR spectroscopy. In the analysis of concatenated data over all the residues in all the NMR ensembles, we observe that the overall nature of inherent local structural variations in NMR ensembles is similar to the nature of local structural differences in homologous proteins with a high correlation coefficient of .94. High correlation at the alignment positions corresponding to helical and β-sheet regions is only expected. However, the correlation coefficient by considering only the loop regions is also quite high (.91). Surprisingly, segregated position-wise analysis shows that this high correlation does not hold true to loop regions at the structurally equivalent positions in NMR ensembles and their homologues of known structure. This suggests that the general nature of local structural changes is unique; however most of the local structural variations in loop regions of NMR ensembles do not correlate to their local structural differences at structurally equivalent positions in homologues.