Structural Dynamics and Conformational Equilibria of SERCA Regulatory Proteins in Membranes by Solid-State NMR Restrained Simulations

Solid-state NMR spectroscopy is emerging as a powerful approach to determine structure, topology, and conformational dynamics of membrane proteins at the atomic level. Conformational dynamics are often inferred and quantified from the motional averaging of the NMR parameters. However, the nature of these motions is difficult to envision based only on spectroscopic data. Here, we utilized restrained molecular dynamics simulations to probe the structural dynamics, topology and conformational transitions of regulatory membrane proteins of the calcium ATPase SERCA, namely sarcolipin and phospholamban, in explicit lipid bilayers. Specifically, we employed oriented solid-state NMR data, such as dipolar couplings and chemical shift anisotropy measured in lipid bicelles, to refine the conformational ensemble of these proteins in lipid membranes. The samplings accurately reproduced the orientations of transmembrane helices and showed a significant degree of convergence with all of the NMR parameters. Unlike the unrestrained simulations, the resulting sarcolipin structures are in agreement with distances and angles for hydrogen bonds in ideal helices. In the case of phospholamban, the restrained ensemble sampled the conformational interconversion between T (helical) and R (unfolded) states for the cytoplasmic region that could not be observed using standard structural refinements with the same experimental data set. This study underscores the importance of implementing NMR data in molecular dynamics protocols to better describe the conformational landscapes of membrane proteins embedded in realistic lipid membranes.     

Conformational Analysis of Neuropeptide Y Segments by CD,NMR Spectroscopy and Restrained Molecular Dynamics

Neuropeptide Y (NPY), a peptide amide comprising 36 residue has been shown to act as a potent vasoconstrictor. In order to shed light on the structural requirements for the biological activities with respect to the different prerequisites for affinity to the NPY receptor subtypes Y₁ and Y₂, in the present study the syntheses and conformational analyses of two C-terminal segments, NPY(18–36) and NPY(13–36), are described. The results obtained by CD measurements, two-dimensional NMR spectros copy and a conformational refinement of the NMR-derived structure by molecular mechanics simulations support the findings of previously published structure –activity relationship studies for biologically active and selective compounds. In particular, the α-helical conformation as well as an appropriate exposure of the side chains of the critical C-terminal dipeptide within NPY(18–36) are in agreement with the prerequisites proposed for Y₂ receptor binding of that segment.     

Conformational Studies of Two Bradykinin Antagonists by Using Two-dimensional NMR Techniques and Molecular Dynamics Simulations

Abstract

The solution conformations of two potent antagonists of bradykinin (Arg¹-Pro²-Pro³-Gly⁴- Phe⁵-Ser⁶-Pro⁷-Phe⁸-Arg⁹), [Aca-¹, DArg⁰, Hyp³, Thi⁵, DPhe⁷,(N-Bzl)Gly⁸]BK (1) and [Aaa- ¹, DArg⁰, Hyp³, Thi⁵,(2-DNal)⁷, Thi⁸]BK (2), were studied by using 2D NMR spectroscopy in DMSO-dg and molecular dynamics simulations. The NMR spectra of peptide 1 reveals the existence of at least two isomers arising from isomerization across the DPhe⁷-(N-Bzl)Gly⁸peptide bond. The more populated isomer possesses the cis peptide bond at this position. The ratio of cis/trans isomers amounted to 7:3. With both antagonists, the NMR data indicate a β-turn structure for the Hyp³-Gly⁴ residues. In addition, for peptide 2, position 2,3 is likely to be occupied by turn-like structures. The cis peptide bond between DPhe⁷ and (N- Bzl)Gly⁸ in analogue 1 suggests type VI β-turn at position 7,8. The molecular dynamics runs were performed on both peptides in DMSO solution. The results indicate that the structure of peptide 1 is characterized by type VIb β-turn comprising residues Ser-Arg⁹ and the βI or βII-turn involving the Pro²-Thi⁵ fragment, whereas peptide 2 shows the tendency towards the formation of type I β-turn at position 2,3. The structures of both antagonists are stabilized by a salt bridge between the guanidine moiety of Arg¹ and the carboxyl group of Arg⁹. Moreover, the side chain of DArg⁰ is apart of the rest of molecule and is not involved in structural elements except for a few calculated structures.     

Light Activation of Rhodopsin: Insights from Molecular Dynamics Simulations Guided by Solid-State NMR Distance Restraints

Structural restraints provided by solid-state NMR measurements of the metarhodopsin II intermediate are combined with molecular dynamics simulations to help visualize structural changes in the light activation of rhodopsin. Since the timescale for the formation of the metarhodopsin II intermediate (> 1 ms) is beyond that readily accessible by molecular dynamics, we use NMR distance restraints derived from ¹³C dipolar recoupling measurements to guide the simulations. The simulations yield a working model for how photoisomerization of the 11-cis retinylidene chromophore bound within the interior of rhodopsin is coupled to transmembrane helix motion and receptor activation. The mechanism of activation that emerges is that multiple switches on the extracellular (or intradiscal) side of rhodopsin trigger structural changes that converge to disrupt the ionic lock between helices H3 and H6 on the intracellular side of the receptor.     

Correlation Spectroscopy and Molecular Dynamics Simulations to Study the Structural Features of Proteins

In this work, we used a combination of fluorescence correlation spectroscopy (FCS) and molecular dynamics (MD) simulation methodologies to acquire structural information on pH-induced unfolding of the maltotriose-binding protein from Thermus thermophilus (MalE2). FCS has emerged as a powerful technique for characterizing the dynamics of molecules and it is, in fact, used to study molecular diffusion on timescale of microsecond and longer. Our results showed that keeping temperature constant, the protein diffusion coefficient decreased from 84±4 µm²/s to 44±3 µm²/s when pH was changed from 7.0 to 4.0. An even more marked decrease of the MalE2 diffusion coefficient (31±3 µm²/s) was registered when pH was raised from 7.0 to 10.0. According to the size of MalE2 (a monomeric protein with a molecular weight of 43 kDa) as well as of its globular native shape, the values of 44 µm²/s and 31 µm²/s could be ascribed to deformations of the protein structure, which enhances its propensity to form aggregates at extreme pH values. The obtained fluorescence correlation data, corroborated by circular dichroism, fluorescence emission and light-scattering experiments, are discussed together with the MD simulations results.     

Conformational Dynamics of Subtilisin-Chymotrypsin Inhibitor 2 Complex by Coarse-Grained Simulations

Abstract

An off-lattice dynamic Monte Carlo (MC) method is used to investigate the conformational dynamics of chymotrypsin inhibitor 2 (CI2) and subtilisin in both free and complex forms over two time windows, referring to short and long time scales. The conformational dynamics of backbone bonds analysed from several independent trajectories reveal that: Both the inhibitor and the enzyme are restricted in their bond rotations, excluding a few bonds, upon binding; the effect being greatest for the loop regions, and for the inhibitor. A cooperativity in the near-neighbor bond rotations are observed on both time scales, whereas the cooperative rotations of the bonds far along the sequence appear only in the long time window, and the latter time window is where most of the interactions between the inhibitor and the enzyme are observed. Upon binding, the cooperatively rotating parts of the inhibitor and the enzyme are readjusted compared to their free forms, and new correlations appear. The binding loop, although it is the closest contact region, is not the only part of the inhibitor involved in the interactions with the enzyme. Loops 3 and 8 and the helices F and G in bound enzyme and the binding loop of the inhibitor contribute at the most to the collective motions of whole structure on the slow time scale and are apparently important for enzyme-inhibitor interactions and function. The results in general provide evidence for the contribution of the loops with cooperative motions to the extensive communication network of the complex.     

Structural Characterization of Polyglutamine Fibrils by Solid-State NMR Spectroscopy

Protein aggregation via polyglutamine stretches occurs in a number of severe neurodegenerative diseases such as Huntington's disease. We have investigated fibrillar aggregates of polyglutamine peptides below, at, and above the toxicity limit of around 37 glutamine residues using solid-state NMR and electron microscopy. Experimental data are consistent with a dry fibril core of at least 70-80 Å in width for all constructs. Solid-state NMR dipolar correlation experiments reveal a largely β-strand character of all samples and point to tight interdigitation of hydrogen-bonded glutamine side chains from different sheets. Two approximately equally frequent populations of glutamine residues with distinct sets of chemical shifts are found, consistent with local backbone dihedral angles compensating for β-strand twist or with two distinct sets of side-chain conformations. Peptides comprising 15 glutamine residues are present as single extended β-strands. Data obtained for longer constructs are most compatible with a superpleated arrangement with individual molecules contributing β-strands to more than one sheet and an antiparallel assembly of strands within β-sheets.     

Active-Site Structure of the Thermophilic Foc-Subunit Ring in Membranes Elucidated by Solid-State NMR

F_oF₁-ATP synthase uses the electrochemical potential across membranes or ATP hydrolysis to rotate the F_oc-subunit ring. To elucidate the underlying mechanism, we carried out a structural analysis focused on the active site of the thermophilic c-subunit (TF_oc) ring in membranes with a solid-state NMR method developed for this purpose. We used stereo-array isotope labeling (SAIL) with a cell-free system to highlight the target. TF_oc oligomers were purified using a virtual ring His tag. The membrane-reconstituted TF_oc oligomer was confirmed to be a ring indistinguishable from that expressed in E. coli on the basis of the H⁺-translocation activity and high-speed atomic force microscopic images. For the analysis of the active site, 2D ¹³C-¹³C correlation spectra of TF_oc rings labeled with SAIL-Glu and -Asn were recorded. Complete signal assignment could be performed with the aid of the C^α_i+1-C^α_i correlation spectrum of specifically ¹³C,¹⁵N-labeled TF_oc rings. The C^δ chemical shift of Glu-56, which is essential for H⁺ translocation, and related crosspeaks revealed that its carboxyl group is protonated in the membrane, forming the H⁺-locked conformation with Asn-23. The chemical shift of Asp-61 C^γ of the E. coli c ring indicated an involvement of a water molecule in the H⁺ locking, in contrast to the involvement of Asn-23 in the TF_oc ring, suggesting two different means of proton storage in the c rings.     

Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes

Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys⁴⁰, residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.     

Structural Dynamics of Human Telomeric G-Quadruplex Loops Studied by Molecular Dynamics Simulations

Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.     

Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.     

Characterization of Stratum Corneum Molecular Dynamics by Natural-Abundance 13C Solid-State NMR

Despite the enormous potential for pharmaceutical applications, there is still a lack of understanding of the molecular details that can contribute to increased permeability of the stratum corneum (SC). To investigate the influence of hydration and heating on the SC, we record the natural-abundance ¹³C signal of SC using polarization transfer solid-state NMR methods. Resonance lines from all major SC components are assigned. Comparison of the signal intensities obtained with the INEPT and CP pulse sequences gives information on the molecular dynamics of SC components. The majority of the lipids are rigid at 32°C, and those lipids co-exist with a small pool of mobile lipids. The ratio between mobile and rigid lipids increases with hydration. An abrupt change of keratin filament dynamics occurs at RH = 80–85%, from completely rigid to a structure with rigid backbone and mobile protruding terminals. Heating has a strong effect on the lipid mobility, but only a weak influence on the keratin filaments. The results provide novel molecular insight into how the SC constituents are affected by hydration and heating, and improve the understanding of enhanced SC permeability, which is associated with elevated temperatures and SC hydration.     

壳聚糖与细菌细胞膜相互作用的分子动力学研究

目的 壳聚糖（chitosan,CS）是一种天然的广谱抗菌活性物质。现有研究表明,壳聚糖与细菌细胞膜的相互作用是其发挥抗菌功能的关键。受限于传统实验技术的表征能力,壳聚糖与细菌细胞膜相互作用的具体机制仍有待研究。本文旨在研究壳聚糖与细菌细胞膜相互作用的分子机制。方法 本研究利用全原子分子动力学模拟技术主要探究了完全脱乙酰化的不同聚合度壳聚糖（八聚糖、十二聚糖和十六聚糖）与革兰氏阴性菌外膜（outer membrane,OM）和革兰氏阳性菌质膜（cytoplasmic membrane,CM）相互作用的动态过程。结果 壳聚糖主要依靠其氨基、碳6位羟基和碳3位羟基与OM和CM的头部极性区发生快速结合。随后壳聚糖末端糖基单元倾向于插入OM内部,深度约1 nm,并与脂质分子脂肪酸链上的羰基形成稳定的氢键相互作用。与之相比,壳聚糖分子难以稳定地插入CM内部。壳聚糖结合对膜结构性质产生影响,主要表现在降低OM和CM的单分子脂质面积,显著减少OM和CM极性区的Ca2+和Na+数目,破坏阳离子介导的脂质间相互作用。结论 本研究证明,壳聚糖带正电的氨基基团是介导其与膜相互作用的关键,并破环脂质间的相互作...     

Nanotube Array Method for Studying Lipid-Induced Conformational Changes of a Membrane Protein by Solid-State NMR

Anodic aluminum oxide substrates with macroscopically aligned homogeneous nanopores of 80 nm in diameter enable two-dimensional, solid-state nuclear magnetic resonance studies of lipid-induced conformational changes of uniformly ¹⁵N-labeled Pf1 coat protein in native-like bilayers. The Pf1 helix tilt angles in bilayers composed of two different lipids are not entirely governed by the membrane thickness but could be rationalized by hydrophobic interactions of lysines at the bilayer interface. The anodic aluminum oxide alignment method is applicable to a broader repertoire of lipids versus bicelle bilayer mimetics currently employed in solid-state nuclear magnetic resonance of oriented samples, thus allowing for elucidation of the role played by lipids in shaping membrane proteins.