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.     

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.     

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.     

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.     

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.     

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.     

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.     

Molecular Dynamics of a Schiff Base Tetramacrocycle Studied by NMR and MD Simulations

Schiff base condensation of the pentadentate tetrapodal amine 1 with 2,6-diformyl-4-methylphenol 2 leads in a [3 + 6] condensation to a previously not described macrocycle 3 which shows C₃-symmetry. X-ray analysis shows a truncated cone shape for 3. At T S 100°C, the ¹H-ROESY/EXSY spectrum reveals pairwise exchange of corresponding sites, indicative of inversion of the whole molecule in an umbrella-like fashion. Molecular dynamics simulations support this hypothesis.     

Capturing Conformational States in Proteins Using Sparse Paramagnetic NMR Data

Capturing conformational changes in proteins or protein-protein complexes is a challenge for both experimentalists and computational biologists. Solution nuclear magnetic resonance (NMR) is unique in that it permits structural studies of proteins under greatly varying conditions, and thus allows us to monitor induced structural changes. Paramagnetic effects are increasingly used to study protein structures as they give ready access to rich structural information of orientation and long-range distance restraints from the NMR signals of backbone amides, and reliable methods have become available to tag proteins with paramagnetic metal ions site-specifically and at multiple sites. In this study, we show how sparse pseudocontact shift (PCS) data can be used to computationally model conformational states in a protein system, by first identifying core structural elements that are not affected by the environmental change, and then computationally completing the remaining structure based on experimental restraints from PCS. The approach is demonstrated on a 27 kDa two-domain NS2B-NS3 protease system of the dengue virus serotype 2, for which distinct closed and open conformational states have been observed in crystal structures. By changing the input PCS data, the observed conformational states in the dengue virus protease are reproduced without modifying the computational procedure. This data driven Rosetta protocol enables identification of conformational states of a protein system, which are otherwise difficult to obtain either experimentally or computationally.     

A H Solid-State NMR Study of Lipid Clustering by Cationic Antimicrobial and Cell-Penetrating Peptides in Model Bacterial Membranes

Domain formation in bacteria-mimetic membranes due to cationic peptide binding was recently proposed based on calorimetric data. We now use ²H solid-state NMR to critically examine the presence and absence of domains in bacterial membranes containing zwitterionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylethanolamine (POPE) and anionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG) lipids. Chain-perdeuterated POPE and POPG are used in single-component membranes, binary POPE/POPG (3:1) membranes, and membranes containing one of four cationic peptides: two antimicrobial peptides (AMPs) of the β-hairpin family of protegrin-1 (PG-1), and two cell-penetrating peptides (CPPs), HIV TAT and penetratin. ²H quadrupolar couplings were measured to determine the motional amplitudes of POPE and POPG acyl chains as a function of temperature. Homogeneously mixed POPE/POPG membranes should give the same quadrupolar couplings for the two lipids, whereas the presence of membrane domains enriched in one of the two lipids should cause distinct ²H quadrupolar couplings that reflect different chain disorder. At physiological temperature (308 K), we observed no or only small coupling differences between POPE and POPG in the presence of any of the cationic peptides. However, around ambient temperature (293 K), at which gel- and liquid-crystalline phases coexist in the peptide-free POPE/POPG membrane, the peptides caused distinct quadrupolar couplings for the two lipids, indicating domain formation. The broad-spectrum antimicrobial peptide PG-1 ordered ∼40% of the POPE lipids while disordering POPG. The Gram-negative selective PG-1 mutant, IB549, caused even larger differences in the POPE and POPG disorder: ∼80% of POPE partitioned into the ordered phase, whereas all of the POPG remained in the disordered phase. In comparison, TAT rigidified POPE and POPG similarly in the binary membrane at ambient temperature, indicating that TAT does not cause dynamic heterogeneity but interacts with the membrane with a different mechanism. Penetratin maintained the POPE order but disordered POPG, suggesting moderate domain separation. These results provide insight into the extent of domain formation in bacterial membranes and the possible peptide structural requirements for this phenomenon.     

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.     

Analysis of Structural Variability in Pharmaceutical Excipients Using Solid-State NMR Spectroscopy

Polysaccharide-based excipients comprise the majority of most solid dosage forms and can vary dramatically in terms of structural and functionally related properties. Analytical methods for characterizing these important formulation components are crucial. Solid-state NMR spectroscopy (SSNMR) can provide a wealth of information on these materials while offering the advantages of non-destructive sample preparation and selectivity. The overall objective of this work is to identify SSNMR parameters that can be used to detect differences among these excipients. Excipients were obtained from a wide range of suppliers and analyzed as received; ¹³C SSNMR spectra were acquired using a Chemagnetics CMX-300 spectrometer operating at approximately 75 MHz. The resolution of SSNMR signals of many excipients allows for positive identification of the major form present. Alginic acid and sodium alginate can be differentiated based on carbonyl peak position. Analysis of relative peak intensities provides insight into the purity of a carrageenan sample compared to known standards. The SSNMR spectrum of starch can be used to identify the source and to quantitate the amorphous and crystalline content. Relaxation values and peak areas of starch derivatives can be related to the degree of hydrolysis, providing an alternative method for determining dextrose equivalent. Differences in peak intensities and relaxation time values of HPMC samples can be correlated to the amount of methoxy subsituent groups. Important characteristics of excipients such as form identification, structural differences, crystalline and amorphous content, and water content variations can be detected using SSNMR spectroscopy.     

Ligand Induced Conformational Changes of the Human Serotonin Transporter Revealed by Molecular Dynamics Simulations

The competitive inhibitor cocaine and the non-competitive inhibitor ibogaine induce different conformational states of the human serotonin transporter. It has been shown from accessibility experiments that cocaine mainly induces an outward-facing conformation, while the non-competitive inhibitor ibogaine, and its active metabolite noribogaine, have been proposed to induce an inward-facing conformation of the human serotonin transporter similar to what has been observed for the endogenous substrate, serotonin. The ligand induced conformational changes within the human serotonin transporter caused by these three different types of ligands, substrate, non-competitive and competitive inhibitors, are studied from multiple atomistic molecular dynamics simulations initiated from a homology model of the human serotonin transporter. The results reveal that diverse conformations of the human serotonin transporter are captured from the molecular dynamics simulations depending on the type of the ligand bound. The inward-facing conformation of the human serotonin transporter is reached with noribogaine bound, and this state resembles a previously identified inward-facing conformation of the human serotonin transporter obtained from molecular dynamics simulation with bound substrate, but also a recently published inward-facing conformation of a bacterial homolog, the leucine transporter from Aquifex Aoelicus. The differences observed in ligand induced behavior are found to originate from different interaction patterns between the ligands and the protein. Such atomic-level understanding of how an inhibitor can dictate the conformational response of a transporter by ligand binding may be of great importance for future drug design.