Document Detail

Visualization of early events in acetic acid denaturation of HIV-1 protease: a molecular dynamics study.
Jump to Full Text
MedLine Citation:
PMID:  21738569     Owner:  NLM     Status:  MEDLINE    
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.
Aditi Narendra Borkar; Manoj Kumar Rout; Ramakrishna V Hosur
Related Documents :
22548209 - Changes in major peanut allergens under different ph conditions.
6738319 - Composition and metabolism of fatty acids in phospholipids of density-separated red cel...
10591409 - Regulation of brain glycosylphosphatidylinositol-specific phospholipase d by natural am...
1182589 - Inhibition of cardiolipin synthesis by end-products and other complex lipids in membran...
4014459 - Energy dependence of lung phosphatidylcholine biosynthesis studied with co hypoxia.
8452869 - Phosphatidic acid stimulates the rolipram-sensitive cyclic nucleotide phosphodiesterase...
10516619 - Quantitative 19f nmr study of trifluorothymidine metabolism in rat brain.
21635219 - Mechanism of acid tolerance in a rhizobium strain isolated from pueraria lobata (willd....
16360039 - Mechanism of lysine 48-linked ubiquitin-chain synthesis by the cullin-ring ubiquitin-li...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2011-06-29
Journal Detail:
Title:  PloS one     Volume:  6     ISSN:  1932-6203     ISO Abbreviation:  PLoS ONE     Publication Date:  2011  
Date Detail:
Created Date:  2011-07-08     Completed Date:  2011-12-01     Revised Date:  2013-06-28    
Medline Journal Info:
Nlm Unique ID:  101285081     Medline TA:  PLoS One     Country:  United States    
Other Details:
Languages:  eng     Pagination:  e19830     Citation Subset:  IM    
Institute of Bioinformatics and Biotechnology, University of Pune, Ganeshkhind, Pune, India.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Acetic Acid / chemistry*
HIV Protease / chemistry*
Molecular Dynamics Simulation*
Protein Denaturation
Reg. No./Substance:
64-19-7/Acetic Acid; EC 3.4.23.-/HIV Protease

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): PLoS One
Journal ID (publisher-id): plos
Journal ID (pmc): plosone
ISSN: 1932-6203
Publisher: Public Library of Science, San Francisco, USA
Article Information
Download PDF
Borkar et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Received Day: 8 Month: 11 Year: 2010
Accepted Day: 18 Month: 4 Year: 2011
collection publication date: Year: 2011
Electronic publication date: Day: 29 Month: 6 Year: 2011
Volume: 6 Issue: 6
E-location ID: e19830
ID: 3126794
PubMed Id: 21738569
Publisher Id: 10-PONE-RA-20796
DOI: 10.1371/journal.pone.0019830

Visualization of Early Events in Acetic Acid Denaturation of HIV-1 Protease: A Molecular Dynamics Study Alternate Title:MD of HIV-PR in Acetic Acid
Aditi Narendra Borkar12¤
Manoj Kumar Rout2
Ramakrishna V. Hosur23*
Darren R. Floweredit1 Role: Editor
1Institute of Bioinformatics and Biotechnology, University of Pune, Ganeshkhind, Pune, India
2Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai, India
3UM-DAE Centre for Excellence in Basic Sciences, Mumbai University Campus, Kalina, Santa Cruz Mumbai, India
Aston University, United Kingdom
Correspondence: * E-mail:
Contributed by footnote: Conceived and designed the experiments: ANB RVH MKR. Performed the experiments: ANB. Analyzed the data: ANB. Contributed reagents/materials/analysis tools: RVH. Wrote the paper: ANB RVH.
¤Current address: Department of Chemistry, University of Cambridge, Cambridge, United Kingdom


Protein folding/unfolding, dynamics and denaturation play crucial role in cellular processes and thus have been the subject of extensive investigation for the last several decades. A variety of experimental and theoretical approaches have been employed to explore various aspects of the processes. While these have led to some generalized concepts, the mechanistic details remain largely unclear. They can vary from system to system - which is perhaps not too surprising – and residue level details are very seldom available.

The HIV-1 protease (PR) is a homodimeric aspartate protease that plays sensitive function in HIV-1 maturation. Hence, it is the subject of extensive pharmaceutical research with strategies encompassing blocking of both PR function and structure dynamics [1], [2]. Thus, both the native state and the denaturation dynamics of PR become crucial topics of study to answer the structure-to-function relation in PR. While there have been many experimental studies on the structure and dynamics of PR complexed to many inhibitors, reports on free protease are rather few because of the autolytic nature of the protease, and the number of reports on the denatured states is even less. Previous studies by NMR [3], [4], [5] and fluorescence experiments [6] have shown that the denatured states of PR created by different denaturants are far from random coils and contain elements of both native and non-native structures. Other experimental analyses of the ‘unfolded’ state of PR even under strong denaturing conditions have also revealed presence of transient folding nuclei (FN) that include residues spanning the active site, the hinge region, and the dimerization domain [7]. Using NMR investigations, acid denatured states [8] have been shown to be significantly different from the guanidine denatured tethered dimer of PR [9] or urea denatured PR precursor having the TFR extension at the N-terminal [10].

However, all these experimental studies can only allude to the regions of the PR that are structured or denatured but cannot provide direct structural visualization of the dynamics of the protein as a whole in the native or denaturing environments. Such visualization is however possible with theoretical techniques like all-atom molecular dynamics simulations (MD) and plethora of other secondary structure calculation algorithms. This provides the motivation for taking up extensive theoretical studies for understanding protein structure dynamics and denaturation.

Many MD studies on PR mechanism [11][13], drug resistance [14][16] and unfolding [17][21] have already been reported. While the mechanism and drug resistance simulations generally employ liganded PR as the starting structure, the unfolding simulations were performed with higher temperature as the denaturing condition. Extensive all-atom 100 ns simulation of both the monomeric and dimeric PR has shown greater flexibility of the termini in the monomer [17] and decoupling between monomer folding and dimerization [18]. Simulation of the 99-residue monomer in water at different temperatures revealed early assembly of the N and C termini into stable β-sheet structures [19]. These studies emphasize that PR is not a two-state dimer, as indicated by equilibrium denaturation experiments, but a three-state dimer with a marginally stable monomeric intermediate which involves the swapping of the flexible termini across the two chains to form the dimer interface. Apart from the termini, the ‘flaps’ too have been shown to be highly flexible regions of the PR [20], [21]. The flap dynamics is considered to be essential in target entry and exit from the PR active site and corresponding open and closed conformations are observed in the crystal structures. MD simulation has also revealed that the unliganded protease predominantly populates the semi open conformation, with closed and fully open structures being a minor component of the overall ensemble and also provides a model for such flap opening and closing [20].

All these all-atom MD simulations of PR dealt with gross structure dynamics only and did not present a consolidated picture of the denaturation process, sequential loss of residue level interactions etc. Thus, in this study, we have derived this information using all-atom molecular dynamics simulation of PR dimer and monomer in acetic acid and water environments separately. For illustration we have chosen here acetic acid denaturation, but such studies can be extended to denaturation by other agents as well.

Results and Discussion
1. Denaturation of the dimeric PR

Due to the presence of high concentration of acetic acid molecules, the simulation is expected to be a non-equilibration run and the root mean square deviations (RMSD) of the backbone atoms of the dimeric PR trajectory with respect to the initial structure [Figure 1(a)] appropriately shows large continuous fluctuations of more than 0.1 nm and do not converge. The RMSD value shows a steep immediate ascent to 0.35 nm from 0.2 nm at around 10 ns. Following this, it rises further up to 0.6 nm and fluctuates substantially but never falls below the 0.3 nm mark. Contrary to this, in the simulation of the dimer in water, soon after 1 ns, the RMSD fluctuates only between 0.25–0.3 nm and thus the structure can be said to have quickly converged in the water simulation [Figure S1]. Thus, at around 10 ns during the acetic acid simulation, some sudden structural change in the PR conformation shifts it away from the native state and then the PR fluctuates between different non-native states.

To locate this structural transition, we first looked at the values of radius of gyration (RG) of the dimeric PR about its x, y and z Cartesian coordinate axes separately. These three different values give a global indication of the shape of the molecule and their fluctuations represent the extent of flexibility along the corresponding axis of the molecule. Figure 1(b) very clearly reveals that during the first 25 ns of simulation in 9 M AcOH the dimeric PR steadily elongates along the z axis [Figure S2] and so the two monomer units in the PR dimer should move away from each other during this initial phase of dimer denaturation in 9 M AcOH.

To confirm this movement, we next looked at the tertiary contact profile of the trajectory. A tertiary contact is defined between two residues when the non-bonded heavy atoms of the residues lie within 0.6 nm of each other [Figure S3]. Figures 2(a) and 2(b) show the tertiary contact profile of the initial structure and 25 ns snapshot from the 9 M AcOH trajectory of the dimeric PR. The latter figure has clear absence of any tertiary contact between the flaps (residues 45 to 55 in both monomeric units) of the dimeric PR. It also reveals reduced contacts between the catalytic triad (residues 25 to 27) of both the monomeric units; whereas the β-sheet formed by the N and C terminal of the two monomers seems to remain intact. Thus, the dimer seems to elongate in such a fashion that the monomer units separate at the flaps but remain tethered by the N and C terminal β-sheet [Figures 2(c) and 2(d)]. Such movement can be clearly visualized in the movie file [Video S1] provided as supporting video with this article.

The hydrogen bond network (Fireman's Grip) between the catalytic triad of both the monomers in the PR is crucial for its function and so its disruption could lead to loss of PR function. We looked at this hydrogen bond network during the simulation of dimeric PR in 9 M AcOH [Figure S4] and this data indicates that the hydrogen bond network completely breaks by 25 ns.

Thus, by summing up these observations from the simulation of dimeric PR in 9 M AcOH, it can be concluded that as early as within the first 25 ns of the simulation, the monomeric units of the PR start separating at the flaps but remain tethered at the terminal β-sheet. This separation leads to disruption of the Fireman's grip and the PR becomes dysfunctional. Even up till further 150 ns, neither the monomers move closer nor the Fireman's grip gets reconstructed. Hence, the dimeric PR can be considered to be denatured and it is only a matter of a little more time when the terminal β-sheet will break and the dimer will truly separate into two monomers. Since the monomer folding and dimerization is reported to be decoupled [18], it can also be assumed that during the time required for the dimer to separate into monomers, the individual monomers will not denature. Consequently, we proceeded with simulation of the native PR monomer under similar conditions to probe further into the denaturation process.

2. Denaturation of the monomeric PR

The RMSD for simulation of monomeric PR in 9 M AcOH and water is given in Figure 1(c). For the initial 4 ns, both simulations show similar RMSD trend; but after this time point, the RMSD of the trajectory in acetic acid shows large fluctuations whereas the RMSD of the trajectory in water converges. Thus, again some conformational transition around 4 ns seems to prevent the structure from equilibrating in the acetic acid simulation. This transition was deciphered from the tertiary contact profiles at this time point [Figure 3]. Although the initial tertiary contact profiles at 0 ns for both the simulations in water [Figure 3(a)] and 9 M AcOH [Figure 3(d)] are similar; as we reach 4 ns, the N and C termini of the PR monomer in water start forming contacts [Figure 3(b)] whereas in the profile for acetic acid simulation, these contacts are distinctly absent [Figure 3(e)]. Thus, the β-sheet structure formed by the N and C termini of the PR must confer some structural stability to the protein and in the absence of these interactions the monomer trajectory gets propelled into a very random non-equilibrium path on the conformation landscape.

Further, secondary structure calculation on the two trajectories [Figure 4(a) and 4(b)] reveals that even in the denaturing 9 M AcOH conditions, much residual secondary structure mainly as the β-sheet remains in PR. However, the C-terminal α-helix transits to a turn within 20 ns simulation of the PR in 9 M AcOH [Figure 3(d) to 3(i) and Figure 4(a)] but not in the water simulation [Figure 3(a) to 3(c) and Figure 4(b)] and neither during the full 150 ns of the dimer simulation in acetic acid [refer to the movie Video S1 and S2 provided as supporting videos with this article]. Gō model [2] simulations have indicated that the (post critical) folding nuclei of PR are or folding nucleus of PR is formed as a spatially closed unit of α helix (83–92) with sheet (74–78) above and another β-strand with sheet (24–34) perpendicular to these elements which essentially coincides with the highly protected structural units and the stabilization core defined by Wallqvist et al. [23]. Thus, the disruption of the α-helical structure and the consequent and persistent loss of the FN interactions cause the PR to denature. It is also interesting to note that lattice model calculations [24], [25] have shown a hierarchical model of PR monomer folding with docking of fragments (83–92) and (24–34) being the first step followed by relaxation of the remaining amino acids in the native structure shortly afterwards. In the unfolding process seen in our simulation, these fragments are also the first to dissociate and start the denaturation of the PR.

3. Insights into the mechanism of denaturation of PR by acetic acid

We next speculated about the driving force for such major structural transitions in the simulation of PR in acetic acid in such a short time period. Figure 5 shows the contact autocorrelation function of acetic acid and water with the PR during the simulation and Table 1 details the integrated correlation time calculated from these graphs. It is clear that acetate ions spend a lot of time on average around the PR, followed by neutral acetic acid molecules. Thus, the acetic acid species interacts with the PR more than water in the simulation.

Further, Figure 6 clearly illustrates the regions in PR where these molecules specifically reside. It clearly shows a continuous envelope of the solvent around the protein stretch from residues 85 to 95 which contain the alpha helix. This stretch also contains only three hydrophilic amino acids (R87, N88 and Q92) and one polar residue (T91) while all the remaining seven residues are hydrophobic. The acetic acid environment (with dielectric constant of 6.82) is amenable to these hydrophobic residues as well and they easily make contacts with the solvent envelope around them. Such new interactions overcome the scanty α-helical interactions easily, and residues 85 to 95 adopt random coil structure. In the highly polar environment of water (dielectric constant 80), this stretch of residues tries to minimize its interaction with the solvent envelope and thus maintains a compact alpha helix.

A similar argument could be extended to explain the stability of the β sheet core domain. The solvent molecules do not penetrate the core domain [Figure S5] to disrupt the hydrophobic interactions in this region and the solvent envelope outside the core domain does not furnish enough interactions to overcome both the extensive hydrogen bonding network and hydrophobic interactions in the β sheets of the core domain.

Thus, we have presented here the first direct visualization of the initial events in the denaturation of mature and functional PR in acetic acid using all atom Molecular Dynamics simulation in explicit solvent. Novel observations from these simulations have helped us to obtain a picture of the initiation of acetic acid denaturation of PR. The dimeric PR begins to denature by separating into monomers. Though the dimer interface and the critical hydrogen bond network of the catalytic triad of PR disrupt within first 25 ns of the simulation, the monomers remain tethered together by the β-sheet formed by the N and C termini for a longer duration of time. However, after separation into monomers, the folding nucleus formed by the α-helix and the surrounding β-strands is the first to disrupt during the denaturation of the monomer. Following the loss of these critical interactions, the PR follows a random and non-equilibrating path on the conformation landscape. Apart from this, we have also been able to look into the mechanism of denaturation by acetic acid and have attributed it to the direct interactions of acetic acid molecules with some sensitive regions of the protein such as the α-helical residues.

However, during the time duration monitored in this study, only some early events could be observed. Nevertheless, these provide satisfactory evidence of initiation of denaturation. Simulations up to microseconds may be required to visualize collapse of the robust β-sheet structure and complete denaturation of the protein. Such studies can also be extended to denaturation by other agents to obtain insights into the overall denaturation landscape of the protein. We believe that the observations reported in this study pinpoint the early and crucial events that can disrupt both the function and structure of the PR and hence can provide a basis for the design of newer targets for blocking PR function.

1. Construction of 9 M acetic acid model

Both uncharged and charged acetic acid molecules were constructed using ArgusLab software [] and their parameters were obtained from PRODRG server at Dundee [26] []. The uncharged molecule was duplicated and a box [5 Å×5 Å×5 Å] containing 334 molecules of acetic acid was constructed using Vega ZZ 2.3.2 [27] []. All the further steps in setting up of the model were then carried by GROMACS 4.0.5 [28], [29] tool and all MD simulations were performed in periodic boundary conditions using GROMACS 4.0.5 tool with GROMOS96 force field [30]. The acetic acid box was solvated with SPC216 water molecules [31]. The ratio of number of acetic acid molecules to that of water in the solvated box was adjusted to 1∶3 and the density was adjusted to 1.02 kg l−1 such that this box models 9 M AcOH. Since the force field parameters for acetic acid molecules obtained in the GROMACS topology from PRODRG server and those for the SPC216 water molecules are already validated, we used these same validated force field parameters for the simulation of the 9 M AcOH solvent. After an initial energy minimization by steepest descent algorithm with tolerance of 100 kJ mol−1 nm−1 to remove any bad contacts, an MD run was set up for 1 ns to completely mix and equilibrate the contents of the box. Complete mixing and equilibration of the contents was confirmed by convergence of potential energy of the system and by convergence of the number of hydrogen bonds between acetic acid and water molecules using GROMACS 4.0.5 analysis tools.

2. Model preparation and MD setup for studying denaturation of PR

The crystal structure of mature PR dimer reported by Pillai et al. [32] (PDB ID 1G6L) and the solution structure of PR monomer reported by Ishima et al. [33] (PDB ID 1Q9P) was adopted as the starting structure. These were separately placed in a cubic box and solvated with the equilibrated 9 M AcOH solution and SPC216 water. The peripheral aspartates, glutamates, histidines and C-terminal of the PR were protonated to model the effect of protonation of the amino acids in low pH conditions (∼1.9) in 9 M AcOH and thirteen random acetic acid molecules were replaced with acetate ions (charged acetic acid molecules) to neutralize the additional charge gained by the system.

All models were energy minimized using steepest descent until they converged to 10 kJ mol−1 nm−1. A 50 ps MD dynamics run was then performed at 300 K with position restraint on the PR atoms. Finally a 150 ns MD was performed on the system of dimeric PR in 9 M AcOH box and a 30 ns MD was performed on the system of monomeric PR in 9 M AcOH and on monomeric and dimeric PR in water. All trajectory analyses were carried out by the analysis tools in GROMACS 4.0.5 package and VMD [34] [] was used for visualization purposes.

Supporting Information Figure S1

The root mean square deviation (RMSD) of the backbone atoms in the trajectory of the dimer in water from the backbone atoms of the NMR structure of PR.


Click here for additional data file (pone.0019830.s001.png)

Figure S2

The definition of the x, y and z axis of the dimeric PR.


Click here for additional data file (pone.0019830.s002.png)

Figure S3

The tertiary contact profiles for a PR dimer. The whole profile can be broken into four quadrants; the lower left and upper right quadrants represent the tertiary contacts within each monomer. The remaining two quadrants contain the interactions between the monomers. The boxed-off values are interactions between the different β-strands; circled values are the α-helical interactions; the values enclosed by a triangle represent interactions between the loops and β-strands; those enclosed by a box with rounded edge represent interactions between the loops and residues of the α-helix and arrows indicate the interactions at the dimer interface.


Click here for additional data file (pone.0019830.s003.png)

Figure S4

The hydrogen bond profile for the catalytic triad during the simulation of dimeric PR in 9 M AcOH. In the panels for monomer units, the first and second H-bonds form a part of the Fireman's grip whereas in the panel for dimer, the 2nd, 5th, 11th and 12th H-bonds form a part of the Fireman's grip. This profile has been generated by keeping H-bond donor-hydrogen distance cutoff of 0.3 nm and angle cutoff of 120°.. Thus, even with these lenient cutoff values, it is clear that only 2 of the eight H-bonds persist continuously throughout the trajectory. The 4 H-bonds in the monomer units exist intermittently whereas the 5th and the 12th H-bond in the dimer panel completely vanishes by 20 ns. Thus, by this time, the Fireman's grip can be considered to be broken and it doesn't get restored in the remainder of the simulation.


Click here for additional data file (pone.0019830.s004.png)

Figure S5

Envelope of solvent 9 M AcOH surrounding the core domain of PR. The envelope is drawn as a surface comprising of all solvent atoms that are within 3.5 Å of the residues.


Click here for additional data file (pone.0019830.s005.png)

Video S1

(9MAcOH-150 ns-PRdimer.wmv): The movie file for the 150 ns trajectory of dimeric PR in 9 M AcOH. The movie contains snapshot of every 300 ps of the simulation.


Click here for additional data file (pone.0019830.s006.wmv)

Video S2

(9MAcOH-30 ns-PRmonomer.wmv): The movie file for the 30 ns trajectory of monomeric PR in 9 M AcOH. The movie contains snapshot of every 100 ps of the simulation.


Click here for additional data file (pone.0019830.s007.wmv)


Competing Interests: The authors have declared that no competing interests exist.

Funding: Department of Atomic Energy (DAE) and Department of Science and Technology (DST), Government of India. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

1. Tiana G,Broglia RA. Year: 2002Folding and design of dimeric proteins.Proteins: Struct Funct Bioinf498294
2. Broglia RA,Tiana G,Sutto L,Provasi D,Simona F. Year: 2005Design of HIV-1 PR inhibitors that do not create resistance: blocking the folding of single monomers.Protein Sci142668268116195553
3. Yi Q,Scalley-Kim ML,Alm E J,Baker D. Year: 2000NMR characterization of residual structure in the denatured state of protein L.J Mol Biol2991341135110873457
4. Mok YK,Kay CM,Kay LE,Forman-Kay J. Year: 1999NOE data demonstrating a compact unfolded state for an SH3 domain under non-denaturing conditions.J Mol Biol28961963810356333
5. Zhang O,Forman-Kay JD. Year: 1995Structural characterization of folded and unfolded states of an SH3 domain in equilibrium in aqueous buffer.Biochemistry34678467947756310
6. Navon A,Ittah V,Landsman P,Scheraga HA,Haas E. Year: 2001Distributions of intramolecular distances in the reduced and denatured states of bovine pancreatic ribonuclease A. Folding initiation structures in the C-terminal portions of the reduced proteins.Biochemistry4010511811141061
7. Bhavesh NS,Sinha R,Mohan PM,Hosur RV. Year: 2003NMR elucidation of early folding hierarchy in HIV-1 protease.J Biol Chem27822199801998512644464
8. Rout MK,Hosur RV. Year: 2009Fluctuating partially native-like topologies in the acid denatured ensemble of autolysis resistant HIV-1 protease.Arch Biochem Biophys482334119100236
9. Chatterjee A,Mridula P,Mishra R,Mittal R,Hosur RV. Year: 2005Folding regulates the autoprocessing of HIV-1 protease precursor.J Biol Chem28012113691137815632156
10. Bhavesh NS,Panchal SC,Mittal R,Hosur RV. Year: 2001NMR identification of local structural preferences in HIV-1 protease tethered heterodimer in 6 M guanidine hydrochloride.FEBS Lett509221822411741592
11. Harrison RW,Weber IT. Year: 1994Molecular dynamics simulations of HIV-1 protease with peptide substrate.Protein Eng711135313637700867
12. Trylska J,Tozzini V,Chang CA,McCammon JA. Year: 2007HIV-1 Protease Substrate Binding and Product Release Pathways Explored with Coarse-Grained Molecular Dynamics.Biophys J92124179418717384072
13. Zhu Z,Schuster DI,Tuckerman ME. Year: 2003Molecular Dynamics Study of the Connection between Flap Closing and Binding of Fullerene-Based Inhibitors of the HIV-1 Protease.Biochemistry4251326133312564936
14. Chen X,Weber IT,Harrison RW. Year: 2004Molecular dynamics simulations of 14 HIV protease mutants in complexes with indinavir.J Mol Model105–637338115597206
15. Wartha F,Horn AHC,Meiselbach H,Sticht H. Year: 2005Molecular Dynamics Simulations of HIV-1 Protease Suggest Different Mechanisms Contributing to Drug Resistance.J Chem Theory Comput12315324
16. Perryman AL,Lin JH,McCammon JA. Year: 2004HIV-1 protease molecular dynamics of a wild type and of the V82F/I84V mutant: possible contributions to drug resistance and a potential new target site for drugs.Protein Sci13411082315044738
17. Levy Y,Caflisch A. Year: 2003Flexibility of Monomeric and Dimeric HIV-1 Protease.J Phys Chem10730683079
18. Levy Y,Caflisch A,Onuchic JN,Wolynes PG. Year: 2004The folding and dimerization of HIV-1 Protease: Evidence for a stable monomer from simulations.J Mol Biol340677915184023
19. Yan M,Sha Y,Wang J,Xiong X,Ren J,Cheng M. Year: 2008Molecular dynamics simulations of HIV-1 protease monomer: Assembly of N-terminus and C-terminus into β-sheet in water solution.Proteins: Struct Funct Bioinf70731738
20. Hornak V,Okur A,Rizzo R,Simmerling C. Year: 2006HIV-1 protease flaps spontaneously open and reclose in molecular dynamics simulations.Proc Natl Acad Sci USA10391592016418268
21. Collins JR,Burt SK,Erickson JW. Year: 1995Flap opening in HIV-1 protease simulated by ‘activated’ molecular dynamics.Nat Struct Biol23343387796268
22. Kabsch W,Sander C. Year: 1983Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features.Biopolymers221225776376667333
23. Wallqvist A,Smythers GW,Covell DG. Year: 1998A cooperative folding unit in HIV-1 protease, Implications for protein stability and occurrence of drug-induced mutations.Protein Eng1199910059876920
24. Broglia RA,Tiana G. Year: 2001aHierarchy of events in the folding of model proteins.J Chem Phys11472677273
25. Broglia RA,Tiana G. Year: 2001bReading the three-dimensional structure of latticed model-designed proteins from their amino acid sequence.Proteins: Struct Funct Bioinf45421427
26. Schuettelkopf AW,van Aalten DMF. Year: 2004PRODRG - a tool for high-throughput crystallography of protein-ligand complexes.Acta CrystallographicaD6013551363
27. Pedretti A,Villa L,Vistoli G. Year: 2004VEGA - An open platform to develop chemo-bio-informatics applications using plug-in architecture and script programming.J Comp Aid Mol Des183167173
28. Berendsen HJC,van der Spoel D,van Drunen R. Year: 1995GROMACS: a message-passing parallel molecular dynamics implementation.Comp Phys Comm914356
29. Lindahl E,Hess B,van der Spoel D. Year: 2001GROMACS 3.0: a package for molecular simulation and trajectory analysis.J Mol Mod7306317
30. Schuler LD,Daura X,van Gunsteren WF. Year: 2001An improved GROMOS96 force field for aliphatic hydrocarbons in the condensed phase.J Comput Chem2212051218
31. Berendsen HJC,Postma JPM,van Gunsteren WF,Hermans J. Year: 1981Interaction models for water in relation to protein hydration in Intermolecular forces,Pullman B331342 (Dordrecht The Netherlands: Reidel).
32. Pillai B,Kannan KK,Hosur MV. Year: 20011.9 Å x-ray study shows closed flap conformation in crystals of tethered HIV-1 PR.Proteins: Struct Funct Bioinf4357564
33. Ishima R,Torchia DA,Lynch SM,Gronenborn AM,Louis JM. Year: 2003Solution structure of the mature HIV-1 protease monomer: Insight into the tertiary fold and stability of a precursor.J Biol Chem278433114331912933791
34. Humphrey W,Dalke A,Schulten K. Year: 1996VMD - Visual Molecular Dynamics.J Molec Graphics143338


[Figure ID: pone-0019830-g001]
doi: 10.1371/journal.pone.0019830.g001.
Figure 1  The root mean square deviation (RMSD) of the backbone atoms in the trajectory of the (a) dimer and (c) monomer in 9 M AcOH from the backbone atoms of the NMR structure of PR and (b) the radii of gyration of PR separately about the x, y and z axes during the 150 ns simulation of dimeric PR in 9 M AcOH.

The RMSD in 9 M AcOH does not equilibrate and fluctuates more than 10 Å indicating a non-equilibrium run.

[Figure ID: pone-0019830-g002]
doi: 10.1371/journal.pone.0019830.g002.
Figure 2  Tertiary contacts at (a) 0 ns and (b) 25 ns and snapshot structures at the same time points i.e. at (c) 0 ns and (d) 25 ns of the simulation of dimeric PR in 9 M AcOH.

The circled region denote the tertiary contacts between the protein flaps in the 0 ns structure and complete absence of these contacts in the 25 ns structure. Corresponding snapshot figures illustrate the marked separation of the protein flaps and the catalytic triads of the monomeric units. For detailed description of the tertiary contact profiles, refer to supplementary material S4.

[Figure ID: pone-0019830-g003]
doi: 10.1371/journal.pone.0019830.g003.
Figure 3  Tertiary contacts between the residues of dimeric PR simulated at (a) 0 ns, (b) 4 ns and (c) 25 ns in water and the corresponding tertiary contacts (d), (e) and (f) and trajectory snapshots (g), (h) and (i) in 9 M AcOH.

Arrows indicate the interactions of the N and C termini that are seen in the simulation in water but not in the simulation in 9 M AcOH. The circled region indicates the α-helical interactions that completely vanish by 25 ns in the monomer simulation in 9 M AcOH.

[Figure ID: pone-0019830-g004]
doi: 10.1371/journal.pone.0019830.g004.
Figure 4  The secondary structure profile calculated using DSSP [22] for the initial 30 ns of the two simulations of monomeric PR a) in 9 M AcOH and b) in water.

[Figure ID: pone-0019830-g005]
doi: 10.1371/journal.pone.0019830.g005.
Figure 5  Contact autocorrelation function of acetate ions, acetic acid and water during the simulation of the PR monomer.

[Figure ID: pone-0019830-g006]
doi: 10.1371/journal.pone.0019830.g006.
Figure 6  Envelope of solvent (a) water and (b) 9 M AcOH surrounding residues 85 to 95 of PR.

This stretch contains the α-helix seen as purple ribbon in the figure. The envelope is drawn as a surface comprising of all solvent atoms that are within 3.5 Å of the residues 85 to 95.

[TableWrap ID: pone-0019830-t001] doi: 10.1371/journal.pone.0019830.t001.
Table 1  Mean residence time of acetate ions, acetic acid and water around the PR calculated by integrating the contact autocorrelation function for these chemical species during the simulation of the PR monomer in 9 M AcOH.
S. No. Chemical Species Mean Residence Time
1. Acetate Ions 2158 ps
2. Neutral Acetate Acid Molecules 373 ps
3. Water Molecules 156 ps

Article Categories:
  • Research Article
Article Categories:
  • Biology
    • Computational Biology
      • Macromolecular Structure Analysis
        • Protein Structure

Previous Document:  Comparative structural analysis of lipid binding START domains.
Next Document:  From SNPs to genes: disease association at the gene level.