Non-native Contacts Research Articles

Analysis of the dynamics of a complex, multipathway reaction: Insulin dimer dissociation.

The protein hormone insulin forms a homodimer that must dissociate to bind to its receptor. Understanding the kinetics and mechanism of dissociation is essential for rational design of therapeutic analogs. In addition to its physiological importance, this dissociation process serves as a paradigm for coupled (un)folding and (un)binding. Based on previous free energy simulations, insulin dissociation is thought to involve multiple pathways with comparable free energy barriers. Here, we analyze the mechanism of insulin dimer dissociation using a recently developed computational framework for estimating kinetic statistics from short-trajectory data. These statistics indicate that the likelihood of dissociation (the committor) closely tracks the decrease in the number of (native and nonnative) intermonomer contacts and the increase in the number of water contacts at the dimer interface; the transition state with equal likelihood of association and dissociation corresponds to an encounter complex with relatively few native contacts and many nonnative contacts. We identify four pathways out of the dimer state and quantify their contributions to the rate, as well as their exchange, by computing reactive fluxes. We show that both the pathways and their extents of exchange can be understood in terms of rotations around three axes of the dimer structure. Our results provide insights into the kinetics of insulin analogues and, more generally, how to characterize complex, multipathway processes.

On the role of native contact cooperativity in protein folding

The consistency of energy landscape theory predictions with available experimental data, as well as direct evidence from molecular simulations, have shown that protein folding mechanisms are largely determined by the contacts present in the native structure. As expected, native contacts are generally energetically favorable. However, there are usually at least as many energetically favorable nonnative pairs owing to the greater number of possible nonnative interactions. This apparent frustration must therefore be reduced by the greater cooperativity of native interactions. In this work, we analyze the statistics of contacts in the unbiased all-atom folding trajectories obtained by Shaw and coworkers, focusing on the unfolded state. By computing mutual cooperativities between contacts formed in the unfolded state, we show that native contacts form the most cooperative pairs, while cooperativities among nonnative or between native and nonnative contacts are typically much less favorable or even anticooperative. Furthermore, we show that the largest network of cooperative interactions observed in the unfolded state consists mainly of native contacts, suggesting that this set of mutually reinforcing interactions has evolved to stabilize the native state.

Impact of Hydrodynamic Interactions on the Kinetic Pathway of Protein Folding.

Protein folding is a fundamental process critical to cellular function and human health, but it remains a grand challenge in biophysics. Hydrodynamic interaction (HI) plays a vital role in the self-organization of soft and biological materials, yet its role in protein folding is not fully understood despite folding occurring in a fluid environment. Here, we use the fluid particle dynamics method to investigate many-body hydrodynamic couplings between amino acid residues and fluid motion in the folding kinetics of a coarse-grained four-α-helices bundle protein. Our results reveal that HI helps select fast folding pathways to the native state without being kinetically trapped, significantly speeding up the folding kinetics compared to its absence. First, the directional flow along the protein backbone expedites protein collapse. Then, the incompressibility-induced squeezing flow effects retard the accumulation of non-native hydrophobic contacts, thus preventing the protein from being trapped in local energy minima during the conformational search of the native structure. We also find that the significance of HI in folding kinetics depends on temperature, with a pronounced effect under biologically relevant conditions. Our findings suggest that HI, particularly the short-range squeezing effect, may be crucial in avoiding protein misfolding.

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water.

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

Nanomechanical footprint of SARS-CoV-2 variants in complex with a potent nanobody by molecular simulations.

Rational design of novel antibody therapeutics against viral infections such as coronavirus relies on surface complementarity and high affinity for their effectiveness. Here, we explore an additional property of protein complexes, the intrinsic mechanical stability, in SARS-CoV-2 variants when complexed with a potent antibody. In this study, we utilized a recent implementation of the GōMartini 3 approach to investigate large conformational changes in protein complexes with a focus on the mechanostability of the receptor-binding domain (RBD) from WT, Alpha, Delta, and XBB.1.5 variants in complex with the H11-H4 nanobody. The analysis revealed moderate differences in mechanical stability among these variants. Also, we identified crucial residues in both the RBD and certain protein segments in the nanobody that contribute to this property. By performing pulling simulations and monitoring the presence of specific native and non-native contacts across the protein complex interface, we provided mechanistic insights into the dissociation process. Force-displacement profiles indicate a tensile force clamp mechanism associated with the type of protein complex. Our computational approach not only highlights the key mechanostable interactions that are necessary to maintain overall stability, but it also paves the way for the rational design of potent antibodies that are mechanostable and effective against emergent SARS-CoV-2 variants.

Benchmarking of force fields to characterize the intrinsically disordered R2-FUS-LC region

Intrinsically Disordered Proteins (IDPs) play crucial roles in numerous diseases like Alzheimer's and ALS by forming irreversible amyloid fibrils. The effectiveness of force fields (FFs) developed for globular proteins and their modified versions for IDPs varies depending on the specific protein. This study assesses 13 FFs, including AMBER and CHARMM, by simulating the R2 region of the FUS-LC domain (R2-FUS-LC region), an IDP implicated in ALS. Due to the flexibility of the region, we show that utilizing multiple measures, which evaluate the local and global conformations, and combining them together into a final score are important for a comprehensive evaluation of force fields. The results suggest c36m2021s3p with mTIP3p water model is the most balanced FF, capable of generating various conformations compatible with known ones. In addition, the mTIP3P water model is computationally more efficient than those of top-ranked AMBER FFs with four-site water models. The evaluation also reveals that AMBER FFs tend to generate more compact conformations compared to CHARMM FFs but also more non-native contacts. The top-ranking AMBER and CHARMM FFs can reproduce intra-peptide contacts but underperform for inter-peptide contacts, indicating there is room for improvement.

Structural Dynamics of Lys11-Selective Deubiquitinylase Cezanne-1 during the Catalytic Cycle.

Deubiquitinylating enzymes (DUBs) regulate the deubiquitinylation process of post-translationally modified proteins and thus control protein signaling in various cellular processes. The DUB Cezanne-1 catalyzes the cleavage of the iso-peptide bond of Lys11-linked polyubiquitin chains with high selectivity. Crystal structures of Cezanne-1 in different states provide important insight regarding the complex formation and global changes during the catalytic cycle but are lacking details of dynamics and control of activation. Activity-based probes are used to isolate intermediate states upon forming covalent bonds with the DUB active site. Those, however, may lead to structures that are non-native. Conformational changes of Cezanne-1, during its process of activation and proteolytic activity, are investigated using all-atom molecular dynamics (MD) simulations of the ubiquitin-free, diubiquitin-bound, and monoubiquitin-bound Cezanne-1 DUB for a total of ∼18 μs. Our results show that ubiquitin-free Cezanne-1 dynamically shuttles between catalytically competent and incompetent states which suggests that its activation is independent of substrate binding. The catalytically competent substrate-free Cezanne-1 promotes distal ubiquitin substrate access to the catalytic center. The subsequent binding of the proximal ubiquitin shifts the equilibrium toward the catalytically competent state of the dyad, thereby promoting proteolysis of the iso-peptide bond. After cleavage of the scissile bond, sequential dissociation of first the proximal ubiquitin induces the inactivation of Cezanne-1. The subsequent release of the distal ubiquitin fully reconstitutes the inactive substrate-free state of Cezanne-1. The process of activation and catalytic turnover of DUB Cezanne-1 is a multistage cycle with several critical dynamic transitions that cannot be characterized based on protein structures alone. Activity-based probes of cysteine proteases lead to non-native protein-protein contacts, which need to be resolved in order to be able to issue statements about physiological states and substrate binding.

On the role of cooperativity in funneled energy landscapes.

The consistency of energy landscape theory predictions with available experimental data, as well as direct evidence from molecular simulations, have shown that protein folding mechanisms are largely determined by the contacts present in the native structure. As expected, native contacts are generally energetically favorable, however there are usually at least as many energetically favorable non-native pairs owing to the greater number of possible non-native interactions. This apparent frustration must therefore be reduced by the greater cooperativity of native interactions. In this work, we analyze the unbiased all-atom folding trajectories obtained by Shaw and co-workers, focusing on the unfolded state. By computing mutual cooperativities between contacts formed in the unfolded state, we are able to show that native contacts are largely more cooperative with each other than non-native contacts. Furthermore, we show that our analysis can determine the cooperative structures in the unfolded state and that native contacts can be more accurately predicted from the unfolded state trajectories than from residue contact frequencies in the unfolded state or from statistical pair potentials.

Beneficial and detrimental effects of non-specific binding during DNA hybridization

DNA strands have to sample numerous states to find the alignment that maximizes Watson-Crick-Franklin base pairing. This process depends strongly on sequence, which affects the stability of the native duplex as well as the prevalence of non-native inter- and intramolecular helices. We present a theory that describes DNA hybridization as a three-stage process: diffusion, registry search, and zipping. We find that non-specific binding affects each of these stages in different ways. Mis-registered intermolecular binding in the registry search stage helps DNA strands sample different alignments and accelerates the hybridization rate. Non-native intramolecular structure affects all three stages by rendering portions of the molecule inert to intermolecular association, limiting mis-registered alignments to be sampled, and impeding the zipping process. Once in-register base pairs are formed, the stability of the native structure is important to hold the molecules together long enough for non-native contacts to break.

Quinoxalinones as A Novel Inhibitor Scaffold for EGFR (L858R/T790M/C797S) Tyrosine Kinase: Molecular Docking, Biological Evaluations, and Computational Insights.

Combating acquired drug resistance of EGFR tyrosine kinase (TK) is a great challenge and an urgent necessity in the management of non-small cell lung cancers. The advanced EGFR (L858R/T790M/C797S) triple mutation has been recently reported, and there have been no specific drugs approved for this strain. Therefore, our research aimed to search for effective agents that could impede the function of EGFR (L858R/T790M/C797S) TK by the integration of in silico and in vitro approaches. Our in-house quinoxalinone-containing compounds were screened through molecular docking and their biological activity was then verified by enzyme- and cell-based assay. We found that the four quinoxalinone-containing compounds including CPD4, CPD15, CPD16, and CPD21 were promising to be novel EGFR (L858R/T790M/C797S) TK inhibitors. The IC50 values measured by the enzyme-based assay were 3.04 ± 1.24 nM; 6.50 ± 3.02 nM,10.50 ± 1.10 nM; and 3.81 ± 1.80 nM, respectively, which are at a similar level to a reference drug; osimertinib (8.93 ± 3.01 nM). Besides that, they displayed cytotoxic effects on a lung cancer cell line (H1975) with IC50 values in the range of 3.47 to 79.43 μM. In this proposed study, we found that all screened compounds could interact with M793 at the hinge regions and two mutated residues including M790 and S797; which may be the main reason supporting the inhibitory activity in vitro. The structural dynamics revealed that the screened compounds have sufficient non-native contacts with surrounding amino acids and could be well-buried in the binding site's cleft. In addition, all predicted physicochemical parameters were favorable to be drug-like based on Lipinski's rule of five, and no extreme violation of toxicity features was found. Altogether, this study proposes a novel EGFR (L858R/T790M/C797S) TK inhibitor scaffold and provides a detailed understanding of compounds' recognition and susceptibility at the molecular level.

Effects of Surface Tethering on the Thermodynamics and Kinetics of Frustrated Protein Folding.

The interaction between the protein and surface plays an important role in biology and biotechnology. To understand how surface tethering influences the folding behavior of frustrated proteins, in this work, we systematically study the thermodynamics and folding kinetics of the bacterial immunity protein Im7 and Fyn SH3 domain tethered to a surface using Langevin dynamics simulations. Upon surface tethering, the stabilization often results from the entropic effect, whereas the destabilization is usually caused by either an energetic or entropic effect. For the Fyn SH3 domain with a two-state folding manner, the influence of nonnative interactions on thermodynamic stability is not significant, while nonnative interactions can weaken the effect of surface tethering on the change in the folding rate. By contrast, for the frustrated protein Im7, depending on where the protein is tethered, the surface tethering can promote or suppress misfolding by modulating specific nonnative contacts, thereby altering the folding rate and folding mechanism. Because surface tethering can change the intrachain diffusivity of unfolding, the kinetic stability cannot be well captured by the thermodynamic stability at some tether points. This study should be helpful in general to understand how surface tethering affects the folding energy landscape of frustrated proteins.

Coevolution-derived native and non-native contacts determine the emergence of a novel fold in a universally conserved family of transcription factors.

The NusG protein family is structurally and functionally conserved in all domains of life. Its members directly bind RNA polymerases and regulate transcription processivity and termination. RfaH, a divergent sub-family in its evolutionary history, is known for displaying distinct features than those in NusG proteins, which allows them to regulate the expression of virulence factors in enterobacteria in a DNA sequence-dependent manner. A striking feature is its structural interconversion between an active fold, which is the canonical NusG three-dimensional structure, and an autoinhibited fold, which is distinctively novel. How this novel fold is encoded within RfaH sequence to encode a metamorphic protein remains elusive. In this work, we used publicly available genomic RfaH protein sequences to construct a complete multiple sequence alignment, which was further augmented with metagenomic sequences and curated by predicting their secondary structure propensities using JPred. Coevolving pairs of residues were calculated from these sequences using plmDCA and GREMLIN, which allowed us to detect the enrichment of key metamorphic contacts after sequence filtering. Finally, we combined our coevolutionary predictions with molecular dynamics to demonstrate that these interactions are sufficient to predict the structures of both native folds, where coevolutionary-derived non-native contacts may play a key role in achieving the compact RfaH novel fold. All in all, emergent coevolutionary signals found within RfaH sequences encode the autoinhibited and active folds of this protein, shedding light on the key interactions responsible for the action of this metamorphic protein.

Role of cotranslational folding for β-sheet-enriched proteins: A perspective from molecular dynamics simulations.

The formations of correct three-dimensional structures of proteins are essential to their functions. Cotranslational folding is vital for proteins to form correct structures in vivo. Although some experiments have shown that cotranslational folding can improve the efficiency of folding, its microscopic mechanism is not yet clear. Previously, we built a model of the ribosomal exit tunnel and investigated the cotranslational folding of a three-helix protein by using all-atom molecular dynamics simulations. Here we study the cotranslational folding of three β-sheet-enriched proteins using the same method. The results show that cotranslational folding can enhance the helical population in most cases and reduce non-native long-range contacts before emerging from the ribosomal exit tunnel. After exiting the tunnel, all proteins fall into local minimal states and the structural ensembles of cotranslational folding show more helical conformations than those of free folding. In particular, for one of the three proteins, the GTT WW domain, we find that one local minimum state of the cotranslational folding is the known folding intermediate, which is not found in free folding. This result suggests that the cotranslational folding may increase the folding efficiency by accelerating the sampling more than by avoiding the misfolded state, which is presently a mainstream viewpoint.

In silico method for identifying the key residues in a knotted protein: with MJ0366 as an example.

The knotted proteins are a class of rare but biologically important proteins, due to the special topology of their native structure. Here we present a simple in silico method to identify the key residues for knotting and unknotting in a knotted protein, using the trefoil protein MJ0366 as an example. We first simulate the folding process via the annealing molecular dynamics (AMD) simulations in the coarse-grained "Go"-like model. From the folding trajectories, we monitor the knotting process using the quantity "length of knot tails". In the meantime, we analyze the evolution of the local geometry of the Cα trace with the help of the Discrete Frenet Frame (DFF). We identify the key residues by correlating the local geometry at each residue with the variable "length of knot tails" in the folding process, where a higher correlation coefficient indicates that the residue is more important for knotting. We validate our method by comparing with the experimental results in the literature. With the same method, we further predict the key residues for unknotting MJ0366 using the AMD simulations in both the coarse-grained "Go"-like model and all-atom (AA) force field model, respectively. We find that the key residues for unknotting are partially overlapped with those for knotting, indicating that the pathways for unknotting and knotting are generally similar except for the existence of some non-native contact interactions in the unknotting process. This in silico method can provide a new insight for understanding the knotting and unknotting processes of a knotted protein.

Small Neutral Crowding Solute Effects on Protein Folding Thermodynamic Stability and Kinetics.

Molecular crowding is a ubiquitous phenomenon in biological systems, with significant consequences on protein folding and stability. Small compounds, such as the osmolyte trimethylamine N-oxide (TMAO), can also present similar effects. To analyze the effects arising from these solute-like molecules, we performed a series of crowded coarse-grained folding simulations. Two well-known proteins were chosen, CI2 and SH3, modeled by the alpha-carbon-structure-based model. In the simulations, the crowding molecules were represented by low-sized neutral atom beads in different concentrations. The results show that a low level of the volume fraction occupied by neutral agents can change protein stability and folding kinetics for the two systems. However, the kinetics were shown to be unaffected in their respective folding temperatures. The faster kinetics correlates with changes in the folding route for systems at the same temperature, where non-native contacts were shown to be relevant. The transition states of the two systems with and without crowders are similar. In the case of SH3, there are differences in the structuring of two strands, which may be associated with the increase in its folding rate, in addition to the destabilization of the denatured ensemble. The present study also detected a crossover in the thermodynamic stability behavior, previously observed experimentally and theoretically. As the temperature increases, crowders change from destabilizing to stabilizing agents.

Frustrated peptide chains at the fibril tip control the kinetics of growth of amyloid-β fibrils

Amyloid fibrillization is an exceedingly complex process in which incoming peptide chains bind to the fibril while concertedly folding. The coupling between folding and binding is not fully understood. We explore the molecular pathways of association of Aβ40 monomers to fibril tips by combining time-resolved in situ scanning probe microscopy with molecular modeling. The comparison between experimental and simulation results shows that a complex supported by nonnative contacts is present in the equilibrium structure of the fibril tip and impedes fibril growth in a supersaturated solution. The unraveling of this frustrated state determines the rate of fibril growth. The kinetics of growth of freshly cut fibrils, in which the bulk fibril structure persists at the tip, complemented by molecular simulations, indicate that this frustrated complex comprises three or four monomers in nonnative conformations and likely is contained on the top of a single stack of peptide chains in the fibril structure. This pathway of fibril growth strongly deviates from the common view that the conformational transformation of each captured peptide chain is templated by the previously arrived peptide. The insights into the ensemble structure of the frustrated complex may guide the search for suppressors of Aβ fibrillization. The uncovered dynamics of coupled structuring and assembly during fibril growth are more complex than during the folding of most globular proteins, as they involve the collective motions of several peptide chains that are not guided by a funneled energy landscape.

Selection for cooperativity causes epistasis predominately between native contacts and enables epistasis-based structure reconstruction

Epistasis and cooperativity of folding both result from networks of energetic interactions in proteins. Epistasis results from energetic interactions among mutants, whereas cooperativity results from energetic interactions during folding that reduce the presence of intermediate states. The two concepts seem intuitively related, but it is unknown how they are related, particularly in terms of selection. To investigate their relationship, we simulated protein evolution under selection for cooperativity and separately under selection for epistasis. Strong selection for cooperativity created strong epistasis between contacts in the native structure but weakened epistasis between nonnative contacts. In contrast, selection for epistasis increased epistasis in both native and nonnative contacts and reduced cooperativity. Because epistasis can be used to predict protein structure only if it preferentially occurs in native contacts, this result indicates that selection for cooperativity may be key for predicting structure using epistasis. To evaluate this inference, we simulated the evolution of guanine nucleotide-binding protein (GB1) with and without cooperativity. With cooperativity, strong epistatic interactions clearly map out the native GB1 structure, while allowing the presence of intermediate states (low cooperativity) obscured the structure. This indicates that using epistasis measurements to reconstruct protein structure may be inappropriate for proteins with stable intermediates.

Converging experimental and computational views of the knotting mechanism of a small knotted protein

MJ0366 from Methanocaldococcus jannaschii is the smallest topologically knotted protein known to date. 92 residues in length, MJ0366 ties a trefoil (31) knot by threading its C-terminal helix through a buttonhole formed by the remainder of the secondary structure elements. By generating a library of point mutations at positions pertinent to the knot formation, we systematically evaluated the contributions of individual residues to the folding stability and kinetics of MJ0366. The experimental Φ-values were used as restraints to computationally generate an ensemble of conformations that correspond to the transition state of MJ0366, which revealed several nonnative contacts. The importance of these nonnative contacts in stabilizing the transition state of MJ0366 was confirmed by a second round of mutagenesis, which also established the pivotal role of F15 in stapling the network of hydrophobic interactions around the threading C-terminal helix. Our converging experimental and computational results show that, despite the small size, the transition state of MJ0366 is formed at a very late stage of the folding reaction coordinate, following a polarized pathway. Eventually, the formation of extensive native contacts, as well as a number of nonnative ones, leads to the threading of the C-terminal helix that defines the topological knot.

Energy landscape remodeling mechanism of Hsp70-chaperone-accelerated protein folding

Hsp70 chaperone is one of the key protein machines responsible for the quality control of protein production in cells. Facilitating in vivo protein folding by counteracting misfolding and aggregation is the essence of its biological function. Although the allosteric cycle during its functional actions has been well characterized both experimentally and computationally, the mechanism by which Hsp70 assists protein folding is still not fully understood. In this work, we studied the Hsp70-mediated folding of model proteins with rugged energy landscape by using molecular simulations. Different from the canonical scenario of Hsp70 functioning, which assumes that folding of substrate proteins occurs spontaneously after releasing from chaperones, our results showed that the substrate protein remains in contacts with the chaperone during its folding process. The direct chaperone-substrate interactions in the open conformation of Hsp70 tend to shield the substrate sites prone to form non-native contacts, which therefore avoids the frustrated folding pathway, leading to a higher folding rate and less probability of misfolding. Our results suggest that in addition to the unfoldase and holdase functions widely addressed in previous studies, Hsp70 can facilitate the folding of its substrate proteins by remodeling the folding energy landscape and directing the folding processes, demonstrating the foldase scenario. These findings add new, to our knowledge, insights into the general molecular mechanisms of chaperone-mediated protein folding.