Reading the three-dimensional structure of lattice model-designed proteins from their amino acid sequence

Visualization of early events in acetic acid denaturation of HIV-1 protease: a molecular dynamics study

Protein denaturation plays a crucial role in cellular processes. In this study, denaturation of HIV-1 Protease (PR) was investigated by all-atom MD simulations in explicit solvent. The PR dimer and monomer were simulated separately in 9 M acetic acid (9 M AcOH) solution and water to study the denaturation process of PR in acetic acid environment. Direct visualization of the denaturation dynamics that is readily available from such simulations has been presented. Our simulations in 9 M AcOH reveal that the PR denaturation begins by separation of dimer into intact monomers and it is only after this separation that the monomer units start denaturing. The denaturation of the monomers is flagged off by the loss of crucial interactions between the α-helix at C-terminal and surrounding β-strands. This causes the structure to transit from the equilibrium dynamics to random non-equilibrating dynamics. Residence time calculations indicate that denaturation occurs via direct interaction of the acetic acid molecules with certain regions of the protein in 9 M AcOH. All these observations have helped to decipher a picture of the early events in acetic acid denaturation of PR and have illustrated that the α-helix and the β-sheet at the C-terminus of a native and functional PR dimer should maintain both the stability and the function of the enzyme and thus present newer targets for blocking PR function.

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.

Sequence of events in folding mechanism: beyond the Gō model

Simplified Gō 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 Gō 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.

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.