Membrane curvature and surface area per lipid affect the conformation and oligomeric state of HIV-1 fusion peptide: a combined FTIR and MD simulation study

Fourier-transformed infrared spectroscopy (FTIR) and molecular dynamics (MD) simulation results are presented to support our hypothesis that the conformation and the oligomeric state of the HIV-1 gp41 fusion domain or fusion peptide (gp41-FP) are determined by the membrane surface area per lipid (APL), which is affected by the membrane curvature. FTIR of the gp41-FP in the Aerosol-OT (AOT) reversed micellar system showed that as APL decreases from approximately 50 to 35 A2 by varying the AOT/water ratio, the FP changes from the monomeric alpha-helical to the oligomeric beta-sheet structure. MD simulations in POPE lipid bilayer systems showed that as the APL decreases by applying a negative surface tension, helical monomers start to unfold into turn-like structures. Furthermore, an increase in the applied lateral pressure during nonequilibrium MD simulations favored the formation of beta-sheet structure. These results provide better insight into the relationship between the structures of the gp41-FP and the membrane, which is essential in understanding the membrane fusion process. The implication of the results of this work on what is the fusogenic structure of the HIV-1 FP is discussed.     

Simulation of the N-terminus of HIV-1 glycoprotein 41000 fusion peptide in micelles.

In this paper, the N-terminus of glycoprotein-41, the HIV-1 fusion peptide, was studied by molecular dynamics simulations in an explicit sodium dodecyl sulfate micelle. The simulation provides a detailed picture of the equilibrium structure and peptide stability as it interacts with the micelle. The equilibrium location of the peptide shows the peptide at the surface of the micelle with hydrophobic residues interacting with the micelle's core. At equilibrium, the peptide adopts an alpha-helical structure from residues 5-16 and a type-1 beta-turn from 17-20 with the other residues exhibiting more flexible conformations. The primary hydrophobic interactions with the micelle are from the leucine and phenylalanine residues (Leu-7, Phe-8, Leu-9, Phe-11, Leu-12) while the alanine and glycine residues (Ala-1, Gly-3, Gly-5, Ala-6, Gly-10, Gly-13, Ala-14, Ala-15, Gly-16, Gly-10, Ala-21) interact favorably with water molecules. The results suggest that Phe-8, part of the highly conserved FLG motif of the fusion peptide, plays a key role in the interaction of the peptide with membranes. Our simulations corroborate experimental investigations of the fusion peptide in SDS micelles, providing a high-resolution picture that explains the experimental findings.     

Conformational study of fragments of envelope proteins (gp120: 254-274 and gp41: 519-541) of HIV-1 by NMR and MD simulations.

The envelope proteins, gp 120 and gp41 of HIV-1, play a crucial role in receptor (CD4+ lymphocytes) binding and membrane fusion. The fragment 254-274 of gp120 is conserved in all strains of HIV and, as a part of the full gp120 protein, behaves as 'immunosilent', but as an individual fragment it is 'immunoreactive'. When this fragment binds to its receptor, it activates the fusion domain of gp41 allowing viral entry into the host CD4+ cells. The conformation of fragment 254-274 of the gp120 domain and fragment 519-541 of the gp41 domain was studied by NMR and MD simulations. The studies were carried out in three varied media--water, DMSO-d6 and hexafluoroacetone (HFA). The fusogenic nature of the gp41 domain peptide was investigated by 31P NMR experiments with model bilayers prepared from dimyristoyl-L-alpha-phosphatidylcholine (DMPC). The solvent was seen to exert a major effect on the structure of the two peptides. Fragment (254-274) of gp120 in DMSO-d6 had a type I beta-turn around the tetrad Val9-Ser10-Thr11-Gln12 while in HFA a helical structure spanning the region Ile5 to Gln12 was seen with the remaining part of the peptide in a random coil structure. It is possible that the beta-turn may constitute an initiation site for the formation of the helix. In water at pH 4.5, the peptide adopted a beta-sheet. The NMR results for fragment 519-541 of gp41 are conclusive of a beta-sheet structure in DMSO-d6, a conformation which may help in insertion into the membrane, a notion also put forward by others. The 31P NMR studies of DMPC vesicles with this fragment show its fusogenic nature, promoting fusion of unilamellar vesicles to larger agglomerates like multilamellar ones.     

Enzymatic surface modification and functionalization of PET: A water contact angle,FTIR, and fluorescence spectroscopy study

The purpose of this study was to investigate the changes induced by a lypolytic enzyme on the surface properties of polyethylene terephthalate (PET). Changes in surface hydrophilicity were monitored by means of water contact angle (WCA) measurements. Fourier Transform Infrared spectroscopy (FTIR) in the Attenuated Total Reflectance mode (ATR) was used to investigate the structural and conformational changes of the ethylene glycol and benzene moieties of PET. Amorphous and crystalline PET membranes were used as substrate. The lipolytic enzyme displayed higher hydrolytic activity towards the amorphous PET substrate, as demonstrated by the decrease of the WCA values. Minor changes were observed on the crystalline PET membrane. The effect of enzyme adhesion was addressed by applying a protease after‐treatment which was able to remove the residual enzyme protein adhering to the surface of PET, as demonstrated by the behavior of WCA values. Significant spectral changes were observed by FTIR–ATR analysis in the spectral regions characteristic of the crystalline and amorphous PET domains. The intensity of the crystalline marker bands increased while that of the amorphous ones decreased. Accordingly, the crystallinity indexes calculated as band intensity ratios (1,341/1,410 cm^?1 and 1,120/1,100 cm^?1) increased. Finally, the free carboxyl groups formed at the surface of PET by enzyme hydrolysis were esterified with a fluorescent alkyl bromide, 2‐(bromomethyl)naphthalene (BrNP). WCA measurements confirmed that the reaction proceeded effectively. The fluorescence results indicate that the enzymatically treated PET films are more reactive towards BrNP. FTIR analysis showed that the surface of BrNP‐modified PET acquired a more crystalline character. Biotechnol. Bioeng. 2009;103: 845–856. © 2009 Wiley Periodicals, Inc.     

1-Alkanols and membranes: A story of attraction

Although 1-alkanols have long been known to act as penetration enhancers and anesthetics, the mode of operation is not yet understood. In this study, long-time molecular dynamics simulations have been performed to investigate the effect of 1-alkanols of various carbon chain lengths onto the structure and dynamics of dimyristoylphosphatidylcholine bilayers. The simulations were complemented by microcalorimetry, continuous bleaching and film balance experiments. In the simulations, all investigated 1-alkanols assembled inside the lipid bilayer within tens of nanoseconds. Their hydroxyl groups bound preferentially to the lipid carbonyl group and the hydrocarbon chains stretched into the hydrophobic core of the bilayer. Both molecular dynamics simulations and experiments showed that all 1-alkanols drastically affected the bilayer properties. Insertion of long-chain 1-alkanols decreased the area per lipid while increasing the thickness of the bilayer and the order of the lipids. The bilayer elasticity was reduced and the diffusive motion of the lipids within the bilayer plane was suppressed. On the other hand, integration of ethanol into the bilayer enlarged the area per lipid. The bilayer became softer and lipid diffusion was enhanced.     

Influence of Ionization State on the Activation of Temocapril by hCES1: A Molecular‐Dynamics Study

Temocapril is a prodrug whose hydrolysis by carboxylesterase 1 (CES1) yields the active ACE inhibitor temocaprilat. This molecular‐dynamics (MD) study uses a resolved structure of the human CES1 (hCES1) to investigate some mechanistic details of temocapril hydrolysis. The ionization constants of temocapril (pK1 and pK3) and temocaprilat (pK1, pK2, and pK3) were determined experimentally and computationally using commercial algorithms. The constants so obtained were in good agreement and revealed that temocapril exists mainly in three ionic forms (a cation, a zwitterion, and an anion), whereas temocaprilat exists in four major ionic forms (a cation, a zwitterion, an anion, and a dianion). All these ionic forms were used as ligands in 5‐ns MS simulations. While the cationic and zwitterionic forms of temocapril were involved in an ion‐pair bond with Glu255 suggestive of an inhibitor behavior, the anionic form remained in a productive interaction with the catalytic center. As for temocaprilat, its cation appeared trapped by Glu255, while its zwitterion and anion made a slow departure from the catalytic site and a partial egress from the protein. Only its dianion was effectively removed from the catalytic site and attracted to the protein surface by Lys residues. A detailed mechanism of product egress emerges from the simulations.     

Effective prediction model and determination of binding residues influential for inhibitors targeting HIV-1 integrase-LEDGF/p75 interface by employing solvent accessible surface area energy as key determinant

Abstract

Development of a highly accurate prediction model for protein–ligand inhibition has been a major challenge in drug discovery. Herein, we describe a novel predictive model for the inhibition of HIV-1 integrase (IN)-LEDGF/p75 protein-protein interaction. The model was constructed using energy parameters approximated from molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations. Chemometric analysis using partial least squares (PLS) regression revealed that solvent accessible surface area energy (ΔG_SASA) is the major determinant parameter contributing greatly to the prediction accuracy. PLS prediction model on the ΔG_SASA values collected from 41 complexes yielded a strong correlation between the predicted and the actual inhibitory activities (R² = 0.9666, RMSEC of pIC₅₀ values = 0.0890). Additionally, for the test set of 14 complexes, the model performed satisfactorily with very low pIC₅₀ errors (Q² = 0.5168, RMSEP = 0.3325). A strong correlation between the buried surface areas on the IN protein, when bound with IN-LEDGF/p75 inhibitors, and the respective ΔG_SASA values was also obtained. Furthermore, the current method could identify ‘hot spots’of amino acid residues highly influential to the inhibitory activity prediction. This could present fruitful implications in binding site determination and future inhibitor developments targeting protein-protein interactions.

Communicated by Ramaswamy H. Sarma     

Investigations of peptide hydration using NMR and molecular dynamics simulations: A study of effects of water on the conformation and dynamics of antamanide

Summary The influence of water binding on the conformational dynamics of the cyclic decapeptide antamanide dissolved in the model lipophilic environment chloroform is investigated by NMR relaxation measurements. The water-peptide complex has a lifetime of 35 s at 250 K, which is longer than typical lifetimes of water-peptide complexes reported in aqueous solution. In addition, there is a rapid intracomplex mobility that probably involves librational motions of the bound water or water molecules hopping between different binding sites. Water binding restricts the flexibility of antamanide. The experimental findings are compared with GROMOS molecular dynamics simulations of antamanide with up to eight bound water molecules. Within the simulation time of 600 ps, no water molecule leaves the complex. Additionally, the simulations show a reduced flexibility for the complex in comparison with uncomplexed antamanide. Thus, there is a qualitative agreement between the experimental NMR results and the computer simulations.     

Functional and structural characterization of HIV-1 gp41 ectodomain regions in phospholipid membranes suggests that the fusion-active conformation is extended

HIV-1 entry into its host cell involves a sequential interaction whereby gp41 is in direct contact with the plasma membrane. Understanding the effect of membrane composition on the fusion mechanism can shed light on the unsolved phases of this complex mechanism. Here, we studied N36, a peptide derived from the N-heptad-repeat (NHR) of the gp41 ectodomain, its six helix bundle (SHB) forming counterpart C34, together with the N-terminal 70-mer wild-type peptide (N70), and additional gp41 ectodomain-derived peptides in the presence of two membranes, modeling inner and outer leaflets of the plasma membrane. Information on the structure of these peptides, their affinity towards phospholipids and their ability to induce vesicle fusion was gathered by a variety of fluorescence, spectroscopic and microscopy methods. We found that N36, having strong affinity towards phospholipids, prominently shifts conformation from alpha-helix in an outer leaflet-like zwitterionic membrane to beta-sheet in a membrane mimicking the negatively charged inner leaflet environment, leading to pronounced fusion-activity. Real-time atomic force microscopy (AFM) was used to study the peptides' effect on the membrane morphology, revealing severe bilayer perturbation and extensive pore formation.We also found, that the N36/C34 core is destabilized by electronegative, but not zwitterionic phospholipids. Taken together, our data suggest that the fusion-active pore forming conformation of gp41 is extended, upstream of the SHB. In this manner, folding of the ectodomain into a SHB might also serve as a negative regulator of fusion by impeding gp41 fusion-active surfaces, thus preventing irreversible damage to the cell membrane. This assumption is supported by the finding that pre-incubation of large unilamellar vesicles (LUV) with C-heptad repeat (CHR)-derived fusion inhibitors reduces the fusogenic activity of N-terminal peptides in a dose-dependant manner, and suggests that CHR-derived fusion inhibitors inhibit HIV entry in an analogous mechanism.     

Deprotonation states of the two active site water molecules regulate the binding of protein phosphatase 5 with its substrate: A molecular dynamics study

Protein phosphatase 5 (PP5), mainly localized in human brain, can dephosphorylate tau protein whose high level of phosphorylation is related to Alzheimer's disease. Similar to other protein phosphatases, PP5 has a conserved motif in the catalytic domain that contains two binding sites for manganese (Mn²⁺) ions. Structural data indicate that two active site water molecules, one bridging the two Mn²⁺ ions and the other terminally coordinated with one of the Mn²⁺ ions (Mn1), are involved in catalysis. Recently, a density functional theory study revealed that the two water molecules can be both deprotonated to keep a neutral active site for catalysis. The theoretical study gives us an insight into the catalytic mechanism of PP5, but the knowledge of how the deprotonation states of the two water molecules affect the binding of PP5 with its substrate is still lacking. To approach this problem, molecular dynamics simulations were performed to model the four possible deprotonation states. Through structural, dynamical and energetic analyses, the results demonstrate that the deprotonation states of the two water molecules affect the structure of the active site including the distance between the two Mn²⁺ ions and their coordination, impact the interaction energy of residues R275, R400 and H304 which directly interact with the substrate phosphoserine, and mediate the dynamics of helix αJ which is involved in regulation of the enzyme's activity. Furthermore, the deprotonation state that is preferable for PP5 binding of its substrate has been identified. These findings could provide new design strategy for PP5 inhibitor.     

A spectroscopic and voltammetric study of the pH-dependent Cu(II) coordination to the peptide GGGTH: relevance to the fifth Cu(II) site in the prion protein

The GGGTH sequence has been proposed to be the minimal sequence involved in the binding of a fifth Cu(II) ion in addition to the octarepeat region of the prion protein (PrP) which binds four Cu(II) ions. Coordination of Cu(II) by the N- and C-protected Ac-GGGTH-NH(2) pentapeptide (P(5)) was investigated by using potentiometric titration, electrospray ionization mass spectrometry, UV-vis spectroscopy, electron paramagnetic resonance (EPR) spectroscopy and cyclic voltammetry experiments. Four different Cu(II) complexes were identified and characterized as a function of pH. The Cu(II) binding mode switches from NO(3) to N(4) for pH values ranging from 6.0 to 10.0. Quasi-reversible reduction of the [Cu(II)(P(5))H(-2)] complex formed at pH 6.7 occurs at E (1/2)=0.04 V versus Ag/AgCl, whereas reversible oxidation of the [Cu(II)(P(5))H(-3)](-) complex formed at pH 10.0 occurs at E (1/2)=0.66 V versus Ag/AgCl. Comparison of our EPR data with those of the rSHaPrP(90-231) (Burns et al. in Biochemistry 42:6794-6803, 2003) strongly suggests an N(3)O binding mode at physiological pH for the fifth Cu(II) site in the protein.     

Structure and plasticity of the human immunodeficiency virus gp41 fusion domain in lipid micelles and bilayers

A thorough understanding of the structure of fusion domains of enveloped viruses in changing lipid environments helps us to formulate mechanistic models on how they might function in mediating viral entry by membrane fusion. We have expressed the N-terminal fusion domain of HIV-1 gp41 as a construct that is water-soluble in the absence of membranes, but that also binds with high affinity to lipid micelles and bilayers in their presence. We have solved the structure and studied the dynamics of this domain bound to dodecylphosphocholine micelles by homo- and heteronuclear NMR spectroscopy. The fusion peptide forms a stable hydrophobic helix from Ile(4) to Ala(14), but is increasingly more disordered and dynamic in a segment of intermediate polarity that stretches from Ala(15) to Ser(23). When bound to lipid bilayers at low concentration, the HIV fusion domain is also largely alpha-helical, as determined by CD and FTIR spectroscopy. However, at higher protein/lipid ratios, the domain is partially converted to form beta-structures in lipid bilayers. Controlled lipid mixing occurs at concentrations that support the alpha-helical, but not the beta-strand conformation.     

Properties and structures of the influenza and HIV fusion peptides on lipid membranes: implications for a role in fusion

The fusion peptides of HIV and influenza virus are crucial for viral entry into a host cell. We report the membrane-perturbing and structural properties of fusion peptides from the HA fusion protein of influenza virus and the gp41 fusion protein of HIV. Our goals were to determine: 1), how fusion peptides alter structure within the bilayers of fusogenic and nonfusogenic lipid vesicles and 2), how fusion peptide structure is related to the ability to promote fusion. Fluorescent probes revealed that neither peptide had a significant effect on bilayer packing at the water-membrane interface, but both increased acyl chain order in both fusogenic and nonfusogenic vesicles. Both also reduced free volume within the bilayer as indicated by partitioning of a lipophilic fluorophore into membranes. These membrane ordering effects were smaller for the gp41 peptide than for the HA peptide at low peptide/lipid ratio, suggesting that the two peptides assume different structures on membranes. The influenza peptide was predominantly helical, and the gp41 peptide was predominantly antiparallel beta-sheet when membrane bound, however, the depths of penetration of Trps of both peptides into neutral membranes were similar and independent of membrane composition. We previously demonstrated: 1), the abilities of both peptides to promote fusion but not initial intermediate formation during PEG-mediated fusion and 2), the ability of hexadecane to compete with this effect of the fusion peptides. Taken together, our current and past results suggest a hypothesis for a common mechanism by which these two viral fusion peptides promote fusion.     

Structural and functional properties of peptides based on the N-terminus of HIV-1 gp41 and the C-terminus of the amyloid-beta protein

Given their high alanine and glycine levels, plaque formation, α-helix to β-sheet interconversion and fusogenicity, FP (i.e., the N-terminal fusion peptide of HIV-1 gp41; 23 residues) and amyloids were proposed as belonging to the same protein superfamily. Here, we further test whether FP may exhibit ‘amyloid-like’ characteristics, by contrasting its structural and functional properties with those of Aβ(26-42), a 17-residue peptide from the C-terminus of the amyloid-beta protein responsible for Alzheimer's. FTIR spectroscopy, electron microscopy, light scattering and predicted amyloid structure aggregation (PASTA) indicated that aqueous FP and Aβ(26-42) formed similar networked β-sheet fibrils, although the FP fibril interactions were weaker. FP and Aβ(26-42) both lysed and aggregated human erythrocytes, with the hemolysis-onsets correlated with the conversion of α-helix to β-sheet for each peptide in liposomes. Congo red (CR), a marker of amyloid plaques in situ, similarly inhibited either FP- or Aβ(26-42)-induced hemolysis, and surface plasmon resonance indicated that this may be due to direct CR-peptide binding. These findings suggest that membrane-bound β-sheets of FP may contribute to the cytopathicity of HIV in vivo through an amyloid-type mechanism, and support the classification of HIV-1 FP as an ‘amyloid homolog’ (or ‘amylog’).     

Complementary biophysical tools to investigate lipid specificity in the interaction between bioactive molecules and the plasma membrane: A review

Plasma membranes are complex entities common to all living cells. The basic principle of their organization appears very simple, but they are actually of high complexity and represent very dynamic structures. The interactions between bioactive molecules and lipids are important for numerous processes, from drug bioavailability to viral fusion. The cell membrane is a carefully balanced environment and any change inflicted upon its structure by a bioactive molecule must be considered in conjunction with the overall effect that this may have on the function and integrity of the membrane. Conceptually, understanding the molecular mechanisms by which bioactive molecules interact with cell membranes is of fundamental importance.     

Is the Sudlow site I of human serum albumin more generous to adopt prospective anti-cancer bioorganic compound than that of bovine: A combined spectroscopic and docking simulation approach

A comparative biophysical study on the individual conformational adaptation embraced by two homologous serum albumins (SA) (bovine and human) towards a potential anticancer bioorganic compound 2-(6-chlorobenzo[d] thiazol-2-yl)-1H-benzo[de] isoquinoline-1,3(2H)- dione (CBIQD) is apparent from the discrimination in binding behavior and the ensuing consequences accomplished by combined in vitro optical spectroscopy, in silico molecular docking and molecular dynamics (MD) simulation. The Sudlow site I of HSA although anion receptive, harbors neutral CBIQD in Sudlow site I (subdomain IIA, close to Trp) of HSA, while in BSA its prefers to snugly fit into Sudlow site II (subdomain IIIA, close to Tyr). Based on discernable diminution of HSA mean fluorescence lifetime as a function of biluminophore concentration, facile occurrence of fluorescence resonance energy transfer (FRET) is substantiated as the probable quenching mechanism accompanied by structural deformations in the protein ensemble. CBIQD establishes itself within HSA close to Trp214, and consequently reduces the micropolarity of the cybotactic environment that is predominantly constituted by hydrophobic amino acid residues. The stronger association of CBIQD with HSA encourages an allosteric modulation leading to slight deformation in its secondary structure whereas for BSA the association is comparatively weaker. Sudlow site I of HSA is capable to embrace a favorable conformation like malleable gold to provide room for incoming CBIQD, whereas for BSA it behaves more like rigid cast-iron which does not admit any change thus forcing CBIQD to occupy an altogether different binding location i.e. the Sudlow site II. The anticancer CBIQD is found to be stable within the HSA scaffold as vindicated by root mean square deviation (RMSD) and root mean square fluctuation (RMSF) obtained by MD simulation. A competitively inhibited esterase-like activity of HSA upon CBIQD binding to Lys199 and Arg257 residues, plausibly envisions that similar naphthalimide based prodrugs, bearing ester functionality, can be particularly activated by Sudlow site I of HSA. The consolidated spectroscopic research described herein may encourage design of naphthalimide based pro-drugs for effective in vivo biodistribution using HSA-based drug delivery systems.     

Conformational analysis by NMR and molecular dynamics of adamantane-doxorubicin prodrugs and their assemblies with β-cyclodextrin: A focus on the design of platforms for controlled drug delivery

Nanoscale design and construction of affinity-based drug delivery systems (ADDS) is an active research area with enormous potential for the improvement of cancer treatment. For the therapeutic load of these ADDS, a promising strategy is the design of pH-sensitive prodrugs based on the construction of conjugates between adamantane and doxorubicin (Ad-Dox), which stands out as an excellent model system to obtain novel supramolecular materials. Construction of these prodrugs involves a modification of three zones of doxorubicin which in principle does not affect the action mechanism: the carbonyl group C13 (hydrazone linker), the primary alcohol neighboring the carbonyl (ester linker) and the 3′ amino group of daunosamine sugar (amide linker). These modifications are aimed to improve the efficacy and reduce the systemic toxicity of the drug chemotherapy by controlling its release in cancer cells. In this work, we performed 2D NMR experiments and molecular dynamics simulations to characterize the conformational changes of three constructed prodrugs. Our results demonstrated that ring A and the daunsamine sugar of the hydrazone and amide linkers conserve the half-chair state ⁹H₈, while the ester linker disrupts this conformation. Our study also showed that the hydrazone-linked compound (Ad-h-Dox) does not modify the conformation of the original drug and maintains cytotoxic activity. Moreover, the inclusion complex (IC) of Ad-h-Dox with β-cyclodextrin (βCD) generated a highly soluble platform in water, whereas the ester-linked compound (Ad-e-Dox) causes the loss of biological activity. This study proves that Ad-h-Dox prodrug can be an optimum prodrug and act as a building block for a more complex drug transport system.