Molecular Dynamics Study of the Solvation of an α-Helical Transmembrane Peptide by DMSO

A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent

Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 μs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 μs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.

In silico folding of hydrophobic peptides that form β-hairpin structures in solution

Peptides designed with residues that have a high propensity to occur in β-turns form β-hairpin structures in apolar as well as in polar organic solvents such as dimethyl sulfoxide (DMSO). Due to limited solubility, their conformations have not been investigated experimentally in water. We have examined the conformations of four of such designed peptides that fold into well-defined β-hairpin structures facilitated by β-turns, in the crystalline state and in solution, by molecular dynamics simulations (MDS). The peptides folded into β-hairpin structures in water, starting from the fully extended conformation. However, in DMSO, neither folding nor unfolding was observed during MDS, when the starting structures were unfolded and folded, respectively. The lack of folding in DMSO was investigated by constructing folding free energy landscapes by umbrella sampling. The folding free energy landscape is smooth in water while in DMSO folded and unfolded structures are separated by high-energy barriers. The folding free energy is less in DMSO compared to water due to a more stable unfolded structure in DMSO compared to water, which in turn is due to stabilization of the unfolded state by hydrophobic interactions in DMSO. This finding will be helpful to researchers to accurately model and/or design small peptide structures in water and organic solvents.

Using Molecular Simulation to Guide Protein Engineering for Biocatalysis in Organic Solvents

Biocatalysis in organic solvents (OSs) is very appealing for the industry in producing bulk and/or fine chemicals, such as pharmaceuticals, biodiesel, and fragrances. The poor performance of enzymes in OSs (e.g., reduced activity, insufficient stability, and deactivation) negates OSs' excellent solvent properties. Molecular dynamics (MD) simulations provide a complementary method to study the relationship between enzymes dynamics and the stability in OSs. Here we describe computational procedure for MD simulation of enzymes in OSs with an example of Bacillus subtilis lipase A (BSLA) in dimethyl sulfoxide (DMSO) cosolvent with software GROMACS. We discuss main essential practical issues considered (such as choice of force field, parameterization, simulation setup, and trajectory analysis). The core part of this protocol (enzyme-OS system setup, analysis of structural-based and solvation-based observables) is transferable to other enzymes and any OS systems. Combining with experimental studies, the obtained molecular knowledge is most likely to guide researchers to access rational protein engineering approaches to tailor OS resistant enzymes and expand the scope of biocatalysis in OS media. Finally, we discuss potential solutions to overcome the remaining challenges of computational biocatalysis in OSs and briefly draw future directions for further improvement in this field.

Is Dissociation of HCl in DMSO Clusters Bistable?

Dissociation of HCl embedded in dimethyl sulfoxide (DMSO) clusters was investigated by projecting the solvent electric field along the HCl bond using B3LYP-D3/6-31+G(d) and MP2/6-31+G(d,p) levels of theory. A large number of distinct structures (about 1500) consisting of up to five DMSO molecules were considered in the present work for statistical reliability. The B3LYP-D3 calculations reveal that the dissociation of HCl embedded in DMSO clusters requires a critical electric field of 138 MV cm^-1 along the H-Cl bond. However, a large number of exceptions wherein the electric field values much higher than the critical electric field of 138 MV cm^-1 did not result in dissociation of HCl were observed, in addition to several cases wherein the HCl dissociates with an electric field less than the critical electric field. On the other hand, the MP2 level calculations reveal that the critical electric field for HCl dissociation is about 181 MV cm^-1 with almost no exceptions. A comparison of calculations carried out using the MP2 and the B3LYP-D3 levels suggests that the dissociation of HCl embedded in DMSO clusters is bistable at the B3LYP-D3 level, which is an artifact, suggesting that care must be exercised in interpreting the processes of proton transfer. The answer to the question raised as the title of this paper is NO.

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala₂₅ ). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala₂₅ . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅ . This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala₂₅ . Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.

Rational Design of Antigen Incorporation Into Subunit Vaccine Biomaterials Can Enhance Antigen-Specific Immune Responses

Peptide subunit vaccines increase safety by reducing the risk of off-target responses and improving the specificity of the induced adaptive immune response. The immunogenicity of most soluble peptides, however, is often insufficient to produce robust and lasting immunity. Many biomaterials and delivery vehicles have been developed for peptide antigens to improve immune response while maintaining specificity. Peptide nanoclusters (PNC) are a subunit peptide vaccine material that has shown potential to increase immunogenicity of peptide antigens. PNC are comprised only of crosslinked peptide antigen and have been synthesized from several peptide antigens as small as 8 amino acids in length. However, as with many peptide vaccine biomaterials, synthesis requires adding residues to the peptide and/or engaging amino acids within the antigen epitope covalently to form a stable material. The impact of antigen modifications made to enable biomaterial incorporation or formation is rarely investigated, since the goal of most studies is to compare the soluble antigen with biomaterial form of antigen. This study investigates PNC as a platform vaccine biomaterial to evaluate how peptide modification and biomaterial formation with different crosslinking chemistries affect epitope-specific immune cell presentation and activation. Several types of PNC were synthesized by desolvation from the model peptide epitope SIINFEKL, which is derived from the immunogenic protein ovalbumin. SIINFEKL was altered to include extra residues on each end, strategically chosen to enable multiple conjugation chemistry options for incorporation into PNC. Several crosslinking methods were used to control which functional groups were used to stabilize the PNC, as well as the reducibility of the crosslinking. These variations were evaluated for immune responses and biodistribution following in vivo immunization. All modified antigen formulations still induced comparable immune responses when incorporated into PNC compared to unmodified soluble antigen alone. However, some crosslinking methods led to a significant increase in desirable immune responses while others did not, suggesting that not all PNC were processed the same. These results help guide future peptide vaccine biomaterial design, including PNC and a wide variety of conjugated and self-assembled peptide antigen materials, to maximize and tune the desired immune response.

The role of the jaw subdomain of peptidoglycan glycosyltransferases for lipid II polymerization

Bacterial peptidoglycan glycosyltransferases (PGT) catalyse the essential polymerization of lipid II into linear glycan chains required for peptidoglycan biosynthesis. The PGT domain is composed of a large head subdomain and a smaller jaw subdomain and can be potently inhibited by the antibiotic moenomycin A (MoeA). We present an X-ray structure of the MoeA-bound Staphylococcus aureus monofunctional PGT enzyme, revealing electron density for a second MoeA bound to the jaw subdomain as well as the PGT donor site. Isothermal titration calorimetry confirms two drug-binding sites with markedly different affinities and positive cooperativity. Hydrophobic cluster analysis suggests that the membrane-interacting surface of the jaw subdomain has structural and physicochemical properties similar to amphipathic cationic α -helical antimicrobial peptides for lipid II recognition and binding. Furthermore, molecular dynamics simulations of the drug-free and -bound forms of the enzyme demonstrate the importance of the jaw subdomain movement for lipid II selection and polymerization process and provide molecular-level insights into the mechanism of peptidoglycan biosynthesis by PGTs.

Scrutiny of electrostatic-driven conformational ordering of polypeptide chains in DMSO: a study with a model oligopeptide

The molecular mechanism of DMSO-induced stabilisation of β-sheets is attributed to the combination of polar electrostatic interactions among side chains, and backbone desolvation through bulky side chains which promotes backbone hydrogen bonding.

Unscrambling micro-solvation of –COOH and –NH groups in neat dimethyl sulfoxide: insights from 1H-NMR spectroscopy and computational studies

This study investigates the interactions of –COOH and –NH groups in neat DMSO solutions, with special focus on their thermodynamics and kinetics.

Combined effects of solvation and aggregation propensity on the final supramolecular structures adopted by hydrophobic, glycine-rich, elastin-like polypeptides

Previous work on elastin-like polypeptides (ELPs) made of hydrophobic amino acids of the type XxxGlyGlyZzzGly (Xxx, Zzz = Val, Leu) has consistently shown that differing dominant supramolecular structures were formed when the suspending media were varied: helical, amyloid-like fibers when suspended in water and globules evolving into "string of bead" structures, poly(ValGlyGlyValGly), or cigar-like bundles, poly(ValGlyGlyLeuGly), when suspended in methyl alcohol. Comparative experiments with poly(LeuGlyGlyValGly) have further indicated that the interface energy plays a significant role and that solvation effects act in concomitance with the intrinsic aggregation propensity of the repeat sequence. Continuing our investigation on ELPs using surface (X-ray photoelectron spectroscopy, atomic force microscopy) and bulk (circular dichroism, Fourier transform infrared spectroscopy) techniques for their characterization, here we have compared the effect of suspending solvents (H(2)O, dimethylsulfoxide, ethylene glycol, and MeOH) on poly(ValGlyGlyValGly), the polypeptide most inclined to form long and well-refined helical fibers in water, searching for the signature of intermolecular interactions occurring between the polypeptide chains in the given suspension. The influence of sequence specificities has been studied by comparing poly(ValGlyGlyValGly) and poly(LeuGlyGlyValGly) with a similar degree of polymerization. Deposits on substrates of the polypeptides were characterized taking into account the differing evaporation rate of solvents, and tests on their stability in ultra high vacuum were performed. Finally, combining experimental and computational studies, we have revaluated the three-dimensional modeling previously proposed for the supramolecular assembly in water of poly(ValGlyGlyValGly). The results were discussed and rationalized also in the light of published data.

Conformational studies of peptides representing a segment of TM7 from H+-VO-ATPase in SDS micelles

The conformation of a transmembrane peptide, sMTM7, encompassing the cytoplasmic hemi-channel domain of the seventh transmembrane section of subunit a from V-ATPase from Saccharomyces cerevisiae solubilized in SDS solutions was studied by circular dichroism (CD) spectroscopy and fluorescence spectroscopy of the single tryptophan residue of this peptide. The results show that the peptide adopts an alpha-helical conformation or aggregated beta-sheet depending on the peptide-to-SDS ratio used. The results are compared with published data about a longer version of the peptide (i.e., MTM7). It is concluded that the bulky, positively charged arginine residue located in the center of both peptides has a destabilizing effect on the helical conformation of the SDS-solubilized peptides, leading to beta-sheet formation and subsequent aggregation.

Predictions suggesting a participation of beta-sheet configuration in the M2 domain of the P2X(7) receptor: a novel conformation?

Scanning experiments have shown that the putative TM2 domain of the P2X(7) receptor (P2X(7)R) lines the ionic pore. However, none has identified an alpha-helix structure, the paradigmatic secondary structure of ion channels in mammalian cells. In addition, some researchers have suggested a beta-sheet conformation in the TM2 domain of P2X(2). These data led us to investigate a new architecture within the P2X receptor family. P2X(7)R is considered an intriguing receptor because its activation induces nonselective large pore formation, in contrast to the majority of other ionic channel proteins in mammals. This receptor has two states: a low-conductance channel (approximately 10 pS) and a large pore (> 400 pS). To our knowledge, one fundamental question remains unanswered: Are the P2X(7)R channel and the pore itself the same entity or are they different structures? There are no structural data to help solve this question. Thus, we investigated the hydrophobic M2 domain with the aim of predicting the fitted position and the secondary structure of the TM2 segment from human P2X(7)R (hP2X(7)R). We provide evidence for a beta-sheet conformation, using bioinformatics algorithms and molecular-dynamics simulation in conjunction with circular dichroism in different environments and Fourier transform infrared spectroscopy. In summary, our study suggests the possibility that a segment composed of residues from part of the M2 domain and part of the putative TM2 segment of P2X(7)R is partially folded in a beta-sheet conformation, and may play an important role in channel/pore formation associated with P2X(7)R activation. It is important to note that most nonselective large pores have a transmembrane beta-sheet conformation. Thus, this study may lead to a paradigmatic change in the P2X(7)R field and/or raise new questions about this issue.