Upgraded molecular models of the human KCNQ1 potassium channel.

The voltage-gated potassium channel KCNQ1 (KV7.1) assembles with the KCNE1 accessory protein to generate the slow delayed rectifier current, IKS, which is critical for membrane repolarization as part of the cardiac action potential. Loss-of-function (LOF) mutations in KCNQ1 are the most common cause of congenital long QT syndrome (LQTS), type 1 LQTS, an inherited genetic predisposition to cardiac arrhythmia and sudden cardiac death. A detailed structural understanding of KCNQ1 is needed to elucidate the molecular basis for KCNQ1 LOF in disease and to enable structure-guided design of new anti-arrhythmic drugs. In this work, advanced structural models of human KCNQ1 in the resting/closed and activated/open states were developed by Rosetta homology modeling guided by newly available experimentally-based templates: X. leavis KCNQ1 and various resting voltage sensor structures. Using molecular dynamics (MD) simulations, the capacity of the models to describe experimentally established channel properties including state-dependent voltage sensor gating charge interactions and pore conformations, PIP2 binding sites, and voltage sensor-pore domain interactions were validated. Rosetta energy calculations were applied to assess the utility of each model in interpreting mutation-evoked KCNQ1 dysfunction by predicting the change in protein thermodynamic stability for 50 experimentally characterized KCNQ1 variants with mutations located in the voltage-sensing domain. Energetic destabilization was successfully predicted for folding-defective KCNQ1 LOF mutants whereas wild type-like mutants exhibited no significant energetic frustrations, which supports growing evidence that mutation-induced protein destabilization is an especially common cause of KCNQ1 dysfunction. The new KCNQ1 Rosetta models provide helpful tools in the study of the structural basis for KCNQ1 function and can be used to generate hypotheses to explain KCNQ1 dysfunction.

Structural, functional, and behavioral insights of dopamine dysfunction revealed by a deletion in SLC6A3.

The human dopamine (DA) transporter (hDAT) mediates clearance of DA. Genetic variants in hDAT have been associated with DA dysfunction, a complication associated with several brain disorders, including autism spectrum disorder (ASD). Here, we investigated the structural and behavioral bases of an ASD-associated in-frame deletion in hDAT at N336 (∆N336). We uncovered that the deletion promoted a previously unobserved conformation of the intracellular gate of the transporter, likely representing the rate-limiting step of the transport process. It is defined by a "half-open and inward-facing" state (HOIF) of the intracellular gate that is stabilized by a network of interactions conserved phylogenetically, as we demonstrated in hDAT by Rosetta molecular modeling and fine-grained simulations, as well as in its bacterial homolog leucine transporter by electron paramagnetic resonance analysis and X-ray crystallography. The stabilization of the HOIF state is associated both with DA dysfunctions demonstrated in isolated brains of Drosophila melanogaster expressing hDAT ∆N336 and with abnormal behaviors observed at high-time resolution. These flies display increased fear, impaired social interactions, and locomotion traits we associate with DA dysfunction and the HOIF state. Together, our results describe how a genetic variation causes DA dysfunction and abnormal behaviors by stabilizing a HOIF state of the transporter.

Mechanisms of KCNQ1 channel dysfunction in long QT syndrome involving voltage sensor domain mutations.

Mutations that induce loss of function (LOF) or dysfunction of the human KCNQ1 channel are responsible for susceptibility to a life-threatening heart rhythm disorder, the congenital long QT syndrome (LQTS). Hundreds of KCNQ1 mutations have been identified, but the molecular mechanisms responsible for impaired function are poorly understood. We investigated the impact of 51 KCNQ1 variants with mutations located within the voltage sensor domain (VSD), with an emphasis on elucidating effects on cell surface expression, protein folding, and structure. For each variant, the efficiency of trafficking to the plasma membrane, the impact of proteasome inhibition, and protein stability were assayed. The results of these experiments combined with channel functional data provided the basis for classifying each mutation into one of six mechanistic categories, highlighting heterogeneity in the mechanisms resulting in channel dysfunction or LOF. More than half of the KCNQ1 LOF mutations examined were seen to destabilize the structure of the VSD, generally accompanied by mistrafficking and degradation by the proteasome, an observation that underscores the growing appreciation that mutation-induced destabilization of membrane proteins may be a common human disease mechanism. Finally, we observed that five of the folding-defective LQTS mutant sites are located in the VSD S0 helix, where they interact with a number of other LOF mutation sites in other segments of the VSD. These observations reveal a critical role for the S0 helix as a central scaffold to help organize and stabilize the KCNQ1 VSD and, most likely, the corresponding domain of many other ion channels.

Computational design of membrane proteins using RosettaMembrane.

Computational membrane protein design is challenging due to the small number of high-resolution structures available to elucidate the physical basis of membrane protein structure, multiple functionally important conformational states, and a limited number of high-throughput biophysical assays to monitor function. However, structural determination of membrane proteins has made tremendous progress in the past years. Concurrently the field of soluble computational design has made impressive inroads. These developments allow us to tackle the formidable challenge of designing functional membrane proteins. Herein, Rosetta is benchmarked for membrane protein design. We evaluate strategies to cope with the often reduced quality of experimental membrane protein structures. Further, we test the usage of symmetry in design protocols, which is particularly important as many membrane proteins exist as homo-oligomers. We compare a soluble scoring function with a scoring function optimized for membrane proteins, RosettaMembrane. Both scoring functions recovered around half of the native sequence when completely redesigning membrane proteins. However, RosettaMembrane recovered the most native-like amino acid property composition. While leucine was overrepresented in the inner and outer-hydrophobic regions of RosettaMembrane designs, it resulted in a native-like surface hydrophobicity indicating that it is currently the best option for designing membrane proteins with Rosetta.

Documentation of an Imperative To Improve Methods for Predicting Membrane Protein Stability.

There is a compelling and growing need to accurately predict the impact of amino acid mutations on protein stability for problems in personalized medicine and other applications. Here the ability of 10 computational tools to accurately predict mutation-induced perturbation of folding stability (ΔΔG) for membrane proteins of known structure was assessed. All methods for predicting ΔΔG values performed significantly worse when applied to membrane proteins than when applied to soluble proteins, yielding estimated concordance, Pearson, and Spearman correlation coefficients of <0.4 for membrane proteins. Rosetta and PROVEAN showed a modest ability to classify mutations as destabilizing (ΔΔG < -0.5 kcal/mol), with a 7 in 10 chance of correctly discriminating a randomly chosen destabilizing variant from a randomly chosen stabilizing variant. However, even this performance is significantly worse than for soluble proteins. This study highlights the need for further development of reliable and reproducible methods for predicting thermodynamic folding stability in membrane proteins.

Protocols for Molecular Modeling with Rosetta3 and RosettaScripts.

Previously, we published an article providing an overview of the Rosetta suite of biomacromolecular modeling software and a series of step-by-step tutorials [Kaufmann, K. W., et al. (2010) Biochemistry 49, 2987-2998]. The overwhelming positive response to this publication we received motivates us to here share the next iteration of these tutorials that feature de novo folding, comparative modeling, loop construction, protein docking, small molecule docking, and protein design. This updated and expanded set of tutorials is needed, as since 2010 Rosetta has been fully redesigned into an object-oriented protein modeling program Rosetta3. Notable improvements include a substantially improved energy function, an XML-like language termed "RosettaScripts" for flexibly specifying modeling task, new analysis tools, the addition of the TopologyBroker to control conformational sampling, and support for multiple templates in comparative modeling. Rosetta's ability to model systems with symmetric proteins, membrane proteins, noncanonical amino acids, and RNA has also been greatly expanded and improved.

An Integrated Framework Advancing Membrane Protein Modeling and Design.

Membrane proteins are critical functional molecules in the human body, constituting more than 30% of open reading frames in the human genome. Unfortunately, a myriad of difficulties in overexpression and reconstitution into membrane mimetics severely limit our ability to determine their structures. Computational tools are therefore instrumental to membrane protein structure prediction, consequently increasing our understanding of membrane protein function and their role in disease. Here, we describe a general framework facilitating membrane protein modeling and design that combines the scientific principles for membrane protein modeling with the flexible software architecture of Rosetta3. This new framework, called RosettaMP, provides a general membrane representation that interfaces with scoring, conformational sampling, and mutation routines that can be easily combined to create new protocols. To demonstrate the capabilities of this implementation, we developed four proof-of-concept applications for (1) prediction of free energy changes upon mutation; (2) high-resolution structural refinement; (3) protein-protein docking; and (4) assembly of symmetric protein complexes, all in the membrane environment. Preliminary data show that these algorithms can produce meaningful scores and structures. The data also suggest needed improvements to both sampling routines and score functions. Importantly, the applications collectively demonstrate the potential of combining the flexible nature of RosettaMP with the power of Rosetta algorithms to facilitate membrane protein modeling and design.

Inverted topologies in membrane proteins: a mini-review.

Helical membrane proteins such as transporters, receptors, or channels often exhibit structural symmetry. Symmetry is perfect in homo-oligomers consisting of two or more copies of the same protein chain. Intriguingly, in single chain membrane proteins, often internal pseudo-symmetry is observed, in particular in transporters and channels. In several cases single chain proteins with pseudo-symmetry exist, that share the fold with homo-oligomers suggesting evolutionary pathways that involve gene duplication and fusion. It has been hypothesized that such evolutionary pathways allow for the rapid development of large proteins with novel functionality. At the same time symmetry can be leveraged to recognize highly symmetric substrates such as ions. Here we review helical transporter proteins with an inverted two-fold pseudo-symmetry. In this special scenario the symmetry axis lies in the membrane plane. As a result, the putative ancestral monomeric protein would insert in both directions into the membrane and its open-to-the-inside and open-to-the-outside conformations would be structurally identical and iso-energetic, giving a possible evolutionary pathway to create a transporter protein that needs to flip between the two states.