ATP dependent NS3 helicase interaction with RNA: insights from molecular simulations.

Non-structural protein 3 (NS3) helicase from hepatitis C virus is an enzyme that unwinds and translocates along nucleic acids with an ATP-dependent mechanism and has a key role in the replication of the viral RNA. An inchworm-like mechanism for translocation has been proposed based on crystal structures and single molecule experiments. We here perform atomistic molecular dynamics in explicit solvent on the microsecond time scale of the available experimental structures. We also construct and simulate putative intermediates for the translocation process, and we perform non-equilibrium targeted simulations to estimate their relative stability. For each of the simulated structures we carefully characterize the available conformational space, the ligand binding pocket, and the RNA binding cleft. The analysis of the hydrogen bond network and of the non-equilibrium trajectories indicates an ATP-dependent stabilization of one of the protein conformers. Additionally, enthalpy calculations suggest that entropic effects might be crucial for the stabilization of the experimentally observed structures.

Prebiotic chemistry and origins of life research with atomistic computer simulations.

Research in origins of life is an intrinsically multi-disciplinary field, aimed at finding answers to the formidably complex problem of understanding the emergence of life from the modern versions of Charles Darwin's celebrated "primordial soup". In the last few years, thanks to the increasing computational power and the development of sophisticated theoretical and numerical methods, several computational chemistry and physics groups have invested this field, providing new microscopic insights on fundamental prebiotic chemistry phenomena possibly occurring in the early Earth and outer space. This review presents the most successful and powerful approaches in computational chemistry, and the main results thus obtained in prebiotic chemistry and origins of life. The aim of this work is both to describe the state-of-the-art in computational prebiotic chemistry, possibly useful both to theorists and experimentalists in origins of life research, and to suggest future directions and new perspectives offered by modern simulation tools.

Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds.

The organic compounds found in carbonaceous chondrite meteorites provide insight into primordial solar system chemistry. Evaluating the formation and decomposition mechanisms of meteoritic amino acids may aid our understanding of the origins of life and homochirality on Earth. The amino acid glycine is widespread in meteorites and other extraterrestrial environments; other amino acids, such as isovaline, are found with enantiomeric excesses in some meteorites. The relationship between meteoritic amino acids and other compounds with similar molecular structures, such as aliphatic monoamines and monocarboxylic acids is unclear; experimental results evaluating the decomposition of amino acids have produced inconclusive results about the preferred pathways, reaction intermediates, and if the conditions applied may be compatible with those occurring inside meteoritic parent bodies. In this work, we performed extensive tandem metadynamics, umbrella sampling, and committor analysis to simulate the neutral mild hydrothermal decomposition mechanisms of glycine and isovaline and put them into context for the origins of meteoritic organic compounds. Our ab initio simulations aimed to determine free energy profiles and decomposition pathways for glycine and isovaline. We found that under our modeled conditions, methylammonium, glycolic acid, and sec-butylamine are the most likely decomposition products. These results suggest that meteoritic aliphatic monocarboxylic acids are not produced from decomposition of meteoritic amino acids. Our results also indicate that the decomposition of L-isovaline prefers an enantioselective pathway resulting in the production of (S)-sec-butylamine.

Synthesis of RNA Nucleotides in Plausible Prebiotic Conditions from ab Initio Computer Simulations.

Understanding the mechanism of spontaneous formation of ribonucleotides under realistic prebiotic conditions is a key open issue of origins-of-life research. In cells, de novo and salvage nucleotide enzymatic synthesis combines 5-phospho-α-d-ribose-1-diphosphate (α-PRPP) and nucleobases. Interestingly, these reactants are also known as prebiotically plausible compounds. Combining ab initio molecular dynamics simulations with recently developed reaction exploration and enhanced sampling methods, we show that nucleobases and α-PRPP should spontaneously combine, under mild hydrothermal conditions, with an exothermic reaction and a facile mechanism, forming both purine and pyrimidine ribonucleotides. Surprisingly, this mechanism is very similar to the biological one and yields ribonucleotides with the same anomeric carbon chirality as in biological systems. Mass spectrometry experiments performed on solutions of adenine and PRPP in similar conditions support the formation of AMP. These results suggest that natural selection might have optimized, through enzymes, a pre-existing ribonucleotide formation mechanism, carrying it forward to modern life forms.

Accurate multiple time step in biased molecular simulations.

Many recently introduced enhanced sampling techniques are based on biasing coarse descriptors (collective variables) of a molecular system on the fly. Sometimes the calculation of such collective variables is expensive and becomes a bottleneck in molecular dynamics simulations. An algorithm to treat smooth biasing forces within a multiple time step framework is here discussed. The implementation is simple and allows a speed up when expensive collective variables are employed. The gain can be substantial when using massively parallel or GPU-based molecular dynamics software. Moreover, a theoretical framework to assess the sampling accuracy is introduced, which can be used to assess the choice of the integration time step in both single and multiple time step biased simulations.