An Extended Model for the Evolution of Prebiotic Homochirality: A Bottom-Up Approach to the Origin of Life

A generalized autocatalytic model for chiral polymerization is investigated in detail. Apart from enantiomeric cross-inhibition, the model allows for the autogenic (non-catalytic) formation of left and right-handed monomers from a substrate with reaction rates $\epsilon_L$ and $\epsilon_R$, respectively. The spatiotemporal evolution of the net chiral asymmetry is studied for models with several values of the maximum polymer length, N. For N=2, we study the validity of the adiabatic approximation often cited in the literature. We show that the approximation obtains the correct equilibrium values of the net chirality, but fails to reproduce the short time behavior. We show also that the autogenic term in the full N=2 model behaves as a control parameter in a chiral symmetry- breaking phase transition leading to full homochirality from racemic initial conditions. We study the dynamics of the N -> infinity model with symmetric ($\epsilon_L = \epsilon_R$) autogenic formation, showing that it only achieves homochirality for $\epsilon < \epsilon_c$, where $\epsilon_c$ is an N-dependent critical value. For $\epsilon \leq \epsilon_c$ we investigate the behavior of models with several values of N, showing that the net chiral asymmetry grows as tanh(N). We show that for a given symmetric autogenic reaction rate, the net chirality and the concentrations of chirally pure polymers increase with the maximum polymer length in the model. We briefly discuss the consequences of our results for the development of homochirality in prebiotic Earth and possible experimental verification of our findings

Chiral Polymerization in Open Systems From Chiral-Selective Reaction Rates

We investigate the possibility that prebiotic homochirality can be achieved exclusively through chiral-selective reaction rate parameters without any other explicit mechanism for chiral bias. Specifically, we examine an open network of polymerization reactions, where the reaction rates can have chiral-selective values. The reactions are neither autocatalytic nor do they contain explicit enantiomeric cross-inhibition terms. We are thus investigating how rare a set of chiral-selective reaction rates needs to be in order to generate a reasonable amount of chiral bias. We quantify our results adopting a statistical approach: varying both the mean value and the rms dispersion of the relevant reaction rates, we show that moderate to high levels of chiral excess can be achieved with fairly small chiral bias, below 10%. Considering the various unknowns related to prebiotic chemical networks in early Earth and the dependence of reaction rates to environmental properties such as temperature and pressure variations, we argue that homochirality could have been achieved from moderate amounts of chiral selectivity in the reaction rates.Comment: 15 pages, 6 figures, accepted for publication in Origins of Life and Evolution of Biosphere

The past, present and future of Scientific discourse.

Published versio

Solid-phase molecular recognition of cytosine based on proton-transfer reaction. Part II. supramolecular architecture in the cocrystals of cytosine and its 5-Fluoroderivative with 5-Nitrouracil

<p>Abstract</p> <p>Background</p> <p>Cytosine is a biologically important compound owing to its natural occurrence as a component of nucleic acids. Cytosine plays a crucial role in DNA/RNA base pairing, through several hydrogen-bonding patterns, and controls the essential features of life as it is involved in genetic codon of 17 amino acids. The molecular recognition among cytosines, and the molecular heterosynthons of molecular salts fabricated through proton-transfer reactions, might be used to investigate the theoretical sites of cytosine-specific DNA-binding proteins and the design for molecular imprint.</p> <p>Results</p> <p>Reaction of cytosine (Cyt) and 5-fluorocytosine (5Fcyt) with 5-nitrouracil (Nit) in aqueous solution yielded two new products, which have been characterized by single-crystal X-ray diffraction. The products include a dihydrated molecular salt (CytNit) having both ionic and neutral hydrogen-bonded species, and a dihydrated cocrystal of neutral species (5FcytNit). In CytNit a protonated and an unprotonated cytosine form a triply hydrogen-bonded aggregate in a self-recognition ion-pair complex, and this dimer is then hydrogen bonded to one neutral and one anionic 5-nitrouracil molecule. In 5FcytNit the two neutral nucleobase derivatives are hydrogen bonded in pairs. In both structures conventional N-H^...O, O-H^...O, N-H^+...N and N-H^...N^-intermolecular interactions are most significant in the structural assembly.</p> <p>Conclusion</p> <p>The supramolecular structure of the molecular adducts formed by cytosine and 5-fluorocytosine with 5-nitrouracil, CytNit and 5FcytNit, respectively, have been investigated in detail. CytNit and 5FcytNit exhibit widely differing hydrogen-bonding patterns, though both possess layered structures. The crystal structures of CytNit (D<it>p</it>k_a= -0.7, molecular salt) and 5FcytNit (D<it>p</it>k_a= -2.0, cocrystal) confirm that, at the present level of knowledge about the nature of proton-transfer process, there is not a strict correlation between the D<it>p</it>k_avalues and the proton transfer, in that the acid/base <it>p</it>k_astrength is not a definite guide to predict the location of H atoms in the solid state. Eventually, the absence in 5FcytNit of hydrogen bonds involving fluorine is in agreement with findings that covalently bound fluorine hardly ever acts as acceptor for available Brønsted acidic sites in the presence of competing heteroatom acceptors.</p

Toward homochiral protocells in noncatalytic peptide systems

The activation-polymerization-epimerization-depolymerization (APED) model of Plasson et al. has recently been proposed as a mechanism for the evolution of homochirality on prebiotic Earth. The dynamics of the APED model in two-dimensional spatially-extended systems is investigated for various realistic reaction parameters. It is found that the APED system allows for the formation of isolated homochiral proto-domains surrounded by a racemate. A diffusive slowdown of the APED network such as induced through tidal motion or evaporating pools and lagoons leads to the stabilization of homochiral bounded structures as expected in the first self-assembled protocells.Comment: 10 pages, 5 figure

Hierarchically Structured Lipid Systems (Part of 2nd Edition)

No abstract available

C-Terminal Substitution of MDM2 Interacting Peptides Modulates Binding Affinity by Distinctive Mechanisms

The complex between the proteins MDM2 and p53 is a promising drug target for cancer therapy. The residues 19–26 of p53 have been biochemically and structurally demonstrated to be a most critical region to maintain the association of MDM2 and p53. Variation of the amino acid sequence in this range obviously alters the binding affinity. Surprisingly, suitable substitutions contiguous to this region of the p53 peptides can yield tightly binding peptides. The peptide variants may differ by a single residue that vary little in their structural conformations and yet are characterized by large differences in their binding affinities. In this study a systematic analysis into the role of single C-terminal mutations of a 12 residue fragment of the p53 transactivation domain (TD) and an equivalent phage optimized peptide (12/1) were undertaken to elucidate their mechanistic and thermodynamic differences in interacting with the N-terminal of MDM2. The experimental results together with atomistically detailed dynamics simulations provide insight into the principles that govern peptide design protocols with regard to protein-protein interactions and peptidomimetic design

Synthetic Biology of Proteins: Tuning GFPs Folding and Stability with Fluoroproline

Proline residues affect protein folding and stability via cis/trans isomerization of peptide bonds and by the C(gamma)-exo or -endo puckering of their pyrrolidine rings. Peptide bond conformation as well as puckering propensity can be manipulated by proper choice of ring substituents, e.g. C(gamma)-fluorination. Synthetic chemistry has routinely exploited ring-substituted proline analogs in order to change, modulate or control folding and stability of peptides.In order to transmit this synthetic strategy to complex proteins, the ten proline residues of enhanced green fluorescent protein (EGFP) were globally replaced by (4R)- and (4S)-fluoroprolines (FPro). By this approach, we expected to affect the cis/trans peptidyl-proline bond isomerization and pyrrolidine ring puckering, which are responsible for the slow folding of this protein. Expression of both protein variants occurred at levels comparable to the parent protein, but the (4R)-FPro-EGFP resulted in irreversibly unfolded inclusion bodies, whereas the (4S)-FPro-EGFP led to a soluble fluorescent protein. Upon thermal denaturation, refolding of this variant occurs at significantly higher rates than the parent EGFP. Comparative inspection of the X-ray structures of EGFP and (4S)-FPro-EGFP allowed to correlate the significantly improved refolding with the C(gamma)-endo puckering of the pyrrolidine rings, which is favored by 4S-fluorination, and to lesser extents with the cis/trans isomerization of the prolines.We discovered that the folding rates and stability of GFP are affected to a lesser extent by cis/trans isomerization of the proline bonds than by the puckering of pyrrolidine rings. In the C(gamma)-endo conformation the fluorine atoms are positioned in the structural context of the GFP such that a network of favorable local interactions is established. From these results the combined use of synthetic amino acids along with detailed structural knowledge and existing protein engineering methods can be envisioned as a promising strategy for the design of complex tailor-made proteins and even cellular structures of superior properties compared to the native forms

Stability mechanisms of a thermophilic laccase probed by molecular dynamics.

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T opt (∼350 K) from 300 K, this change correlated with a beginning "unzipping" of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F(-) (but not Cl(-)) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F(-) intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes

Peptide Bond Distortions from Planarity: New Insights from Quantum Mechanical Calculations and Peptide/Protein Crystal Structures

By combining quantum-mechanical analysis and statistical survey of peptide/protein structure databases we here report a thorough investigation of the conformational dependence of the geometry of peptide bond, the basic element of protein structures. Different peptide model systems have been studied by an integrated quantum mechanical approach, employing DFT, MP2 and CCSD(T) calculations, both in aqueous solution and in the gas phase. Also in absence of inter-residue interactions, small distortions from the planarity are more a rule than an exception, and they are mainly determined by the backbone ψ dihedral angle. These indications are fully corroborated by a statistical survey of accurate protein/peptide structures. Orbital analysis shows that orbital interactions between the σ system of Cα substituents and the π system of the amide bond are crucial for the modulation of peptide bond distortions. Our study thus indicates that, although long-range inter-molecular interactions can obviously affect the peptide planarity, their influence is statistically averaged. Therefore, the variability of peptide bond geometry in proteins is remarkably reproduced by extremely simplified systems since local factors are the main driving force of these observed trends. The implications of the present findings for protein structure determination, validation and prediction are also discussed