Calculation of the solution properties of flexible macromolecules: methods and applications

While the prediction of hydrodynamic properties of rigid particles is nowadays feasible using simple and efficient computer programs, the calculation of such properties and, in general, the dynamic behavior of flexible macromolecules has not reached a similar situation. Although the theories are available, usually the computational work is done using solutions specific for each problem. We intend to develop computer programs that would greatly facilitate the task of predicting solution behavior of flexible macromolecules. In this paper, we first present an overview of the two approaches that are most practical: the Monte Carlo rigid-body treatment, and the Brownian dynamics simulation technique. The Monte Carlo procedure is based on the calculation of properties for instantaneous conformations of the macromolecule that are regarded as if they were instantaneously rigid. We describe how a Monte Carlo program can be interfaced to the programs in the HYDRO suite for rigid particles, and provide an example of such calculation, for a hypothetical particle: a protein with two domains connected by a flexible linker. We also describe briefly the essentials of Brownian dynamics, and propose a general mechanical model that includes several kinds of intramolecular interactions, such as bending, internal rotation, excluded volume effects, etc. We provide an example of the application of this methodology to the dynamics of a semiflexible, wormlike DNA.     

Representation of a protein molecule as a tree and application to modular computer programs which calculate and modify atomic coordinates

The concept and representation of a logical tree as defined in computer science is applied to obtain a suitable representation of protein molecules in computer programs which handle or calculate atomic coordinates of protein molecules. On the basis of this analysis and of the analysis of the calculation and modification of the structure of a protein from bond lengths, bond angles, and dihedral angles, which is reproduced in an appendix, program modules which accomplish the various required computations are described. Three such modules are given in Iverson notation; in fact, it is hoped that this article will serve as a reasonably complete basis for the preparation of machine programs by moderately proficient programmers.     

Melting of Cross-Linked DNA: I. Model and Theoretical Methods

Abstract

Covalent and strong coordination binding to DNA of a large number of antitumour drugs and other compounds leads to interstrand cross-link formation. To investigate cross-link influence on double helix stability, two methods are developed for the calculation of melting curves. The first method is based on Poland's approach. It requires computer time proportional to u-N, where u is the average distance (in base pairs) between neighboring crosslinks, and N is the number of base pairs in the DNA chain. The method is more suitable when u is not large, and small loops formed by interstrand cross-links in melted regions strongly affect DNA melting. The computer time for the second method, based on the Fixman-Freire approach, does not depend on the number of cross-links and is proportional to I N (I is the number of exponential functions used for a decomposition of the loop entropy factor). It is more appropriate when N and u are large, and therefore particular values of the entropy factors of small loops do not influence DNA melting behavior.     

Accessible and total molecular surfaces of unsmooth protein globules

A method for the calculation of accessible and total molecular surfaces of protein globules is proposed. It is based on consideration of structuring of protein globule surface. As a whole, the protein globule is approximated by ellipsoids. By computer imitations, the areas of accessible and total molecular surfaces as well as the volumes for 58 globular proteins were calculated. The calculation showed that the best agreement between the calculated areas and volumes and the experimental results is obtained at the radii of protein surface structures of 1-2 A.     

A stereochemical model for DNA-Na+ (H20) complex]

By computer calculation with the use of atom-atom potentials a stereochemical model of the DNA-hydrate Na+ complex was obtained according to which the Na+(H2O)6 octahedrons are localized in the narrow groove of DNA with formation of a great number of van der vaals contacts and hydrogen bonds. The model explains why the C-form of DNA is stabilized in the Na-DNA complex which is formed in the presence of 80% methanol (v/v). A possibility of the existence of such complexes in vivo at the DNA regions surrounded with media of diminished water activity (the DNA-membrane complex, chromosomes, phage heads) is discussed.     

The Mechanism of Inactivation of T4 Bacteriophage by Tritium Decay

Coliphage T4 was used as a model system to study the mechanism of biological inactivation produced by tritium decay. Experimentally, tritiated precursors were incorporated into phage DNA (thymidine-³H) or into phage protein (³H-amino acids). The ratio of killing efficiencies for decays originating in phage DNA to those originating in phage protein was 2.6. Inactivation by decays from labeled amino acids was assumed to occur exclusively from β-particle irradiation of phage DNA. If decays originating in DNA are due solely to irradiation of DNA, then the killing efficiencies reflect the energy transfer paths in phage DNA for decays originating in phage DNA and in the protein coat. The energy transfer paths were determined for the two cases with the help of a computer and found to be very nearly equal to the experimentally determined ratio (2.6). The killing efficiencies for decays originating in phage DNA were 0.12 and for decays originating in protein 0.046.     

Determination of the molecular weight of DNA-bound protein(s) responsible for gel electrophoretic mobility shift of linear DNA fragments examplified with purified viral myb protein.

A protein-DNA complex has less gel electrophoretic mobility than the free DNA fragment. One parameter for the degree of retardation of a linear DNA fragment in a protein-DNA complex is the molecular weight of the bound protein(s). The quotient of the migration distances of free DNA (m) and protein-DNA complex (m') is a function of the molecular weight (MW) of the bound protein(s). Based on the evaluation of the lac repressor induced mobility shift of a 203 bp DNA fragment containing the lac operator in a 5% non-denaturating polyacrylamide gel a direct proportionality could be shown between (m/m'-1) and MW with the proportionality factor K = 215 kDa. The factor K depends on the acrylamide concentration in the gel, getting lower values with increasing acrylamide concentrations. A calculation is given to determine the molecular weight of DNA-binding factors responsible for the decreased electrophoretic mobility of a linear DNA fragment. As an example this calculation was used in order to analyse DNA-binding of the isolated viral myb protein. It could be demonstrated that the viral myb protein binds to DNA as a monomer and as a dimer.     

Analog-digital molecular cell computer.

The living cell is controlled by a molecular stochastic computer of parallel-successive action. MCC may be equivalent to the universal computer. MCC operates with molecule-words (DNA, RNA, proteins) according to the programme recorded in DNA and RNA. Operations are produced by molecular devices (RNA- and DNA-polymerases, ligases, proteinases and so on). Molecular devices operating with molecule-words are recorded on molecules themselves, and they are read off by ribosomes. Therefore the programme of the reorganization of the programme itself may be recorded on the molecule-words. MCC operates with molecular words having definite addresses. The words and the operators collide by Brownian movement and combine if the molecular surface of address segment is complementary and properly oriented. It is possible to reproduce not only the programmes but also the operators of MCC. The molecular computer operates with word-molecules according to the programme, recorded in DNA, with the aim of predicting an outer situation in the next time-moment and selecting of a correct answer by synthesis of suitable proteins and other substances and also by macroscopic motion. Each step of directed calculation is needed of the consumption of minimally necessary portion of free energy and search is due to the Brownian movement without free energy loss. Cyclic nucleotides are intraneuron membrane potential controlling systems which can be described as molecular diffusional analog computers, well fitted to solve mathematical physics equations if there are high frequency generators and regulators of cyclic nucleotides sources and sinks. It is suggested that molecular proton channels in an electric field are such generators of 10(11)-10(12) Hz. Biophysics cannot use the ordinary laws of physics and must take into account the influence on the phenomena to be studied, not only of a measurement but also of a calculation process in the real device predicting the future.     

锌指蛋白的核酸识别特异性研究进展

锌指蛋白是最大的DNA结合蛋白,它能和DNA进行特异性识别,是研究蛋白—DNA相互作用的理想对象。改变锌指元件上的几个保守的氨基酸位点可设计筛选出序列特异的全新锌指蛋白,计算机在锌指蛋白设计方面的应用,使得全新的锌指蛋白识别特异性明显增强。这在基因治疗等方面,具有广阔的应用前景。     

Modeling H3 histone N-terminal tail and linker DNA interactions

Molecular dynamics computer simulations were performed for the 25-residue N-terminal tail of the H3 histone protein in the proximity of a DNA segment of 10 base pairs (bp), representing a model for the linker DNA in chromatin. Several least biased configurations were used as initial configurations. The secondary structure content of the protein was increased by the presence of DNA close to it, but the locations of the secondary motifs were different for different initial orientations of the DNA grooves with respect to the protein. As a common feature to all simulations, the electrostatic attraction between negatively charged DNA and positively charged protein was screened by the water solvent and counterbalanced by the intrinsic compaction of the protein due to hydrophobic effects. The protein secondary structure limited the covering of DNA by the protein to 4-5 bp. The degree of compaction and charge density of the bound protein suggests a possible role of H3 tail in a nonspecific bending and plasticity of the linker DNA when the protein is located in the crowded dense chromatin.     

Discover 1: a new program to search for unusually represented DNA motifs.

DISCOVER1 (DIStribution COunter VERsion 1) is a new program that can identify DNA motifs occurring with a high deviation from the expected frequency. The program generates families of patterns, each family having a common set of defined bases. Undefined bases are inserted amongst the defined bases in different ways, thus generating the diverse patterns of each family. The occurrences of the different patterns are then compared and analysed within each family, assuming that all patterns should have the same probability of occurrence. An extensive use of computer memory, combined with the immediate sorting of counts by address calculation allow a complete counting of all DNA motifs on a single pass on the DNA sequence. This approach offers a very fast way to search for unusually distributed patterns and can identify inexact patterns as well as exact patterns.     

KorB protein of promiscuous plasmid RP4 recognizes inverted sequence repetitions in regions essential for conjugative plasmid transfer.

We have constructed a RP4 KorB overproducing strain and purified the protein to near homogeneity. KorB is a DNA binding protein recognizing defined palindromic 13-bp sequences (TTTAGCSGCTAAA). Inverted sequence repetitions of this type, designated OB, are present on RP4 12 times. OB-sequences are localized in replication and maintenance regions as well as in the regions Tra1 and Tra2 essential for conjugative transfer. All sites found in Tra regions by computer search act as targets for specific binding of KorB protein. KorB-DNA complexes were detected by DNA fragment retardation assay using polyacrylamide gels. The 13-bp symmetric arrangement of the consensus OB-sequence constitutes the core for binding KorB protein since any truncation of this sequence prevents complex assembly or leads to a considerable destabilization of the KorB-DNA complexes. A hydroxyl radical footprint analysis demonstrated complex formation of KorB with the OB-sequence directly and suggests the presence of an unusual DNA structure within the nucleoprotein complex.     

Long-range interactions between ligands bound to a DNA molecule give rise to adsorption with the character of phase transition of the first kind

Influence of long-range interactions between ligands bound to DNA molecule on the character of their adsorption is studied using computer modeling. For this investigation, two calculation procedures are developed. They are based upon the method of the free energy minimum and on the partition function method. The both procedures demonstrate that in the case of a strong enough attraction between all the bound ligands their binding to DNA has the character of phase transition of the first kind. There is a break in the binding curve c(c0) where c - relative concentration of bound ligands, c0 - molar concentration of free ligands. The break occurs because there is an interval of central degrees of binding (approximately 50% of the maximum c value) that are prohibited for individual DNA molecules. Such a transition might be caused by some types of DNA condensation. Attraction between the neighboring ligands only, adjacent or/and separated by double helix regions, does not cause this effect.     

A program for selecting DNA fragments to detect mutations by denaturing gel electrophoresis methods.

A computer program was developed to automate the selection of DNA fragments for detecting mutations within a long DNA sequence by denaturing gel electrophoresis methods. The program, MELTSCAN, scans through a user specified DNA sequence calculating the melting behavior of overlapping DNA fragments covering the sequence. Melting characteristics of the fragments are analyzed to determine the best fragment for detecting mutations at each base pair position in the sequence. The calculation also determines the optimal fragment for detecting mutations within a user specified mutational hot spot region. The program is built around the statistical mechanical model of the DNA melting transition. The optimal fragment for a given position is selected using the criteria that its melting curve has at least two steps, the base pair position is in the fragment's lowest melting domain, and the melting domain has the smallest number of base pairs among fragments that meet the first two criteria. The program predicted fragments for detecting mutations in the cDNA and genomic DNA of the human p53 gene.     

GENEUS, a computer system for DNA and protein sequence analysis containing an information retrieval system for the EMBL data library.

We describe a comprehensive computer system, GENEUS, for extensive DNA, RNA and protein sequence analysis. The analysis system is developed for the DEC VAX/VMS computer and uses the EMBL Nucleic Acid Sequence Data Library. Help information is available on-line on terminal screen. To speed up system handling, a qualifier oriented user communication is employed. All results are stored on files making them accessible to the computer editor. An information retrieval system for the EMBL Nucleotide Sequence Data Library is also described. A defined data-base interface allows connection to other analysis programs.+     

人线粒体DNA的信息结构

本文报道了运用FORTRAN-77语言,在SIRIUS-1微机上计算遗传信息的冗余结构D_1、D_2、D_3的程序。计算出人线粒体DNA(16569个核苷酸残基)的H_1=1.930554,H_2=3.849254,H_3=5.760944,D_1=0.069446,D_2=0.011853,D_3=0.007011。 D_1、D_2的结果表明,人线粒体DNA的信息结构远比脊椎动物DNA的低级,这支持线粒体的共生起源学说。并对D_3的结果进行了分析,对其意义作了初步探讨。     

DNA melting profiles from a matrix method

In this article we give a new method for the calculation of DNA melting profiles. Based on the matrix formulation of the DNA partition function, the method relies for its efficiency on the fact that the required matrices are very sparse, essentially reducing matrix multiplication to vector multiplication and thus making the computer time required to treat a DNA molecule containing N base pairs proportional to N(2). A key ingredient in the method is the result that multiplication by the inverse matrix can also be reduced to vector multiplication. The task of calculating the melting profile for the entire genome is further reduced by treating regions of the molecule between helix-plateaus, thus breaking the molecule up into independent parts that can each be treated individually. The method is easily modified to incorporate changes in the assignment of statistical weights to the different structural features of DNA. We illustrate the method using the genome of Haemophilus influenzae.     

Calculation of binding isotherms for heterogenous polymers

The matrix method of statistical mechanics is used to calculate equilibria for the binding of small molecules to polymers. When there is only one kind of binding site the problem is simple; some examples are given for illustrative purposes. If, however, the binding sites are not all equivalent and the bound molecules interact or interfere with each other, the problem is no longer trivial, being formally analogous with calculation of the helix–coil transition equilibrium in a heterogeneous polypeptide. Particular difficulties arise when the sequence of binding sites is aperiodic; most naturally occurring materials fall in this class. The purpose of this paper is to point out that problems of this type are readily solved with good accuracy by use of random-number methods on a high-speed digital computer. One such calculation is presented for illustration. The methods developed are applicable to such systems as the binding of actinomycin, Hg^–, and acridine dyes to DNA.