Heiko Lammert

Robustness and generalization of structure‐based models for protein folding and function

Functional dynamics of native proteins share the energy landscape that guides folding into the native state. Folding simulations of structure‐based protein models, using an minimally frustrated energy landscape dominated by native interactions, can describe the geometrical aspects of the folding mechanism. Technical limitations imposed by the fixed shape of conventional contact potentials are a key obstacle toward advanced applications of structure‐based models like allostery or ligand binding, which require multiple stable conformations. Generalizations of existing models, commonly using Lennard‐Jones‐like potentials, lead to inevitable clashes between their repulsive branches. To resolve these challenges, a new contact potential is developed that combines an attractive part based on Gaussians with a separate repulsive term allowing flexibility for adjustments of the potential shape. With this new model multiple minima for studies of functional transitions can be introduced easily and consistently. A sensitivity analysis for five small proteins confirms the robust behavior of structure‐based models with our adaptable potential and explores their capacity for quantitative adjustment of the folding thermodynamics. We demonstrate its ability to incorporate alternative contact distances in simulations of structural transitions for the well‐studied ROP dimer. Individual contact pairs can switch between distinct states to match the competing syn and anti structures. The flexibility of the new potential facilitates advanced uses of structure‐based models. Depending on the application, features can be chosen from physical considerations or to match experiments. Generalized models can be built from multiple structures to study structural transitions or effects of disorder. Proteins 2009.

SMOG 2: A Versatile Software Package for Generating Structure-Based Models

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.

Allostery in the ferredoxin protein motif does not involve a conformational switch

Regulation of protein function via cracking, or local unfolding and refolding of substructures, is becoming a widely recognized mechanism of functional control. Oftentimes, cracking events are localized to secondary and tertiary structure interactions between domains that control the optimal position for catalysis and/or the formation of protein complexes. Small changes in free energy associated with ligand binding, phosphorylation, etc., can tip the balance and provide a regulatory functional switch. However, understanding the factors controlling function in single-domain proteins is still a significant challenge to structural biologists. We investigated the functional landscape of a single-domain plant-type ferredoxin protein and the effect of a distal loop on the electron-transfer center. We find the global stability and structure are minimally perturbed with mutation, whereas the functional properties are altered. Specifically, truncating the L1,2 loop does not lead to large-scale changes in the structure, determined via X-ray crystallography. Further, the overall thermal stability of the protein is only marginally perturbed by the mutation. However, even though the mutation is distal to the iron–sulfur cluster (∼20 Å), it leads to a significant change in the redox potential of the iron–sulfur cluster (57 mV). Structure-based all-atom simulations indicate correlated dynamical changes between the surface-exposed loop and the iron–sulfur cluster-binding region. Our results suggest intrinsic communication channels within the ferredoxin fold, composed of many short-range interactions, lead to the propagation of long-range signals. Accordingly, protein interface interactions that involve L1,2 could potentially signal functional changes in distal regions, similar to what is observed in other allosteric systems.

Pierced Lasso Bundles Are a New Class of Knot-like Motifs

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.

The Dominant Folding Route Minimizes Backbone Distortion in SH3

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.

The role of atomic level steric effects and attractive forces in protein folding

Protein folding into tertiary structures is controlled by an interplay of attractive contact interactions and steric effects. We investigate the balance between these contributions using structure‐based models using an all‐atom representation of the structure combined with a coarse‐grained contact potential. Tertiary contact interactions between atoms are collected into a single broad attractive well between the Cβ atoms between each residue pair in a native contact. Through the width of these contact potentials we control their tolerance for deviations from the ideal structure and the spatial range of attractive interactions. In the compact native state dominant packing constraints limit the effects of a coarse‐grained contact potential. During folding, however, the broad attractive potentials allow an early collapse that starts before the native local structure is completely adopted. As a consequence the folding transition is broadened and the free energy barrier is decreased. Eventually two‐state folding behavior is lost completely for systems with very broad attractive potentials. The stabilization of native‐like residue interactions in non‐perfect geometries early in the folding process frequently leads to structural traps. Global mirror images are a notable example. These traps are penalized by the details of the repulsive interactions only after further collapse. Successful folding to the native state requires simultaneous guidance from both attractive and repulsive interactions. Proteins 2012.

How Relevant is the Network Dynamics for the Alkali Diffusion in Alkali Silicate Glasses

Abstract The relevance of network dynamics in a disordered alkali silicate ion conductor is studied via computer simulations. It is shown that the network dynamics below the glass transition is very localized and does not give rise to major fluctuations of ionic sites. These local network fluctuations are, however, very important to understand the alkali dynamics. After immobilising the network the activation energy roughly doubles. Furthermore dynamic heterogeneities as well as backdragging forces are much stronger. Qualitatively, these observations can be rationalized by a scenario similar to the Anderson-Stuart model of ion dynamics.

Structure of the human monomeric NEET protein MiNT and its role in regulating iron and reactive oxygen species in cancer cells

Significance NEET proteins belong to a unique family of iron-sulfur (Fe-S) proteins that regulate iron and reactive oxygen homeostasis and are involved in the progression of cancer, diabetes, neurodegeneration, and aging. Of the three human NEET proteins, the structure and function of the mitochondrial matrix-localized mitochondrial inner NEET protein (MiNT) are unknown. Here, we show that MiNT is a pseudosymmetrical monomeric protein that contains two distinct Fe-S cluster-binding motifs. MiNT transfers its clusters to the human mitochondrial ferredoxins FDX1/adrenodoxin and FDX2 and is required for regulating iron and reactive oxygen species levels in the mitochondria. Our study demonstrates that although MiNT differs in structure and localization from its homodimeric NEET counterparts, it nevertheless cooperates with them in the same important pathway. The NEET family is a relatively new class of three related [2Fe-2S] proteins (CISD1–3), important in human health and disease. While there has been growing interest in the homodimeric gene products of CISD1 (mitoNEET) and CISD2 (NAF-1), the importance of the inner mitochondrial CISD3 protein has only recently been recognized in cancer. The CISD3 gene encodes for a monomeric protein that contains two [2Fe-2S] CDGSH motifs, which we term mitochondrial inner NEET protein (MiNT). It folds with a pseudosymmetrical fold that provides a hydrophobic motif on one side and a relatively hydrophilic surface on the diametrically opposed surface. Interestingly, as shown by molecular dynamics simulation, the protein displays distinct asymmetrical backbone motions, unlike its homodimeric counterparts that face the cytosolic side of the outer mitochondrial membrane/endoplasmic reticulum (ER). However, like its counterparts, our biological studies indicate that knockdown of MiNT leads to increased accumulation of mitochondrial labile iron, as well as increased mitochondrial reactive oxygen production. Taken together, our study suggests that the MiNT protein functions in the same pathway as its homodimeric counterparts (mitoNEET and NAF-1), and could be a key player in this pathway within the mitochondria. As such, it represents a target for anticancer or antidiabetic drug development.

Constructing a folding model for protein S6 guided by native fluctuations deduced from NMR structures

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 proteins 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.

A magnesium-induced triplex pre-organizes the SAM-II riboswitch

Our 13C- and 1H-chemical exchange saturation transfer (CEST) experiments previously revealed a dynamic exchange between partially closed and open conformations of the SAM-II riboswitch in the absence of ligand. Here, all-atom structure-based molecular simulations, with the electrostatic effects of Manning counter-ion condensation and explicit magnesium ions are employed to calculate the folding free energy landscape of the SAM-II riboswitch. We use this analysis to predict that magnesium ions remodel the landscape, shifting the equilibrium away from the extended, partially unfolded state towards a compact, pre-organized conformation that resembles the ligand-bound state. Our CEST and SAXS experiments, at different magnesium ion concentrations, quantitatively confirm our simulation results, demonstrating that magnesium ions induce collapse and pre-organization. Agreement between theory and experiment bolsters microscopic interpretation of our simulations, which shows that triplex formation between helix P2b and loop L1 is highly sensitive to magnesium and plays a key role in pre-organization. Pre-organization of the SAM-II riboswitch allows rapid detection of ligand with high selectivity, which is important for biological function.