Structural refinement of membrane proteins by restrained molecular dynamics and solvent accessibility data

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

A novel approach to modeling side chain ensembles of the bifunctional spin label RX

We introduce a novel approach to modeling side chain ensembles of bifunctional spin labels. This approach utilizes rotamer libraries to generate side chain conformational ensembles. Because the bifunctional label is constrained by two attachment sites, the label is split into two monofunctional rotamers which are first attached to their respective sites, then rejoined by a local optimization in dihedral space. We validate this method against a set of previously published experimental data using the bifunctional spin label, RX. This method is relatively fast and can readily be used for both experimental analysis and protein modeling, providing significant advantages over modeling bifunctional labels with molecular dynamics simulations. Use of bifunctional labels for site directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy dramatically reduces label mobility, which can significantly improve resolution of small changes in protein backbone structure and dynamics. Coupling the use of bifunctional labels with side chain modeling methods allows for improved quantitative application of experimental SDSL EPR data to protein modeling.

Hybrid Refinement of Heterogeneous Conformational Ensembles Using Spectroscopic Data

Multistructured biomolecular systems play crucial roles in a wide variety of cellular processes but have resisted traditional methods of structure determination, which often resolve only a few low-energy states. High-resolution structure determination using experimental methods that yield distributional data remains extremely difficult, especially when the underlying conformational ensembles are quite heterogeneous. We have therefore developed a method to integrate sparse, multimultimodal spectroscopic data to obtain high-resolution estimates of conformational ensembles. We have tested our method by incorporating double electron-electron resonance data on the soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) protein syntaxin-1a into biased molecular dynamics simulations. We find that our method substantially outperforms existing state-of-the-art methods in capturing syntaxin's open-closed conformational equilibrium and further yields new conformational states that are consistent with experimental data and may help in understanding syntaxin's function. Our improved methods for refining heterogeneous conformational ensembles from spectroscopic data will greatly accelerate the structural understanding of such systems.

Resting state of the human proton channel dimer in a lipid bilayer

The voltage-gated proton channel Hv1 plays a critical role in the fast proton translocation that underlies a wide range of physiological functions, including the phagocytic respiratory burst, sperm motility, apoptosis, and metastatic cancer. Both voltage activation and proton conduction are carried out by a voltage-sensing domain (VSD) with strong similarity to canonical VSDs in voltage-dependent cation channels and enzymes. We set out to determine the structural properties of membrane-reconstituted human proton channel (hHv1) in its resting conformation using electron paramagnetic resonance spectroscopy together with biochemical and computational methods. We evaluated existing structural templates and generated a spectroscopically constrained model of the hHv1 dimer based on the Ci-VSD structure at resting state. Mapped accessibility data revealed deep water penetration through hHv1, suggesting a highly focused electric field, comprising two turns of helix along the fourth transmembrane segment. This region likely contains the H(+) selectivity filter and the conduction pore. Our 3D model offers plausible explanations for existing electrophysiological and biochemical data, offering an explicit mechanism for voltage activation based on a one-click sliding helix conformational rearrangement.

Structural Refinement of Proteins by Restrained Molecular Dynamics Simulations with Non-interacting Molecular Fragments

The knowledge of multiple conformational states is a prerequisite to understand the function of membrane transport proteins. Unfortunately, the determination of detailed atomic structures for all these functionally important conformational states with conventional high-resolution approaches is often difficult and unsuccessful. In some cases, biophysical and biochemical approaches can provide important complementary structural information that can be exploited with the help of advanced computational methods to derive structural models of specific conformational states. In particular, functional and spectroscopic measurements in combination with site-directed mutations constitute one important source of information to obtain these mixed-resolution structural models. A very common problem with this strategy, however, is the difficulty to simultaneously integrate all the information from multiple independent experiments involving different mutations or chemical labels to derive a unique structural model consistent with the data. To resolve this issue, a novel restrained molecular dynamics structural refinement method is developed to simultaneously incorporate multiple experimentally determined constraints (e.g., engineered metal bridges or spin-labels), each treated as an individual molecular fragment with all atomic details. The internal structure of each of the molecular fragments is treated realistically, while there is no interaction between different molecular fragments to avoid unphysical steric clashes. The information from all the molecular fragments is exploited simultaneously to constrain the backbone to refine a three-dimensional model of the conformational state of the protein. The method is illustrated by refining the structure of the voltage-sensing domain (VSD) of the Kv1.2 potassium channel in the resting state and by exploring the distance histograms between spin-labels attached to T4 lysozyme. The resulting VSD structures are in good agreement with the consensus model of the resting state VSD and the spin-spin distance histograms from ESR/DEER experiments on T4 lysozyme are accurately reproduced.

Knowledge of multiple conformational states of membrane transport proteins is a prerequisite to understand their function. However, the determination of atomic structures for all these states with conventional high-resolution approaches can be very challenging due to inherent difficulties in high yield purification of functional membrane transport proteins. Various complementary structural information of proteins in their native states can be obtained by a variety of biophysical and biochemical methods with site-directed mutations. Here, a novel restrained molecular dynamics structural refinement method is developed to help derive a structural model that is consistent with experimental data by incorporating all the experimental constraints simultaneously through the use of non-interacting all-atom molecular fragments. The method can be easily and effectively extended to incorporate many kinds of structural constraints from a variety of biophysical and biochemical experiments, and should be very useful in generating and refining models of proteins in specific functional states.

Molecular mechanism of Mg2+-dependent gating in CorA

CorA is the major transport system responsible for Mg²⁺ uptake in bacteria and can functionally substitute for its homologue Mrs2p in the yeast inner mitochondrial membrane. Although several CorA crystal structures are available, the molecular mechanism of Mg²⁺ uptake remains to be established. Here we use EPR spectroscopy, electrophysiology and molecular dynamic simulations to show that CorA is regulated by cytoplasmic Mg²⁺ acting as a ligand and elucidate the basic conformational rearrangements responsible for Mg²⁺-dependent gating. Mg²⁺ unbinding at the divalent cation sensor triggers a conformational change that leads to the inward motion of the stalk helix, which propagates to the pore forming transmembrane helix TM1. Helical tilting and rotation in TM1 generates an iris-like motion that increases the diameter of the permeation pathway, triggering ion conduction. This work establishes the molecular basis of a Mg²⁺-driven negative feedback loop in CorA as the key physiological event controlling Mg²⁺ uptake and homeostasis in prokaryotes.

Structural basis of lipid-driven conformational transitions in the KvAP voltage-sensing domain

Voltage-gated ion channels respond to transmembrane electric fields through reorientations of the positively charged S4 helix within the voltage-sensing domain (VSD). Despite a wealth of structural and functional data, the details of this conformational change remain controversial. Recent electrophysiological evidence showed that equilibrium between the resting ('down') and activated ('up') conformations of the KvAP VSD from Aeropyrum pernix can be biased through reconstitution in lipids with or without phosphate groups. We investigated the structural transition between these functional states, using site-directed spin-labeling and EPR spectroscopic methods. Solvent accessibility and interhelical distance determinations suggest that KvAP gates through S4 movements involving an ∼3-Å upward tilt and simultaneous ∼2-Å axial shift. This motion leads to large accessibly changes in the intracellular water-filled crevice and supports a new model of gating that combines structural rearrangements and electric-field remodeling.

The voltage sensor module in sodium channels

The mechanism by which voltage-gated ion channels respond to changes in membrane polarization during action potential signaling in excitable cells has been the subject of research attention since the original description of voltage-dependent sodium and potassium flux in the squid giant axon. The cloning of ion channel genes and the identification of point mutations associated with channelopathy diseases in muscle and brain has facilitated an electrophysiological approach to the study of ion channels. Experimental approaches to the study of voltage gating have incorporated the use of thiosulfonate reagents to test accessibility, fluorescent probes, and toxins to define domain-specific roles of voltage-sensing S4 segments. Crystallography, structural and homology modeling, and molecular dynamics simulations have added computational approaches to study the relationship of channel structure to function. These approaches have tested models of voltage sensor translocation in response to membrane depolarization and incorporate the role of negative countercharges in the S1 to S3 segments to define our present understanding of the mechanism by which the voltage sensor module dictates gating particle permissiveness in excitable cells.

Influences of membrane mimetic environments on membrane protein structures

The number of membrane protein structures in the Protein Data Bank is becoming significant and growing. Here, the transmembrane domain structures of the helical membrane proteins are evaluated to assess the influences of the membrane mimetic environments. Toward this goal, many of the biophysical properties of membranes are discussed and contrasted with those of the membrane mimetics commonly used for structure determination. Although the mimetic environments can perturb the protein structures to an extent that potentially gives rise to misinterpretation of functional mechanisms, there are also many structures that have a native-like appearance. From this assessment, an initial set of guidelines is proposed for distinguishing native-like from nonnative-like membrane protein structures. With experimental techniques for validation and computational methods for refinement and quality assessment and enhancement, there are good prospects for achieving native-like structures for these very important proteins.

Structural investigation of MscL gating using experimental data and coarse grained MD simulations

The mechanosensitive channel of large conductance (MscL) has become a model system in which to understand mechanosensation, a process involved in osmoregulation and many other physiological functions. While a high resolution closed state structure is available, details of the open structure and the gating mechanism remain unknown. In this study we combine coarse grained simulations with restraints from EPR and FRET experiments to study the structural changes involved in gating with much greater level of conformational sampling than has previously been possible. We generated a set of plausible open pore structures that agree well with existing open pore structures and gating models. Most interestingly, we found that membrane thinning induces a kink in the upper part of TM1 that causes an outward motion of the periplasmic loop away from the pore centre. This previously unobserved structural change might present a new mechanism of tension sensing and might be related to a functional role in osmoregulation.

Cells in biological organisms have to be able to respond to mechanical forces during processes such as touch, hearing, pain sensation and tissue growth. One way this is achieved is through mechanosensitive ion channels, membrane embedded proteins that initiate electrical signalling upon tension within the cell or cell membrane. The malfunction of such channels is also associated with a range of diseases including muscular dystrophy and cardiac arrhythmia. In this manuscript, we study in detail the mechanosensitive channel of large conductance (MscL) from bacteria, a model system in which to understand the principles of mechanosensation. Despite many years of investigative work the details of how the protein senses tension in the surrounding membrane remain unknown. By combining structural data from experiments with computer simulation we are able to model the open channel structure of the protein and report previously unobserved structural changes that might present a new mechanism of sensing tension. The methods developed in this paper are not limited to the study of mechanosensitive ion channels and may be useful in understanding the structure and function of other membrane proteins.

Toward the fourth dimension of membrane protein structure: insight into dynamics from spin-labeling EPR spectroscopy

Trapping membrane proteins in the confines of a crystal lattice obscures dynamic modes essential for interconversion between multiple conformations in the functional cycle. Moreover, lattice forces could conspire with detergent solubilization to stabilize a minor conformer in an ensemble thus confounding mechanistic interpretation. Spin labeling in conjunction with electron paramagnetic resonance (EPR) spectroscopy offers an exquisite window into membrane protein dynamics in the native-like environment of a lipid bilayer. Systematic application of spin labeling and EPR identifies sequence-specific secondary structures, defines their topology and their packing in the tertiary fold. Long range distance measurements (60 Å-80 Å) between pairs of spin labels enable quantitative analysis of equilibrium dynamics and triggered conformational changes. This review highlights the contribution of spin labeling to bridging structure and mechanism. Efforts to develop methods for determining structures from EPR restraints and to increase sensitivity and throughput promise to expand spin labeling applications in membrane protein structural biology.

Structural dynamics of the magnesium-bound conformation of CorA in a lipid bilayer

The transmembrane conformation of Thermotoga maritima CorA, a magnesium transport system, has been studied in its Mg(2+)-bound form by site-directed spin labeling and electron paramagnetic resonance spectroscopy. Probe mobility together with accessibility data were used to evaluate the overall dynamics and relative arrangement of individual transmembrane segments TM1 and TM2. TM1 extends toward the cytoplasmic side creating a water-filled cavity, while TM2 is located in the periphery of the oligomer, contacting the lipid bilayer. A structural model for the conserved extracellular loop was generated based on EPR data and MD simulations, in which residue E316 is located toward the five-fold symmetry axis in position to electrostatically influence divalent ion translocation. Electrostatic analysis of our model suggest that, in agreement with the crystal structure, Mg(2+) -bound CorA is in a closed conformation. The present results suggest that long-range structural rearrangements are necessary to allow Mg(2+) translocation.

The activated state of a sodium channel voltage sensor in a membrane environment

Direct structural insights on the fundamental mechanisms of permeation, selectivity, and gating remain unavailable for the Na(+) and Ca(2+) channel families. Here, we report the spectroscopic structural characterization of the isolated Voltage-Sensor Domain (VSD) of the prokaryotic Na(+) channel NaChBac in a lipid bilayer. Site-directed spin-labeling and EPR spectroscopy were carried out for 118 mutants covering all of the VSD. EPR environmental data were used to unambiguously assign the secondary structure elements, define membrane insertion limits, and evaluate the activated conformation of the isolated-VSD in the membrane using restrain-driven molecular dynamics simulations. The overall three-dimensional fold of the NaChBac-VSD closely mirrors those seen in KvAP, Kv1.2, Kv1.2-2.1 chimera, and MlotiK1. However, in comparison to the membrane-embedded KvAP-VSD, the structural dynamics of the NaChBac-VSD reveals a much tighter helix packing, with subtle differences in the local environment of the gating charges and their interaction with the rest of the protein. Using cell complementation assays we show that the NaChBac-VSD can provide a conduit to the transport of ions in the resting or "down" conformation, a feature consistent with our EPR water accessibility measurements in the activated or "up" conformation. These results suggest that the overall architecture of VSD's is remarkably conserved among K(+) and Na(+) channels and that pathways for gating-pore currents may be intrinsic to most voltage-sensors. Cell complementation assays also provide information about the putative location of the gating charges in the "down/resting" state and hence a glimpse of the extent of conformational changes during activation.

Structural studies of the Enterococcus faecalis SufU [Fe-S] cluster protein

Background

Iron-sulfur clusters are ubiquitous and evolutionarily ancient inorganic prosthetic groups, the biosynthesis of which depends on complex protein machineries. Three distinct assembly systems involved in the maturation of cellular Fe-S proteins have been determined, designated the NIF, ISC and SUF systems. Although well described in several organisms, these machineries are poorly understood in Gram-positive bacteria. Within the Firmicutes phylum, the Enterococcus spp. genus have recently assumed importance in clinical microbiology being considered as emerging pathogens for humans, wherein Enterococcus faecalis represents the major species associated with nosocomial infections. The aim of this study was to carry out a phylogenetic analysis in Enterococcus faecalis V583 and a structural and conformational characterisation of it SufU protein.

Results

BLAST searches of the Enterococcus genome revealed a series of genes with sequence similarity to the Escherichia coli SUF machinery of [Fe-S] cluster biosynthesis, namely sufB, sufC, sufD and SufS. In addition, the E. coli IscU ortholog SufU was found to be the scaffold protein of Enterococcus spp., containing all features considered essential for its biological activity, including conserved amino acid residues involved in substrate and/or co-factor binding (Cys^50,76,138 and Asp⁵²) and, phylogenetic analyses showed a close relationship with orthologues from other Gram-positive bacteria. Molecular dynamics for structural determinations and molecular modeling using E. faecalis SufU primary sequence protein over the PDB:1su0 crystallographic model from Streptococcus pyogenes were carried out with a subsequent 50 ns molecular dynamic trajectory. This presented a stable model, showing secondary structure modifications near the active site and conserved cysteine residues. Molecular modeling using Haemophilus influenzae IscU primary sequence over the PDB:1su0 crystal followed by a MD trajectory was performed to analyse differences in the C-terminus region of Gram-positive SufU and Gram-negative orthologous proteins, in which several modifications in secondary structure were observed.

Conclusion

The data describe the identification of the SUF machinery for [Fe-S] cluster biosynthesis present in the Firmicutes genome, showing conserved sufB, sufC, sufD and sufS genes and the presence of the sufU gene coding for scaffold protein, instead of sufA; neither sufE nor sufR are present. Primary sequences and structural analysis of the SufU protein demonstrated its structural-like pattern to the scaffold protein IscU nearby on the ISC machinery. E. faecalis SufU molecular modeling showed high flexibility over the active site regions, and demonstrated the existence of a specific region in Firmicutes denoting the Gram positive region (GPR), suggested as a possible candidate for interaction with other factors and/or regulators.

A structural mechanism for MscS gating in lipid bilayers

The mechanosensitive channel of small conductance (MscS) is a key determinant in the prokaryotic response to osmotic challenges. We determined the structural rearrangements associated with MscS activation in membranes, using functorial measurements, electron paramagnetic resonance spectroscopy, and computational analyses. MscS was trapped in its open conformation after the transbilayer pressure profile was modified through the asymmetric incorporation of lysophospholipids. The transition from the closed to the open state is accompanied by the downward tilting of the transmembrane TM1-TM2 hairpin and by the expansion, tilt, and rotation of the TM3 helices. These movements expand the permeation pathway, leading to an increase in accessibility to water around TM3. Our open MscS model is compatible with single-channel conductance measurements and supports the notion that helix tilting is associated with efficient pore widening in mechanosensitive channels.