Lorenzo Bongini

Atomic force microscopy images suggest aggregation mechanism in cerato-platanin.

Cerato-platanin (CP), the first member of the “cerato-platanin family”, is a moderately hydrophobic protein produced by Ceratocystis fimbriata, the causal agent of a severe plant disease called “canker stain”. The protein is localized in the cell wall of the fungus and it seems to be involved in the host-plane interaction and induces both cell necrosis and phytoalexin synthesis (one of the first plant defence-related events). Recently, it has been determined that CP, like other fungal surface protein, is able to self assemble in vitro. In this paper we characterize the aggregates of CP by Atomic Force Microscopy (AFM) images. We observe that CP tends to form early annular-shaped oligomers that seem to constitute the fundamental bricks of a hierarchical aggregation process, eventually resulting in large macrofibrillar assemblies. A simple model, based on the hypothesis that the aggregation is energetically favourable when the exposed surface is reduced, is compatible with the measured aggregates’ shape and size. The proposed model can help to understand the mechanism by which CP and many other fungal surface proteins exert their effects.

PicoNewton-Millisecond Force Steps Reveal the Transition Kinetics and Mechanism of the Double-Stranded DNA Elongation

We study the kinetics of the overstretching transition in λ-phage double-stranded (ds) DNA from the basic conformation (B state) to the 1.7-times longer and partially unwound conformation (S state), using the dual-laser optical tweezers under force-clamp conditions at 25°C. The unprecedented resolution of our piezo servo-system, which can impose millisecond force steps of 0.5-2 pN, reveals the exponential character of the elongation kinetics and allows us to test the two-state nature of the B-S transition mechanism. By analyzing the load-dependence of the rate constant of the elongation, we find that the elementary elongation step is 5.85 nm, indicating a cooperativity of ~25 basepairs. This mechanism increases the free energy for the elementary reaction to ~94 k(B)T, accounting for the stability of the basic conformation of DNA, and explains why ds-DNA can remain in equilibrium as it overstretches.

A dynamical study of antibody-antigen encounter reactions

The effects of internal dynamics in diffusion-driven encounters between macro-molecules represent a problem of broad relevance in molecular biology. In this view, we investigate a typical antigen-antibody reaction chain, based on a coarse-grained mechanical model parameterized directly upon results from single-molecule experiments. We demonstrate that the internal dynamics is a crucial factor in the encounter process. To describe our numerical results, we formulate a simple, intuitive theoretical framework, and we develop it analytically. This enables us to show that the inner dynamics of antibody molecules results in a cooperative behavior of their individual sub-units. Along the same lines, we also investigate the case of double binding to multi-valent antigens. Our results quantify the enhancement of avidity afforded by the double binding in excellent agreement with the available experimental data.

Transient kinetics measured with force steps discriminate between double-stranded DNA elongation and melting and define the reaction energetics

Under a tension of ∼65 pN, double-stranded DNA undergoes an overstretching transition from its basic (B-form) conformation to a 1.7 times longer conformation whose nature is only recently starting to be understood. Here we provide a structural and thermodynamic characterization of the transition by recording the length transient following force steps imposed on the λ-phage DNA with different melting degrees and temperatures (10–25°C). The shortening transient following a 20–35 pN force drop from the overstretching force shows a sequence of fast shortenings of double-stranded extended (S-form) segments and pauses owing to reannealing of melted segments. The lengthening transients following a 2–35 pN stretch to the overstretching force show the kinetics of a two-state reaction and indicate that the whole 70% extension is a B-S transition that precedes and is independent of melting. The temperature dependence of the lengthening transient shows that the entropic contribution to the B-S transition is one-third of the entropy change of thermal melting, reinforcing the evidence for a double-stranded S-form that maintains a significant fraction of the interstrand bonds. The cooperativity of the unitary elongation (22 bp) is independent of temperature, suggesting that structural factors, such as the nucleic acid sequence, control the transition.

The transition mechanism of DNA overstretching: a microscopic view using molecular dynamics

The overstretching transition in torsionally unconstrained DNA is studied by means of atomistic molecular dynamics simulations. The free-energy profile as a function of the length of the molecule is determined through the umbrella sampling technique providing both a thermodynamic and a structural characterization of the transition pathway. The zero-force free-energy profile is monotonic but, in accordance with recent experimental evidence, becomes two-state at high forces. A number of experimental results are satisfactorily predicted: (i) the entropic and enthalpic contributions to the free-energy difference between the basic (B) state and the extended (S) state; (ii) the longitudinal extension of the transition state and (iii) the enthalpic contribution to the transition barrier. A structural explanation of the experimental finding that overstretching is a cooperative reaction characterized by elementary units of approximately 22 base pairs is found in the average distance between adenine/thymine-rich regions along the molecule. The overstretched DNA adopts a highly dynamical and structurally disordered double-stranded conformation which is characterized by residual base pairing, formation of non-native intra-strand hydrogen bonds and effective hydrophobic screening of apolar regions.

Thermally activated processes in polymer dynamics

Jumps between neighboring minima in the energy landscape of both homopolymeric and heteropolymeric chains are numerically investigated by determining the average escape time from different valleys. The numerical results are compared to the theoretical expression derived by Langer [J.S. Langer, Ann. Phys. (N.Y.) 54, 258 (1969)] with reference to a 2N-dimensional space. Our simulations indicate that the dynamics within the native valley is well described by a sequence of thermally activated process up to temperatures well above the folding temperature. At larger temperatures, systematic deviations from the Langers estimate are instead observed. Several sources for such discrepancies are thoroughly discussed.

Vibrational entropy and the structural organization of proteins

Abstract.In this paper we analyze the vibrational spectra of a large ensemble of non-homologous protein structures by means of a novel tool, that we coin Hierarchical Network Model (HNM). Our coarse-grained scheme accounts for the intrinsic heterogeneity of force constants displayed by protein arrangements and also incorporates side chain degrees of freedom. Our analysis shows that vibrational entropy per unit residue correlates with the content of secondary structure. Furthermore, we assess the individual contribution to vibrational entropy of the novel features of our scheme as compared with the predictions of state-of-the-art network models. This analysis highlights the importance of properly accounting for the intrinsic hierarchy in force strengths typical of the different atomic bonds that build up and stabilize protein scaffolds. Finally, we discuss possible implications of our findings in the context of protein aggregation phenomena.

Resolving the geometry of biomolecules imaged by cryo electron tomography

In this paper, we describe two methods for computerized analysis of cryo electron tomography reconstructions of biomolecules. Both methods aim at quantifying the degree of structural flexibility of macromolecules and eventually resolving the inner dynamics through automatized protocols. The first method performs a Brownian dynamics evolution of a simplified molecular model into a fictitious force field generated by the tomograms. This procedure enables us to dock the simplified model into the experimental profiles. The second uses a fuzzy framework to delineate the subparts of the proteins and subsequently determine their interdomain relations. The two methods are discussed and their complementarities highlighted with reference to the case of the immonoglobulin antibody. Both artificial maps, constructed from immunoglobulin G entries in the Protein Data Bank and real tomograms are analyzed. Robustness issues and agreement with previously reported measurements are discussed.

Exploring the energy landscape of model proteins: a metric criterion for the determination of dynamical connectivity.

A method to reconstruct the energy landscape of small peptides is presented with reference to a two-dimensional off-lattice model. The starting point is a statistical analysis of the configurational distances between generic minima and directly connected pairs (DCP). As the mutual distance of DCP is typically much smaller than that of generic pairs, a metric criterion can be established to identify the great majority of DCP. Advantages and limits of this approach are thoroughly analyzed for three different heteropolymeric chains. A funnel-like structure of the energy landscape is found in all of the three cases, but the escape rates clearly reveal that the native configuration is more easily accessible (and is significantly more stable) for the sequence that is expected to behave as a real protein.

Graph theoretical analysis of the energy landscape of model polymers

In systems characterized by a rough potential-energy landscape, local energetic minima and saddles define a network of metastable states whose topology strongly influences the dynamics. Changes in temperature, causing the merging and splitting of metastable states, have nontrivial effects on such networks and must be taken into account. We do this by means of a recently proposed renormalization procedure. This method is applied to analyze the topology of the network of metastable states for different polypeptidic sequences in a minimalistic polymer model. A smaller spectral dimension emerges as a hallmark of stability of the global energy minimum and highlights a nonobvious link between dynamic and thermodynamic properties.