An AT-barrier mechanically controls DNA reannealing under tension.

Regulation of genomic activity occurs through the manipulation of DNA by competent mechanoenzymes. Force-clamp optical tweezers that allow the structural dynamics of the DNA molecule to be measured were used here to investigate the kinetics of mechanically-driven strand reannealing. When the force on the torsionally unconstrained λ-phage DNA is decreased stepwise from above to below the overstretching transition, reannealing occurs via discrete shortening steps separated by exponentially distributed time intervals. Kinetic analysis reveals a transition barrier 0.58 nm along the reaction coordinate and an average reannealing-step size of ∼750 bp, consistent with the average bp interval separating segments of more than 10 consecutive AT bases. In an AT-rich DNA construct, in which the distance between segments of more than 10 consecutive AT is reduced to ∼210 bps, the reannealing step reduces accordingly without changes in the position of the transition barrier. Thus, the transition barrier for reannealing is determined by the presence of segments of more than 10 consecutive AT bps independent of changes in sequence composition, while the length of the reannealing strand changes according to the distance between poly-AT segments at least 10 bps long.

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.

Vibrational entropy and the structural organization of proteins.

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 immunoglobulin 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.

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.

Anti-cooperativity in diffusion-controlled reactions with pairs of anisotropic domains: a model for the antigen-antibody encounter.

The encounter between anisotropic agents in diffusion-controlled reactions is a topic of very general relevance in chemistry and biology. Here we introduce a simplified model of encounter of an isotropic molecule with a pair of partially reacting agents and apply it to the encounter reaction between an antibody and its antigen. We reduce the problem to the solution of dual series relations, which can be solved iteratively, yielding the exact solution for the encounter rate constant at any desired order of accuracy. We quantify the encounter effectiveness by means of a simple indicator and show that the two binding centers systematically behave in an anti-cooperative fashion. However, we demonstrate that a reduction of the binding active sites allows the composite molecule to recover binding effectiveness, in spite of the overall reduction of the rate constant. In addition, we provide a simple formula that enables one to calculate the anti-cooperativity as a function of the size of the binding site for any values of the separation between the two active lobes and of the antigen size. Finally, some biological implications of our results are discussed.

Topological properties of the energy landscape of small peptides.

The topology of the potential energy landscape (PEL) underlying the dynamics of a two dimensional off-lattice model for a heteropolymer is analyzed for different sequences of amino-acids. A statistical characterization of the metastable minima and first-order saddles of the PEL highlights structural differences in the landscape of good and bad folding sequences and provides insight on the chain dynamics during folding.

Dynamics of antibodies from cryo-electron tomography.

The issue of protein dynamics and its implications in the biological function of proteins are arousing greater and greater interest in molecular biology. In cryo-electron tomography experiments one takes several snapshots of a given biological macromolecule. In principle, a large enough collection of snapshots may then be used to calculate its equilibrium configuration in terms of the experimentally accessible degrees of freedom, and hence estimate its potential energy. Consequently, one could analyze the biological functions of biomolecules by directly accessing their dynamics. In this work, we analyze the results of cryo-electron tomography experiments on monoclonal murine IgG2a antibodies. With the aid of a novel software for image processing, we measure the equilibrium distribution of the angles which describe the configuration of the molecule. This helps us shed some critical light on recent results from X-ray crystallography. We then build a model of the antibody dynamics, which enables us to use the measured angular distribution in order to derive an explicit expression of the IgG potential energy. Finally, as a preliminary application of our results, we investigate the dynamical effects in the rate of formation of the antigen-antibody encounter complex. In particular, we suggest that the dynamics of antibodies operates in the direction of decreasing anticooperativity of the two antigen binding arms.

Freezing immunoglobulins to see them move.

The issue of protein dynamics and its implications in the biological function of proteins are arousing greater and greater interest in different fields of molecular biology. In cryo-electron tomography experiments one may take several snapshots of a given biological macromolecule. In principle, a large enough collection of snapshots of the molecule may then be used to calculate its equilibrium configuration in terms of the experimentally accessible degrees of freedom and, hence, to estimate its potential energy. This information would be crucial in order to analyze the biological functions of biomolecules by directly accessing the relevant dynamical indicators. In this article, we analyze the results of cryo-electron tomography experiments performed on monoclonal murine IgG2a antibodies. We measure the equilibrium distribution of the molecule in terms of the relevant angular coordinates and build a mechanical model of the antibody dynamics. This approach enables us to derive an explicit expression of the IgG potential energy. Furthermore, we discuss the configuration space at equilibrium in relation to results from other techniques, and we set our discussion in the context of the current debate regarding conformation and flexibility of antibodies.