The denatured state of HIV-1 protease under native conditions.

The denatured state of several proteins has been shown to display transient structures that are relevant for folding, stability, and aggregation. To detect them by nuclear magnetic resonance (NMR) spectroscopy, the denatured state must be stabilized by chemical agents or changes in temperature. This makes the environment different from that experienced in biologically relevant processes. Using high-resolution heteronuclear NMR spectroscopy, we have characterized several denatured states of a monomeric variant of HIV-1 protease, which is natively structured in water, induced by different concentrations of urea, guanidinium chloride, and acetic acid. We have extrapolated the chemical shifts and the relaxation parameters to the denaturant-free denatured state at native conditions, showing that they converge to the same values. Subsequently, we characterized the conformational properties of this biologically relevant denatured state under native conditions by advanced molecular dynamics simulations and validated the results by comparison to experimental data. We show that the denatured state of HIV-1 protease under native conditions displays rich patterns of transient native and non-native structures, which could be of relevance to its guidance through a complex folding process.

Cold Denaturation of the HIV-1 Protease Monomer.

The human immunodeficiency virus-1 (HIV-1) protease is a complex protein that in its active form adopts a homodimer dominated by ß-sheet structures. We have discovered a cold-denatured state of the monomeric subunit of HIV-1 protease that is populated above 0 °C and therefore directly accessible to various spectroscopic approaches. Using nuclear magnetic resonance secondary chemical shifts, temperature coefficients, and protein dynamics, we suggest that the cold-denatured state populates a compact wet globule containing transient non-native-like α-helical elements. From the linearity of the temperature coefficients and the hydrodynamic radii, we propose that the overall architecture of the cold-denatured state is maintained over the temperature range studied.

The maturation of HIV-1 protease precursor studied by discrete molecular dynamics.

The equilibrium properties of a HIV-1-protease precursor are studied by means of an efficient molecular dynamics scheme, which allows for the simulation of the folding of the protein monomers and their dimerization into an active form and compare them with those of the mature protein. The results of the model provide, with atomic detail, an overall account of several experimental findings, including the NMR conformation of the mature dimer, the calorimetric properties of the system, the effects of the precursor tail on the dimerization constant, the secondary chemical shifts of the monomer, and the paramagnetic relaxation enhancement data associated with the conformations of the precursor. It is found that although the mature protein can dimerize in a unique, single way, the precursor populates several dimeric conformations in which monomers are always native-like, but their binding can be non-native.

Hierarchy of folding and unfolding events of protein G, CI2, and ACBP from explicit-solvent simulations.

The study of the mechanism which is at the basis of the phenomenon of protein folding requires the knowledge of multiple folding trajectories under biological conditions. Using a biasing molecular-dynamics algorithm based on the physics of the ratchet-and-pawl system, we carry out all-atom, explicit solvent simulations of the sequence of folding events which proteins G, CI2, and ACBP undergo in evolving from the denatured to the folded state. Starting from highly disordered conformations, the algorithm allows the proteins to reach, at the price of a modest computational effort, nativelike conformations, within a root mean square deviation (RMSD) of approximately 1 Å. A scheme is developed to extract, from the myriad of events, information concerning the sequence of native contact formation and of their eventual correlation. Such an analysis indicates that all the studied proteins fold hierarchically, through pathways which, although not deterministic, are well-defined with respect to the order of contact formation. The algorithm also allows one to study unfolding, a process which looks, to a large extent, like the reverse of the major folding pathway. This is also true in situations in which many pathways contribute to the folding process, like in the case of protein G.

Identification of the folding inhibitors of hen-egg lysozyme: gathering the right tools.

The unfolded state of proteins displays a surprisingly rich amount of local native structure, which appears to be critical for driving the protein to its native state. Peptides with the same sequence of the corresponding structured segments can be used to interfere with the correct folding of the protein. Using model simulations, we investigate the folding of hen-egg lysozyme, identifying its key segments. Activity assays, NMR and circular dichroism experiments are used to screen the peptides which are able to inhibit the folding of lysozyme. Few peptides, corresponding to the segments of the protein which are structured in the unfolded state, are identified to have significant inhibitory effects.

Atomistic simulations of the HIV-1 protease folding inhibition.

Biochemical experiments have recently revealed that the p-S8 peptide, with an amino-acid sequence identical to the conserved fragment 83-93 (S8) of the HIV-1 protease, can inhibit catalytic activity of the enzyme by interfering with protease folding and dimerization. In this study, we introduce a hierarchical modeling approach for understanding the molecular basis of the HIV-1 protease folding inhibition. Coarse-grained molecular docking simulations of the flexible p-S8 peptide with the ensembles of HIV-1 protease monomers have revealed structurally different complexes of the p-S8 peptide, which can be formed by targeting the conserved segment 24-34 (S2) of the folding nucleus (folding inhibition) and by interacting with the antiparallel termini beta-sheet region (dimerization inhibition). All-atom molecular dynamics simulations of the inhibitor complexes with the HIV-1 PR monomer have been independently carried out for the predicted folding and dimerization binding modes of the p-S8 peptide, confirming the thermodynamic stability of these complexes. Binding free-energy calculations of the p-S8 peptide and its active analogs are then performed using molecular dynamics trajectories of the peptide complexes with the HIV-1 PR monomers. The results of this study have provided a plausible molecular model for the inhibitor intervention with the HIV-1 PR folding and dimerization and have accurately reproduced the experimental inhibition profiles of the active folding inhibitors.

Exploring the protein G helix free-energy surface by solute tempering metadynamics.

The free-energy landscape of the alpha-helix of protein G is studied by means of metadynamics coupled with a solute tempering algorithm. Metadynamics allows to overcome large energy barriers, whereas solute tempering improves the sampling with an affordable computational effort. From the sampled free-energy surface we are able to reproduce a number of experimental observations, such as the fact that the lowest minimum corresponds to a globular conformation displaying some degree of beta-structure, that the helical state is metastable and involves only 65% of the chain. The calculations also show that the system populates consistently a pi-helix state and that the hydrophobic staple motif is present only in the free-energy minimum associated with the helices, and contributes to their stabilization. The use of metadynamics coupled with solute tempering results then particularly suitable to provide the thermodynamics of a short peptide, and its computational efficiency is promising to deal with larger proteins.

Optical absorption of a green fluorescent protein variant: environment effects in a density functional study.

We present an ab initio study of the optical absorption properties of a particularly interesting fluorescent protein (E2GFP), whose complex photophysics still escapes elucidation. In particular, we focus on the role of the protein environment, showing that the effects of both nearby residues and the external field due to residues not accounted for explicitly are needed to properly reproduce the experimental data. The spectra calculated taking such contributions into account provide for the first time a robust identification of the states relevant for the photophysics of this system.

Sequence of events in folding mechanism: beyond the Go model.

Simplified Go models, where only native contacts interact favorably, have proven useful to characterize some aspects of the folding of small proteins. The success of these models is limited by the fact that all residues interact in the same way so that the folding features of a protein are determined only by the geometry of its native conformation. We present an extended version of a Calpha-based Go model where different residues interact with different energies. The model is used to calculate the thermodynamics of three small proteins (Protein G, Src-SH3, and CI2) and the effect of mutations (DeltaDeltaGU-N, DeltaDeltaGdouble dagger-N, DeltaDeltaGdouble dagger-U, and phi-values) on the wild-type sequence. The model allows us to investigate some of the most controversial areas in protein folding, such as its earliest stages and the nature of the unfolded state, subjects that have lately received particular attention.

The determinants of stability in the human prion protein: insights into folding and misfolding from the analysis of the change in the stabilization energy distribution in different conditions.

The dynamic evolution of the PrP(C) from its NMR-derived conformation to a beta-sheet-rich, aggregation-prone conformation is studied through all-atom, explicit solvent molecular dynamics in different temperature and pH conditions. The trajectories are analyzed by means of a recently introduced energy decomposition approach aimed at identifying the key residues for the stabilization and folding of the protein. It is shown that under native conditions the stabilization energy is concentrated in regions of the helices H1 and H3, whereas under misfolding conditions (low pH, high temperature, or mutations in selected sites) it is spread out over helix H2. Misfolding appears to be a rearrangement of the chain that disrupts most of the native secondary structure of the protein, producing some beta-rich conformations with an energy distribution similar to that of the native state.

Design of HIV-1-PR inhibitors that do not create resistance: blocking the folding of single monomers.

The main problems found in designing drugs are those of optimizing the drug-target interaction and of avoiding the insurgence of resistance. We suggest a scheme for the design of inhibitors that can be used as leads for the development of a drug and that do not face either of these problems, and then apply it to the case of HIV-1-PR. It is based on the knowledge that the folding of single-domain proteins, such as each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES), stabilized by local contacts among hydrophobic, strongly interacting, and highly conserved amino acids that play a central role in the folding process. Because LES have evolved over many generations to recognize and strongly interact with each other so as to make the protein fold fast and avoid aggregation with other proteins, highly specific (and thus little toxic) as well as effective folding-inhibitor molecules suggest themselves: short peptides (or eventually their mimetic molecules) displaying the same amino acid sequence of that of LES (p-LES). Aside from being specific and efficient, these inhibitors are expected not to induce resistance; in fact, mutations in HIV-1-PR that successfully avoid the action of p-LES imply the destabilization of one or more LES and thus should lead to protein denaturation. Making use of Monte Carlo simulations, we first identify the LES of the HIV-1-PR and then show that the corresponding p-LES peptides act as effective inhibitors of the folding of the protease.

Understanding the determinants of stability and folding of small globular proteins from their energetics.

The results of minimal model calculations indicate that the stability and the kinetic accessibility of the native state of small globular proteins are controlled by few "hot" sites. By means of molecular dynamics simulations around the native conformation, which describe the protein and the surrounding solvent at the all-atom level, an accurate and compact energetic map of the native state of the protein is generated. This map is further simplified by means of an eigenvalue decomposition. The components of the eigenvector associated with the lowest eigenvalue indicate which hot sites are likely to be responsible for the stability and for the rapid folding of the protein. The comparison of the results of the model with the findings of mutagenesis experiments performed for four small proteins show that the eigenvalue decomposition method is able to identify between 60% and 80% of these (hot) sites.

Design and folding of dimeric proteins.

In a similar way in which the folding of single-domain proteins provides an important test in the study of self-organization, the folding of homodimers constitutes a basic challenge in the quest for the mechanisms that are the basis of biological recognition. Dimerization is studied by following the evolution of two identical 20-letter amino acid chains within the framework of a lattice model and using Monte Carlo simulations. It is found that when design (evolution pressure) selects few, strongly interacting (conserved) amino acids to control the process, a three-state folding scenario follows, where the monomers first fold forming the halves of the eventual dimeric interface independently of each other, and then dimerize ("lock and key" kind of association). On the other hand, if design distributes the control of the folding process on a large number of (conserved) amino acids, a two-state folding scenario ensues, where dimerization takes place at the beginning of the process, resulting in an "induced type" of association. Making use of conservation patterns of families of analogous dimers, it is possible to compare the model predictions with the behavior of real proteins. It is found that theory provides an overall account of the experimental findings.

Beta-hairpin conformation of fibrillogenic peptides: structure and alpha-beta transition mechanism revealed by molecular dynamics simulations.

Understanding the conformational transitions that trigger the aggregation and amyloidogenesis of otherwise soluble peptides at atomic resolution is of fundamental relevance for the design of effective therapeutic agents against amyloid-related disorders. In the present study the transition from ideal alpha-helical to beta-hairpin conformations is revealed by long timescale molecular dynamics simulations in explicit water solvent, for two well-known amyloidogenic peptides: the H1 peptide from prion protein and the Abeta(12-28) fragment from the Abeta(1-42) peptide responsible for Alzheimer's disease. The simulations highlight the unfolding of alpha-helices, followed by the formation of bent conformations and a final convergence to ordered in register beta-hairpin conformations. The beta-hairpins observed, despite different sequences, exhibit a common dynamic behavior and the presence of a peculiar pattern of the hydrophobic side-chains, in particular in the region of the turns. These observations hint at a possible common aggregation mechanism for the onset of different amyloid diseases and a common mechanism in the transition to the beta-hairpin structures. Furthermore the simulations presented herein evidence the stabilization of the alpha-helical conformations induced by the presence of an organic fluorinated cosolvent. The results of MD simulation in 2,2,2-trifluoroethanol (TFE)/water mixture provide further evidence that the peptide coating effect of TFE molecules is responsible for the stabilization of the soluble helical conformation.

Modeling the alpha-helix to beta-hairpin transition mechanism and the formation of oligomeric aggregates of the fibrillogenic peptide Abeta(12-28): insights from all-atom molecular dynamics simulations.

In this paper, all-atom molecular dynamics simulations in explicit solvent are used to investigate the structural and dynamical determinants of the alpha-helical to beta-hairpin conformational transition of the 12-28 fragment from the full length Abeta Alzheimer's peptide. The transition from alpha-helical to beta-structure requires the peptide to populate intermediate beta-bend geometries in which several mainly hydrophobic interactions are partially formed. This is followed by the sudden collapse to ordered beta-hairpin structures and the simultaneous disruption of the hydrophobic side-chain interactions with a consequent increase in solvent exposure. The solvent exposure of hydrophobic side-chains belonging to a sequence of five consecutive residues in the beta-hairpin defines a possible starting point for the onset of the aggregation mechanisms. Several different conformations of model oligomeric (dimeric and tetrameric) aggregates are then investigated. These simulations show that while hydrophobic contacts are important to bring together different monomers with a beta-hairpin like conformation, more specific interactions such as hydrogen-bonding and coulombic interactions, should be considered necessary to provide further stabilization and ordering to the nascent fibrillar aggregates.

Early events in protein folding: Is there something more than hydrophobic burst?

The presence of native contacts in the denatured state of many proteins suggests that elements of the biologically active structure of these molecules are formed during the initial stage of the folding process. The rapidity with which these events take place makes it difficult to study them in vitro, but, by the same token, suitable for studies in silico. With the help of all-atom, explicit solvent, molecular dynamics simulations we have followed in time, starting from elongated structureless conformations, the early events in the folding of src-SH3 domain and of proteins G, L, and CI2. It is observed that within the first 50 ns two important events take place, essentially independent of each other: hydrophobic collapse and formation of a few selected native contacts. The same contacts are also found in simulations carried out in the presence of guanidinium chloride in order to reproduce the conditions used to characterize experimentally the denatured state and testify to the fact that these contacts are to be considered a resilient characterizing property of the denaturated state.

Insight into the folding inhibition of the HIV-1 protease by a small peptide.

It has recently been shown that the highly protected segments 24-34 (S2) and 83-93 (S8) of each of the two 99-mers of human immunodeficiency virus type 1 protease play an essential role in the folding of the monomers, giving rise to the so-called (postcritical) folding nucleus ((FN) minimum condensation unit ensuring folding) when they dock. This scenario received further support from model calculations that demonstrated that the peptide p-S8, displaying an amino acid sequence identical to the corresponding (83-93) segment of the monomer, can be used to interfere with the formation of the FN and eventually to inhibit folding by docking the fragment 24-34. Experiments in vitro and in cells infected with ex vivo wild-type and multiresistant HIV isolates confirm that the inhibition power of p-S8 is robust. On the other hand, there is no direct evidence demonstrating the validity of the proposed mechanism of inhibition associated with p-S8. To shed light on this question and to provide the basis for the design of a molecule mimetic to p-S8, to be used as lead of an eventual drug against AIDS, we study, in this paper, with the help of all-atom simulations in explicit solvent and the novel method of metadynamics combined with parallel tempering: a), the free energy and the equilibrium structure of each of the peptides p-S2 and p-S8; b), the details of the docking mechanism of the two peptides and the free energy associated with this process. Whereas p-S8 is found to be well structured, p-S2 is rather flexible, wrapping itself around p-S8 to give rise to the FN, which is stabilized by three particular hydrogen bonds.

The evolution dynamics of model proteins.

Explicit simulations of protein evolution, where protein chains are described at a molecular, although simplified, level provide important information to understand the similarities found to exist between known proteins. The results of such simulations suggest that a number of evolutionary-related quantities, such as the distribution of sequence similarity for structurally similar proteins, are controlled by evolutionary kinetics and do not reflect an equilibrium state. An important result for phylogeny is that a subset of the residues of each protein evolve on a much larger time scale than the other residues.