Denaturation of HIV-1 protease (PR) monomer by acetic acid: mechanistic and trajectory insights from molecular dynamics simulations and NMR

Inside a living cell there can be a variety of interactions for any given protein, which serve to regulate denaturation and renaturation processes. Insights into some of them can be obtained by in vitro studies using various denaturing agents. In this study, all-atom MD simulations in explicit solvent and NMR relaxation studies were performed on HIV-1 Protease (PR) in 9 M acetic acid (AcOH) (the commonly used denaturant during PR preparation). Following previous reports that denaturation proceeds via dissociation of the dimer into monomers, unfolding of the monomer by acetic acid has been explicitly investigated here. Direct visualization of the denaturation process and evidence for the mechanism of denaturation have been presented. Our simulations reveal that the denaturation of the PR monomer is caused due to direct interaction between acetic acid molecules and PR. Autocorrelation of N-H vectors calculated from the simulations have revealed that the α-helix and the surrounding β-strands represent the sensitive regions of the PR that respond maximally to the change in the solvent environment around the PR and are prone to disruption by acetic acid. This disruption is caused due to increased penetration of the acetic acid molecules into the PR structure by formation of preferred tertiary contacts and hydrogen bonds between the PR and acetic acid molecules. Following the loss of these critical interactions, the PR follows a random and non-equilibrating path on the conformation landscape and cycles between different denatured extended and compact states.     

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.     

The Sequential Mechanism of Guanidine Hydrochloride-Induced Denaturation of cAMP Receptor Protein from Escherichia coli. A Fluorescent Study Using 8-Anilino-1-Naphthalenesulfonic Acid

cAMP receptor protein (CRP) regulates expression of a number of genes in Escherichia coli. The protein is a homodimer and each monomer is folded into two structural domains. The biological activation of CRP upon cAMP binding may involve the subunit realignment as well as reorientation between the domains within each subunit. In order to study the interactions between the subunits or domains, we performed stopped-flow measurements of the guanidine hydrochloride (GuHCl)-induced denaturation of CRP. The changes in CRP structure induced by GuHCl were monitored using both intrinsic Trp fluorescence as well as the fluorescence of an extrinsic probe, 8-anilino-1-Naphthalenesulfonic acid (ANS). Results of CRP denaturation using Trp fluorescence detection are consistent with a two-step model [Malecki, and Wasylewski, (1997), Eur. J. Biochem. 243, 660], where the dissociation of dimer into subunits is followed by the monomer unfolding. The denaturation of CRP monitored by ANS fluorescence reveals the existence of two additional processes. One occurs before the dissociation of CRP into subunits, whereas the second takes place after the dissociation, but prior to proper subunit unfolding. These additional processes suggest that CRP denaturation is described by a more complicated mechanism than a simple three-state equilibrium and may involve additional changes in both inter- and intrasubunit interactions. We also report the effect of cAMP on the kinetics of CRP subunit unfolding and refolding.     

The Sequential Mechanism of Guanidine Hydrochloride-Induced Denaturation of cAMP Receptor Protein from Escherichia coli. A Fluorescent Study Using 8-Anilino-1-Naphthalenesulfonic Acid

cAMP receptor protein (CRP) regulates expression of a number of genes in Escherichia coli. The protein is a homodimer and each monomer is folded into two structural domains. The biological activation of CRP upon cAMP binding may involve the subunit realignment as well as reorientation between the domains within each subunit. In order to study the interactions between the subunits or domains, we performed stopped-flow measurements of the guanidine hydrochloride (GuHCl)-induced denaturation of CRP. The changes in CRP structure induced by GuHCl were monitored using both intrinsic Trp fluorescence as well as the fluorescence of an extrinsic probe, 8-anilino-1-Naphthalenesulfonic acid (ANS). Results of CRP denaturation using Trp fluorescence detection are consistent with a two-step model [Malecki, and Wasylewski, (1997), Eur. J. Biochem. 243, 660], where the dissociation of dimer into subunits is followed by the monomer unfolding. The denaturation of CRP monitored by ANS fluorescence reveals the existence of two additional processes. One occurs before the dissociation of CRP into subunits, whereas the second takes place after the dissociation, but prior to proper subunit unfolding. These additional processes suggest that CRP denaturation is described by a more complicated mechanism than a simple three-state equilibrium and may involve additional changes in both inter- and intrasubunit interactions. We also report the effect of cAMP on the kinetics of CRP subunit unfolding and refolding.     

Soni-removal of nucleic acids from inclusion bodies

Inclusion bodies (IBs) are commonly formed in Escherichiacoli due to over expression of recombinant proteins in non-native state. Isolation, denaturation and refolding of these IBs is generally performed to obtain functional protein. However, during this process IBs tend to form non-specific interactions with sheared nucleic acids from the genome, thus getting carried over into downstream processes. This may hinder the refolding of IBs into their native state. To circumvent this, we demonstrate a methodology termed soni-removal which involves disruption of nucleic acid–inclusion body interaction using sonication; followed by solvent based separation. As opposed to conventional techniques that use enzymes and column-based separations, soni-removal is a cost effective alternative for complete elimination of buried and/or strongly bound short nucleic acid contaminants from IBs.     

Synthesis of A New Intercalating Nucleic Acid 6H-INDOLO[2,3-b] Quinoxaline Oligonucleotides to Improve Thermal Stability Of Hoogsteen-Type Triplexes

A new intercalating nucleic acid monomer X was obtained in high yield starting from alkylation of 4-iodophenol with (S)-(+)-2-(2,2-dimethyl-1,3-dioxolan-4-yl)ethanol under Mitsunobu conditions followed by hydrolysis with 80% aqueous acetic acid to give a diol which was coupled under Sonogashira conditions with trimethylsilylacetylene (TMSA) to achieve the TMS protected (S)-4-(4-((trimethylsilyl)ethynyl)phenoxy)butane-1,2-diol. Tetrabutylammonium flouride was used to remove the silyl protecting group to obtain (S)-4-(4-ethynylphenoxy)butane-1,2-diol which was coupled under Sonogashira conditions with 2-(9-bromo-6H-indolo[2,3-b]quinoxalin-6-yl)-N,N-dimethylethanamine to achieve (S)-4-(4-((6-(2-(dimethylamino)ethyl)-6H-indolo[2,3-b]quinoxalin-9-yl)ethynyl)phenoxy)butane-1,2-diol. This compound was tritylated with 4,4′-dimethoxytrityl chloride followed by treatment with 2-cyanoethyltetraisopropylphosphordiamidite in the presence of N,N′-diisopropyl ammonium tetrazolide to afford the corresponding phosphoramidite. This phosphoramidite was used to insert the monomer X into an oligonucleotide which was used for thermal denaturation studies of a corresponding parallel triplex.     

Effects of hesperidin,a flavanone glycoside interaction on the conformation,stability, and aggregation of lysozyme: multispectroscopic and molecular dynamic simulation studies?

Hesperidin (HESP), a flavanone glycoside, shows high antioxidant properties and posses ability to go through the blood–brain barrier. Therefore, it could be a potential drug molecule against aggregation based diseases such as Alzheimer’s, Parkinson’s, and systemic amyloidoses. In this work, we investigated the potential of HESP to interact with hen egg-white lysozyme (HEWL) monomer and prevent its aggregation. The HESP–HEWL binding studies were performed using a fluorescence quenching technique, molecular docking and molecular dynamics simulations. We found a strong interaction of HESP with the lysozyme monomer (K_a, ~ 5 × 10⁴ M^?1) mainly through hydrogen bonding, water bridges, and hydrophobic interactions. We showed that HESP molecule spanned the highly aggregation prone region (amino acid residues 48-101) of HEWL and prevented its fibrillar aggregation. Further, we found that HESP binding completely inhibited amorphous aggregation of the protein induced by disulfide-reducing agent tries-(2-carboxyethyl) phosphine. Conformational and stability studies as followed by various tertiary and secondary structure probes revealed that HESP binding only marginally affected the lysozyme monomer conformation and increased both stability and reversibility of the protein against thermal denaturation. Future studies should investigate detail effects of HESP on solvent dynamics, structure, and toxicity of various aggregates. The answers to these questions will not only target the basic sciences, but also have application in biomedical and biotechnological sciences.     

The role of ion-lipid interactions and lipid packing in transient defects caused by phenolic compounds

The transient disruption of membranes for the passive permeation of ions or small molecules is a complex process relevant to understanding physiological processes and biotechnology applications. Phenolic compounds are widely studied for their antioxidant and antimicrobial properties, and some of these activities are based on the interactions of the phenolic compound with membranes. Ions are ubiquitous in cells and are known to alter the structure of phospholipid bilayers. Yet, ion-lipid interactions are usually ignored when studying the membrane-altering properties of phenolic compounds. This study aims to assess the role of Ca²⁺ ions on the membrane-disrupting activity of two phenolic acids and to highlight the role of local changes in lipid packing in forming transient defects or pores. Results from tethered bilayer lipid membrane electrical impedance spectroscopy experiments showed that Ca²⁺ significantly reduces membrane disruption by caffeic acid methyl ester and caffeic acid. As phenolic acids are known metal chelators, we used UV-vis and fluorescence spectroscopy to exclude the possibility that Ca²⁺ interferes with membrane disruption by binding to the phenolic compound and subsequently preventing membrane binding. Molecular dynamics simulations showed that Ca²⁺ but not caffeic acid methyl ester or caffeic acid increases lipid packing in POPC bilayers. The combined data confirm that Ca²⁺ reduces the membrane-disrupting activity of the phenolic compounds, and that Ca²⁺-induced changes to lipid packing govern this effect. We discuss our data in the context of ion-induced pores and transient defects and how lipid packing affects membrane disruption by small molecules.     

The Effect of Acetic Acid on the Association and the State of the Tyrosine Residue of an Active Fragment of Bovine Growth Hormone

Self-association and the state of a tyrosine residue of the active fragment of bovine growth hormone were investigated by gel filtration and difference spectroscopy with various concentrations of acetic acid. In acetic acid solutions more concentrated than 1.5 m, the active fragment was found to exist as a monomeric form. Decrease in the concentration of acetic acid promoted self-association of the active fragment. Dimeric and trimeric forms were observed in 0.1 m acetic acid. By difference spectroscopy, the tyrosine residue in the active fragment showed a typical blue shift in 2.5 M acetic acid compared to 0.5 M acetic acid. The transition concentration, 1.25 M of acetic acid for the association-dissociation of the fragment was consistent with that for the appearance of the blue shift difference spectrum. It was concluded that a tyrosine residue in the active fragment was exposed on the surface of the fragment molecule in the monomer form and was held between the fragment molecules under conditions in which the active fragment forms dimer and trimer.     

Quantum-chemical and empirical calculations on phospholipids IV. Glycero-phosphorylethanolamine in a quasihexagonal planar lattice

The energetically most stable head group conformations of a racemic mixture of diacyl-glycero-phosphorylethanolamine in a planar quasihexagonal lattice were calculated using empirical 1-6-12 atom-atom potential functions for intra- and intermolecular interactions. The results demonstrate that the conformation of phospholipid head groups in bilayer systems is determined by intramolecular interactions as well as by intermolecular interactions with neighbouring phospholipid molecules and with solvent molecules. The most stable conformers are that with a φ2 = guache^? conformation of the phosphodiester group. All conformers with a φ2 = gauche⁽⁺⁾ or trans conformation have energies more than 15 kcal ☆ mol^?1 above that of the global minimum. The calculated torsional angles ?1 and φ1 are in very good agreement with the results of the X-ray diffraction analysis of 1,2-dilauroyl-DL-phosphatidylethanolamine (DLPE) acetic acid single crystals.     

Monte-Carlo Calculations for Methane and Argon over a Wide Range of Density and Temperature,Including the Two-Phase Vapor-Liquid Region

Abstract

The Gibbs-ensemble simulation technique provides a powerful method to calculate vapor-liquid phase behavior [1]. To evaluate the configurational energy of a system of molecules, commonly used experessions describe the interaction between two molecules. Contributions from higher-body forces are usually implicitly taken into account by adjusting two-body potential parameters to give agreement with experimental data. Explicit expressions for higher-body potentials are not commonly used in simulations [8]. The work by Smit et al. [9] gives the appropriate expressions to evaluate the pressure as well as the chemical potential from a density-dependent two-body potential in an NVT ensemble.

In the present work, contributions to the potential from two-body interactions are separated from those due to higher-body interactions; to take higher-body forces into account, a mean-field term, proportional to (density)^0.9, is added to the two-body potential. NPT-simulations over a wide range of temperature and density, as well as Gibbs-ensemble simulations, are used to evaluate phase behavior of argon and of methane. The results indicate that a simple mean-field correction to the “true” two-body Kihara potential provides good agreement between experiment and simulation.     

X-ray Studies on Crystalline Complexes Involving Amino Acids and Peptides. XXVII. Effect of Chirality,Specific Interactions and Characteristic Aggregation Patterns in the Structures of Arginine and its Complexes with Formic Acid

Abstract

The crystal structures of DL- arginine dihydrate, DL-arginine formate dihydrate and L-arginine formate have been determined and refined using X-ray crystallographic techniques. The three structures, along with other related ones, demonstrate the conformational variability of arginine. The amino acid molecules aggregate essentially in a similar manner in DL- arginine dihydrate and in the known structure of L-arginine dihydrate; the effects arising out of the reversal of the chirality of half the amino acid molecules are absorbed by small local adjustments. However, such a reversal leads to profound differences in aggregation in DL- arginine and L- arginine formates, in contrast to the situation in the corresponding acetates. Thus the effect of chirality on biomolecular aggregation cannot be easily predicted or even rationalized. Arginine-carboxylate interactions in the complexes primarily involve the guanidyl groups and contain specific interactions. Indeed the primary mode of arginine-carboxylic acid aggregation is substantially invariant in the arginine complexes of succinic, acetic and formic acids.     

A CFF 91-based Continuum Solvation Model: Solvation Free Energies of Small Organic Molecules and Conformations of the Alanine Dipeptide in Solution

Abstract

For molecular mechanics simulations of solvated molecules, it is important to use a consistent approach for calculating both the force field energy and the solvation free energy. A continuum solvation model based upon the atomic charges provided with the CFF91 force field is derived. The electrostatic component of the solvation free energy is described by the Poisson-Bolzmann equation while the nonpolar comonent of the solvation energy is assumed to be proportional to the solvent accessible surface area of the solute. Solute atomic radii used to describe the interface between the solute and solvent are fitted to reproduce the energies of small organic molecules. Data for 140 compounds are presented and compared to experiment and to the results from the well-characterized quantum mechanical solvation model AM1-SM2. In particular, accurate results are obtained for amino acid neutral analogues (mean unsigned error of 0.3 kcal/mol). The conformational energetics of the solvated alanine dipeptide is discussed.     

Computer Aided Study of Ligand Binding With Catalytic Domain of Avian Sarcoma Virus Integrase and Its Ligand Binding Loops

Abstract

We have carried out 1 nanosecond (ns) Molecular Dynamics (MD) simulations of the drug Y3 (4-acetylamino-5-hydroxynaphthalene-2, 7-disulfonic acid) complexed with catalytic domain of Avian sarcoma virus Integrase (ASV-IN), both in vacuum and in the presence of explicit solvent. Starting models were obtained on the basis of PDB co-ordinates (1A5X) of ASV-IN-Y3 complex, by Lubkowski et al [1]. Mn²⁺ cation was present in the active site. To neutralize the positive charge in the presence of explicit solvent, eight Cl^? anions were added. Energy Minimization (EM) and MD simulations, for both the systems, were carried out using Sander's module of AMBER5.0 [2] with all atom force field. Analysis of ligand- protein interaction in both environments is discussed in the paper. We also carried out 1 ns MD simulation on two flexible loops—L1 (Gly54-Gln62) and L2 (Trp138-Met155) playing crucial role in interaction of IN with the drug [3], under differing environmental conditions (vacuum, aqueous and organic solvent methanol). Comparison of the conformational changes in the loops, monomer and dimer is presented in the paper. Our results showed that the conformation of the loop region was closest to crystallographic data in case of monomer and constrained loops in aqueous environment. However, the dimer in vacuum was more stable than monomer. The β sheet structure of the monomer in aqueous environment was unstable. Latter also took long time for equilibration. The box formed by loops L1 and L2 from two sub units IINA and INB) of the dimer satisfies prerequisites for ligand recognition site and seems to be the functional biological unit.     

Theory of cooperative transitions in protein molecules. I. Why denaturation of globular protein is a first-order phase transition

A theory of equilibrium denaturation of proteins is suggested. According to this theory, a cornerstone of protein denaturation is disruption of tight packing of side chains in protein core. Investigation of this disruption is the object of this paper. It is shown that this disruption is an "all-or-none" transition (independent of how compact is the denatured state of a protein and independent of the protein-solvent interactions) because expansion of a globule must exceed some threshold to release rotational isomerization of side chains. Smaller expansion cannot produce entropy compensation of nonbonded energy loss; this is the origin of a free-energy barrier (transition state) between the native and denatured states. The density of the transition state is so high that the solvent cannot penetrate into protein in this state. The results obtained in this paper make it possible to present in the following paper a general phase diagram of protein molecule in solution.     

Molecular dynamic simulation and DFT study on the Drug-DNA interaction; Crocetin as an anti-cancer and DNA nanostructure model

In this research, the interaction of Crocetin as an anti-cancer drug and a Dickerson DNA has been investigated. 25 ns molecular dynamic simulations of Crocetin and DNA composed of 12 base pairs and a sequence of d(CGCGAATTCGCG)₂ were done in water. Three definite parts of the B-DNA were considered in analyzing the best interactive site from the thermodynamic point of view. Binding energy analysis showed that van der Waals interaction is the most important part related to the reciprocal O and H atoms of the Crocetin and DNA. Stabilizing interactions, obtained by ΔG calculations, showed that maximum and minimum interactions are related to the S1 and S3 regions, respectively. This means that the most probable van der Waals interaction site of the Dickerson B-DNA and Crocetin is located in the minor groove of DNA. Two sharp peaks at 2.55 and 1.75 Å in radial distribution functions of the PO?HO and NH?OC parts are related to new hydrogen bonds between the Crocetin and DNA in the complex which can be considered as the driving force of the anti-cancer mechanism of the Crocetin. Average values of 0.3 au and zero for the electron densities of the H?O bonds for DNA and complex, obtained by Quantum theory of atoms in molecules (QTAIM), means that the origin of DNA instability after complexation may be related to the H-bond denaturation by Crocetin. Finally, the evaluation of the dispersion interactions using the dispersion functional, -148.76 kcal.mol^?1, confirmed the importance of the dispersion interaction in drug-DNA complex.     

Characterization of the denaturation of human alpha-lactalbumin in urea by molecular dynamics simulations

Molecular dynamics (MD) simulations were used to characterize the non-cooperative denaturation of the molten globule A-state of human alpha-lactalbumin by urea. A solvent of explicit urea and water molecules was used, corresponding to a urea concentration of approximately 6M. Three simulations were performed at temperatures of 293K, 360K and 400K, with lengths of 2 ns, 8 ns and 8 ns respectively. The results of the simulations were compared with experimental data from NMR studies of human alpha-lactalbumin and related peptides. During the simulations, hydrogen bonds were formed from the protein to both urea and water molecules as intra-protein hydrogen bonds were lost. Urea was shown to compete efficiently with water as both a hydrogen bond donor and acceptor. Radial distribution functions of water and urea around hydrophobic side chain atoms showed a significant increase in urea molecules in the solvation shell as the side chains became exposed during denaturation. A considerable portion of the native-like secondary structure persisted throughout the simulations. However, in the simulations at 360K and 400K, there were substantial changes in the packing of aromatic and other hydrophobic side chains in the protein, and many native contacts were lost. The results suggest that during the non-cooperative denaturation of the molten globule, secondary structure elements are stabilized by non-specific, non-native interactions.     

Molecular Simulations of Membrane Based Separations of Supercritical Electrolyte Solutions

Abstract

Computer simulations of solutions of electrolytes (NaCl and KCl) in supercritical water undergoing membrane based separations have been carried out. These studies used a technique developed recently, in which the system is maintained at steady state by periodically recycling the solvent molecules that permeated the membrane. Our results showed that ionic clusters, formed as a result of water molecules surrounding the ions, play a significant role in these separations. The effect of the main osmotic driving forces, such as pressure, temperature, concentration, and electric fields on the rate of permeation across the membrane was studied. In addition, we also looked at the effect of changes in the pore size and the attractive force between the membrane and solvent/solute. Finally, we examined the stability of the ionic clusters in these simulations.     

Solvent Interactions and Protein Dynamics in Spin-labeled T4 Lysozyme

Abstract

Aspects of T4 lysozyme dynamics and solvent interaction are investigated using atomically detailed Molecular Dynamics (MD) simulations. Two spin-labeled mutants of T4 lysozyme are analyzed (T4L-N40C and T4L-K48C), which have been found from electronic paramagnetic resonance (EPR) experiments to exhibit different mobilities at the site of spin probe attachment (N- and C-terminus of helix B, respectively). Similarities and differences in solvent distribution and diffusion around the spin label, as well as around exposed and buried residues within the protein, are discussed. The purpose is to capture possible strong interactions between the spin label (ring) and solvent molecules, which may affect EPR lineshapes. The effect of backbone motions on the water density profiles is also investigated. The focus is on the domain closure associated with the T4 lysozyme hinge-bending motion, which is analyzed by Essential Dynamics (ED). The N-terminus of helix B is found to be a “hinge” residue, which explains the high degree of flexibility and motional freedom at this site.