ABSTRACT
Membrane remodeling and repair are essential for all cells. Proteins that perform these functions include Vipp1/IM30 in photosynthetic plastids, PspA in bacteria, and ESCRT-III in eukaryotes. Here, using a combination of evolutionary and structural analyses, we show that these protein families are homologous and share a common ancient evolutionary origin that likely predates the last universal common ancestor. This homology is evident in cryo-electron microscopy structures of Vipp1 rings from the cyanobacterium Nostoc punctiforme presented over a range of symmetries. Each ring is assembled from rungs that stack and progressively tilt to form dome-shaped curvature. Assembly is facilitated by hinges in the Vipp1 monomer, similar to those in ESCRT-III proteins, which allow the formation of flexible polymers. Rings have an inner lumen that is able to bind and deform membranes. Collectively, these data suggest conserved mechanistic principles that underlie Vipp1, PspA, and ESCRT-III-dependent membrane remodeling across all domains of life.
Subject(s)
Bacterial Proteins/metabolism , Cell Membrane/metabolism , Endosomal Sorting Complexes Required for Transport/metabolism , Heat-Shock Proteins/metabolism , Multigene Family , Nostoc/metabolism , Amino Acid Sequence , Animals , Bacterial Proteins/chemistry , Bacterial Proteins/isolation & purification , Bacterial Proteins/ultrastructure , Chickens , Cryoelectron Microscopy , Endosomal Sorting Complexes Required for Transport/chemistry , Evolution, Molecular , Heat-Shock Proteins/chemistry , Heat-Shock Proteins/ultrastructure , Humans , Models, Molecular , Protein Structure, Secondary , Sequence Homology, Amino Acid , ThermodynamicsABSTRACT
Balanced fusion and fission are key for the proper function and physiology of mitochondria1,2. Remodelling of the mitochondrial inner membrane is mediated by the dynamin-like protein mitochondrial genome maintenance 1 (Mgm1) in fungi or the related protein optic atrophy 1 (OPA1) in animals3-5. Mgm1 is required for the preservation of mitochondrial DNA in yeast6, whereas mutations in the OPA1 gene in humans are a common cause of autosomal dominant optic atrophy-a genetic disorder that affects the optic nerve7,8. Mgm1 and OPA1 are present in mitochondria as a membrane-integral long form and a short form that is soluble in the intermembrane space. Yeast strains that express temperature-sensitive mutants of Mgm19,10 or mammalian cells that lack OPA1 display fragmented mitochondria11,12, which suggests that Mgm1 and OPA1 have an important role in inner-membrane fusion. Consistently, only the mitochondrial outer membrane-not the inner membrane-fuses in the absence of functional Mgm113. Mgm1 and OPA1 have also been shown to maintain proper cristae architecture10,14; for example, OPA1 prevents the release of pro-apoptotic factors by tightening crista junctions15. Finally, the short form of OPA1 localizes to mitochondrial constriction sites, where it presumably promotes mitochondrial fission16. How Mgm1 and OPA1 perform their diverse functions in membrane fusion, scission and cristae organization is at present unknown. Here we present crystal and electron cryo-tomography structures of Mgm1 from Chaetomium thermophilum. Mgm1 consists of a GTPase (G) domain, a bundle signalling element domain, a stalk, and a paddle domain that contains a membrane-binding site. Biochemical and cell-based experiments demonstrate that the Mgm1 stalk mediates the assembly of bent tetramers into helical filaments. Electron cryo-tomography studies of Mgm1-decorated lipid tubes and fluorescence microscopy experiments on reconstituted membrane tubes indicate how the tetramers assemble on positively or negatively curved membranes. Our findings convey how Mgm1 and OPA1 filaments dynamically remodel the mitochondrial inner membrane.
Subject(s)
Chaetomium/chemistry , Cryoelectron Microscopy , Fungal Proteins/chemistry , Fungal Proteins/metabolism , GTP-Binding Proteins/chemistry , Mitochondrial Membranes/metabolism , Mitochondrial Proteins/chemistry , Crystallography, X-Ray , Fungal Proteins/ultrastructure , GTP-Binding Proteins/metabolism , GTP-Binding Proteins/ultrastructure , Galactosylceramides/metabolism , Mitochondrial Proteins/metabolism , Mitochondrial Proteins/ultrastructure , Models, Molecular , Protein Domains , Protein MultimerizationABSTRACT
Dynamin oligomerizes into helical filaments on tubular membrane templates and, through constriction, cleaves them in a GTPase-driven way. Structural observations of GTP-dependent cross-bridges between neighboring filament turns have led to the suggestion that dynamin operates as a molecular ratchet motor. However, the proof of such mechanism remains absent. Particularly, it is not known whether a powerful enough stroke is produced and how the motor modules would cooperate in the constriction process. Here, we characterized the dynamin motor modules by single-molecule Förster resonance energy transfer (smFRET) and found strong nucleotide-dependent conformational preferences. Integrating smFRET with molecular dynamics simulations allowed us to estimate the forces generated in a power stroke. Subsequently, the quantitative force data and the measured kinetics of the GTPase cycle were incorporated into a model including both a dynamin filament, with explicit motor cross-bridges, and a realistic deformable membrane template. In our simulations, collective constriction of the membrane by dynamin motor modules, based on the ratchet mechanism, is directly reproduced and analyzed. Functional parallels between the dynamin system and actomyosin in the muscle are seen. Through concerted action of the motors, tight membrane constriction to the hemifission radius can be reached. Our experimental and computational study provides an example of how collective motor action in megadalton molecular assemblies can be approached and explicitly resolved.
Subject(s)
Dynamins/metabolism , Models, Biological , Biomechanical Phenomena , Dynamins/chemistry , Fluorescence Resonance Energy Transfer , Kinetics , Molecular Motor Proteins/chemistry , Molecular Motor Proteins/metabolism , Nucleotides/metabolism , Protein Domains , Protein Multimerization , SolutionsABSTRACT
Influenza hemagglutinin (HA) mediates viral entry into host cells through a large-scale conformational rearrangement at low pH that leads to fusion of the viral and endosomal membranes. Crystallographic and biochemical data suggest that a loop-to-coiled-coil transition of the B-loop region of HA is important for driving this structural rearrangement. However, the microscopic picture for this proposed "spring-loaded" movement is missing. In this study, we focus on understanding the transition of the B loop and perform a set of all-atom molecular dynamics simulations of the full B-loop trimeric structure with the CHARMM36 force field. The free-energy profile constructed from our simulations describes a B loop that stably folds half of the postfusion coiled coil in tens of microseconds, but the full coiled coil is unfavorable. A buried hydrophilic residue, Thr59, is implicated in destabilizing the coiled coil. Interestingly, this conserved threonine is the only residue in the B loop that strictly differentiates between the group 1 and 2 HA molecules. Microsecond-scale constant temperature simulations revealed that kinetic traps in the structural switch of the B loop can be caused by nonnative, intramonomer, or intermonomer ß-sheets. The addition of the A helix stabilized the postfusion state of the B loop, but introduced the possibility for further ß-sheet structures. Overall, our results do not support a description of the B loop in group 2 HAs as a stiff spring, but, rather, it allows for more structural heterogeneity in the placement of the fusion peptides during the fusion process.
Subject(s)
Hemagglutinin Glycoproteins, Influenza Virus/chemistry , Influenza A virus/chemistry , Molecular Dynamics Simulation , Hemagglutinin Glycoproteins, Influenza Virus/metabolism , Influenza A virus/metabolism , Protein Structure, Quaternary , Protein Structure, SecondaryABSTRACT
With the broad accessibility of high-performance computing resources, the significance of a molecular dynamics simulation is now rarely limited by hardware and/or software availability. Rather, the scientific value of each calculation is determined by the principles that underlie the theoretical model. The current review addresses this topic in the context of simplified models applied to large-scale (â¼20-100â¯Å) dynamics in the ribosome. Specifically, we focus on applications of the "SMOG" class of structure-based models, which can be used to simulate spontaneous (i.e. non-targeted) conformational rearrangements in complex molecular assemblies. Here, we aim to provide an entry-level assessment of the methods, which can help bridge conceptual and communication gaps between the experimental and computational communities. In addition, inspecting the strategies that have been deployed previously can provide guidelines for future computational investigations into the relationship between structure, energetics and dynamics in other assemblies.
Subject(s)
Computational Biology/methods , Molecular Dynamics Simulation , RNA, Transfer/metabolism , Ribosomes/metabolism , Molecular Conformation , Peptide Chain Elongation, Translational , SoftwareABSTRACT
Peripheral membrane proteins with intrinsic curvature can act both as sensors of membrane curvature and shape modulators of the underlying membranes. A well-studied example of such proteins is the mechanochemical GTPase dynamin, which assembles into helical filaments around membrane tubes and catalyzes their scission in a GTPase-dependent manner. It is known that the dynamin coat alone, without GTP, can constrict membrane tubes to radii of â¼10 nm, indicating that the intrinsic shape and elasticity of dynamin filaments should play an important role in membrane remodeling. However, molecular and dynamic understanding of the process is lacking. Here, we develop a dynamical polymer-chain model for a helical elastic filament bound on a deformable membrane tube of conserved mass, accounting for thermal fluctuations in the filament and lipid flows in the membrane. The model is based on the locally cylindrical helix approximation for dynamin. We obtain the elastic parameters of the dynamin filament by molecular dynamics simulations of its tetrameric building block and also from coarse-grained structure-based simulations of a 17-dimer filament. The results show that the stiffness of dynamin is comparable to that of the membrane. We determine equilibrium shapes of the filament and the membrane and find that mostly the pitch of the filament, not its radius, is sensitive to variations in membrane tension and stiffness. The close correspondence between experimental estimates of the inner tube radius and those predicted by the model suggests that dynamin's "stalk" region is responsible for its GTP-independent membrane-shaping ability. The model paves the way for future mesoscopic modeling of dynamin with explicit motor function.
Subject(s)
Cell Membrane/metabolism , Dynamins/metabolism , Models, Biological , Polymers/metabolism , Biomechanical Phenomena , Elasticity , Molecular Dynamics Simulation , Temperature , Time FactorsABSTRACT
Focal adhesions are dynamic constructs at the leading edge of migrating cells, linking them to the extracellular matrix and enabling force sensing and transmission. The lifecycle of a focal adhesion is a highly coordinated process involving spatial and temporal variations of protein composition, interaction, and cellular tension. The assembly of focal adhesions requires the recruitment and activation of vinculin. Vinculin is present in the cytoplasm in an autoinhibited conformation in which its tail is held pincerlike by its head domains, further stabilized by two high-affinity head-tail interfaces. Vinculin has binding sites for talin and F-actin, but effective binding requires vinculin activation to release its head-tail associations. In migrating cells, it has been shown that the locations of vinculin activation coincide with areas of high cellular tension, and that the highest recorded tensions across vinculin are associated with adhesion assembly. Here, we use a structure-based model to investigate vinculin activation by talin modulated by tensile force generated by transient associations with F-actin. We show that vinculin activation may proceed from an intermediate state stabilized by partial talin-vinculin association. There is a low-force regime and a high-force regime where vinculin activation is dominated by two different pathways with distinct responses to force. Specifically, at zero or low forces, vinculin activation requires substantial destabilization of the main head-tail interface, which is rigid and undergoes very limited fluctuations, despite the other being relatively flexible. This pathway is not significantly affected by force; instead, higher forces favor an alternative pathway, which seeks to release the vinculin tail from its pincerlike head domains before destabilizing the head-tail interfaces. This pathway has a force-sensitive activation barrier and is significantly accelerated by force. Experimental data of vinculin during various stages of the focal adhesion lifecycle are consistent with the proposed force-regulated activation pathway.
Subject(s)
Focal Adhesions/metabolism , Vinculin/metabolism , Actins/chemistry , Actins/metabolism , Biomechanical Phenomena , Computer Simulation , Cytoplasm/chemistry , Cytoplasm/metabolism , Focal Adhesions/chemistry , Humans , Kinetics , Models, Molecular , Protein Domains , Protein Stability , Protein Structure, Secondary , Talin/chemistry , Talin/metabolism , Vinculin/chemistryABSTRACT
Molecular dynamics simulations with coarse-grained or simplified Hamiltonians have proven to be an effective means of capturing the functionally important long-time and large-length scale motions of proteins and RNAs. Originally developed in the context of protein folding, structure-based models (SBMs) have since been extended to probe a diverse range of biomolecular processes, spanning from protein and RNA folding to functional transitions in molecular machines. The hallmark feature of a structure-based model is that part, or all, of the potential energy function is defined by a known structure. Within this general class of models, there exist many possible variations in resolution and energetic composition. SMOG 2 is a downloadable software package that reads user-designated structural information and user-defined energy definitions, in order to produce the files necessary to use SBMs with high performance molecular dynamics packages: GROMACS and NAMD. SMOG 2 is bundled with XML-formatted template files that define commonly used SBMs, and it can process template files that are altered according to the needs of each user. This computational infrastructure also allows for experimental or bioinformatics-derived restraints or novel structural features to be included, e.g. novel ligands, prosthetic groups and post-translational/transcriptional modifications. The code and user guide can be downloaded at http://smog-server.org/smog2.
Subject(s)
Algorithms , Models, Chemical , Molecular Dynamics Simulation , Proteins/chemistry , Proteins/ultrastructure , Software , Protein Conformation , Software Design , Software ValidationABSTRACT
Influenza hemagglutinin (HA), a homotrimeric glycoprotein crucial for membrane fusion, undergoes a large-scale structural rearrangement during viral invasion. X-ray crystallography has shown that the pre- and postfusion configurations of HA2, the membrane-fusion subunit of HA, have disparate secondary, tertiary, and quaternary structures, where some regions are displaced by more than 100 Å. To explore structural dynamics during the conformational transition, we studied simulations of a minimally frustrated model based on energy landscape theory. The model combines structural information from both the pre- and postfusion crystallographic configurations of HA2. Rather than a downhill drive toward formation of the central coiled-coil, we discovered an order-disorder transition early in the conformational change as the mechanism for the release of the fusion peptides from their burial sites in the prefusion crystal structure. This disorder quickly leads to a metastable intermediate with a broken threefold symmetry. Finally, kinetic competition between the formation of the extended coiled-coil and C-terminal melting results in two routes from this intermediate to the postfusion structure. Our study reiterates the roles that cracking and disorder can play in functional molecular motions, in contrast to the downhill mechanical interpretations of the "spring-loaded" model proposed for the HA2 conformational transition.
Subject(s)
Hemagglutinin Glycoproteins, Influenza Virus/chemistry , Molecular Dynamics SimulationABSTRACT
RfaH is a virulence factor from Escherichia coli whose C-terminal domain (CTD) undergoes a dramatic α-to-ß conformational transformation. The CTD in its α-helical fold is stabilized by interactions with the N-terminal domain (NTD), masking an RNA polymerase binding site until a specific recruitment site is encountered. Domain dissociation is triggered upon binding to DNA, allowing the NTD to interact with RNA polymerase to facilitate transcription while the CTD refolds into the ß-barrel conformation that interacts with the ribosome to activate translation. However, structural details of this transformation process in the context of the full protein remain to be elucidated. Here, we explore the mechanism of the α-to-ß conformational transition of RfaH in the full-length protein using a dual-basin structure-based model. Our simulations capture several features described experimentally, such as the requirement of disruption of interdomain contacts to trigger the α-to-ß transformation, confirms the roles of previously indicated residues E48 and R138, and suggests a new important role for F130, in the stability of the interdomain interaction. These native basins are connected through an intermediate state that builds up upon binding to the NTD and shares features from both folds, in agreement with previous in silico studies of the isolated CTD. We also examine the effect of RNA polymerase binding on the stabilization of the ß fold. Our study shows that native-biased models are appropriate for interrogating the detailed mechanisms of structural rearrangements during the dramatic transformation process of RfaH.
Subject(s)
DNA-Directed RNA Polymerases/chemistry , DNA-Directed RNA Polymerases/ultrastructure , Escherichia coli Proteins/chemistry , Escherichia coli Proteins/ultrastructure , Models, Chemical , Molecular Dynamics Simulation , Peptide Elongation Factors/chemistry , Peptide Elongation Factors/ultrastructure , Trans-Activators/chemistry , Trans-Activators/ultrastructure , Binding Sites , Computer Simulation , Protein Binding , Protein Conformation , Protein Folding , Protein Structure, TertiaryABSTRACT
RNA is highly sensitive to the ionic environment and typically requires Mg(2+) to form compact structures. There is a need for models capable of describing the ion atmosphere surrounding RNA with quantitative accuracy. We present a model of RNA electrostatics and apply it within coarse-grained molecular dynamics simulation. The model treats Mg(2+) ions explicitly to account for ion-ion correlations neglected by mean-field theories. Since mean-field theories capture KCl well, it is treated implicitly by a generalized Manning counterion condensation model. The model extends Manning condensation to deal with arbitrary RNA conformations, nonlimiting KCl concentrations, and the ion inaccessible volume of RNA. The model is tested against experimental measurements of the excess Mg(2+) associated with the RNA, Γ(2+), because Γ(2+) is directly related to the Mg(2+)-RNA interaction free energy. The excellent agreement with experiment demonstrates that the model captures the ionic dependence of the RNA free energy landscape.
Subject(s)
Magnesium/chemistry , Models, Chemical , RNA/chemistry , Cations, Monovalent/chemistry , Nucleic Acid Conformation , Static ElectricityABSTRACT
A four-helix bundle is a well-characterized motif often used as a target for designed pharmaceutical therapeutics and nutritional supplements. Recently, we discovered a new structural complexity within this motif created by a disulphide bridge in the long-chain helical bundle cytokine leptin. When oxidized, leptin contains a disulphide bridge creating a covalent-loop through which part of the polypeptide chain is threaded (as seen in knotted proteins). We explored whether other proteins contain a similar intriguing knot-like structure as in leptin and discovered 11 structurally homologous proteins in the PDB. We call this new helical family class the Pierced Lasso Bundle (PLB) and the knot-like threaded structural motif a Pierced Lasso (PL). In the current study, we use structure-based simulation to investigate the threading/folding mechanisms for all the PLBs along with three unthreaded homologs as the covalent loop (or lasso) in leptin is important in folding dynamics and activity. We find that the presence of a small covalent loop leads to a mechanism where structural elements slipknot to thread through the covalent loop. Larger loops use a piercing mechanism where the free terminal plugs through the covalent loop. Remarkably, the position of the loop as well as its size influences the native state dynamics, which can impact receptor binding and biological activity. This previously unrecognized complexity of knot-like proteins within the helical bundle family comprises a completely new class within the knot family, and the hidden complexity we unraveled in the PLBs is expected to be found in other protein structures outside the four-helix bundles. The insights gained here provide critical new elements for future investigation of this emerging class of proteins, where function and the energetic landscape can be controlled by hidden topology, and should be take into account in ab initio predictions of newly identified protein targets.
Subject(s)
Models, Molecular , Protein Structure, Tertiary , Proteins/chemistry , Amino Acid Motifs , Animals , Computational Biology , Computer Simulation , Cysteine/chemistry , Cytokines/chemistry , Databases, Protein , Drug Design , Humans , Leptin/chemistry , Molecular Dynamics Simulation , Multiprotein Complexes/chemistry , Oxidation-Reduction , Protein Folding , Protein Structure, Secondary , Receptors, Leptin/chemistry , Structural Homology, ProteinABSTRACT
The diversity in a set of protein nuclear magnetic resonance (NMR) structures provides an estimate of native state fluctuations that can be used to refine and enrich structure-based protein models (SBMs). Dynamics are an essential part of a protein's functional native state. The dynamics in the native state are controlled by the same funneled energy landscape that guides the entire folding process. SBMs apply the principle of minimal frustration, drawn from energy landscape theory, to construct a funneled folding landscape for a given protein using only information from the native structure. On an energy landscape smoothed by evolution towards minimal frustration, geometrical constraints, imposed by the native structure, control the folding mechanism and shape the native dynamics revealed by the model. Native-state fluctuations can alternatively be estimated directly from the diversity in the set of NMR structures for a protein. Based on this information, we identify a highly flexible loop in the ribosomal protein S6 and modify the contact map in a SBM to accommodate the inferred dynamics. By taking into account the probable native state dynamics, the experimental transition state is recovered in the model, and the correct order of folding events is restored. Our study highlights how the shared energy landscape connects folding and function by showing that a better description of the native basin improves the prediction of the folding mechanism.
Subject(s)
Nuclear Magnetic Resonance, Biomolecular , Protein Folding , Ribosomal Protein S6/chemistry , Models, Molecular , Protein ConformationABSTRACT
Recent experiments have conclusively shown that proteins are able to fold from an unknotted, denatured polypeptide to the knotted, native state without the aid of chaperones. These experiments are consistent with a growing body of theoretical work showing that a funneled, minimally frustrated energy landscape is sufficient to fold small proteins with complex topologies. Here, we present a theoretical investigation of the folding of a knotted protein, 2ouf, engineered in the laboratory by a domain fusion that mimics an evolutionary pathway for knotted proteins. Unlike a previously studied knotted protein of similar length, we see reversible folding/knotting and a surprising lack of deep topological traps with a coarse-grained structure-based model. Our main interest is to investigate how evolution might further select the geometry and stiffness of the threading region of the newly fused protein. We compare the folding of the wild-type protein to several mutants. Similarly to the wild-type protein, all mutants show robust and reversible folding, and knotting coincides with the transition state ensemble. As observed experimentally, our simulations show that the knotted protein folds about ten times slower than an unknotted construct with an identical contact map. Simulated folding kinetics reflect the experimentally observed rollover in the folding limbs of chevron plots. Successful folding of the knotted protein is restricted to a narrow range of temperature as compared to the unknotted protein and fits of the kinetic folding data below folding temperature suggest slow, nondiffusive dynamics for the knotted protein.
Subject(s)
Protein Folding , Proteins/chemistry , Evolution, Molecular , Kinetics , Models, Molecular , Molecular Dynamics Simulation , Mutation , Proteins/metabolismABSTRACT
To reveal the molecular determinants of biological function, one seeks to characterize the interactions that are formed in conformational and chemical transition states. In other words, what interactions govern the molecule's energy landscape? To accomplish this, it is necessary to determine which degrees of freedom can unambiguously identify each transition state. Here, we perform simulations of large-scale aminoacyl-transfer RNA (aa-tRNA) rearrangements during accommodation on the ribosome and project the dynamics along experimentally accessible atomic distances. From this analysis, we obtain evidence for which coordinates capture the correct number of barrier-crossing events and accurately indicate when the aa-tRNA is on a transition path. Although a commonly used coordinate in single-molecule experiments performs poorly, this study implicates alternative coordinates along which rearrangements are accurately described as diffusive movements across a one-dimensional free-energy profile. From this, we provide the theoretical foundation required for single-molecule techniques to uncover the energy landscape governing aa-tRNA selection by the ribosome.
Subject(s)
Molecular Dynamics Simulation , RNA, Transfer/chemistry , Ribosomes/chemistry , Amino Acid Sequence , Base Sequence , Molecular Sequence Data , Nucleic Acid Conformation , Protein Conformation , RNA, Transfer/metabolism , Ribosomes/metabolismABSTRACT
Molecular dynamics simulations supplement single-molecule pulling experiments by providing the possibility of examining the full free energy landscape using many coordinates. Here, we use an all-atom structure-based model to study the force and temperature dependence of the unfolding of the protein filamin by applying force at both termini. The unfolding time-force relation τ(F) indicates that the force-induced unfolding behavior of filamin can be characterized into three regimes: barrier-limited low- and intermediate-force regimes, and a barrierless high-force regime. Slope changes of τ(F) separate the three regimes. We show that the behavior of τ(F) can be understood from a two-dimensional free energy landscape projected onto the extension X and the fraction of native contacts Q. In the low-force regime, the unfolding rate is roughly force-independent due to the small (even negative) separation in X between the native ensemble and transition state ensemble (TSE). In the intermediate-force regime, force sufficiently separates the TSE from the native ensemble such that τ(F) roughly follows an exponential relation. This regime is typically explored by pulling experiments. While X may fail to resolve the TSE due to overlap with the unfolded ensemble just below the folding temperature, the overlap is minimal at lower temperatures where experiments are likely to be conducted. The TSE becomes increasingly structured with force, whereas the average order of structural events during unfolding remains roughly unchanged. The high-force regime is characterized by barrierless unfolding, and the unfolding time approaches a limit of â¼10 µs for the highest forces we studied. Finally, a combination of X and Q is shown to be a good reaction coordinate for almost the entire force range.
Subject(s)
Filamins/chemistry , Protein Folding , Transition Temperature , Amino Acid Sequence , Molecular Dynamics Simulation , Molecular Sequence DataABSTRACT
The stability of RNA tertiary structures depends heavily on Mg(2+). The Mg(2+)-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ2+, the number of excess Mg(2+) associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg(2+) ions interact closely and strongly with the RNA, unlike monovalent ions such as K(+), suggesting that an explicit treatment of Mg(2+) is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg(2+)-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg(2+) through the explicit representation of Mg(2+) ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K(+) and models its competition with condensed Mg(2+). The model contains one fitted parameter, the number of condensed K(+) ions in the absence of Mg(2+). Values of Γ2+ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ2+ to the Mg(2+)-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg(2+) and RNA dynamics.
Subject(s)
Magnesium/chemistry , Models, Molecular , Nucleic Acid Conformation , Riboswitch , Base Sequence , Calibration , Molecular Dynamics Simulation , Molecular Sequence Data , Nonlinear Dynamics , Potassium/chemistry , Solvents/chemistry , ThermodynamicsABSTRACT
Proteins with a pierced lasso topology (PLT) have a covalent loop created by a disulfide bond, and the backbone circles back to thread the loop. This threaded topology has unique features compared to knotted topologies; notably, the topology is controlled by the chemical environment and the covalent loop remains intact even when denatured. In this work, we use the hormone leptin as our model PLT system and study its folding using molecular dynamics simulations that employ a structure-based (GoÌ -like) model. We find that the reduced protein has a two-state folding mechanism with a transition state ensemble (TSE) that can be characterized by the reaction coordinate Q, the fraction of native contacts formed. In contrast, the oxidized protein, which must thread part of the polypeptide chain through a covalent loop, has a folding process that is poorly characterized by Q. Instead, we find that a topological coordinate that monitors the residue crossing the loop can identify the TSE of oxidized leptin. By precisely identifying the predicted TSE, one may now reliably calculate theoretical phi-values for the PLT protein, thereby enabling a comparison with experimental measurements. We find the loop-threading constraint leads to noncanonical phi-values that are uniformly small because this PLT protein has a flat energy landscape through the TSE.
Subject(s)
Leptin , Protein Folding , Leptin/chemistry , Molecular Dynamics Simulation , Software , ThermodynamicsABSTRACT
Most proteins, in order to perform their biological function, have to fold to a compact native state. The increasing number of knotted and slipknotted proteins identified suggests that proteins are able to manoeuvre around topological barriers during folding. In the present article, we review the current progress in elucidating the knotting process in proteins. Although we concentrate on theoretical approaches, where a knotted topology can be unambiguously detected, comparison with experiments is also reviewed. Numerical simulations suggest that the folding process for small knotted proteins is composed of twisted loop formation and then threading by either slipknot geometries or flipping. As the size of the knotted proteins increases, particularly for more deeply threaded termini, the prevalence of traps in the free energy landscape also increases. Thus, in the case of longer knotted and slipknotted proteins, the folding mechanism is probably supported by chaperones. Overall, results imply that knotted proteins can be folded efficiently and survive evolutionary pressure in order to perform their biological functions.
Subject(s)
Proteins/chemistry , Animals , Humans , Protein Conformation , Protein Engineering , Protein Folding , Proteins/metabolism , ThermodynamicsABSTRACT
Energetic frustration in protein folding is minimized by evolution to create a smooth and robust energy landscape. As a result the geometry of the native structure provides key constraints that shape protein folding mechanisms. Chain connectivity in particular has been identified as an essential component for realistic behavior of protein folding models. We study the quantitative balance of energetic and geometrical influences on the folding of SH3 in a structure-based model with minimal energetic frustration. A decomposition of the two-dimensional free energy landscape for the folding reaction into relevant energy and entropy contributions reveals that the entropy of the chain is not responsible for the folding mechanism. Instead the preferred folding route through the transition state arises from a cooperative energetic effect. Off-pathway structures are penalized by excess distortion in local backbone configurations and contact pair distances. This energy cost is a new ingredient in the malleable balance of interactions that controls the choice of routes during protein folding.