Search

Biological systems have been shown to have quantum-like behaviors by applying the adaptive dynamics view on their interaction networks. In particular, in the process of lactose–glucose metabolism, cells generate probabilistic interference patterns similarly to photons in the two-slit experiment. Such quantum-like interference patterns can be found in biological data, on all scales, from proteins to cognitive, ecological, and social systems. The adaptive dynamics approach covers both biological and physical phenomena, including the ones which are typically associated with quantum physics. We guess that the adaptive dynamics can be used for the clarification of quantum foundations, and the present paper is the first step in this direction. We suggest the use of an algorithm for the numerical simulation of the behavior of a billiard ball-like particle passing through two slits by explicitly considering the influence of the two-slit environment (experimental context). Our simulation successfully mimics the interference pattern obtained experimentally in quantum physics. The interference of photons or electrons by two slits is known as a typical quantum mechanical effect. We do not claim that the adaptive dynamics can reproduce the whole body of quantum mechanics, but we hope that this numerical simulation example will stimulate further extensive studies in this direction—the representation of quantum physical phenomena in an adaptive dynamical framework.

The class I ligase ribozyme consists of 121 nucleotides and shows a high catalytic rate comparable to that found in natural proteinaceous polymerases. In this study, we aimed to identify the smaller active unit of the class I ligase ribozyme comprising ~50 nucleotides, comparable to the estimated length of prebiotically synthesized RNA. Based on the three-dimensional structure of the class I ligase ribozyme, mutants were prepared and their ligation activities were analyzed. Sufficient ligation activity was maintained even when shortening to 94 nucleotides. However, because it would be difficult to approach the target of ~50 nucleotides by removing only the partial structure, the class I ligase ribozyme was then split into two molecules. The ligation activity was maintained even when splitting into two molecules of 55 and 39 nucleotides. Using a system with similar split ribozymes, we analyzed the ligation activity of mutants C30, C47, and A71, which have been previously identified as the positions that contribute to catalytic activity, and discussed the structural basis of the activity of these bases. Our findings suggest the rationale for the class I ligase ribozyme’s assembling from multiple fragments that would be achievable with prebiotic synthesis.

Aminoacylation of a primordial RNA minihelix composed of D-ribose shows L-amino acid preference over D-amino acid without any ribozymes or enzymes. This preference in the amino acylation reaction likely plays an important role in the establishment of homochirality in L-amino acid in modern proteins. However, molecular mechanisms of the chiral selective reaction remain unsolved mainly because of difficulty in direct observation of the reaction at the molecular scale by experiments. For seeking a possible mechanism of the chiral selectivity, quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations of the aminoacylation reactions in a modeled RNA were performed to investigate differences in their free-energy profiles along the reactions for L- and D-alanine and its physicochemical origin. The reaction is initiated by approaching a 3′-oxygen of the RNA minihelix to the carbonyl carbon of an aminoacyl phosphate oligonucleotide. The QM/MM umbrella sampling MD calculations showed that the height of the free-energy barrier for L-alanine aminoacylation reaction was 17 kcal/mol, which was 9 kcal/mol lower than that for the D-alanine system. At the transition state, the distance between the negatively charged 3′-oxygen and the positively charged amino group of L-alanine was shorter than that of D-alanine, which was caused by the chirality difference of the amino acid. These results indicate that the transition state for L-alanine is more electrostatically stabilized than that for D-alanine, which would be a plausible mechanism previously unexplained for chiral selectivity in the RNA minihelix aminoacylation.

In this study, we performed molecular dynamics (MD) simulations of the filling process of few-nanometer-wide trenches with various resist materials in ultraviolet nanoimprint lithography (UV-NIL) to identify the main molecular features necessary for a successful filling process. The 2- and 3-nm wide trenches were successfully filled with the resist materials that had (experimentally determined) viscosities less than 10 mPa·s. The resist composed of a three-armed bulky and highly viscous molecule could not fill the trenches. The radius of gyration of this molecule was smaller than half of the distance in which the first peak of its radial distribution function occurred. The available shapes of 1,6-hexanediol diacrylate (HDDA) and tri(propylene glycol) diacrylate (TPGDA), which are linear photopolymers, were compared to reveal that TPGDA is more flexible and adopts more conformations than HDDA. The terminal functional groups of TPGDA can be close due to its flexibility, which would increase the probability of intramolecular crosslinking of the molecule. This simulation result could explain the difference in hardness between the UV-cured HDDA and TPGDA based materials observed by experiments. The findings revealed by our MD simulations provide useful information for selecting and designing resists for fine patterning by UV-NIL.

We report a patient with a 6q22.1 deletion, who presented with a rare syndrome of generalized epilepsy, myoclonic tremor, and intellectual disability. There was no clinical progression after follow-up for more than 10 years. Our report presents the genetic basis for a phenotype involving a non-progressive generalized epilepsy with tremor. The efficacy of valproic acid for seizure control and the partial efficacy of deep brain stimulation with propranolol for myoclonic tremor is detailed.

The peptidyl transferase center (PTC) in the ribosome is composed of two symmetrically arranged tRNA-like units that contribute to peptide bond formation. We prepared units of the PTC components with putative tRNA-like structure and attempted to obtain peptide bond formation between aminoacyl-minihelices (primordial tRNAs, the structures composed of a coaxial stack of the acceptor stem on the T-stem of tRNA). One of the components of the PTC, P1c2UGGU (74-mer), formed a dimer and a peptide bond was formed between two aminoacyl-minihelices tethered by the dimeric P1c2UGGU. Peptide synthesis depended on both the existence of the dimeric P1c2UGGU and the sequence complementarity between the ACCA-3′ sequence of the minihelix. Thus, the tRNA-like structures derived from the PTC could have originated as a scaffold of aminoacyl-minihelices for peptide bond formation through an interaction of the CCA sequence of minihelices. Moreover, with the same origin, some would have evolved to constitute the present PTC of the ribosome, and others to function as present tRNAs.

In the present in vitro selection study, we isolated and characterized RNA aptamers for a tRNA-binding protein (Trbp) from an extremophile archaeon Aeropyrum pernix. Trbp-like structures are frequently found not only in aminoacyl-tRNA synthetases but also in diverse types of proteins from different organisms. They likely arose early in evolution and have played important roles in evolution through interactions with key RNA structures. RNA aptamers specific for A. pernix Trbp were successfully selected from a pool of RNAs composed of 60 nucleotides, including a random 30-nucleotide region. From the secondary structures, we obtained a shortened sequence composed of 21 nucleotides, of which the 3′-terminal single stranded CA nucleotides are essential for binding. This may be related to the initial evolutionary role of the universal CCA-3′ terminus of tRNA in the interaction with Trbp-like structures.

Biological macromolecules function in highly crowded cellular environments. The structure and dynamics of proteins and nucleic acids are well characterized in vitro, but in vivo crowding effects remain unclear. Using molecular dynamics simulations of a comprehensive atomistic model cytoplasm we found that protein-protein interactions may destabilize native protein structures, whereas metabolite interactions may induce more compact states due to electrostatic screening. Protein-protein interactions also resulted in significant variations in reduced macromolecular diffusion under crowded conditions, while metabolites exhibited significant two-dimensional surface diffusion and altered protein-ligand binding that may reduce the effective concentration of metabolites and ligands in vivo. Metabolic enzymes showed weak non-specific association in cellular environments attributed to solvation and entropic effects. These effects are expected to have broad implications for the in vivo functioning of biomolecules. This work is a first step towards physically realistic in silico whole-cell models that connect molecular with cellular biology.

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.