Anharmonicity and self-similarity of the free energy landscape of protein G

The near-native free energy landscape of protein G is investigated through 0.4 microseconds-long atomistic molecular dynamics simulations in explicit solvent. A theoretical and computational framework is used to assess the time-dependence of salient thermodynamical features. While the quasi-harmonic character of the free energy is found to degrade in a few ns, the slow modes display a very mild dependence on the trajectory duration. This property originates from a striking self-similarity of the free energy landscape embodied by the consistency of the principal directions of the local minima, where the system dwells for several ns, and of the virtual jumps connecting them.Comment: revtex, 6 pages, 5 figure

Role of Secondary Motifs in Fast Folding Polymers: A Dynamical Variational Principle

A fascinating and open question challenging biochemistry, physics and even geometry is the presence of highly regular motifs such as alpha-helices in the folded state of biopolymers and proteins. Stimulating explanations ranging from chemical propensity to simple geometrical reasoning have been invoked to rationalize the existence of such secondary structures. We formulate a dynamical variational principle for selection in conformation space based on the requirement that the backbone of the native state of biologically viable polymers be rapidly accessible from the denatured state. The variational principle is shown to result in the emergence of helical order in compact structures.Comment: 4 pages, RevTex, 4 eps figure

Farmland Use Transitions After the CAP Greening: a Preliminary Analysis Using Markov Chains Approach

This paper represents a preliminary attempt to evaluate ex-post impact of the CAP greening payment on farmland use changes, testing by a Markov Chain approach whether farmland use transitions dynamics changed after the introduction of this new policy instrument. Unlike previous contributions, relying on ex-ante simulations, this analysis is based on the actual behaviour of farmers over the period immediately after the last CAP reform. Such ex-post assessment was based on real georeferenced data on farmland allocation, collected in the Lombardy Region, in Northern Italy, over the period 2011-2016. As the current CAP has recently entered in force (in 2015), the present analysis covers the \ufb01rst two years of implementation of the new rules along with the previous four years. Results are in line with previous ex-ante simulations in the same region, detecting a deep discontinuity for those farmland uses characterised by monoculture before the introduction of the greening. They show a signi\ufb01cant discontinuity of farmland use transitions in the reference area after the introduction of greening rules, pointing to a decrease in maize monoculture, in favour of other cereals and legume crops like soybean and alfalfa. Unlike some critical opinions that see current greening rules as a \u201clow pro\ufb01le\u201d compromise, the present analysis points to a strong e\ufb00ect of such rules on regions with high-intensity agriculture

Optimal shapes of compact strings

Optimal geometrical arrangements, such as the stacking of atoms, are of relevance in diverse disciplines. A classic problem is the determination of the optimal arrangement of spheres in three dimensions in order to achieve the highest packing fraction; only recently has it been proved that the answer for infinite systems is a face-centred-cubic lattice. This simply stated problem has had a profound impact in many areas, ranging from the crystallization and melting of atomic systems, to optimal packing of objects and subdivision of space. Here we study an analogous problem--that of determining the optimal shapes of closely packed compact strings. This problem is a mathematical idealization of situations commonly encountered in biology, chemistry and physics, involving the optimal structure of folded polymeric chains. We find that, in cases where boundary effects are not dominant, helices with a particular pitch-radius ratio are selected. Interestingly, the same geometry is observed in helices in naturally-occurring proteins.Comment: 8 pages, 3 composite ps figure

What thermodynamic features characterize good and bad folders? Results from a simplified off-lattice protein model

The thermodynamics of the small SH3 protein domain is studied by means of a simplified model where each bead-like amino acid interacts with the others through a contact potential controlled by a 20x20 random matrix. Good folding sequences, characterized by a low native energy, display three main thermodynamical phases, namely a coil-like phase, an unfolded globule and a folded phase (plus other two phases, namely frozen and random coil, populated only at extremes temperatures). Interestingly, the unfolded globule has some regions already structured. Poorly designed sequences, on the other hand, display a wide transition from the random coil to a frozen state. The comparison with the analytic theory of heteropolymers is discussed

Ruthenium-thymine acetate binding modes: Experimental and theoretical studies

Ruthenium complexes have proved to exhibit antineoplastic activity, related to the interaction of the metal ion with DNA. In this context, synthetic and theoretical studies on ruthenium binding modes of thymine acetate (THAc) have been focused to shed light on the structure-activity relationship. This report deals with the reaction between dihydride ruthenium mer-[Ru(H)2(CO)(PPh3)3], 1 and the thymine acetic acid (THAcOH) selected as model for nucleobase derivatives. The reaction in refluxing toluene between 1 and THAcOH excess, by H2 release affords the double coordinating species k1-(O)THAc-, k2-(O,O)THAc-[Ru(CO)(PPh3)2], 2. The X-ray crystal structure confirms a simultaneous monohapto, dihapto- THAc coordination in a reciprocal facial disposition. Stepwise additions of THAcOH allowed to intercept the monohapto mer-k1(O)THAc-Ru(CO)H(PPh3)3] 3 and dihapto trans(P,P)-k2(O,O)THAc-[Ru(CO)H(PPh3)2] 4 species. Nuclear magnetic resonance (NMR) studies, associated with DFT (Density Function Theory)-calculations energies and analogous reactions with acetic acid, supported the proposed reaction path. As evidenced by the crystal supramolecular hydrogen-binding packing and 1H NMR spectra, metal coordination seems to play a pivotal role in stabilizing the minor [(N=C(OH)] lactim tautomers, which may promote mismatching to DNA nucleobase pairs as a clue for its anticancer activity

Exact Solution of the Munoz-Eaton Model for Protein Folding

A transfer-matrix formalism is introduced to evaluate exactly the partition function of the Munoz-Eaton model, relating the folding kinetics of proteins of known structure to their thermodynamics and topology. This technique can be used for a generic protein, for any choice of the energy and entropy parameters, and in principle allows the model to be used as a first tool to characterize the dynamics of a protein of known native state and equilibrium population. Applications to a $\beta$-hairpin and to protein CI-2, with comparisons to previous results, are also shown.Comment: 4 pages, 5 figures, RevTeX 4. To be published in Phys. Rev. Let

Glassy transition in a disordered model for the RNA secondary structure

We numerically study a disordered model for the RNA secondary structure and we find that it undergoes a phase transition, with a breaking of the replica symmetry in the low temperature region (like in spin glasses). Our results are based on the exact evaluation of the partition function.Comment: 4 pages, 3 figure

Influence of conformational fluctuations on enzymatic activity: modelling the functional motion of beta-secretase

Considerable insight into the functional activity of proteins and enzymes can be obtained by studying the low-energy conformational distortions that the biopolymer can sustain. We carry out the characterization of these large scale structural changes for a protein of considerable pharmaceutical interest, the human $\beta$-secretase. Starting from the crystallographic structure of the protein, we use the recently introduced beta-Gaussian model to identify, with negligible computational expenditure, the most significant distortion occurring in thermal equilibrium and the associated time scales. The application of this strategy allows to gain considerable insight into the putative functional movements and, furthermore, helps to identify a handful of key regions in the protein which have an important mechanical influence on the enzymatic activity despite being spatially distant from the active site. The results obtained within the Gaussian model are validated through an extensive comparison against an all-atom Molecular Dynamics simulation.Comment: To be published in a special issue of J. Phys.: Cond. Mat. (Bedlewo Workshop

Conformations of Proteins in Equilibrium

We introduce a simple theoretical approach for an equilibrium study of proteins with known native state structures. We test our approach with results on well-studied globular proteins, Chymotrypsin Inhibitor (2ci2), Barnase and the alpha spectrin SH3 domain and present evidence for a hierarchical onset of order on lowering the temperature with significant organization at the local level even at high temperatures. A further application to the folding process of HIV-1 protease shows that the model can be reliably used to identify key folding sites that are responsible for the development of drug resistance .Comment: 6 pages, 3 eps figure