Exhaustive conformational search and simulated annealing for models of lattice peptides

We consider simple lattice models for short peptide chains whose states can be exhaustively enumerated to find the lowest energy conformation. Using these exact results and numerical simulations, we compute the distributions for the mean time t_N, required to find the global minimum energy state by simulated annealing (SA), as a function of N, the number of units in the chain. On the basis of scaling arguments, the time t_N, to find the global minimum energy of longer chains, beyond the range covered by exhaustive enumeration, can be estimated. On the basis of the observed exponential increase in folding time of the standard SA algorithms, it is imperative that better algorithms be found for minimizing longer chains. © 1993 John Wiley & Sons, Inc.     

Structural Basis for the Binding Affinity of a Homologous Series of Synthetic Phenoxazone Drugs with DNA: NMR and Molecular Mechanics Analysis

Abstract

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5′-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear ¹H-¹H 2D-NOESY and heteronuclear ¹H-³¹P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH₂ groups in the dimethylami- noalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH₃)₂ groups of the side chains.     

How human IgGs against DNA recognize oligonucleotides and DNA

In the literature, there are no available data on how anti‐DNA antibodies recognize DNA. In the present work, to study the molecular mechanism of DNA recognition by antibodies, we have used anti‐DNA IgGs from blood sera of patients with multiple sclerosis. A stepwise increase in ligand complexity approach was used to estimate the relative contributions of virtually every nucleotide unit of different single‐ (ss) and double‐stranded (ds) oligonucleotides to their affinity for IgG fraction having high affinity to DNA‐cellulose. DNA‐binding site disposed on the heavy chain demonstrates higher affinity to different dNMPs (K_d = 0.63μM‐3.8μM) than the site located on the light chain (28μM‐170μM). The heavy and light chains interact independently forming relatively strong contacts with 2 to 4 nucleotides of short homo‐ and hetero‐d(pN)_2‐9. Then the increase in the affinity of different d(pN)_n became minimal, and at n ≥ 8 to 9, all dependencies reached plateaus: approximately 3.2nM to 20nM and approximately 200nM to 460nM for the heavy and light chains, respectively. A similar situation was observed for different ribooligonucleotides, in which their affinity is 6‐fold to 100‐fold lower than that for d(pN)_n. Transition from ss to ds d(pN)_n leads to a moderate increase in affinity of ligands to DNA‐binding site of heavy chains, while light chains demonstrate the same affinity for ss and ds d(pN)_n. Long supercoiled DNA interacts with both heavy and light chains with affinity of approximately 10‐fold higher than that for short oligonucleotides. The thermodynamic models were constructed to describe the interactions of IgGs light and heavy chains with DNA.     

Transfection of Mycobacterium smegmatis SN2 with mycobacteriophage I3 DNA

Mycobacterium smegmatis SN2 does not exhibit natural competence for the uptake of phage I3 DNA. Competence can artificially be induced by treatment with glycine or CaCl₂, and the combination of both is even more effective. The efficiency of transfection can be improved by inclusion of protamine sulphate and heterologous RNA in the system. From ³²P DNA uptake studies the major barrier for the entry of DNA has been found to be the complex cell wall. The efficiency of transfection calculated on the basis of fraction of DNA which has entered the cell is comparable to that of other bacterial systems. The phage development takes a longer time (7 h for one cycle) after transfection, as compared to infection (4 h).     

Fracture of DNA in transient extensional flow. A numerical simulation study

Using the Brownian dynamics simulation technique, we studied the fracture process of DNA chains subjected to transient extensional flow, letting the solution with DNA molecules pass through a very small orifice (radius = 0.0065 cm), thus experiencing extensional flow of the convergent (sink) type. The DNA molecules were modeled as FENE bead-spring chains with the springs obeying a modified Warner force law, as proposed by Reese and Zimm. The fracture yield was strongly dependent on flow rate and molecular weight, reaching, in our setup, a level of 100% at a flow rate of around 0.001 cm³/s for DNA with molecular weight 26 × 10⁶ (T7 DNA). There was found to exist a critical flow rate (Q_crit) below which fracture did not occur, in accordance with what was observed in studies on polystyrene in transient extensional flow. We found that for DNA, the critical flow rate depended on the molecular weight as Q_crit ∼ M⁻¹⁴ when the hydrodynamic interaction effect (HI) was not included in the simulations. When HI was accounted for, the relation was found to be Q_crit ∼ M^−1.1, close to the theoretical prediction for fracture of partly extended chains in transient extensional flow. © 1996 John Wiley & Sons, Inc.     

Isolation of globin messenger RNA of Xenopus laevis

The polypeptide chains of Xenopus laevis hemoglobin have been analyzed by sodium dodecyl sulfate (SDS) and acid-urea gel electrophoresis. Four components can be distinguished, each having an approximate molecular weight of 13,000 daltons. Messenger RNA coding for the globin chains has been isolated and characterized. In a denaturing acrylamide gel the mRNA has an approximate molecular weight of 250,000 daltons. The complexity of the RNA is consistent with the presence of four different mRNA molecules, each of this molecular weight. When the mRNA is assayed in a wheat germ in vitro translation system, four polypeptides are synthesized corresponding to the four globin subunits. The relative proportion of the four synthesized polypeptides appears to vary according to the developmental stage of the red blood cells used for mRNA isolation. Hybridization of a complementary DNA (cDNA) copy of the globin mRNA to Xenopus laevis DNA in DNA excess indicates that each of the globin genes is present in one to three copies per haploid genome.     

Total intensity light scattering from short,optically anisotropic wormlike chains

Simple exact equations are derived for intensity light scattering from optically anisotropic wormlike chains in the absence of excluded volume. The results are valid at low scattering angles (q² 〈R²_G〉 ? 1) for all sormilke chains from rigid rods to random couils. The present work and an earlier theory [Nagai, K. (1972) Polym. J. 3 , 67–83] appear to be equivalent, although they were both derived using different methods. The present work is primarily concerned with short wormlike chains, since intensity light scattering from short fragments may provide valuable information about DNA flexibility. By using the results of this work to reanalyze some older light-scattering studies [Godfrey, J. E. & Eisenberg, H. (1976) Biophys. Chem. 5 , 301–318], it is shown that anisotropy corrections to polarized light-scattering measurements have been overcorrected in the past. It can be anticipated that future light-scattering experiments will determine the base-pair anisotropy.     

Proton magnetic resonance studies on conformation and association of oligo-γ-benzyl-L-glutamates in solution

The proton magnetic resonance spectra of o-nitrophenylthio-tetra- and hexa-γ-benzyl-L -glutamate ethylamides have been measured at different concentrations in CDCl₃ and CD₂₂C1. The NH and α-CH resonances of the tetrapeptide show downfield shifts with increasing concentration, accompanying disappearance of their fine structure and line broadening. The apparent feature of chemical shifts against concentration is sigmoidal, and it can be interpreted by assuming the presence of a step or more of association–dissociation equilibria of tetrapeptide. With increasing concentration, small aggregates are formed by intermolecular hydrogen bonding, the size of which is not sufficiently large to exhibit critical micelle concentrations. In contrast to the tetrapeptide, the hexapeptide has constant chemical shifts of the NH and α-CH resonances, independent of concentration, which implies that only the unassociated molecules show observable sharp resonances. In the hexapeptide, the phenyl CH and benzyl CH₂ groups of the side chains exhibit new resonances above certain critical concentrations, indicating the restriction of rotational freedom of the side chains in the aggregated states.     

Genetic transformation of the C3–C4 intermediate plant, Flaveria pubescens (Asteraceae)

The dicot genus Flaveria (Asteraceae), besides species with C₃ or C₄ photosynthesis, contains taxa with a broad range of different states of transition between the two major photosynthetic types. We have developed a reproducible and efficient Agrobacterium-mediated method for the stable genetic transformation of the C₃–C₄ intermediate species F. pubescens. Fusion constructs of the reporter gene β-glucuronidase (uidA, GUS) to several plant promoters, mainly derived from genes encoding subunits of the glycine cleavage system (gdcs), have been used to confirm the reproducibility and efficiency of the method. The stable integration of the foreign DNA has been examined by Southern analysis, kanamycin resistance, GUS enzyme activities and histochemical staining. Transformed shoots can be routinely obtained within 8–10 weeks after co-cultivation with A. tumefaciens. Received: 16 April 1996 / Revision received: 12 July 1996 / Accepted: 28 July 1996     

Mode‐coupling Smoluchowski dynamics of a double‐stranded DNA oligomer

The local dynamics of a double‐stranded DNA d(TpCpGpCpG)₂ is obtained to second order in the mode‐coupling expansion of the Smoluchowski diffusion theory. The time correlation functions of bond variables are derived and the ¹³C‐nmr spin–lattice relaxation times T₁ of different ¹³C along the chains are calculated and compared to experimental data from the literature at three frequencies. The DNA is considered as a fluctuating three‐dimensional structure undergoing rotational diffusion. The fluctuations are evaluated using molecular dynamics simulations, with the ensemble averages approximated by time averages along a trajectory of length 1 ns. Any technique for sampling the configurational space can be used as an alternative. For a fluctuating three‐dimensional (3D) structure using the three first‐order vector modes of lower rates, higher order basis sets of second‐rank tensor are built to give the required mode coupling dynamics. Second‐ and even first‐order theories are found to be in close agreement with the experimental results, especially at high frequency, where the differences in T₁ for ¹³C in the base pairs, sugar, and backbone are well described. These atomistic calculations are of general application for studying, on a molecular basis, the local dynamics of fluctuating 3D structures such as double‐helix DNA fragments, proteins, and protein–DNA complexes. © 1999 John Wiley & Sons, Inc. Biopoly 50: 613–629, 1999     

Interpretation of energy-transfer experiments by theoretical studies of model compounds using semiempirical potential functions. II. Monte Carlo calculations on oligopeptides

The properties relevant to nonradiative energy transfer have been computed in the unperturbed chain model for oligopeptides composed of from 4 to 21 residues of the formula Tyr-(Ala)_n-Tyr and Trp-(X)_n-Tyr, X being either Ala or Gly. A Monte Carlo method has been used for the generation of the chains. The relation between the distribution functions of the distances between the luminophores and the various properties in energy transfer has been examined for chains of different lengths and compositions. The averge of the orientation factor κ² has been computed as a functions of chain length both for molecules in a randomly coiled state and for molecules with backbones in a well-defined three-dimensional structure. The various averaging regimes of energy-transfer efficiency and of fluorescence decay are compared. Theoretical curves relating experimental efficiencies to the mean distance between the luminophores are proposed.     

Molecular analysis of the mitochondrial genome of Helianthus annuus in relation to cytoplasmic male sterility and phylogeny

Summary A circular supercoiled mitochondrial DNA plasmid P₁ (1.45 kb) is shown in both normal fertile plants of Helianthus annuus, and some cytoplasmic male sterile lines (CMS A and CMS P). In contrast, no plasmid is found in some other types of CMS C, I, B and K. A circular supercoiled DNA (P₂) of higher molecular weight (1.8 kb) is observed in CMS F. The mitochondrial plasmid P₁ was cloned, nick-translated and hybridized with native mitochondrial DNA from different lines of male fertile, CMS or wild Helianthus. No sequence homology has been detected between plasmid DNA P₁ and high molecular weight mitochondrial DNA in any line examined. A slight hybridization occurs between plasmids P₁ and P₂. Thus, there is no apparent relationship between mitochondrial plasmid DNA and CMS or Helianthus species. On the contrary, each Helianthus CMS and male fertile strain can be characterized by digestion fragment patterns (Sal I and Bgl I). Analysis of mitochondrial DNA from wild Helianthus strains indicated a relation between some CMS and the strain from which they were maternally derived, as for example CMS I and H. annuus ssp lenticularis and CMS F and H. petiolaris fallax. On the basis of restriction endonuclease patterns, a CMS phylogenic tree is proposed which illustrates a molecular polymorphism in the mitochondrial genome of Helianthus.     

Structural Elucidation of N-Linked Sugar Chains of Storage Glycoproteins in Mature Pea (Pisum sativum) Seeds by Ion-spray Tandem Mass Spectrometry (IS-MS/MS)

The structures of N-linked sugar chains of the storage glycoproteins in mature pea seeds have been estimated. Nine pyridylaminated (PA-) N-linked sugar chains were derived and purified from the hydrazinolysate of the storage glycoproteins by reversed-phase HPLC and size-fractionation HPLC. The structures of the PA-sugar chains purified were first identified by two-dimensional PA-sugar chain mapping, considering the data of sugar composition analysis or sequential exoglycosidase digestions. The deduced structures were further analyzed by IS-MS/MS analysis. Every relevant fragment ion derived from all PA-sugar chains could be assigned on the basis of deduced structures. The estimated nine structures fell into two categories; the first was a typical oligomannose type (Man_8-3GlcNAc₂; 77.7%) which can be hydrolyzed by endo-β-N-acetylglucosaminidase PS [Y. Kimura et al., Biosci. Biotech. Biochem., 60, 228–232 (1996)], the second was a xylose-containing type (Man_4-3Xyl₁GlcNAc₂, Man₃Fuc₁Xyl₁GlcNAc₂; 22.3%). Among these structures, Man₈GlcNAc₂ (19.7%), Man₆GlcNAc₂ (24.7%), and Man₃Fuc₁Xyl₁GlcNAc₂ (18.8%) were the dominant structures.     

The structure of native cellulose

Native cellulose has been shown to consist of a crystalline array of parallel chains, based on the X-ray diffraction data for specimens from the sea alga Valonia ventricosa. The unit cell is monoclinic with dimensions a = 16.34 Å, b = 15.72 Å, c = 10.38 Å (fiber axis), and β = 97.0°. The space group is P2₁ and the cell contains disaccharide segments of eight chains. Models containing chains with the same sense (parallel) or alternating sense (antiparallel) were refined against the intensity data using rigidbody least squares procedures. The results show a preference for a parallel chain structure with specific chain polarity with respect to the c axis. The refinement places the ? CH₂OH side chains approximately 20′ from the so-called tg conformation, with a result that an 02′? H…06 intramolecular bond is formed. The structure also contains an 03? H…05′ intramolecular bond and an 06? H…03 intermolecular bond along the a axis. All these bonds lie in the 020 planes, and the structure is an array of hydrogen-bonded sheets. A major consequence of this work is that regular chain folding can be ruled out and cellulose is seen as extended chain polymer single crystals.     

Studies on the Interaction of [SnMe2Cl2(bu2bpy)] Complex with ct-DNA Using Multispectroscopic,Atomic Force Microscopy (AFM) and Molecular Docking

The interaction of SnMe₂Cl₂(bu₂bpy)complex with calf thymus DNA (ct-DNA) has been explored following, using spectroscopic methods, viscosity measurements, Atomic force microscopy, Thermal denaturation and Molecular docking. It was found that Sn(IV) complex could bind with DNA via intercalation mode as evidenced by hyperchromism and bathochromic in UV–Vis spectrum; these spectral characteristics suggest that the Sn(IV) complex interacts with DNA most likely through a mode that involves a stacking interaction between the aromatic chromophore and the base pairs of DNA. In addition, the fluorescence emission spectra of intercalated methylene blue (MB) with increasing concentrations of SnMe₂Cl₂(bu₂bpy) represented a significant increase of MB intensity as to release MB from MB-DNA system. Positive values of ΔH and ΔS imply that the complex is bound to ct-DNA mainly via the hydrophobic attraction. Large complexes contain the DNA chains with an average size of 859?nm were observed by using AFM for Sn(IV) Complex–DNA. The Fourier transform infrared study showed a major interaction of Sn(IV) complex with G-C and A-T base pairs and a minor perturbation of the backbone PO₂ group. Addition of the Sn(IV)complex results in a noticeable rise in the Tm of DNA. In addition, the results of viscosity measurements suggest that SnMe₂Cl₂(bu₂bpy) complex may bind with the classical intercalative mode. From spectroscopic and hydrodynamic studies, it has been found that Sn(IV)complex interacts with DNA by intercalation mode. Optimized docked model of DNA–complex mixture confirmed the experimental results.     

A human homolog of the mouse CD8 molecule,Lyt-3: genomic sequence and expression

The CD8 antigen is a marker for those cytotoxic T cells that recognize antigen in the context of class I major histocompatibility antigens (MHC) and has now been identified in many species. In rodents the CD8 antigen is a heterodimer of two distinct chains, Lyt-2 and Lyt-3 in the mouse and OX-8 M _r 32 000 and 37 000 chains in the rat. Human CD8 has consistently been described as a homodimer/homomultimer on mature T cells made up of one chain homologous to the Lyt-2 and OX-8 M _r 32 000 chains. This paper identifies a human equivalent of the second rodent CD8 chain (Lyt-3 and OX-8 M _r 37 000 chains) at the genomic level and shows that this gene is transcribed in human thymocytes and in some acute leukemic T-cell lines. The existence of a human Lyt-3 homolog raises the possibility that human CD8, like mouse CD8, may exist as a heterodimer.     

Specificity of DNA-basic polypeptide interactions. III. Sequential and random copolymers of lysine and aromatic amino acids.

Sequential and random lysine copolymers containing various amounts of different aromatic amino acids were synthesized. The sequential copolypeptides exhibited strong dependence of yield and degree of polymerization on the amino acid sequence of the repeating unit. To elucidate the specific contributions of aromatic side chains to the interaction of these copolymers with DNA, direct-mixed complexes were studied by thermal denaturation and CD. The melting behaviour of peptide-bound DNA was found to be strongly affected by amino acid composition and sequence. The contribution of the different aromatic amino acids to thermal stability decreased in the order: polylysine > [Lys, Tyr]_n > [Lys,Phe]_n > [Lys,(OMe)Tyr]_n. The CD spectrum of DNA was altered by random copolymers, whereas sequential copolymers exhibited no changes. The influence of the random copolymers on the CD spectrum of DNA decreased in the series: polylysine > [Lys,Phe]_n > [Lys,(OMe)Tyr]_n > [Lys,Tyr]_n. The contribution of the different aromatic amino acids to thermal stability is interpreted as stacking tendencies toward denatured and, in the case of Tyr, H-bond formation with native DNA. The differences found for the random and the sequential polypeptides can best be explained by assuming a cooperative action of rather small peptide segments.     

The mouse C1q genes are clustered on chromosome 4 and show conservation of gene organization

Mouse complement component C1q is a serum glycoprotein which consists of six A chains, six B chains and six C chains. The three polypeptides are 223, 228, and 217 residues long, respectively, and are encoded by three genes. DNA probes for mouse C1q A, B, and C chains were hybridized to Southern blots of DNA obtained from various inbred mouse strains. On the basis of fragment length polymorphisms, two different alleles of each of the genes could be identified. The distribution of these alleles was determined in the BXD and LXPL recombinant inbred strain series. Comparison with previously reported strain distribution patterns shows that the genes encoding mouseClq map to the same locus on distal chromosome 4. Overlapping clones spanning the entire gene cluster ofClq were isolated from genomic libraries using specific cDNA probes. The three genesClqA, ClqB, andClqC are closely arranged on a 19 kilobase stretch of DNA in the 5 to 3 orientation A-C-B. Each gene consists of two exons separated by one intron. Sequence comparison of Clq from three different species have shown that the B chains have the strongest similarity. Southern blot analysis of chromosomal DNA from 14 vertebrate species demonstrated highest similarity between theClqB genes, followed byClqC and finallyClqA.The nucleotide sequence data reported in this paper have been submitted to the EMBL, GenBank, and DDBJ nucleotide sequence databases and have been assigned the accession numbers X92958 (ClqA), X92959 (ClqB), and X92960 (ClqC)     

Excluded volume of an intermediate-molecular-weight DNA. A Monte Carlo analysis

A Monte Carlo procedure wasused to determine the effect of excluded volume on the dimensions of an intermediate-molecular-weight DNA for different Na⁺ concentrations. The calculation of α, the parameter for the linear expansion due to excluded volume, was accomplished by generating sets of chains and, for each set, comparing the average radius of gyration for the set of chains that do not overlap to that averaged over the entire set of chains. Each chain was defined by cylinders linked with free rotation and with bend angles generated according to a weighted Gaussian distribution. The chain parameters—contour lenght, cylinder lenght and diameter—were fixed in order to resemble published light-scattering experiments on Col E₁ DNA. Values for α were less than 1.08. for Na⁺ concentrations between 0.007 and 1.0M. A previously reported analytical calculation of the excluded-volume correction of intermediate-sized DNA gave results that are closely similar to those from the Monte Carlo analysis.