Compositional properties of human cDNA libraries: practical implications

The strikingly wide and bimodal gene distribution exhibited by the human genome has prompted us to study the correlations between EST-counts (expression levels) and base composition of genes, especially since existing data are contradictory. Here we investigate how cDNA library preparation affects the GC distributions of ESTs and/or genes found in the library, and address consequences for expression studies. We observe that strongly anomalous GC distributions often indicate experimental biases or deficits during their preparation. We propose the use of compositional distributions of raw ESTs from a cDNA library, and/or of the genes they represent, as a simple and effective tool for quality control.     

The structure of lombricine kinase: implications for phosphagen kinase conformational changes

Lombricine kinase is a member of the phosphagen kinase family and a homolog of creatine and arginine kinases, enzymes responsible for buffering cellular ATP levels. Structures of lombricine kinase from the marine worm Urechis caupo were determined by x-ray crystallography. One form was crystallized as a nucleotide complex, and the other was substrate-free. The two structures are similar to each other and more similar to the substrate-free forms of homologs than to the substrate-bound forms of the other phosphagen kinases. Active site specificity loop 309-317, which is disordered in substrate-free structures of homologs and is known from the NMR of arginine kinase to be inherently dynamic, is resolved in both lombricine kinase structures, providing an improved basis for understanding the loop dynamics. Phosphagen kinases undergo a segmented closing on substrate binding, but the lombricine kinase ADP complex is in the open form more typical of substrate-free homologs. Through a comparison with prior complexes of intermediate structure, a correlation was revealed between the overall enzyme conformation and the substrate interactions of His(178). Comparative modeling provides a rationale for the more relaxed specificity of these kinases, of which the natural substrates are among the largest of the phosphagen substrates.     

Representation of DNA sequences in recombinant DNA libraries prepared by restriction enzyme partial digestion

We present a theoretical study of the fraction of sequences incorporated in a recombinant DNA partial digest library as a function of the size of the library. The fraction incorporated depends on the degree of restriction enzyme partial digestion. If all restriction sites in the target DNA can be cleaved with the same rate, optimum incorporation of sequences is observed when the number average length of the digested DNA equals the desired average length of the cloned insert. Overdigestion severely reduces the fraction of sequences present in a sample of clones. Heterogeneity in restriction enzyme cleavage rates also reduces the fraction incorporated, and underdigestion improves sequence representation in the face of cleavage rate heterogeneity. Practical methods for determining the number average length of partially digested DNAs are also presented.     

A small tripeptide AFA undergoes two state cooperative conformational transitions: implications for conformational biases in unfolded states

It is important to understand the conformational biases that are present in unfolded states to understand protein folding. In this context, it is surprising that even a short tripeptide like AFA samples folded/ordered conformation as demonstrated recently by NMR experiments of the peptide in aqueous solution at 280 K. In this paper, we present molecular dynamics simulation of the peptide in explicit water using OPLS-AA/L all-atom force field. The results are in overall agreement with NMR results and provide some further insights. The peptide samples turn and extended conformational forms corresponding to minima in free energy landscape. Frequent transitions between the minima are observed due to modest free energy barriers. The turn conformation seems to be stabilized by hydrophobic interactions and possibly by bridging water molecules between backbone donors and acceptors. Thus the peptide does not sample conformations randomly, but samples well defined conformations. The peptide served as a model for folding-unfolding equilibrium in the context of peptide folding. Further, implications for drug design are also discussed.     

Calcium-induced conformational changes of Thrombospondin-1 signature domain: implications for vascular disease

Context: Thrombospondin1 (TSP1) participates in numerous signaling pathways critical for vascular physiology and disease. The conserved signature domain of thrombospondin 1 (TSP1-Sig1) comprises three epidermal growth factor (EGF), 13 calcium-binding type 3 thrombospondin (T3) repeats, and one lectin-like module arranged in a stalk-wire-globe topology. TSP1 is known to be present in both calcium-replete (Holo-) and calcium-depleted (Apo-) state, each with distinct downstream signaling effects.

Objective: To prepare a homology model of TSP1-Sig1 and investigate the effect of calcium on its dynamic structure and interactions.

Methods: A homology model of Holo-TSP1-Sig1 was prepared with TSP2 as template in Swissmodel workspace. The Apo-form of the model was obtained by omitting the bound calcium ions from the homology model. Molecular dynamics (MD) simulation studies (100?ns) were performed on the Holo- and Apo- forms of TSP1 using Gromacs4.6.5.

Results and discussion: After simulation, Holo-TSP1-Sig1 showed significant reorientation at the interface of the EGF1-2 and EGF2-3 modules. The T3 wire is predicted to show the maximum mobility and deviation from the initial model. In Apo-TSP1-Sig1 model, the T3 repeats unfolded and formed coils with predicted increase in flexibility. Apo-TSP1-Sig1model also predicted the exposure of the binding sites for neutrophil elastase, integrin and fibroblast growth factor 2. We present a structural model and hypothesis for the role of TSP1-Sig1 interactions in the development of vascular disorders.

Conclusion: The simulated model of the fully calcium-loaded and calcium-depleted TSP1-Sig1 may enable the development of its interactions as a novel therapeutic target for the treatment of vascular diseases.     

Distal heme pocket conformers of carbonmonoxy derivatives of Ascaris hemoglobin: evidence of conformational trapping in porous sol-gel matrices

We report the ligand dependence of the conformer distribution in the distal heme pocket of Ascaris suum hemoglobin (Hb) studied by resonance Raman spectroscopy. The heme-bound CO is used as a spectroscopic antenna to probe the original distribution of conformers in the dioxygen derivative of Ascaris Hb, by utilizing sol-gel encapsulation. The first step is to encapsulate the dioxygen derivative in the porous sol-gel and let the gel age, thus trapping the equilibrium conformational distribution of Ascaris dioxygen Hb. In the second step, the dioxygen ligand is replaced by CO. The sol-gel environment impedes any large scale movements, drastically slowing down the conformational relaxation triggered by the ligation change, essentially "locking in" the initial quaternary and even tertiary structure of the protein. Studying the Fe-CO frequencies of the latter sample allows evaluation of the distribution of the distal heme pocket conformers that was originally associated with the dioxygen derivative. Extending the study to the Ascaris mutants allows for examination of the effect of specific residues in the distal pocket on the conformational distribution. The choice of mutants was largely based on the anticipated variation in hydrogen bonding patterns. The results show that the sol-gel encapsulation can slow or prevent re-equilibration within the distal heme pocket of Ascaris Hb and that the distribution of distal heme pocket conformers for the CO derivative of Ascaris Hb in the sol-gel is highly dependent on the history of the sample. Additionally, we report a detailed study of the CO complex of the mutants in solution for assignment of the various heme pocket conformers, and we present a comparison of the sol-gel data with solution data. The results support a picture in which the dioxygen derivative biases the population strongly toward a tightly packed configuration that favors the network of strong hydrogen bonding interactions, and suggest that Ascaris Hb is uniquely designed for dioxygen capture.     

Lysosomal ceroid depletion by drugs: therapeutic implications for a hereditary neurodegenerative disease of childhood

Neuronal ceroid lipofuscinoses (NCLs) are the most common hereditary neurodegenerative diseases of childhood. The infantile form, INCL, is caused by lysosomal palmitoyl-protein thioesterase (PPT) deficiency, which impairs the cleavage of thioester linkages in palmitoylated proteins, preventing their hydrolysis by lysosomal proteinases. Consequent accumulation of these lipid-modified proteins (constituents of ceroid) in lysosomes leads to INCL. Because thioester linkages are susceptible to nucleophilic attack, drugs with this property may have therapeutic potential for INCL. We report here that two such drugs, phosphocysteamine and N-acetylcysteine, disrupt thioester linkages in a model thioester compound, [14C]palmitoyl approximately CoA. Most importantly, in lymphoblasts derived from INCL patients, phosphocysteamine, a known lysosomotrophic drug, mediates the depletion of lysosomal ceroids, prevents their re-accumulation and inhibits apoptosis. Our results define a novel pharmacological approach to lysosomal ceroid depletion and raise the possibility that nucleophilic drugs such as phosphocysteamine hold therapeutic potential for INCL.     

Limited sampling of conformational space by the distance geometry algorithm: implications for structures generated from NMR data

Calculations with a metric matrix distance geometry algorithm were performed that show that the standard implementation of the algorithm generally samples a very limited region of conformational space. This problem is most severe when only a small amount of distance information is used as input for the algorithm. Control calculations were performed on linear peptides, disulfide-linked peptides, and a double-stranded DNA decamer where only distances defining the covalent structures of the molecules (as well as the hydrogen bonds for the base pairs in the DNA) were included as input. Since the distance geometry algorithm is commonly used to generate structures of biopolymers from distance data obtained from NMR experiments, simulations were performed on the small globular protein basic pancreatic trypsin inhibitor (BPTI) that mimic calculations performed with actual NMR data. The results on BPTI and on the control peptides indicate that the standard implementation of the algorithm has two main problems: first, that it generates extended structures; second, that it has a tendency to consistently produce similar structures instead of sampling all structures consistent with the input distance information. These results also show that use of a simple root-mean-square deviation for evaluating the quality of the structures generated from NMR data may not be generally appropriate. The main sources of these problems are identified, and our results indicate that the problems are not a fundamental property of the distance geometry algorithm but arise from the implementations presently used to generate structures from NMR data. Several possible methods for alleviating these problems are discussed.     

An analysis of conformational changes on protein-protein association: implications for predictive docking.

Conformational changes on complex formation have been measured for 39 pairs of structures of complexed proteins and unbound equivalents, averaged over interface and non-interface regions and for individual residues. We evaluate their significance by comparison with the differences seen in 12 pairs of independently solved structures of identical proteins, and find that just over half have some substantial overall movement. Movements involve main chains as well as side chains, and large changes in the interface are closely involved with complex formation, while those of exposed non-interface residues are caused by flexibility and disorder. Interface movements in enzymes are similar in extent to those of inhibitors. All eight of the complexes (six enzyme-inhibitor and two antibody-antigen) that have structures of both components in an unbound form available show some significant interface movement. However, predictive docking is successful even when some of the largest changes occur. We note however that the situation may be different in systems other than the enzyme-inhibitors which dominate this study. Thus the general model is induced fit but, because there is only limited conformational change in many systems, recognition can be treated as lock and key to a first approximation.     

Perceived pattern regularity computed as a summary statistic: implications for camouflage

Why do the equally spaced dots in figure 1 appear regularly spaced? The answer ‘because they are’ is naive and ignores the existence of sensory noise, which is known to limit the accuracy of positional localization. Actually, all the dots in figure 1 have been physically perturbed, but in the case of the apparently regular patterns to an extent that is below threshold for reliable detection. Only when retinal pathology causes severe distortions do regular grids appear perturbed. Here, we present evidence that low-level sensory noise does indeed corrupt the encoding of relative spatial position, and limits the accuracy with which observers can detect real distortions. The noise is equivalent to a Gaussian random variable with a standard deviation of approximately 5 per cent of the inter-element spacing. The just-noticeable difference in positional distortion between two patterns is smallest when neither of them is perfectly regular. The computation of variance is statistically inefficient, typically using only five or six of the available dots.     

The Mad2 conformational dimer: structure and implications for the spindle assembly checkpoint

The 25 kDa Mad2 protein is a key player in the spindle assembly checkpoint, a safeguard against chromosome segregation errors in mitosis. Mad2 combines three unusual properties. First, Mad2 adopts two conformations with distinct topologies, open (O) and closed (C) Mad2. Second, C-Mad2 forms topological links with its two best-characterized protein ligands, Mad1 and Cdc20. Third, O-Mad2 and C-Mad2 engage in a "conformational" dimer that is essential for spindle checkpoint function in different organisms. The crystal structure of the O-Mad2-C-Mad2 conformational dimer, reported here, reveals an asymmetric interface that explains the selective dimerization of the O-Mad2 and C-Mad2 conformers. The structure also identifies several buried hydrophobic residues whose rearrangement correlates with the Mad2 topological change. The structure of the O-Mad2-C-Mad2 conformational dimer is consistent with a catalytic model in which a C-Mad2 template facilitates the binding of O-Mad2 to Cdc20, the target of Mad2 in the spindle checkpoint.     

Novel conformational states of peptide deformylase from pathogenic bacterium Leptospira interrogans: implications for population shift

Peptide deformylase is an attractive target for developing novel antibiotics. Previous studies at pH 3.0 showed peptide deformylase from Leptospira interrogans (LiPDF) exists as a dimer in which one monomer is in a closed form and the other is in an open form, with different conformations of the CD-loop controlling the entrance to the active pocket. Here we present structures of LiPDF at its active pH range. LiPDF forms a similar dimer at pH values 6.5-8.0 as it does at pH 3.0. Interestingly, both of the monomers are almost in the same closed form as that observed at pH 3.0. However, when the enzyme is complexed with the natural inhibitor actinotin, the conformation of the CD-loop is half-open. Two pairs of Arg109-mediated cation-pi interactions, as well as hydrogen bonds, have been identified to stabilize the different CD-loop conformations. These results indicate that LiPDF may be found in different structural states, a feature that has never before been observed in the peptide deformylase family. Based on our results, a novel substrate binding model, featured by an equilibrium between the closed and the open forms, is proposed. Our results present crystallographic evidence supporting population shift theory, which is distinguished from the conventional lock-and-key or induced-fit models. These results not only facilitate the development of peptide deformylase-targeted drugs but also provide structural insights into the mechanism of an unusual type of protein binding event.     

Predominant torsional forms adopted by oligopeptide conformers in solution: parameters for molecular recognition.

In this paper, we describe the predominant conformational forms adopted by tripeptides and higher oligopeptides in aqueous solution. About 50 tripeptides and almost 20 higher oligopeptides (4-6 residues) were subjected to conformational analysis using SYBYL Random Search. As with dipeptides (Grail BM, Payne JW. J. Peptide Sci. 2000; 6: 186-199), both tripeptides and higher oligopeptides were found to occupy relatively few combinations of psi-phi space that were distinct from those associated with predominant protein secondary structures (e.g. helices and beta-sheets). Again, the preferred psi (psi) values for the first residue (i - 1) were in sectors encompassed by the ranges from +150 degrees to +/-180 degrees, +60 degrees to +90 degrees and -60 degrees to -90 degrees, which were combined with preferred phi (phi) values for the second residue (i) in sectors with ranges from -150 degrees to +/-180 degrees, -60 degrees to -90 degrees and +30 degrees to +60 degrees. It was notable that tripeptides and, to a greater extent, higher oligopeptides adopted an initial psi (psi) (Tor2) from +150 degrees to +/-180 degrees. For tripeptides, their N-C distances (distance between N-terminal nitrogen and C-terminal carbon atoms) distribute about 6.5 A to give shorter, 'folded' conformers that are similar in length to dipeptides, and longer, 'extended' conformers that are distinct. Furthermore, for higher oligopeptides, their N-C distances did not increment in relation to their increasing number of residues and short, 'folded' conformers were still present. These findings have a bearing upon the recognition of these molecules as substrates for widely distributed peptidases and peptide transporters.     

Anomalous electrochemical behavior of clostripain: Evidence for disulfide isolation conformers

Summary Purified preparations of clostripain exhibit two distinct components on analytical and preparative acrylamide gel electrophoresis as well as adsorption chromatography on hydroxylapatite. Both components are of identical molecular size and specific activity. By reducing the enzyme for an extended period of time prior to chromatography, the specific activity increases by a factor of four and the enzyme elutes from the hydroxylapatite column as a homogeneous peak. Enzyme labeled at the active site with³H-TLCK exhibits a similar chromatographic behavior to native enzyme on hydroxylapatite.It is inferred that such behavior may be attributed to a two phase disulfide reduction, one involving reduction of a disulfide thereby freeing an active site SH group and a second disulfide reduction resulting in the chromatographic transition.     

Predominant torsional forms adopted by dipeptide conformers in solution: parameters for molecular recognition.

The present paper describes the predominant conformational forms adopted by dipeptides in aqueous solution. More than 50 dipeptides were subjected to conformational analysis using SYBYL Random Search. The resultant collections of conformers for individual dipeptides, for small groups with related side chain residues and for large groups of about 50 dipeptides were visualized graphically and analysed using a novel three-dimensional pseudo-Ramachandran plot. The distribution of conformers, weighted according to the percentage of each in the total conformer pool, was found to be restricted to nine main combinations of backbone psi (psi) and phi (phi) torsion angles. The preferred psi values were in sectors A7 (+150 degrees to +/-180 degrees), A10 (+60 degrees to +90 degrees) and A4 (-60 degrees to -90 degrees), and these were combined with preferred phi values in sectors B12 (-150 degrees to +/-180 degrees), B9 (-60 degrees to -90 degrees) and B2 (+30 degrees to +60 degrees). These combinations of psi and phi values are distinct from those found in common secondary structures of proteins. These results show that although dipeptides can each adopt many conformations in solution, each possesses a profile of common conformers that is quantifiable. A similarly weighted distribution of dipeptide conformers according to distance between amino-terminal nitrogen and carboxyl-terminal carbon shows how the preferred combinations of backbone torsional angles result in particular N-C geometries for the conformers. This approach gives insight into the important conformational parameters of dipeptides that provide the basis for their molecular recognition as substrates by widely distributed peptide transporters. It offers a basis for the rational design of peptide-based bioactive compounds able to exploit these transporters for targeting and delivery.     

Effects of backbone substitutions on the conformational behavior of Shigella flexneri O-antigens: implications for vaccine strategy

The O-antigen (O-Ag), the polysaccharide part of the lipopolysaccharide, is the major target of the serotype-specific protective humoral response elicited upon host infection by Shigella flexneri, the main causal agent of the endemic form of bacillary dysentery. The O-Ag repeat units (RUs) of 12 S. flexneri serotypes share the tetrasaccharide backbone →2)-α-l-Rhap-(1?→?2)-α-l-Rhap-(1?→?3)-α-l-Rhap-(1?→?3)-β-d-GlcpNAc-(1→, with site-selective glucosylation(s) and/or O-acetylation defining the serotypes. To investigate the conformational basis of serotype specificity, we sampled conformational behaviors during 60?ns of molecular dynamic simulations for oligosaccharides representing three RUs of each one of the O-Ags corresponding to serotypes 1a, 1b, 2a, 2b, 3a, 3b, 4a, 4b, 5a, 5b, X and Y, respectively. The calculated trajectories were checked by nuclear magnetic resonance (NMR) for 1a, 2a, 3a and 5a O-Ags. The simulations predict that in all O-Ags, but 1a and 1b, serotype-specific substitutions of the backbone do not induce any new backbone conformations compared with the linear type O-Ag Y, although they restrain locally the accessible conformational space. Moreover, the influence of any given substituent on the backbone is independent of the eventual presence of other substituents. Finally, only slight differences in conformational behavior between terminal and inner RUs were observed. These results suggest that the reported serotype-specificity of the protective immune response is not due to recognition of distinct backbone conformations, but to binding of the serotype-defining substituents in the O-Ag context. The gained knowledge is discussed in terms of impact on the development of a broad-serotype coverage vaccine.     

Sequence context of oligomer tracts in eukaryotic DNA: biological and conformational implications

Recent studies of homooligomer tracts suggest different characteristics from random sequence DNA (dA).(dT) and (dG).(dC) tracts are frequent in upstream regions and in some cases have been shown to be essential for regulation. Here we examine homooligomer occurrences in non-coding and coding eukaryotic sequences, focusing on the context in which the homooligomers occur. This analysis of sequences in the junction areas yields distinct and consistent characteristics. In particular, the nucleotide interrupting a run is most frequently complementary to the run. The base next to it is most frequently identical to the one constituting the run. For A or T runs the least frequent nearest and next to nearest neighbors are G or C. For G or C tracts the least frequent are A or T. Complementary oligomers behave similarly. These and additional trends are strongest for run lengths greater than or equal to 3. The computations are carried out on the whole eukaryotic database of greater than 4 x 10(6) nucleotides, separately for coding and non-coding regions. These same trends are evident for both groups, but are somewhat stronger for the non-coding regions. The context in which the homooligomers occur may yield some clues to DNA conformation and its biological implications.     

Cofactor triggers the conformational change in thymidylate synthase: implications for an ordered binding mechanism.

We have solved crystal structures of two complexes with Escherichia coli thymidylate synthase (TS) bound either to the cofactor analog N10-propargyl-5,8-dideazafolate (CB3717) or to a tighter binding polygutamyl derivative of CB3717. These structures suggest that cofactor binding alone is sufficient to induce the conformational change in TS; dUMP binding is not required. Because polyglutamyl folates are the primary cofactor form in vivo, and because they can bind more tightly than dUMP to TS, these structures may represent a key intermediate along the TS reaction pathway. These structures further suggest that the dUMP binding site is accessible in the TS-cofactor analog binary complexes. Conformational flexibility of the binary complex may permit dUMP to enter the active site of TS while the cofactor is bound. Alternatively, dUMP may enter the active site from the opposite side that the cofactor appears to enter; that is, through a portal flanked by arginines that also coordinate the phosphate group in the active site. Entry of dUMP through this portal may allow dUMP to bind to a TS-cofactor binary complex in which the complex has completed its conformational transition to the catalytically competent structure.