Search

tRNATyr of Nanoarchaeum equitans has a remarkable feature with an extra guanosine residue at the 5ʹ-terminus. However, the N. equitans tRNATyr mutant without extra guanosine at the 5ʹ-end was tyrosylated by tyrosyl-tRNA synthase (TyrRS). We solved the crystal structure of N. equitans TyrRS at 2.80 A resolution. By comparing the present solved structure with the complex structures TyrRS with tRNATyr of Thermus thermophilus and Methanocaldococcus jannaschii, an arginine substitution mutant of N. equitans TyrRS at Ile200 (I200R), which is the putative closest candidate to the 5ʹ-phosphate of C1 of N. equitans tRNATyr, was prepared. The I200R mutant tyrosylated not only wild-type tRNATyr but also the tRNA without the G-1 residue. Further tyrosylation analysis revealed that the second base of the anticodon (U35), discriminator base (A73), and C1:G72 base pair are strong recognition sites.

The nucleosome is the basic organizational unit of the genome. The folding structure of nucleosomes is closely related to genome functions, and has been reported to be in dynamic interplay with binding of various nuclear proteins to genomic loci. Here, we describe our high-throughput chromosome conformation capture with nucleosome orientation (Hi-CO) technology to derive 3D nucleosome positions with their orientations at every genomic locus in the nucleus. This technology consists of an experimental procedure for nucleosome proximity analysis and a computational procedure for 3D modeling. The experimental procedure is based on an improved method of high-throughput chromosome conformation capture (Hi-C) analysis. Whereas conventional Hi-C allows spatial proximity analysis among genomic loci with 1-10 kbp resolution, our Hi-CO allows proximity analysis among DNA entry or exit points at every nucleosome locus. This analysis is realized by carrying out ligations among the entry/exit points in every nucleosome in a micrococcal-nuclease-fragmented genome, and by quantifying frequencies of ligation products with next-generation sequencing. Our protocol has enabled this analysis by cleanly excluding unwanted non-ligation products that are abundant owing to the frequent genome fragmentation by micrococcal nuclease. The computational procedure is based on simulated annealing-molecular dynamics, which allows determination of optimized 3D positions and orientations of every nucleosome that satisfies the proximity ligation data sufficiently well. Typically, examination of the Saccharomyces cerevisiae genome with 130 million sequencing reads facilitates analysis of a total of 66,360 nucleosome loci with 6.8 nm resolution. The technique requires 2-3 weeks for sequencing library preparation and 2 weeks for simulation.

Vasohibin-1 and small vasohibin-binding protein (SVBP) form an intermolecular heterotetramer in the crystal. The heterotetramer was stabilized by exchange of the conserved N-terminal region., Vasohibins regulate angiogenesis, tumor growth, metastasis and neuronal differentiation. They form a complex with small vasohibin-binding protein (SVBP) and show tubulin tyrosine carboxypeptidase activity. Recent crystal structure determinations of vasohibin–SVBP complexes have provided a molecular basis for complex formation, substrate binding and catalytic activity. However, the regulatory mechanism and dynamics of the complex remain elusive. Here, the crystal structure of the VASH1–SVBP complex and a molecular-dynamics simulation study are reported. The overall structure of the complex was similar to previously reported structures. Importantly, however, the structure revealed a domain-swapped heterotetramer that was formed between twofold symmetry-related molecules. This heterotetramerization was stabilized by the mutual exchange of ten conserved N-terminal residues from the VASH1 structural core, which was intramolecular in other structures. Interestingly, a comparison of this region with previously reported structures revealed that the patterns of hydrogen bonding and hydrophobic interactions vary. In the molecular-dynamics simulations, differences were found between the heterotetramer and heterodimer, where the fluctuation of the N-terminal region in the heterotetramer was suppressed. Thus, heterotetramer formation and flexibility of the N-terminal region may be important for enzyme activity and regulation.

Nanoarchaeum equitans is a species of hyperthermophilic archaea with the smallest genome size. Its alanyl-tRNA synthetase genes are split into AlaRS-α and AlaRS-β, encoding the respective subunits. In the current report, we surveyed N. equitans AlaRS-dependent alanylation of RNA minihelices, composed only of the acceptor stem and the T-arm of tRNAAla. Combination of AlaRS-α and AlaRS-β showed a strong alanylation activity specific to a single G3:U70 base pair, known to mark a specific tRNA for charging with alanine. However, AlaRS-α alone had a weak but appreciable alanylation activity that was independent of the G3:U70 base pair. The shorter 16-mer RNA tetraloop substrate mimicking only the first four base pairs of the acceptor stem of tRNAAla, but with C3:G70 base pair, was also successfully aminoacylated by AlaRS-α. The end of the acceptor stem, including CCA-3' terminus and the discriminator A73, was able to function as a minimal structure for the recognition by the enzyme. Our findings imply that aminoacylation by N. equitans AlaRS-α may represent a vestige of a primitive aminoacylation system, before the appearance of the G3:U70 pair as an identity element for alanine.

This study reports the X-ray crystallographic structure of the glycyl-tRNA synthetase (GlyRS) of Nanoarchaeum equitans — a hyperthermophilic archaeal species. This is the first archaeal GlyRS crystal structure elucidated. The GlyRS comprises an N-terminal catalytic domain and a C-terminal anticodon-binding domain with a long β-sheet inserted between these domains. An unmodified transcript of the wild-type N. equitans tRNAGly was successfully glycylated using GlyRS. Substitution of the discriminator base A73 of tRNAGly with any other nucleotide caused a significant decrease in glycylation activity. Mutational analysis of the second base-pair C2G71 of the acceptor stem of tRNAGly elucidated the importance of the base-pair, especially G71, as an identity element for recognition by GlyRS. Glycylation assays using tRNAGly G71 substitution mutants and a GlyRS mutant where Arg223 is mutated to alanine strengthen the possibility that the carbonyl oxygen at position 6 of G71 would hydrogen-bond with the guanidine nitrogen of Arg223 in N. equitans GlyRS.

The origin of homochirality in L-amino acid in proteins is one of the mysteries of the evolution of life. Experimental studies show that a non-enzymatic aminoacylation reaction of an RNA minihelix has a preference for L-amino acid over D-amino acid. The reaction initiates by approaching of a 3′-oxygen of the RNA minihelix to the carbonyl carbon of an aminoacyl phosphate oligonucleotide. Here, employing molecular dynamics simulations, we examined the possible mechanisms that determine this chiral selectivity. The simulation system adopted a geometry required for the chemical reaction to occur more frequently with L-alanine than that with D-alanine. For L-alanine, the structure with this geometry was formed by a combination of stable dihedral angles along alanyl phosphate backbone with a canonical RNA structure, where the methyl group of alanine was placed on the opposite side of the approaching 3′-hydroxyl group with respect to the carbonyl plane. For D-alanine, the methyl group and the 3′-hydroxyl group were placed on the same side with respect to the carbonyl plane, which significantly decreased its ability to approach 3′-oxygen close to the carbonyl carbon compared to L-alanine. The mechanism suggested herein can explain experimentally observed chiral preferences.

The cytoplasm of a cell is crowded with many different kinds of macromolecules. The macromolecular crowding affects the thermodynamics and kinetics of biological reactions in a living cell, such as protein folding, association, and diffusion. Theoretical and simulation studies using simplified models focus on the essential features of the crowding effects and provide a basis for analyzing experimental data. In most of the previous studies on the crowding effects, a uniform crowder size is assumed, which is in contrast to the inhomogeneous size distribution of macromolecules in a living cell. Here, we evaluate the free energy changes upon macromolecular association in a cell-like inhomogeneous crowding system via a theory of hard-sphere fluids and free energy calculations using Brownian dynamics trajectories. The inhomogeneous crowding model based on 41 different types of macromolecules represented by spheres with different radii mimics the physiological concentrations of macromolecules in the cytoplasm of Mycoplasma genitalium. The free energy changes of macromolecular association evaluated by the theory and simulations were in good agreement with each other. The crowder size distribution affects both specific and non-specific molecular associations, suggesting that not only the volume fraction but also the size distribution of macromolecules are important factors for evaluating in vivo crowding effects. This study relates in vitro experiments on macromolecular crowding to in vivo crowding effects by using the theory of hard-sphere fluids with crowder-size heterogeneity.

Hydrodynamic interactions (HI) give rise to collective motions between molecules, which are known to be important in the dynamics of random coil polymers and colloids. However, their role in the biological self-assembly of many molecule systems has not been investigated. Here, using Brownian dynamics simulations, we evaluate the importance of HI on the kinetics of self-assembly of lipid membranes. One-thousand coarse-grained lipid molecules in periodic simulation boxes were allowed to assemble into stable bilayers in the presence and absence of intermolecular HI. Hydrodynamic interactions reduce the monomer-monomer association rate by 50%. In contrast, the rate of association of lipid clusters is much faster in the presence of intermolecular HI. In fact, with intermolecular HI, the membrane self-assembly rate is 3–10 times faster than that without intermolecular HI. We introduce an analytical model to describe the size dependence of the diffusive encounter rate of particle clusters, which can qualitatively explain our simulation results for the early stage of the membrane self-assembly process. These results clearly suggest that HI greatly affects the kinetics of self-assembly and that simulations without HI will significantly underestimate the kinetic parameters of such processes.

To evaluate the enthalpic and entropic contributions of polar residues in alanine-based peptides to α-helical propensities of peptides, we applied a Brownian dynamics simulation, together with the umbrella sampling, to two alanine-based 21-residue peptides: one composed of alanine only (AAA), the other composed of 18 alanines and 3 arginines (ARA). Higher α-helical propensity of ARA than that of AAA was obtained. However, they showed similar conformational stability in enthalpy, considering the contribution of the solvation energy and the potential energy. As an evaluation of entropic effects, the fluctuation of a dihedral angle, ψ, was investigated. Arginine residue showed smaller fluctuation than alanine in elongated states. Higher α-helical propensity might originate from entropic effects. Furthermore, arginine seems to affect the α-helical propensity of alanines interacting with arginine.

To begin to elucidate the principles of intermolecular dynamics in the crowded environment of cells, employing Brownian dynamics (BD) simulations, we examined possible mechanism(s) responsible for the great reduction in diffusion constants of macromolecules in vivo from that at infinite dilution. In an Escherichia coli cytoplasm model comprised of 15 different macromolecule types at physiological concentrations, BD simulations of molecular-shaped and equivalent sphere representations were performed with a soft repulsive potential. At cellular concentrations, the calculated diffusion constant of GFP is much larger than experiment, with no significant shape dependence. Next, using the equivalent sphere system, hydrodynamic interactions (HI) were considered. Without adjustable parameters, the in vivo experimental GFP diffusion constant was reproduced. Finally, the effects of nonspecific attractive interactions were examined. The reduction in diffusivity is very sensitive to macromolecular radius with the motion of the largest macromolecules dramatically slowed down; this is not seen if HI dominate. In addition, long-lived clusters involving the largest macromolecules form if attractions dominate, whereas HI give rise to significant, size independent intermolecular dynamic correlations. These qualitative differences provide a testable means of differentiating the importance of HI vs. nonspecific attractive interactions on macromolecular motion in cells.

Chromatin dynamics and arrangement are involved in many biological processes in nuclei of eukaryotes including plants. Plants have to respond rapidly to various environmental stimuli to achieve growth and development because they cannot move. It is assumed that the alteration of chromatin dynamics and arrangement support the response to these stimuli; however, there is little information in plants. In this study, we investigated the chromatin dynamics and arrangement with DNA damage in Arabidopsis thaliana by live-cell imaging with the lacO/LacI-EGFP system and simulation analysis. It was revealed that homologous loci kept a constant distance in nuclei of A. thaliana roots in general growth. We also found that DNA double-strand breaks (DSBs) induce the approach of the homologous loci with γ-irradiation. Furthermore, AtRAD54, which performs an important role in the homologous recombination repair pathway, was involved in the pairing of homologous loci with γ-irradiation. These results suggest that homologous loci approach each other to repair DSBs and AtRAD54 mediates these phenomena.

DNA binding proteins efficiently search for their cognitive sites on long genomic DNA by combining 3D diffusion and 1D diffusion (sliding) along the DNA. Recent experimental results and theoretical analyses revealed that the proteins show a rotation-coupled sliding along DNA helical pitch. Here, we performed Brownian dynamics simulations using newly developed coarse-grained protein and DNA models for evaluating how hydrodynamic interactions between the protein and DNA molecules, binding affinity of the protein to DNA, and DNA fluctuations affect the one dimensional diffusion of the protein on the DNA. Our results indicate that intermolecular hydrodynamic interactions reduce 1D diffusivity by 30%. On the other hand, structural fluctuations of DNA give rise to steric collisions between the CG-proteins and DNA, resulting in faster 1D sliding of the protein. Proteins with low binding affinities consistent with experimental estimates of non-specific DNA binding show hopping along the CG-DNA. This hopping significantly increases sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome., Author Summary DNA binding proteins efficiently search for their cognitive sites on long genomic DNA in cells to control biological activities. Recent experimental studies have revealed that the proteins use not only three-dimensional diffusion, but also one-dimensional diffusion (sliding) on DNA for this search process. For a better understanding of this biological process, we need to elucidate the mechanism of sliding. We report here molecular simulations using newly developed coarse-grained protein and DNA models for elucidating the nature of the sliding motions. Our simulation results show that: 1) hydrodynamic interactions between protein and DNA reduce sliding rate by 30%, 2) structural fluctuations of DNA give rise to steric collisions between proteins and DNA, which facilitate sliding motions, and 3) proteins with low binding affinities to DNA can hop along the DNA, resulting in a significant increase in sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome.

We applied an atomistic Brownian dynamics (BD) simulation with multiple time step method for the folding simulation of a 13-mer α-helical peptide and a 12-mer β-hairpin peptide, giving successful folding simulations. In this model, the driving energy contribution towards folding came from both electrostatic and van der Waals interactions for the α-helical peptide and from van der Waals interactions for the β-hairpin peptide. Although, many non-native structures having the same or lower energy than that of native structure were observed, the folded states formed the most populated cluster when the structures obtained by the BD simulations were subjected to the cluster analysis based on distance-based root mean square deviation of side-chains between different structures. This result indicates that we can predict the native structures from conformations sampled by BD simulation.

A new simple implicit solvent model, effective charge (EC) model, was introduced into the Brownian dynamics algorithm based on AMBER united-atom force field. In the EC model, an atomic charge was decreased as a function of solvent-accessible surface area of the atom. We carried out the Brownian dynamics simulations of a 28-mer ββα fold peptide using four implicit solvent models: a generalized Born/solvent-accessible surface area (GB/SA) model, a solvent-accessible surface area (SA) based solvent model, a SA in combination with distance-dependent dielectric (DD/SA) and the EC combined with DD/SA (DD/SA/EC) model; and the calculated results on structure and dynamics of the peptide were compared with those of molecular dynamics simulation using explicit solvent model. Several artifacts were observed in the simulation using the GB/SA model. On the other hand, simulation using the DD/SA and DD/SA/EC implicit solvent models were free from such artifacts. Especially BD with the DD/SA/EC model gave the most stable trajectory as judged by root mean square deviations from the initial structure without large computational cost.

We report a multiple time step algorithm applied to an atomistic Brownian dynamics simulation for simulating the long time scale dynamics of biomolecules. The algorithm was based on the original multiple time step method; a short time step was used to keep faster motions in local equilibrium. When applied to a 28-mer # # ! folded peptide, the simulation gave stable trajectories and the computation time was reduced by a factor of 160 compared to a conventional molecular dynamics simulation using explicit water molecules. We applied it for the folding simulation of a 13-mer ! -helical peptide, giving a successful folding simulation. These results indicate that the Brownian dynamics with the multiple time step algorithm is useful for studies of biomolecular motions by long time simulation.

We report the implementation of an atomistic Brownian dynamics simulation of proteins. A protein was described by united-atom model with AMBER91 force field. The solvent was treated by distance-dependent dielectric/surface area model. The computation time of the Brownian dynamics and the calculated results on structure and dynamics of 28-mer ββα fold peptide with the implicit solvent model were compared with those of molecular dynamics simulation with explicit solvent. The Brownian dynamics simulation was 53 times faster than molecular dynamics simulation with explicit solvent. The simulation was stable and the artifacts often observed in simulations at vacuum condition were reduced. The results indicate that the Brownian dynamics with the implicit solvent model can be widely used in studies of protein dynamics by long time simulation in future.

The interiors of all living cells are highly crowded with macromolecules, which results in a considerable difference between the thermodynamics and kinetics of biological reactions in vivo from that in vitro. To begin to elucidate the principles of intermolecular dynamics in the crowded environment of cells, employing Brownian dynamics (BD) simulations, we examined possible mechanism(s) responsible for the great reduction in diffusion constants of macromolecules in vivo from that at infinite dilution. In an E. coli cytoplasm modelcomprised of 15 different macromolecule types at physiological concentrations, where macromolecules were represented by spheres with their Stokes radii, BD simulations were performed with and without hydrodynamic interactions (HI). Without HI, the calculated diffusion constant of green fluorescent protein (GFP) is much larger than experiment. On the other hand, when HI were considered, the in vivo experimental GFP diffusion constant is almost reproduced without adjustable parameters. In addition, HI give rise to significant, size independent intermolecular dynamic correlations. These results suggest that HI play an important role on macromolecular dynamics in vivo.

A new tripodal ligand with an asymmetric carbon atom, Ph2PCH2CH(Me)N(CH2CH2PPh2)2 (S-mnp3) was prepared using L-alanine as an optically active starting material; a series of its complexes, [M(X)(S-mnp3)]BPh4(M = Ni(II) and Co(II); X = NCS−, Cl−, Br− and I−) was prepared. Electronic spectra of these complexes indicated that their coordination geometries are essentially the same as those of corresponding complexes containing the tripodal ligand tris (diphenylphosphinoethyl) amine. Thus, all the present nickel complexes in solution exhibit electronic spectra typical of trigonal-bipyramidal complexes of diamagnetic configuration. CD spectra of these nickel complexes have been discussed on the basis of the selection rules for magnetic dipole transitions under C 3v symmetry. Cobalt(II) complexes are either high-spin or low-spin, depending on the nature of the unidentate ligand. The high-spin cobalt(II) complexes did not show CD spectra typical of trigonal-bipyramidal geometry.

Asp187 and Gln190 were predicted as conserved and closely located at the Na(+) binding site in a topology and homology model structure of Na(+)/proline symporter (PutP) of Escherichia coli. The replacement of Asp187 with Ala or Leu did not affect proline transport activity; whereas, change to Gln abolished the active transport. The binding affinity for Na(+) or proline of these mutants was similar to that of wild-type (WT) PutP. This result indicates Asp187 to be responsible for active transport of proline without affecting the binding. Replacement of Gln190 with Ala, Asn, Asp, Leu and Glu had no effect on transport or binding, suggesting that it may not have a role in the transport. However, in the negative D187Q mutant, a second mutation, of Gln190 to Glu or Leu, restored 46 or 7% of the transport activity of WT, respectively, while mutation to Ala, Asn or Asp had no effect. Thus, side chain at position 190 has a crucial role in suppressing the functional defect of the D187Q mutant. We conclude that Asp187 is responsible for transport activity instead of coupling-ion binding by constituting the translocation pathway of the ion and Gln190 provides a suppressing mutation site to regain PutP functional activity.

The interiors of all living cells are highly crowded with macromolecules, which differs considerably the thermodynamics and kinetics of biological reactions between in vivo and in vitro. For example, the diffusion of green fluorescent protein (GFP) in E. coli is ~10-fold slower than in dilute conditions. In this study, we performed Brownian dynamics (BD) simulations of rigid macromolecules in a crowded environment mimicking the cytosol of E. coli to study the motions of macromolecules. The simulation systems contained 35 70S ribosomes, 750 glycolytic enzymes, 75 GFPs, and 392 tRNAs in a 100 nm × 100 nm × 100 nm simulation box, where the macromolecules were represented by rigid-objects of one bead per amino acid or four beads per nucleotide models. Diffusion tensors of these molecules in dilute solutions were estimated by using a hydrodynamic theory to take into account the diffusion anisotropy of arbitrary shaped objects in the BD simulations. BD simulations of the system where each macromolecule is represented by its Stokes radius were also performed for comparison. Excluded volume effects greatly reduce the mobility of molecules in crowded environments for both molecular-shaped and equivalent sphere systems. Additionally, there were no significant differences in the reduction of diffusivity over the entire range of molecular size between two systems. However, the reduction in diffusion of GFP in these systems was still 4–5 times larger than for the in vivo experiment. We will discuss other plausible factors that might cause the large reduction in diffusion in vivo.

Eight new cobalt(III) complexes of the type, [CoX3(mmtp)]n+ (mmtp = 1,1,1-tris(dimethylphosphinomethyl)ethane, X3 = (Cl)3, (Br)3, (CN)3, (NH3)3, mmtp, 3-azapentane-1,5-diamine (dien), 1,4,7-triazacyclononane (tacn), and 1,4,7-trithiacyclononane (ttcn)), were prepared and characterized. All of the complexes are stable in the solid state, but the halogeno complexes aquate rapidly in water, and the ammonia complex liberates NH3 slowly in 1 mol dm−3 HCl to give [CoCl2(NH3)(mmtp)]+, other complexes including [Co(mmtp)2]3+ being stable in solution. The crystal structure of [Co(mmtp)2][Co(CN)6]·2.25H2O was determined by X-ray analysis. Crystal data: orthorhombic, space group Pca21, a = 31.579(5), b = 10.521(1), c = 23.458(3) A, V = 7794(2) A3, Z= 8, R = 0.084 for 4196 reflections. The Co–P bond distances are 2.291(9)–2.357(7) A (av 2.331 A), which are appreciably longer than those usually found in Co(III)–phosphine complexes. The absorption spectra of the complexes were compared with one another and with those o...

鶏白血病ウイルス(ALV)抗体測定法としての酵素抗体法(ELISA)について抗原抗体反応系の解析を行い,さらに血清中和法(SN),蛍光抗体法(FA)との相関関係と感度の比較を行った。1)A亜群,B亜群2種の抗原に対する抗体価の推移をELISA,SN,FAで確認すると,ELISA,FAでは主として群特異抗体(gs抗体)を,SNは亜群特異抗体を検出していると思われた。2)gs抗体を産生する抗原成分として3種の蛋白が確認された。そして,それぞれに対する抗体の発現順序は,まずp27,つぎにgp85,さらにp15に対する抗体であった。3)1野外鶏群について,3種の測定方法で算出した抗体価の関係を見ると,いずれにも相関が認められた。4)測定方法による検出感度を比較すると,ELISAは,SNの約10倍,FAはSNとほぼ同感度であった。

Self-organization is one of the most fundamental phenomena in biology. Elucidating the mechanisms of the self-assembly of biological molecules, such as protein folding/binding and formation of lipid bilayer membranes, have been a fascinating topic over decades in biological field. Hydrodynamic interactions (HIs) give rise to collective motions between molecules. Therefore, they may play an important role on dynamics of self-assembly. Although there are significant advances in theoretical and computational studies on thermodynamics and kinetics of self-assembling systems recent years, analysis on effects of HIs to the reactions are limited due to their high complexity. Here, by using a Brownian dynamics (BD) simulation technique, we evaluate the importance of HIs on kinetic of self-assembly in a set of examples, formation of lipid membrane and actin polymerization reaction. We first built coarse-grained models of lipid and actin monomer. A cluster of lipid molecules was represented by two particles, one is hydrophilic head group and the other is hydrophobic tail group, and attractive interaction between tail groups was applied. For actin polymerization system, each monomer was represented by two particles and an angle-dependent attractive interaction was applied to form long filaments. In BD simulations, the Rotne-Prager-Yamakawa tensor was used to account for HIs between particles. In the early stage of self-assembly, HIs between inter-monomer particles decelerated the association rate of monomers to form oligomers compared with that in the simulations without inter-monomer HIs in both systems. On the other hand, the HIs accelerate the reaction of oligomers to form a fully-assembled state or much longer polymers. This phenomenon could be explained by the Kirkwood-Riseman theory. Those results clearly suggest that HIs greatly affect kinetics of self-assembly reactions and considering the interactions is crucial for studying dynamics of biological systems.

The presence of chick helper factor (chf) in chick embryo fibroblasts (CEF) which lack the natural group-specific (gs) antigen of avian leukosis was determined in a simple Mendelian manner from the dominant autosomal gene. The presence or absence of chf determined the susceptibility to subgroup E virus in the WL-O line, but not in the WL-B line, which had the fixed C/BE character. Dimethyl sulfoxide (DMSO) together with BH-RSV(RAV-7) promoted the appearance of chf in gs − chf − CEF after long-term culture.