An all-atom structure-based potential for proteins: bridging minimal models with all-atom empirical forcefields

Reduced model captures Mg(2+)-RNA interaction free energy of riboswitches

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.

From structure to function: the convergence of structure based models and co-evolutionary information

Understanding protein folding and function is one of the most important problems in biological research. Energy landscape theory and the folding funnel concept have provided a framework to investigate the mechanisms associated to these processes. Since protein energy landscapes are in most cases minimally frustrated, structure based models (SMBs) have successfully determined the geometrical features associated with folding and functional transitions. However, structural information is limited, particularly with respect to different functional configurations. This is a major limitation for SBMs. Alternatively, statistical methods to study amino acid co-evolution provide information on residue-residue interactions useful for the study of structure and function. Here, we show how the combination of these two methods gives rise to a novel way to investigate the mechanisms associated with folding and function. We use this methodology to explore the mechanistic aspects of protein translocation in the integral membrane protease FtsH. Dual basin-SBM simulations using the open and closed state of this hexameric motor reveals a functionally important paddling motion in the catalytic cycle. We also find that Direct Coupling Analysis (DCA) predicts physical contacts between AAA and peptidase domains of the motor, which are crucial for the open to close transition. Our combined method, which uses structural information from the open state experimental structure and co-evolutionary couplings, suggests that this methodology can be used to explore the functional landscape of complex biological macromolecules previously inaccessible to methods dependent on experimental structural information. This efficient way to sample the conformational space of large systems creates a theoretical/computational framework capable of better characterizing the functional landscape in large biomolecular assemblies.

Toward rationally redesigning bacterial two-component signaling systems using coevolutionary information

A challenge in molecular biology is to distinguish the key subset of residues that allow two-component signaling (TCS) proteins to recognize their correct signaling partner such that they can transiently bind and transfer signal, i.e., phosphoryl group. Detailed knowledge of this information would allow one to search sequence space for mutations that can be used to systematically tune the signal transmission between TCS partners as well as potentially encode a TCS protein to preferentially transfer signals to a nonpartner. Motivated by the notion that this detailed information is found in sequence data, we explore the sequence coevolution between signaling partners to better understand how mutations can positively or negatively alter their ability to transfer signal. Using direct coupling analysis for determining evolutionarily conserved protein-protein interactions, we apply a metric called the direct information score to quantify mutational changes in the interaction between TCS proteins and demonstrate that it accurately correlates with experimental mutagenesis studies probing the mutational change in measured in vitro phosphotransfer. Furthermore, by subtracting from our metric an appropriate null model corresponding to generic, conserved features in TCS signaling pairs, we can isolate the determinants that give rise to interaction specificity and recognition, which are variable among different TCS partners. Our methodology forms a potential framework for the rational design of TCS systems by allowing one to quickly search sequence space for mutations or even entirely new sequences that can increase or decrease our metric, as a proxy for increasing or decreasing phosphotransfer ability between TCS proteins.

Balancing bond, nonbond, and gō-like terms in coarse grain simulations of conformational dynamics

Characterization of the protein conformational landscape remains a challenging problem, whether it concerns elucidating folding mechanisms, predicting native structures or modeling functional transitions. Coarse-grained molecular dynamics simulation methods enable exhaustive sampling of the energetic landscape at resolutions of biological interest. The general utility of structure-based models is reviewed along with their differing levels of approximation. Simple Gō models incorporate attractive native interactions and repulsive nonnative contacts, resulting in an ideal smooth landscape. Non-Gō coarse-grained models reduce the parameter set as needed but do not include bias to any desired native structure. While non-Gō models have achieved limited success in protein coarse-graining, they can be combined with native structured-based potentials to create a balanced and powerful force field. Recent applications of such Gō-like models have yielded insight into complex folding mechanisms and conformational transitions in large macromolecules. The accuracy and usefulness of reduced representations are also revealed to be a function of the mathematical treatment of the intrinsic bonded topology.

COMPARING FOLDING MECHANISMS OF DIFFERENT PRION PROTEINS BY Gō MODEL

Prions are associated with neurodegenerative diseases induced by transmissible spongiform encephalopathies. The infectious scrapie form is referred to as PrP Sc , which has conformational change from normal prion with predominant α-helical conformation to the abnormal PrP Sc that is rich in β-sheet content. Neurodegenerative diseases have been found from both human and bovine sources, but there are no reports about infected by transmissible spongiform encephalopathies from rabbit, canine and horse sources. Here we used coarse-grained Gō model to compare the difference among human, bovine, rabbit, canine, and horse normal (cellular) prion proteins. The denatured state of normal prion has relation with the conversion from normal to abnormal prion protein, so we used all-atom Gō model to investigate the folding pathway and energy landscape for human prion protein. Through using coarse-grained Gō model, the cooperativity of the five prion proteins was characterized in terms of calorimetric criterion, sigmoidal transition, and free-energy profile. The rabbit and horse prion proteins have higher folding free-energy barrier and cooperativity, and canine prion protein has slightly higher folding free-energy barrier comparing with human and bovine prion proteins. The results from all-atom Gō model confirmed the validity of C α-Gō model. The correlations of our results with previous experimental and theoretical researches were discussed.

Differences between cotranscriptional and free riboswitch folding

Riboswitches are part of noncoding regions of messenger RNA (mRNA) that act as RNA sensors regulating gene expression of the downstream gene. Typically, one out of two distinct conformations is formed depending on ligand binding when the transcript leaves RNA polymerase (RNAP). Elongation of the RNA chain by RNAP, folding and binding all occurs simultaneously and interdependently on the seconds' timescale. To investigate the effect of transcript elongation velocity on folding for the S-adenosylmethionine (SAM)-I and adenine riboswitches we employ two complementary coarse-grained in silico techniques. Native structure-based molecular dynamics simulations provide a 3D, atomically resolved model of folding with homogenous energetics. Energetically more detailed kinetic Monte Carlo simulations give access to longer timescale by describing folding on the secondary structure level and feature the incorporation of competing aptamer conformations and a ligand-binding model. Depending on the extrusion scenarios, we observe and quantify different pathways in structure formation with robust agreements between the two techniques. In these scenarios, free-folding riboswitches exhibit different folding characteristics compared with transcription-rate limited folding. The critical transcription rate distinguishing these cases is higher than physiologically relevant rates. This result suggests that in vivo folding of the analyzed SAM-I and adenine riboswitches is transcription-rate limited.

ModBase, a database of annotated comparative protein structure models and associated resources

ModBase (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by ModPipe, an automated modeling pipeline that relies primarily on Modeller for fold assignment, sequence-structure alignment, model building and model assessment (http://salilab.org/modeller/). ModBase currently contains almost 30 million reliable models for domains in 4.7 million unique protein sequences. ModBase allows users to compute or update comparative models on demand, through an interface to the ModWeb modeling server (http://salilab.org/modweb). ModBase models are also available through the Protein Model Portal (http://www.proteinmodelportal.org/). Recently developed associated resources include the AllosMod server for modeling ligand-induced protein dynamics (http://salilab.org/allosmod), the AllosMod-FoXS server for predicting a structural ensemble that fits an SAXS profile (http://salilab.org/allosmod-foxs), the FoXSDock server for protein–protein docking filtered by an SAXS profile (http://salilab.org/foxsdock), the SAXS Merge server for automatic merging of SAXS profiles (http://salilab.org/saxsmerge) and the Pose & Rank server for scoring protein–ligand complexes (http://salilab.org/poseandrank). In this update, we also highlight two applications of ModBase: a PSI:Biology initiative to maximize the structural coverage of the human alpha-helical transmembrane proteome and a determination of structural determinants of human immunodeficiency virus-1 protease specificity.

Unfolding thermodynamics of cysteine-rich proteins and molecular thermal-adaptation of marine ciliates

Euplotes nobilii and Euplotes raikovi are phylogenetically closely allied species of marine ciliates, living in polar and temperate waters, respectively. Their evolutional relation and the sharply different temperatures of their natural environments make them ideal organisms to investigate thermal-adaptation. We perform a comparative study of the thermal unfolding of disulfide-rich protein pheromones produced by these ciliates. Recent circular dichroism (CD) measurements have shown that the two psychrophilic (E. nobilii) and mesophilic (E. raikovi) protein families are characterized by very different melting temperatures, despite their close structural homology. The enhanced thermal stability of the E. raikovi pheromones is realized notwithstanding the fact that these proteins form, as a rule, a smaller number of disulfide bonds. We perform Monte Carlo (MC) simulations in a structure-based coarse-grained (CG) model to show that the higher stability of the E. raikovi pheromones is due to the lower locality of the disulfide bonds, which yields a lower entropy increase in the unfolding process. Our study suggests that the higher stability of the mesophilic E. raikovi phermones is not mainly due to the presence of a strongly hydrophobic core, as it was proposed in the literature. In addition, we argue that the molecular adaptation of these ciliates may have occurred from cold to warm, and not from warm to cold. To provide a testable prediction, we identify a point-mutation of an E. nobilii pheromone that should lead to an unfolding temperature typical of that of E. raikovi pheromones.

Solution and crystallographic structures of the central region of the phosphoprotein from human metapneumovirus

Human metapneumovirus (HMPV) of the family Paramyxoviridae is a major cause of respiratory illness worldwide. Phosphoproteins (P) from Paramyxoviridae are essential co-factors of the viral RNA polymerase that form tetramers and possess long intrinsically disordered regions (IDRs). We located the central region of HMPV P (P(ced)) which is involved in tetramerization using disorder analysis and modeled its 3D structure ab initio using Rosetta fold-and-dock. We characterized the solution-structure of P(ced) using small angle X-ray scattering (SAXS) and carried out direct fitting to the scattering data to filter out incorrect models. Molecular dynamics simulations (MDS) and ensemble optimization were employed to select correct models and capture the dynamic character of P(ced). Our analysis revealed that oligomerization involves a compact central core located between residues 169-194 (P(core)), that is surrounded by flexible regions with α-helical propensity. We crystallized this fragment and solved its structure at 3.1 Å resolution by molecular replacement, using the folded core from our SAXS-validated ab initio model. The RMSD between modeled and experimental tetramers is as low as 0.9 Å, demonstrating the accuracy of the approach. A comparison of the structure of HMPV P to existing mononegavirales P(ced) structures suggests that P(ced) evolved under weak selective pressure. Finally, we discuss the advantages of using SAXS in combination with ab initio modeling and MDS to solve the structure of small, homo-oligomeric protein complexes.

Structural intermediates and folding events in the early assembly of the ribosomal small subunit

Using all-atom explicit solvent molecular dynamics (MD) simulations, we investigated the early structural intermediates of the 5' domain of the 16S rRNA in Escherichia coli upon the removal of the primary binding r-proteins S4, S17, and S20 and the secondary binding r-protein S16. Removal of each r-protein corresponded to the disappearance of subdomains with correlated dynamics. Correlation-based network analysis of the MD trajectories of the naked rRNA showed that the different subdomains are connected via multiple pathways with high betweenness. These pathways cross at the internal loop of helix 17 (h17) in the five-way junction (5WJ). The structure of the internal loop is disrupted by the binding of S17 and rescued by the addition of S16, suggesting an important function of the secondary binding protein in biasing the rRNA folding landscape toward the native basin. Using structure-based Gō simulations, we investigated the folding barriers of the lower four-way junction (4WJ) with h6, which is the primary binding site of S20 and the first to be transcribed. The folding of the 4WJ is consistent with the protection patterns observed in hydroxyl radical footprinting. Results from the all-atom simulations show that the fluctuations in the 5WJ are independent of the fluctuations in the 4WJ, suggesting that the subdomains fold independently and are stabilized by primary r-proteins.

Mapping polymerization and allostery of hemoglobin S using point mutations

Hemoglobin is a complex system that undergoes conformational changes in response to oxygen, allosteric effectors, mutations, and environmental changes. Here, we study allostery and polymerization of hemoglobin and its variants by application of two previously described methods: (i) AllosMod for simulating allostery dynamics given two allosterically related input structures and (ii) a machine-learning method for dynamics- and structure-based prediction of the mutation impact on allostery (Weinkam et al. J. Mol. Biol. 2013, 425, 647-661), now applicable to systems with multiple coupled binding sites, such as hemoglobin. First, we predict the relative stabilities of substates and microstates of hemoglobin, which are determined primarily by entropy within our model. Next, we predict the impact of 866 annotated mutations on hemoglobin's oxygen binding equilibrium. We then discuss a subset of 30 mutations that occur in the presence of the sickle cell mutation and whose effects on polymerization have been measured. Seven of these HbS mutations occur in three predicted druggable binding pockets that might be exploited to directly inhibit polymerization; one of these binding pockets is not apparent in the crystal structure, but only in structures generated by AllosMod. For the 30 mutations, we predict that mutation-induced conformational changes within a single tetramer tend not to significantly impact polymerization; instead, these mutations more likely impact polymerization by directly perturbing a polymerization interface. Finally, our analysis of allostery allows us to hypothesize why hemoglobin evolved to have multiple subunits and a persistent low frequency sickle cell mutation.

Structure based design of protein linkers for zinc finger nuclease

Zinc finger nucleases are a promising tool to edit DNA in many biological applications, in particular for gene knockout. Despite many efforts the number of genes that can be effectively targeted with ZFNs remains severely limited, as available constructs cannot address arbitrary gene sequences. Here, we develop a novel concept to significantly enhance the number of DNA sequences that can be targeted by ZFN. Using an efficient computational model, we provide an extensive library of possible linker molecules between individual zinc finger motifs in the construct that can skip up to 10 base pairs between adjacent zinc finger recognition sites in the DNA sequence, which increases the number of genes that can be efficiently targeted by more than an order of magnitude.

Sampling of Protein Folding Transitions: Multicanonical Versus Replica Exchange Molecular Dynamics

We compare the efficiency of multicanonical and replica exchange molecular dynamics for the sampling of folding/unfolding events in simulations of proteins with end-to-end β-sheet. In Go-model simulations of the 75-residue MNK6, we observe improvement factors of 30 in the number of folding/unfolding events of multicanonical molecular dynamics over replica exchange molecular dynamics. As an application, we use this enhanced sampling to study the folding landscape of the 36-residue DS119 with an all-atom physical force field and implicit solvent. Here, we find that the rate-limiting step is the formation of the central helix that then provides a scaffold for the parallel β-sheet formed by the two chain ends.

Analysis of the interaction interfaces of the N-terminal domain from Pseudomonas aeruginosa MutL

Mismatch Repair System corrects mutations arising from DNA replication that escape from DNA polymerase proofreading activity. This system consists of three main proteins, MutS-L-H, responsible for lesion recognition and repair. MutL is a member of GHKL ATPase family and its ATPase cycle has been proposed to modulate MutL activity during the repair process. Pseudomonas aeruginosa MutL (PaMutL) contains an N-terminal (NTD) ATPase domain connected by a linker to a C-terminal (CTD) dimerization domain that possesses metal ion-dependent endonuclease activity. With the aim to identify characteristics that allow the PaMutL NTD allosteric control of CTD endonuclease activity, we used an in silico and experimental approach to determine the interaction surfaces of P. aeruginosa NTD (PaNTD), and compared it with the well characterized Escherichia coli MutL NTD (EcNTD). Molecular dynamics simulations of PaNTD and EcNTD bound to or free of adenosine nucleotides showed that a significant difference exists between the behavior of the EcNTD and PaNTD dimerization interface, particularly in the ATP lid. Structure based simulations of MutL homologues with endonuclease activity were performed that allowed an insight of the dimerization interface behavior in this family of proteins. Our experimental results show that, unlike EcNTD, PaNTD is dimeric in presence of ADP. Simulations in mixed solvent allowed us to identify the PaNTD putative DNA binding patch and a putative interaction patch located opposite to the dimerization face. Structure based simulations of PaNTD dimer in presence of ADP or ATP suggest that nucleotide binding could differentially modulate PaNTD protein-protein interactions. Far western assays performed in presence of ADP or ATP are in agreement with our in silico analysis.

Analyzing the effect of homogeneous frustration in protein folding

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a Cα structure-based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well-separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non-native contact interactions in different folding scenarios. These findings strongly correlate with the protein free-energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.

Folding and self-assembly of the TatA translocation pore based on a charge zipper mechanism

We propose a concept for the folding and self-assembly of the pore-forming TatA complex from the Twin-arginine translocase and of other membrane proteins based on electrostatic "charge zippers." Each subunit of TatA consists of a transmembrane segment, an amphiphilic helix (APH), and a C-terminal densely charged region (DCR). The sequence of charges in the DCR is complementary to the charge pattern on the APH, suggesting that the protein can be "zipped up" by a ladder of seven salt bridges. The length of the resulting hairpin matches the lipid bilayer thickness, hence a transmembrane pore could self-assemble via intra- and intermolecular salt bridges. The steric feasibility was rationalized by molecular dynamics simulations, and experimental evidence was obtained by monitoring the monomer-oligomer equilibrium of specific charge mutants. Similar "charge zippers" are proposed for other membrane-associated proteins, e.g., the biofilm-inducing peptide TisB, the human antimicrobial peptide dermcidin, and the pestiviral E(RNS) protein.

Native-based simulations of the binding interaction between RAP74 and the disordered FCP1 peptide

By dephosphorylating the C-terminal domain (CTD) of RNA polymerase II (Pol II), the Transcription Factor IIF (TFIIF)-associating CTD phosphatase (FCP1) performs an essential function in recycling Pol II for subsequent rounds of transcription. The interaction between FCP1 and TFIIF is mediated by the disordered C-terminal tail of FCP1, which folds to form an α-helix upon binding the RAP74 subunit of TFIIF. The present work reports a structure-based simulation study of this interaction between the folded winged-helix domain of RAP74 and the disordered C-terminal tail of FCP1. The comparison of measured and simulated chemical shifts suggests that the FCP1 peptide samples 40-60% of its native helical structure in the unbound disordered ensemble. Free energy calculations suggest that productive binding begins when RAP74 makes hydrophobic contacts with the C-terminal region of the FCP1 peptide. The FCP1 peptide then folds into an amphipathic helix by zipping up the binding interface. The relative plasticity of FCP1 results in a more cooperative binding mechanism, allows for a greater diversity of pathways leading to the bound complex, and may also eliminate the need for "backtracking" from contacts that form out of sequence.

SOD1 exhibits allosteric frustration to facilitate metal binding affinity

Superoxide dismutase-1 (SOD1) is a ubiquitous, Cu and Zn binding, free-radical defense enzyme whose misfolding and aggregation play a potential key role in amyotrophic lateral sclerosis, an invariably fatal neurodegenerative disease. Over 150 mutations in SOD1 have been identified with a familial form of the disease, but it is presently not clear what unifying features, if any, these mutants share to make them pathogenic. Here, we develop several unique computational assays for probing the thermo-mechanical properties of both ALS-associated and rationally designed SOD1 variants. Allosteric interaction-free energies between residues and metals are calculated, and a series of atomic force microscopy experiments are simulated with variable tether positions to quantify mechanical rigidity "fingerprints" for SOD1 variants. Mechanical fingerprinting studies of a series of C-terminally truncated mutants, along with an analysis of equilibrium dynamic fluctuations while varying native constraints, potential energy change upon mutation, frustratometer analysis, and analysis of the coupling between local frustration and metal binding interactions for a glycine scan of 90 residues together, reveal that the apo protein is internally frustrated, that these internal stresses are partially relieved by mutation but at the expense of metal-binding affinity, and that the frustration of a residue is directly related to its role in binding metals. This evidence points to apo SOD1 as a strained intermediate with "self-allostery" for high metal-binding affinity. Thus, the prerequisites for the function of SOD1 as an antioxidant compete with apo state thermo-mechanical stability, increasing the susceptibility of the protein to misfold in the apo state.

Molecular model of a soluble guanylyl cyclase fragment determined by small-angle X-ray scattering and chemical cross-linking

Soluble guanylyl/guanylate cyclase (sGC) converts GTP to cGMP after binding nitric oxide, leading to smooth muscle relaxation and vasodilation. Impaired sGC activity is common in cardiovascular disease, and sGC stimulatory compounds are vigorously sought. sGC is a 150 kDa heterodimeric protein with two H-NOX domains (one with heme, one without), two PAS domains, a coiled-coil domain, and two cyclase domains. Binding of NO to the sGC heme leads to proximal histidine release and stimulation of catalytic activity. To begin to understand how binding leads to activation, we examined truncated sGC proteins from Manduca sexta (tobacco hornworm) that bind NO, CO, and stimulatory compound YC-1 but lack the cyclase domains. We determined the overall shape of truncated M. sexta sGC using analytical ultracentrifugation and small-angle X-ray scattering (SAXS), revealing an elongated molecule with dimensions of 115 Å × 90 Å × 75 Å. Binding of NO, CO, or YC-1 had little effect on shape. Using chemical cross-linking and tandem mass spectrometry, we identified 20 intermolecular contacts, allowing us to fit homology models of the individual domains into the SAXS-derived molecular envelope. The resulting model displays a central parallel coiled-coil platform upon which the H-NOX and PAS domains are assembled. The β1 H-NOX and α1 PAS domains are in contact and form the core signaling complex, while the α1 H-NOX domain can be removed without a significant effect on ligand binding or overall shape. Removal of 21 residues from the C-terminus yields a protein with dramatically increased proximal histidine release rates upon NO binding.

Assembling a xylanase-lichenase chimera through all-atom molecular dynamics simulations

Multifunctional enzyme engineering can improve enzyme cocktails for emerging biofuel technology. Molecular dynamics through structure-based models (SB) is an effective tool for assessing the tridimensional arrangement of chimeric enzymes as well as for inferring the functional practicability before experimental validation. This study describes the computational design of a bifunctional xylanase-lichenase chimera (XylLich) using the xynA and bglS genes from Bacillus subtilis. In silico analysis of the average solvent accessible surface area (SAS) and the root mean square fluctuation (RMSF) predicted a fully functional chimera, with minor fluctuations and variations along the polypeptide chains. Afterwards, the chimeric enzyme was built by fusing the xynA and bglS genes. XylLich was evaluated through small-angle X-ray scattering (SAXS) experiments, resulting in scattering curves with a very accurate fit to the theoretical protein model. The chimera preserved the biochemical characteristics of the parental enzymes, with the exception of a slight variation in the temperature of operation and the catalytic efficiency (kcat/Km). The absence of substantial shifts in the catalytic mode of operation was also verified. Furthermore, the production of chimeric enzymes could be more profitable than producing a single enzyme separately, based on comparing the recombinant protein production yield and the hydrolytic activity achieved for XylLich with that of the parental enzymes.

Consensus among multiple approaches as a reliability measure for flexible fitting into cryo-EM data

Cryo-electron microscopy (cryo-EM) can provide low-resolution density maps of large macromolecular assemblies. As the number of structures deposited in the Protein Data Bank by fitting a high-resolution structure into a low-resolution cryo-EM map is increasing, there is a need to revise the protocols and improve the measures for fitting. A recent study suggested using a combination of multiple automated flexible fitting approaches to improve the interpretation of cryo-EM data. The current work further explores the use of multiple approaches by validating this "consensus" fitting approach and deriving a local reliability measure. Here four different flexible fitting approaches are applied for fitting an initial structure into a simulated density map of known target structure from a dataset of proteins. It is found that the models produced from different approaches often have a consensus in conformation and are also near to the target structure, whereas cases not showing consensus are away from the target. A high correlation is also observed between the RMSF profiles calculated with respect to the average and the target structures, which indicates that the relation between consensus and accuracy can also be extended to a per-residue level. Therefore, the RMSF among the fitted models is proposed as a local reliability measure, which can be used to assess the reliability of the fit at specific regions. Hence, we encourage the community to use consensus flexible fitting with different methods to report on local reliability of the resulting models and improve the interpretation of cryo-EM data.

Backtracking due to residual structure in the unfolded state changes the folding of the third fibronectin type III domain from tenascin-C

Residual structure in the unfolded state of a protein may play a crucial role in folding and stability. In the present study, using an all (heavy)-atom structure based model and replica exchange molecular dynamics simulations, we explored the folding landscape of the third fibronectin type III domain from tenascin-C (TNfn3). Specifically, both the wild type (WT) and a variant with two additional amino acids, Gly-Leu (GL), at the C-terminus (WT(+GL)) were studied. We found that, although both domains of TNfn3 are topologically frustrated, the early formation of the native contacts from the C-terminal end of WT(+GL) causes more "backtracking" than in the WT. As a result, the WT exhibits a two-state folding behavior with a broad transition-state ensemble, whereas the WT(+GL) folds through a metastable intermediate state. Furthermore, our study confirmed that the core of both proteins is conformationally heterogeneous and noncompact, and folds late mainly due to backtracking of the part of the core. Finally, in agreement with the previous experimental studies, our results clearly demonstrated distinct thermodynamic behavior of the two proteins with WT(+GL) being more stable.

Folding and association of a homotetrameric protein complex in an all-atom Go model

The 84-residue homotetrameric BBAT1 is one of the smallest stable protein complexes and therefore is a good test system to study the self-assembly of multimeric proteins. We have researched for this protein the interplay between the folding of monomers and their assembly into tetramers. Replica exchange molecular dynamics simulations relying on a Go model are compared with earlier simulations that use the physics-based coarse-grained UNRES model.

Allosteric control in a metalloprotein dramatically alters function

Metalloproteins (MPs) comprise one-third of all known protein structures. This diverse set of proteins contain a plethora of unique inorganic moieties capable of performing chemistry that would otherwise be impossible using only the amino acids found in nature. Most of the well-studied MPs are generally viewed as being very rigid in structure, and it is widely thought that the properties of the metal centers are primarily determined by the small fraction of amino acids that make up the local environment. Here we examine both theoretically and experimentally whether distal regions can influence the metal center in the diabetes drug target mitoNEET. We demonstrate that a loop (L2) 20 Å away from the metal center exerts allosteric control over the cluster binding domain and regulates multiple properties of the metal center. Mutagenesis of L2 results in significant shifts in the redox potential of the [2Fe-2S] cluster and orders of magnitude effects on the rate of [2Fe-2S] cluster transfer to an apo-acceptor protein. These surprising effects occur in the absence of any structural changes. An examination of the native basin dynamics of the protein using all-atom simulations shows that twisting in L2 controls scissoring in the cluster binding domain and results in perturbations to one of the cluster-coordinating histidines. These allosteric effects are in agreement with previous folding simulations that predicted L2 could communicate with residues surrounding the metal center. Our findings suggest that long-range dynamical changes in the protein backbone can have a significant effect on the functional properties of MPs.

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.

Thermodynamics of strongly allosteric inhibition: a model study of HIV-1 protease

Protein inhibitors that shift the thermodynamic equilibrium towards a denatured state escape, in general, the straightforward framework of competitive or allosteric inhibitors. The equilibrium properties of peptides which compete with the folding, or more precisely destabilize the native state, of the human immunodeficiency virus (HIV)-1 protease monomer are studied within a structure-based model. The effect of peptides that disrupt the hydrophobic core of the protein can still be summarized in terms of an inhibition constant, which depends on the thermal stability of the protein. The state of the protein denatured by such a peptide is more structured than its intrinsic denatured state, but displays the same degree of compactness. Peptides that target less buried regions of the protein are less efficient and display a more complex thermodynamics that cannot be captured in a simple way.

Topography of funneled landscapes determines the thermodynamics and kinetics of protein folding

The energy landscape approach has played a fundamental role in advancing our understanding of protein folding. Here, we quantify protein folding energy landscapes by exploring the underlying density of states. We identify three quantities essential for characterizing landscape topography: the stabilizing energy gap between the native and nonnative ensembles δE, the energetic roughness ΔE, and the scale of landscape measured by the entropy S. We show that the dimensionless ratio between the gap, roughness, and entropy of the system Λ=δE/(ΔE√(2S)) accurately predicts the thermodynamics, as well as the kinetics of folding. Large Λ implies that the energy gap (or landscape slope towards the native state) is dominant, leading to more funneled landscapes. We investigate the role of topological and energetic roughness for proteins of different sizes and for proteins of the same size, but with different structural topologies. The landscape topography ratio Λ is shown to be monotonically correlated with the thermodynamic stability against trapping, as characterized by the ratio of folding temperature versus trapping temperature. Furthermore, Λ also monotonically correlates with the folding kinetic rates. These results provide the quantitative bridge between the landscape topography and experimental folding measurements.

The unique cysteine knot regulates the pleotropic hormone leptin

Leptin plays a key role in regulating energy intake/expenditure, metabolism and hypertension. It folds into a four-helix bundle that binds to the extracellular receptor to initiate signaling. Our work on leptin revealed a hidden complexity in the formation of a previously un-described, cysteine-knotted topology in leptin. We hypothesized that this unique topology could offer new mechanisms in regulating the protein activity. A combination of in silico simulation and in vitro experiments was used to probe the role of the knotted topology introduced by the disulphide-bridge on leptin folding and function. Our results surprisingly show that the free energy landscape is conserved between knotted and unknotted protein, however the additional complexity added by the knot formation is structurally important. Native state analyses led to the discovery that the disulphide-bond plays an important role in receptor binding and thus mediate biological activity by local motions on distal receptor-binding sites, far removed from the disulphide-bridge. Thus, the disulphide-bridge appears to function as a point of tension that allows dissipation of stress at a distance in leptin.

Substrate-specific reorganization of the conformational ensemble of CSK implicates novel modes of kinase function

Protein kinases use ATP as a phosphoryl donor for the posttranslational modification of signaling targets. It is generally thought that the binding of this nucleotide induces conformational changes leading to closed, more compact forms of the kinase domain that ideally orient active-site residues for efficient catalysis. The kinase domain is oftentimes flanked by additional ligand binding domains that up- or down-regulate catalytic function. C-terminal Src kinase (Csk) is a multidomain tyrosine kinase that is up-regulated by N-terminal SH2 and SH3 domains. Although the X-ray structure of Csk suggests the enzyme is compact, X-ray scattering studies indicate that the enzyme possesses both compact and open conformational forms in solution. Here, we investigated whether interactions with the ATP analog AMP-PNP and ADP can shift the conformational ensemble of Csk in solution using a combination of small angle x-ray scattering and molecular dynamics simulations. We find that binding of AMP-PNP shifts the ensemble towards more extended rather than more compact conformations. Binding of ADP further shifts the ensemble towards extended conformations, including highly extended conformations not adopted by the apo protein, nor by the AMP-PNP bound protein. These ensembles indicate that any compaction of the kinase domain induced by nucleotide binding does not extend to the overall multi-domain architecture. Instead, assembly of an ATP-bound kinase domain generates further extended forms of Csk that may have relevance for kinase scaffolding and Src regulation in the cell.

The Src protein kinases are integral members of numerous signaling pathways involved in cellular growth and differentiation. The master regulator of the Src family is the protein kinase Csk, which adds a phosphate to the C-terminal tail, inhibiting Src Kinase function. Proper regulation of these signaling pathways by Csk is essential as unregulated activity in these pathways is correlated with the development of various cancers and autoimmune diseases. Understanding the nature of the mechanism and structure of Csk may lead to therapeutics and a better understanding of Src signaling pathways. Conformational changes associated with nucleotide binding and release have been shown to regulate the efficiency of Src down-regulation by Csk. To obtain insights into the nature of these nucleotide-induced structural changes, we examined the conformation of Csk in solution while bound to the ATP analog AMP-PNP and product ADP using a combination of small angle x-ray scattering and molecular dynamics. Surprisingly, both nucleotides induce extended conformations of Csk compared to the apo-enzyme, suggesting a novel mode of function. Further understanding of this mode of function may aid in the design of cancer therapeutics that act by regulating Src signaling pathways by modulating the function of Csk.

Importance of electrostatic interactions in the association of intrinsically disordered histone chaperone Chz1 and histone H2A.Z-H2B

Histone chaperones facilitate assembly and disassembly of nucleosomes. Understanding the process of how histone chaperones associate and dissociate from the histones can help clarify their roles in chromosome metabolism. Some histone chaperones are intrinsically disordered proteins (IDPs). Recent studies of IDPs revealed that the recognition of the biomolecules is realized by the flexibility and dynamics, challenging the century-old structure-function paradigm. Here we investigate the binding between intrinsically disordered chaperone Chz1 and histone variant H2A.Z-H2B by developing a structure-based coarse-grained model, in which Debye-Hückel model is implemented for describing electrostatic interactions due to highly charged characteristic of Chz1 and H2A.Z-H2B. We find that major structural changes of Chz1 only occur after the rate-limiting electrostatic dominant transition state and Chz1 undergoes folding coupled binding through two parallel pathways. Interestingly, although the electrostatic interactions stabilize bound complex and facilitate the recognition at first stage, the rate for formation of the complex is not always accelerated due to slow escape of conformations with non-native electrostatic interactions at low salt concentrations. Our studies provide an ionic-strength-controlled binding/folding mechanism, leading to a cooperative mechanism of “local collapse or trapping” and “fly-casting” together and a new understanding of the roles of electrostatic interactions in IDPs' binding.

Histone chaperones facilitate the assembly and disassembly of nucleosome by interacting with the corresponding histones or histone variants. As the biomolecules in nucleosome are highly charged, electrostatic interactions are particularly important in these processes. The experiments have explored that the histone chaperon Chz1 as an intrinsically disordered protein (IDP) can fold by binding to its histone variants H2A.Z-H2B. Here, we developed a molecular simulation program that treated electrostatic interactions with Debye-Hückel model to study the mechanism of the association. We found that the inter-chain electrostatic interactions facilitate the coupled folding and binding transitions, consistent with the kinetic experiments and microscopic structural perspectives. Furthermore, we show that the intra-chain electrostatic interactions collapse Chz1 and slow the binding rate. The collapsed structure in IDPs caused by intra-chain electrostatic interactions has been widely found in experiments and the effect in binding is well studied in our work. Our theoretical approach shed new light on the role of electrostatics on inter-chain and intra-chain interactions for IDPs' binding and is applicable to the binding/folding of other IDPs to their targets.

Denaturant-dependent folding of GFP

We use molecular simulations using a coarse-grained model to map the folding landscape of Green Fluorescent Protein (GFP), which is extensively used as a marker in cell biology and biotechnology. Thermal and Guanidinium chloride (GdmCl) induced unfolding of a variant of GFP, without the chromophore, occurs in an apparent two-state manner. The calculated midpoint of the equilibrium folding in GdmCl, taken into account using the Molecular Transfer Model (MTM), is in excellent agreement with the experiments. The melting temperatures decrease linearly as the concentrations of GdmCl and urea are increased. The structural features of rarely populated equilibrium intermediates, visible only in free energy profiles projected along a few order parameters, are remarkably similar to those identified in a number of ensemble experiments in GFP with the chromophore. The excellent agreement between simulations and experiments show that the equilibrium intermediates are stabilized by the chromophore. Folding kinetics, upon temperature quench, show that GFP first collapses and populates an ensemble of compact structures. Despite the seeming simplicity of the equilibrium folding, flux to the native state flows through multiple channels and can be described by the kinetic partitioning mechanism. Detailed analysis of the folding trajectories show that both equilibrium and several kinetic intermediates, including misfolded structures, are sampled during folding. Interestingly, the intermediates characterized in the simulations coincide with those identified in single molecule pulling experiments. Our predictions, amenable to experimental tests, show that MTM is a practical way to simulate the effect of denaturants on the folding of large proteins.

Biomolecular dynamics: order-disorder transitions and energy landscapes

While the energy landscape theory of protein folding is now a widely accepted view for understanding how relatively weak molecular interactions lead to rapid and cooperative protein folding, such a framework must be extended to describe the large-scale functional motions observed in molecular machines. In this review, we discuss (1) the development of the energy landscape theory of biomolecular folding, (2) recent advances toward establishing a consistent understanding of folding and function and (3) emerging themes in the functional motions of enzymes, biomolecular motors and other biomolecular machines. Recent theoretical, computational and experimental lines of investigation have provided a very dynamic picture of biomolecular motion. In contrast to earlier ideas, where molecular machines were thought to function similarly to macroscopic machines, with rigid components that move along a few degrees of freedom in a deterministic fashion, biomolecular complexes are only marginally stable. Since the stabilizing contribution of each atomic interaction is on the order of the thermal fluctuations in solution, the rigid body description of molecular function must be revisited. An emerging theme is that functional motions encompass order-disorder transitions and structural flexibility provides significant contributions to the free energy. In this review, we describe the biological importance of order-disorder transitions and discuss the statistical-mechanical foundation of theoretical approaches that can characterize such transitions.

Equilibrium unfolding of the PDZ domain of β2-syntrophin

β2-syntrophin, a dystrophin-associated protein, plays a pivotal role in insulin secretion by pancreatic β-cells. It contains a PDZ domain (β2S-PDZ) that, in complex with protein-tyrosine phosphatase ICA512, anchors the dense insulin granules to actin filaments. The phosphorylation state of β2-syntrophin allosterically regulates the affinity of β2S-PDZ for ICA512, and the disruption of the complex triggers the mobilization of the insulin granule stores. Here, we investigate the thermal unfolding of β2S-PDZ at different pH and urea concentrations. Our results indicate that, unlike other PDZ domains, β2S-PDZ is marginally stable. Thermal denaturation experiments show broad transitions and cold denaturation, and a two-state model fit reveals a significant unfolded fraction under physiological conditions. Furthermore, T(m) and T(max) denaturant-dependent shifts and noncoincidence of melting curves monitored at different wavelengths suggest that two-state and three-state models fail to explain the equilibrium data properly and are in better agreement with a downhill scenario. Its higher stability at pH >9 and the results of molecular dynamics simulations indicate that this behavior of β2S-PDZ might be related to its charge distribution. All together, our results suggest a link between the conformational plasticity of the native ensemble of this PDZ domain and the regulation of insulin secretion.

Structural basis of histidine kinase autophosphorylation deduced by integrating genomics, molecular dynamics, and mutagenesis

Signal transduction proteins such as bacterial sensor histidine kinases, designed to transition between multiple conformations, are often ruled by unstable transient interactions making structural characterization of all functional states difficult. This study explored the inactive and signal-activated conformational states of the two catalytic domains of sensor histidine kinases, HisKA and HATPase. Direct coupling analyses, a global statistical inference approach, was applied to >13,000 such domains from protein databases to identify residue contacts between the two domains. These contacts guided structural assembly of the domains using MAGMA, an advanced molecular dynamics docking method. The active conformation structure generated by MAGMA simultaneously accommodated the sequence derived residue contacts and the ATP-catalytic histidine contact. The validity of this structure was confirmed biologically by mutation of contact positions in the Bacillus subtilis sensor histidine kinase KinA and by restoration of activity in an inactive KinA(HisKA):KinD(HATPase) hybrid protein. These data indicate that signals binding to sensor domains activate sensor histidine kinases by causing localized strain and unwinding at the end of the C-terminal helix of the HisKA domain. This destabilizes the contact positions of the inactive conformation of the two domains, identified by previous crystal structure analyses and by the sequence analysis described here, inducing the formation of the active conformation. This study reveals that structures of unstable transient complexes of interacting proteins and of protein domains are accessible by applying this combination of cross-validating technologies.

Assembly of the five-way junction in the ribosomal small subunit using hybrid MD-Gō simulations

Assembly of the bacterial ribosomal small subunit (SSU) begins with the folding of the five-way junction upon interaction with the primary binding protein S4. This complex contains the largest contiguous molecular signature, which is a conserved feature in all bacterial 16S rRNAs. In a previous study, we used all-atom molecular dynamics simulations to demonstrate that the co-evolving signature in the N-terminus of S4 is intrinsically disordered and capable of accelerating the binding process through a fly casting mechanism. In this paper, comparisons between the all-atom MD simulations and FRET experiments identify multiple metastable conformations of the naked five-way junction without the presence of S4. Furthermore, we capture the simultaneous folding and binding of the five-way junction and r-protein S4 by use of a structure-based Gō potential implemented within the framework of the all-atom molecular dynamics CHARMM force field. Different folding pathways are observed for the refolding of the five-way junction upon partial binding of S4. Our simulations illustrate the complex nature of RNA folding in the presence of a protein binding partner and provide insight into the role of population shift and the induced fit mechanisms in the protein:RNA folding and binding process.

Modeling Structural Flexibility of Proteins with Go-Models

Structure-based models are an efficient tool for folding studies of proteins since by construction their energy landscape is only minimal frustrated. However, their intrinsic drawback is a lack of structural flexibility as usually only one target structure is employed to construct the potentials. Hence, a Go-model may not capture differences in mutation-induced protein dynamics, if - as in the case of the disease-related A629P mutant of the Menkes protein ATP7A - the structural differences between mutant and wild type are small. In this work, we introduced three implementations of Go-models that take into account the flexibility of proteins in the NMR ensemble. Comparing the wild type and the mutant A629P of the 75-residue large 6th domain Menkes protein, we find that these new Go-potentials lead to broader distributions than Go-models relying on a single member of the NMR-ensemble. This allows us to detect the transient unfolding of a loosely formed β1β4-sheet in the mutant protein. Our results are consistent with previous simulations using physical force field and an explicit solvent, and suggests a mechanism by which these mutations cause Menkes disease. In addition, the improved Go-models suggest differences in the folding pathway between wild type and mutant, an observation that was not accessible to simulations of this 75-residue protein with a physical all-atom force field and explicit solvent.

The shadow map: a general contact definition for capturing the dynamics of biomolecular folding and function

Structure-based models (SBMs) are simplified models of the biomolecular dynamics that arise from funneled energy landscapes. We recently introduced an all-atom SBM that explicitly represents the atomic geometry of a biomolecule. While this initial study showed the robustness of the all-atom SBM Hamiltonian to changes in many of the energetic parameters, an important aspect, which has not been explored previously, is the definition of native interactions. In this study, we propose a general definition for generating atomically grained contact maps called "Shadow". The Shadow algorithm initially considers all atoms within a cutoff distance and then, controlled by a screening parameter, discards the occluded contacts. We show that this choice of contact map is not only well behaved for protein folding, since it produces consistently cooperative folding behavior in SBMs but also desirable for exploring the dynamics of macromolecular assemblies since, it distributes energy similarly between RNAs and proteins despite their disparate internal packing. All-atom structure-based models employing Shadow contact maps provide a general framework for exploring the geometrical features of biomolecules, especially the connections between folding and function.

Coarse-grain Protein Models

Coarse-graining is a powerful approach for modeling biomolecules that, over the last few decades, has been extensively applied to proteins. Coarse-grain models offer access to large systems and to slow processes without becoming computationally unmanageable. In addition, they are very versatile, enabling both the protein representation and the energy function to be adapted to the biological problem in hand. This review concentrates on modeling soluble proteins and their assemblies. It presents an overview of the coarse-grain representations, of the associated interaction potentials, and of the optimization procedures used to define them. It then shows how coarse-grain models have been used to understand processes involving proteins, from their initial folding to their functional properties, their binary interactions, and the assembly of large complexes.

Mirror images as naturally competing conformations in protein folding

Evolution has selected a protein's sequence to be consistent with the native state geometry, as this configuration must be both thermodynamically stable and kinetically accessible to prevent misfolding and loss of function. In simple protein geometries, such as coiled-coil helical bundles, symmetry produces a competing, globally different, near mirror image with identical secondary structure and similar native contact interactions. Experimental techniques such as circular dichroism, which rely on probing secondary structure content, cannot readily distinguish these folds. Here, we want to clarify whether the native fold and mirror image are energetically competitive by investigating the free energy landscape of three-helix bundles. To prevent a bias from a specific computational approach, the present study employs the structure prediction forcefield PFF01/02, explicit solvent replica exchange molecular dynamics (REMD) with the Amber94 forcefield, and structure-based simulations based on energy landscape theory. We observe that the native fold and its mirror image have a similar enthalpic stability and are thermodynamically competitive. There is evidence that the mirror fold has faster folding kinetics and could function as a kinetic trap. All together, our simulations suggest that mirror images might not just be a computational annoyance but are competing folds that might switch depending on environmental conditions or functional considerations.

Configuration-dependent diffusion dynamics of downhill and two-state protein folding

Configuration-dependent diffusion (CDD) is important for protein folding kinetics with small thermodynamic barriers. CDD can be even more crucial in downhill folding without thermodynamic barriers. We explored the CDD of a downhill protein (BBL), and a two-state protein (CI2). The hidden kinetic barriers due to CDD were revealed. The increased ~1 k(B)T kinetic barrier is in line with experimental value based on other fast folding proteins. Compared to that of CI2, the effective free-energy profile of BBL is found to be significantly influenced by CDD, and the kinetics are totally determined by diffusion. These findings are consistent with both earlier bulk and single-molecule fluorescence measurements. In addition, we found the temperature dependence of CDD. We also found that the ratio of folding transition temperature against optimal kinetic folding temperature can provide both a quantitative measure for the underlying landscape topography and an indicator for the possible appearance of downhill folding. Our study can help for a better understanding of the role of diffusion in protein folding dynamics.