Contributions of the GABAergic system of the prelimbic cortex and basolateral amygdala to morphine withdrawal-induced contextual fear.

Morphine withdrawal can trigger disruptions in neuronal pathways involved in the modulation and expression of anxiety and fear-related behaviors, particularly those involved in associative learning. When it comes to contextual fear, specific subdivisions of the medial prefrontal cortex (mPFC) regulate the expression of defensive behaviors through projections to specific amygdala (AM) nuclei, such as the prelimbic cortex (PrL). The basolateral nucleus (BLA) of the AM has been shown to be involved in the modulation and expression of associative memories of fear, including those associated with opiate withdrawal-related aversive events. The purpose of this study is to determine the role of GABA mechanisms in the PrL and BLA in startle potentiation and freezing behavior caused by morphine-precipitated withdrawal. Our findings show that morphine withdrawal promotes the emergence of contextual conditioned fear in animals when they are exposed to the same environment where the withdrawal sessions were performed. This suggests that the neural circuits underlying the organism's response to conditioned stressors and the circuits modulating the negative affective states induced by drug withdrawal may overlap. The pharmacological manipulation of GABAergic neurotransmission in the PrL and BLA can reverse contextual fear in morphine-withdrawn rats, an effect that appears to be mediated, at least in part, by GABAA receptors.

Different pulling modes in DNA overstretching: a theoretical analysis.

We study the thermally driven denaturation of a double-stranded polymer in the presence of a stretching force via Monte-Carlo simulations. When one strand only is stretched, the denaturation transition is first order, while when both strands are stretched, melting is second order. By revisiting the Poland-Scheraga model for DNA melting, we show that at room temperature, the most likely scenario is that DNA melts as it overstretches. Our results are in general agreement with the most recent experiments and suggest how varying temperature and stretching mode may help settle the question whether S-DNA exists or not.

Simple solvation potential for coarse-grained models of proteins.

We formulate a simple solvation potential based on a coarsed-grained representation of amino acids with two spheres modeling the C(alpha) atom and an effective side-chain centroid. The potential relies on a new method for estimating the buried area of residues, based on counting the effective number of burying neighbors in a suitable way. This latter quantity shows a good correlation with the buried area of residues computed from all atom crystallographic structures. We check the discriminatory power of the solvation potential alone to identify the native fold of a protein from a set of decoys and show the potential to be considerably selective.

Form of growing strings.

Patterns and forms adopted by nature are often the results of simple dynamical paradigms. Here we show that a growing self-interacting string attached to a tracking origin, modeled to resemble nascent polypeptides in vivo, develops helical structures which are more pronounced at the growing end. We also show that the dynamic growth ensemble shares several features of an equilibrium ensemble in which the growing end of the polymer is under an effective stretching force. A statistical analysis of native states of proteins shows that the signature of this nonequilibrium phenomenon has been fixed by evolution at the C terminus, the growing end of a nascent protein. These findings suggest how evolution may have built on the properties of a generic nonequilibrium growth process in favoring helical structures in nascent chains.

Stepwise unfolding of collapsed polymers.

Motivated by recent experimental data on DNA stretching in presence of polyvalent counterions, we study the force-induced unfolding of a homopolymer on and off lattice. In the fixed force ensemble the globule unravels via a series of steps due to surface effects which play an important role for finite-size chains. This holds both for flexible and stiff polymers. We discuss in a qualitative way how this result may impact on the interpretation of DNA stretching experiments showing peaks in the characteristic curves, by extracting from the raw data the corresponding elongation- versus-force characteristic curves. Furthermore, approximate analytical and numerical calculations, valid in a quasi-equilibrium fixed stretch ensemble, and if the initial low-temperature state is ordered in a spool, show that the average force versus elongation displays peaks related to the geometry of the initial configuration. We finally argue how the proposed mechanisms identified for the arising of peaks may couple in the experiments, and comment on the role of dynamic effects.

Elucidation of the disulfide-folding pathway of hirudin by a topology-based approach.

A theoretical model for the folding of proteins containing disulfide bonds is introduced. The model exploits the knowledge of the native state to favor the progressive establishment of native interactions. At variance with traditional approaches based on native topology, not all native bonds are treated in the same way; in particular, a suitable energy term is introduced to account for the special strength of disulfide bonds, as well as their ability to undergo intramolecular reshuffling. The model thus possesses the minimal ingredients necessary to investigate the much debated issue of whether the refolding process occurs through partially structured intermediates with native or non-native disulfide bonds. This strategy is applied to a context of particular interest, the refolding process of hirudin, a thrombin-specific protease inhibitor, for which conflicting folding pathways have been proposed. We show that the only two parameters in the model (temperature and disulfide strength) can be tuned to reproduce well a set of experimental transitions between species with different number of formed disulfides. This model is then used to provide a characterization of the folding process and a detailed description of the species involved in the rate-limiting step of hirudin refolding.

Mechanical unfolding of directed polymers in a poor solvent: critical exponents.

We study the thermodynamics of an exactly solvable model of a self-interacting, partially directed self-avoiding walk in two dimensions when a force is applied on one end of the chain. The critical force for the unfolding is determined exactly, as a function of the temperature, below the Theta transition. The transition is of second order and is characterized by new critical exponents that are determined by a careful numerical analysis. The usual polymer critical index nu on the critical line, and another one which we call zeta, takes a nontrivial value that is numerically close to 2/3.

Stretching of a polymer below the theta point.

The unfolding of a polymer below the theta point when pulled by an external force is studied both in d=2 on the lattice and in d=3 off the lattice. At T=0 and for finite length chains, it is found that the globule unfolds via multiple steps, corresponding to transitions between different minima, in both cases. In d=3 one of these intermediates is a regular helix. In the infinite length limit, these steps have a qualitative effect only in d=2. The phase diagram in d=2 is determined via the transfer matrix. To rationalize these results, energy-entropy and renormalization group arguments are given.

Optimal potentials for predicting inter-helical packing in transmembrane proteins.

A set of pairwise contact potentials between amino acid residues in transmembrane helices was determined from the known native structure of the transmembrane protein (TMP) bacteriorhodopsin by the method of perceptron learning, using Monte Carlo dynamics to generate suitable "decoy" structures. The procedure of finding these decoys is simpler than for globular proteins, since it is reasonable to assume that helices behave as independent, stable objects and, therefore, the search in the conformational space is greatly reduced. With the learnt potentials, the association of the helices in bacteriorhodopsin was successfully simulated. The folding of a second TMP (the helix-dimer glycophorin A) was then accomplished with only a refinement of the potentials from a small number of decoys.

Dynamical scaling of the DNA unzipping transition.

We report studies of the dynamics of a set of exactly solvable lattice models for the force-induced DNA unzipping transition. Besides yielding the whole equilibrium phase diagram, which reveals a reentrance, these models enable us to characterize the dynamics of the process starting from a nonequilibrium initial condition. The thermal melting of DNA displays a model dependent time evolution. On the contrary, the dynamical mechanism for the unzipping by force is very robust and the scaling behavior is independent of the details of the description and, hence, superuniversal.

A self-consistent knowledge-based approach to protein design.

A simple and very efficient protein design strategy is proposed by developing some recently introduced theoretical tools which have been successfully applied to exactly solvable protein models. The design approach is implemented by using three amino acid classes and it is based on the minimization of an appropriate energy function. For a given native state the results of the design procedure are compared, through a statistical analysis, with the properties of an ensemble of sequences folding in the same conformation. If the success rate is computed on those sites designed with high confidence, it can be as high as 80%. The method is also able to identify key sites for the folding process: results for 2ci2 and barnase are in very good agreement with experimental results.

Learning effective amino acid interactions through iterative stochastic techniques.

The prediction of the three-dimensional structures of the native states of proteins from the sequences of their amino acids is one of the most important challenges in molecular biology. An essential task for solving this problem within coarse-grained models is the deduction of effective interaction potentials between the amino acids. Over the years, several techniques have been developed to extract potentials that are able to discriminate satisfactorily between the native and nonnative folds of a preassigned protein sequence. In general, when these potentials are used in actual dynamical folding simulations, they lead to a drift of the native structure outside the quasinative basin. In this article, we present and validate an approach to overcome this difficulty. By exploiting several numerical and analytical tools, we set up a rigorous iterative scheme to extract potentials satisfying a prerequisite of any viable potential: the stabilization of proteins within their native basin (less than 3-4 A RMSD). The scheme is flexible and is demonstrated to be applicable to a variety of parameterizations of the energy function, and it provides in each case the optimal potentials.

Deciphering the folding kinetics of transmembrane helical proteins.

Nearly a quarter of genomic sequences and almost half of all receptors that are likely to be targets for drug design are integral membrane proteins. Understanding the detailed mechanisms of the folding of membrane proteins is a largely unsolved, key problem in structural biology. Here, we introduce a general model and use computer simulations to study the equilibrium properties and the folding kinetics of a C(alpha)-based two-helix bundle fragment (comprised of 66 aa) of bacteriorhodopsin. Various intermediates are identified and their free energy are calculated together with the free energy barrier between them. In 40% of folding trajectories, the folding rate is considerably increased by the presence of nonobligatory intermediates acting as traps. In all cases, a substantial portion of the helices is rapidly formed. This initial stage is followed by a long period of consolidation of the helices accompanied by their correct packing within the membrane. Our results provide the framework for understanding the variety of folding pathways of helical transmembrane proteins.

Adsorptionlike collapse of diblock copolymers

A linear copolymer made of two reciprocally attracting N-monomer blocks collapses to a compact phase through a novel transition, whose exponents are determined with extensive Monte Carlo simulations in two and three dimensions. In the former case, an identification with the statistical geometry of suitable percolation paths allows one to predict that the number of contacts between the blocks grows like N9/16. In the compact phase the blocks are mixed and, in two dimensions, also zipped, in such a way to form a spiral, double chain structure.

Recurrent oligomers in proteins: an optimal scheme reconciling accurate and concise backbone representations in automated folding and design studies.

A novel scheme is introduced to capture the spatial correlations of consecutive amino acids in naturally occurring proteins. This knowledge-based strategy is able to carry out optimally automated subdivisions of protein fragments into classes of similarity. The goal is to provide the minimal set of protein oligomers (termed "oligons" for brevity) that is able to represent any other fragment. At variance with previous studies in which recurrent local motifs were classified, our concern is to provide simplified protein representations that have been optimised for use in automated folding and/or design attempts. In such contexts, it is paramount to limit the number of degrees of freedom per amino acid without incurring loss of accuracy of structural representations. The suggested method finds, by construction, the optimal compromise between these needs. Several possible oligon lengths are considered. It is shown that meaningful classifications cannot be done for lengths greater than six or smaller than four. Different contexts are considered for which oligons of length five or six are recommendable. With only a few dozen oligons of such length, virtually any protein can be reproduced within typical experimental uncertainties. Structural data for the oligons are made publicly available.

Simple model to study insertion of a protein into a membrane.

A simple coarse grained model on a two-dimensional lattice is presented to elucidate the main effects ruling the insertion of a protein into a polar environment such as a lipidic membrane. The amino acids are divided into two classes (hydrophobic or polar), and they behave differently according to their surroundings. In aqueous solution the hydrophobic amino acids are forced to minimize contacts with water, whereas in the apolar environment all the amino acids try to aggregate regardless to their specificity. The lattice is employed in order to perform exact calculations and to generate a fictitious protein data bank. Despite the simplicity of the model, some morphological features of the protein-like lattice structures obtained by our model are compatible with the observed phenomenology of transmembrane proteins. These results seem to corroborate the hypothesis that the number of classes into which the amino acids can be divided that correctly describe the phenomena may be extremely low.

Design of proteins with hydrophobic and polar amino acids.

A two amino acid (hydrophobic and polar) scheme is used to perform the design on target conformations corresponding to the native states of 20 single chain proteins. Strikingly, the percentage of successful identification of the nature of the residues benchmarked against naturally occurring proteins and their homologues is around 75%, independent of the complexity of the design procedure. Typically, the lowest success rate occurs for residues such as alanine that have a high secondary structure functionality. Using a simple lattice model, we argue that one possible shortcoming of the model studied may involve the coarse-graining of the 20 kinds of amino acids into just two effective types.

Structure-based design of model proteins.

A structure-based, sequence-design procedure is proposed in which one considers a set of decoy structures that compete significantly with the target structure in being low energy conformations. The decoy structures are chosen to have strong overlaps in contacts with the putative native state. The procedure allows the design of sequences with large and small stability gaps in a random-bond heteropolymer model in both two and three dimensions by an appropriate assignment of the contact energies to both the native and nonnative contacts. The design procedure is also successfully applied to the two-dimensional HP model.

Interaction potentials for protein folding.

We outline a general strategy for determining the effective coarse-grained interactions between the amino acids of a protein from the experimentally derived native-state structures. The method is, in principle, free from any adjustable or empirically determined parameters, and it is tested on simple models and compared with other existing approaches.