Folding of Top7 in unbiased all‐atom Monte Carlo simulations

For computational studies of protein folding, proteins with both helical and β‐sheet secondary structure elements are very challenging, as they expose subtle biases of the physical models. Here, we present reproducible folding of a 92 residue α/β protein (residues 3–94 of Top7, PDB ID: 1QYS) in computer simulations starting from random initial conformations using a transferable physical model which has been previously shown to describe the folding and thermodynamic properties of about 20 other smaller proteins of different folds. Top7 is a de novo designed protein with two α‐helices and a five stranded β‐sheet. Experimentally, it is known to be unusually stable for its size, and its folding transition distinctly deviates from the two‐state behavior commonly seen in natural single domain proteins. In our all‐atom implicit solvent parallel tempering Monte Carlo simulations, Top7 shows a rapid transition to a group of states with high native‐like secondary structure, and a much slower subsequent transition to the native state with a root mean square deviation of about 3.5 Å from the experimentally determined structure. Consistent with experiments, we find Top7 to be thermally extremely stable, although the simulations also find a large number of very stable non‐native states with high native‐like secondary structure. Proteins 2013; 81:1446–1456. © 2013 Wiley Periodicals, Inc.     

Folding simulations of Trp‐cage mini protein in explicit solvent using biasing potential replica‐exchange molecular dynamics simulations

Replica exchange molecular dynamics (RexMD) simulations are frequently used for studying structure formation and dynamics of peptides and proteins. A significant drawback of standard temperature RexMD is, however, the rapid increase of the replica number with increasing system size to cover a desired temperature range. A recently developed Hamiltonian RexMD method has been used to study folding of the Trp‐cage protein. It employs a biasing potential that lowers the backbone dihedral barriers and promotes peptide backbone transitions along the replica coordinate. In two independent applications of the biasing potential RexMD method including explicit solvent and starting from a completely unfolded structure the formation of near‐native conformations was observed after 30–40 ns simulation time. The conformation representing the most populated cluster at the final simulation stage had a backbone root mean square deviation of ～1.3 Å from the experimental structure. This was achieved with a very modest number of five replicas making it well suited for peptide and protein folding and refinement studies including explicit solvent. In contrast, during five independent continuous 70 ns molecular dynamics simulations formation of collapsed states but no near native structure formation was observed. The simulations predict a largely collapsed state with a significant helical propensity for the helical domain of the Trp‐cage protein already in the unfolded state. Hydrogen bonded bridging water molecules were identified that could play an active role by stabilizing the arrangement of the helical domain with respect to the rest of the chain already in intermediate states of the protein. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Nonuniform chain collapse during early stages of staphylococcal nuclease folding detected by fluorescence resonance energy transfer and ultrarapid mixing methods

The development of tertiary structure during folding of staphylococcal nuclease (SNase) was studied by time‐resolved fluorescence resonance energy transfer measured using continuous‐ and stopped‐flow techniques. Variants of this two‐domain protein containing intradomain and interdomain fluorescence donor/acceptor pairs (Trp and Cys‐linked fluorophore or quencher) were prepared to probe the intradomain and interdomain structural evolution accompanying SNase folding. The intra‐domain donor/acceptor pairs are within the β‐barrel domain (Trp27/Cys64 and Trp27/Cys97) and the interdomain pair is between the α‐helical domain and the β‐barrel domain (Trp140/Cys64). Time‐resolved energy transfer efficiency accompanying folding and unfolding at different urea concentrations was measured over a time range from 30 μs to ～10 s. Information on average donor/acceptor distances at different stages of the folding process was obtained by using a quantitative kinetic modeling approach. The average distance for the donor/acceptor pairs in the β‐barrel domain decreases to nearly native values whereas that of the interdomain donor/acceptor pairs remains unchanged in the earliest intermediate (<500 μs of refolding). This indicates a rapid nonuniform collapse resulting in an ensemble of heterogeneous conformations in which the central region of the β‐barrel domain is well developed while the C‐terminal α‐helical domain remains disordered. The distance between Trp140 and Cys64 decreases to native values on the 100‐ms time scale, indicating that the α‐helical domain docks onto the preformed β‐barrel at a late stage of the folding. In addition, the unfolded state is found to be more compact under native conditions, suggesting that changes in solvent conditions may induce a nonspecific hydrophobic collapse.     

Constraining local structure can speed up folding by promoting structural polarization of the folding pathway

The pathway which proteins take to fold can be influenced from the earliest events of structure formation. In this light, it was both predicted and confirmed that increasing the stiffness of a beta hairpin turn decreased the size of the transition state ensemble (TSE), while increasing the folding rate. Thus, there appears to be a relationship between conformationally restricting the TSE and increasing the folding rate, at least for beta hairpin turns. In this study, we hypothesize that the enormous sampling necessary to fold even two-state folding proteins in silico could be reduced if local structure constraints were used to restrict structural heterogeneity by polarizing folding pathways or forcing folding into preferred routes. Using a Gō model, we fold Chymotrypsin Inhibitor 2 (CI-2) and the src SH3 domain after constraining local sequence windows to their native structure by rigid body dynamics (RBD). Trajectories were monitored for any changes to the folding pathway and differences in the kinetics compared with unconstrained simulations. Constraining local structure decreases folding time two-fold for 41% of src SH3 windows and 45% of CI-2 windows. For both proteins, folding times are never significantly increased after constraining any window. Structural polarization of the folding pathway appears to explain these rate increases. Folding rate enhancements are consistent with the goal to reduce sampling time necessary to reach native structures during folding simulations. As anticipated, not all constrained windows showed an equal decrease in folding time. We conclude by analyzing these differences and explain why RBD may be the preferred way to constrain structure.     

Information and redundancy in the burial folding code of globular proteins within a wide range of shapes and sizes

Recent ab initio folding simulations for a limited number of small proteins have corroborated a previous suggestion that atomic burial information obtainable from sequence could be sufficient for tertiary structure determination when combined to sequence‐independent geometrical constraints. Here, we use simulations parameterized by native burials to investigate the required amount of information in a diverse set of globular proteins comprising different structural classes and a wide size range. Burial information is provided by a potential term pushing each atom towards one among a small number L of equiprobable concentric layers. An upper bound for the required information is provided by the minimal number of layers L^min still compatible with correct folding behavior. We obtain L^min between 3 and 5 for seven small to medium proteins with 50 ≤ N_r ≤ 110 residues while for a larger protein with N_r = 141 we find that L ≥ 6 is required to maintain native stability. We additionally estimate the usable redundancy for a given L ≥ L^min from the burial entropy associated to the largest folding‐compatible fraction of “superfluous” atoms, for which the burial term can be turned off or target layers can be chosen randomly. The estimated redundancy for small proteins with L = 4 is close to 0.8. Our results are consistent with the above‐average quality of burial predictions used in previous simulations and indicate that the fraction of approachable proteins could increase significantly with even a mild, plausible, improvement on sequence‐dependent burial prediction or on sequence‐independent constraints that augment the detectable redundancy during simulations. Proteins 2016; 84:515–531. © 2016 Wiley Periodicals, Inc.     

Partial unfolding and refolding for structure refinement: A unified approach of geometric simulations and molecular dynamics

The most successful protein structure prediction methods to date have been template‐based modeling (TBM) or homology modeling, which predicts protein structure based on experimental structures. These high accuracy predictions sometimes retain structural errors due to incorrect templates or a lack of accurate templates in the case of low sequence similarity, making these structures inadequate in drug‐design studies or molecular dynamics simulations. We have developed a new physics based approach to the protein refinement problem by mimicking the mechanism of chaperons that rehabilitate misfolded proteins. The template structure is unfolded by selectively (targeted) pulling on different portions of the protein using the geometric based technique FRODA, and then refolded using hierarchically restrained replica exchange molecular dynamics simulations (hr‐REMD). FRODA unfolding is used to create a diverse set of topologies for surveying near native‐like structures from a template and to provide a set of persistent contacts to be employed during re‐folding. We have tested our approach on 13 previous CASP targets and observed that this method of folding an ensemble of partially unfolded structures, through the hierarchical addition of contact restraints (that is, first local and then nonlocal interactions), leads to a refolding of the structure along with refinement in most cases (12/13). Although this approach yields refined models through advancement in sampling, the task of blind selection of the best refined models still needs to be solved. Overall, the method can be useful for improved sampling for low resolution models where certain of the portions of the structure are incorrectly modeled. Proteins 2015; 83:2279–2292. © 2015 Wiley Periodicals, Inc.     

A hierarchical order within protein structures underlies large separations between strands in β‐sheets

Protein β‐sheets often involve nonlocal interactions between parts of the polypeptide chain that are separated by hundreds of residues, raising the question of how these nonlocal contacts form. A recent study of the smallest β‐sheets found that their formation was not driven by signals hidden in the primary sequence. Instead, the strands in these sheets were either local in sequence, or, when separated by large sequential distances, the intervening residues were found to fold into compact modules that anchored distant parts of the chain in close spatial proximity. Here, we examine larger β‐sheets to investigate the extensibility of this principle. From an analysis of the β‐sheets in a nonredundant protein dataset, we find that a highly ordered hierarchical relationship exists in the intervening structure between nonlocal β‐strands. This observation is almost universal: virtually all β‐sheets, no matter their complexity, appear to adopt an antiparallel model to manage the nonlocal aspects of their assembly, one where the chain, having left the vicinity of an unfinished β‐sheet, retraces its steps via the same route to complete the initial sheet. Exceptions typically involve unstructured regions at chain termini. Moreover, an analysis of the residues involved in nonlocal crossstrand interactions did not produce any evidence of a signal hidden in the sequence that might direct long‐range interactions. These results build on those reported for the smallest sheets, suggesting that sheet formation is either local in sequence or local in space following prior folding events that anchor disparate parts of the chain in close proximity. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Topological and sequence information predict that foldons organize a partially overlapped and hierarchical structure

It has been suggested that proteins have substructures, called foldons, which can cooperatively fold into the native structure. However, several prior investigations define foldons in various ways, citing different foldon characteristics, thereby making the concept of a foldon ambiguous. In this study, we perform a Gō model simulation and analyze the characteristics of substructures that cooperatively fold into the native‐like structure. Although some results do not agree well with the experimental evidence due to the simplicity of our coarse‐grained model, our results strongly suggest that cooperatively folding units sometimes organize a partially overlapped and hierarchical structure. This view makes us easy to interpret some different proposal about the foldon as a difference of the hierarchical structure. On the basis of this finding, we present a new method to assign foldons and their hierarchy, using structural and sequence information. The results show that the foldons assigned by our method correspond to the intermediate structures identified by some experimental techniques. The new method makes it easy to predict whether a protein folds sequentially into the native structure or whether some foldons fold into the native structure in parallel. Proteins 2015; 83:1900–1913. © 2015 Wiley Periodicals, Inc.     

Free‐energy function based on an all‐atom model for proteins

We have developed a free‐energy function based on an all‐atom model for proteins. It comprises two components, the hydration entropy (HE) and the total dehydration penalty (TDP). Upon a transition to a more compact structure, the number of accessible configurations arising from the translational displacement of water molecules in the system increases, leading to a water‐entropy gain. To fully account for this effect, the HE is calculated using a statistical‐mechanical theory applied to a molecular model for water. The TDP corresponds to the sum of the hydration energy and the protein intramolecular energy when a fully extended structure, which possesses the maximum number of hydrogen bonds with water molecules and no intramolecular hydrogen bonds, is chosen as the standard one. When a donor and an acceptor (e.g., N and O, respectively) are buried in the interior after the break of hydrogen bonds with water molecules, if they form an intramolecular hydrogen bond, no penalty is imposed. When a donor or an acceptor is buried with no intramolecular hydrogen bond formed, an energetic penalty is imposed. We examine all the donors and acceptors for backbone‐backbone, backbone‐side chain, and side chain‐side chain intramolecular hydrogen bonds and calculate the TDP. Our free‐energy function has been tested for three different decoy sets. It is better than any other physics‐based or knowledge‐based potential function in terms of the accuracy in discriminating the native fold from misfolded decoys and the achievement of high Z‐scores. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Conformational dynamics is more important than helical propensity for the folding of the all α‐helical protein Im7

Im7 folds via an on‐pathway intermediate that contains three of the four native α‐helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94‐residue variant of the 87‐residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model‐free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea‐denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α‐helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α‐helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea‐denatured states promote formation of a three‐helix intermediate in which the residues that form helix III remain non‐helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.     

Roles of non‐native hydrogen‐bonding interaction in helix‐coil transition of a single polypeptide as revealed by comparison between Gō‐like and non‐Gō models

To explore the role of non‐native interactions in the helix‐coil transition, a detailed comparison between a Gō‐like model and a non‐Gō model has been performed via lattice Monte Carlo simulations. Only native hydrogen bonding interactions occur in the Gō‐like model, and the non‐native ones with sequence interval more than 4 is also included into the non‐Gō model. Some significant differences between the results from those two models have been found. The non‐native hydrogen bonds were found most populated at temperature around the helix‐coil transition. The rearrangement of non‐native hydrogen bonds into native ones in the formation of α‐helix leads to the increase of susceptibility of chain conformation, and even two peaks of susceptibility of radius of gyration versus temperature exist in the case of non‐Gō model for a non‐short peptide, while just a single peak exists in the case of Gō model for a single polypeptide chain with various chain lengths. The non‐native hydrogen bonds have complicated the temperature‐dependence of Zimm‐Bragg nucleation constant. The increase of relative probability of non‐native hydrogen bonding for long polypeptide chains leads to non‐monotonous chain length effect on the transition temperature. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Hierarchical protein folding pathways: a computational study of protein fragments

We have previously presented a building block folding model. The model postulates that protein folding is a hierarchical top-down process. The basic unit from which a fold is constructed, referred to as a hydrophobic folding unit, is the outcome of combinatorial assembly of a set of "building blocks." Results obtained by the computational cutting procedure yield fragments that are in agreement with those obtained experimentally by limited proteolysis. Here we show that as expected, proteins from the same family give very similar building blocks. However, different proteins can also give building blocks that are similar in structure. In such cases the building blocks differ in sequence, stability, contacts with other building blocks, and in their 3D locations in the protein structure. This result, which we have repeatedly observed in many cases, leads us to conclude that while a building block is influenced by its environment, nevertheless, it can be viewed as a stand-alone unit. For small-sized building blocks existing in multiple conformations, interactions with sister building blocks in the protein will increase the population time of the native conformer. With this conclusion in hand, it is possible to develop an algorithm that predicts the building block assignment of a protein sequence whose structure is unknown. Toward this goal, we have created sequentially nonredundant databases of building block sequences. A protein sequence can be aligned against these, in order to be matched to a set of potential building blocks.     

Structure refinement of protein model decoys requires accurate side‐chain placement

In this study, the application of temperature‐based replica‐exchange (T‐ReX) simulations for structure refinement of decoys taken from the I‐TASSER dataset was examined. A set of eight nonredundant proteins was investigated using self‐guided Langevin dynamics (SGLD) with a generalized Born implicit solvent model to sample conformational space. For two of the protein test cases, a comparison of the SGLD/T‐ReX method with that of a hybrid explicit/implicit solvent molecular dynamics T‐ReX simulation model is provided. Additionally, the effect of side‐chain placement among the starting decoy structures, using alternative rotamer conformations taken from the SCWRL4 modeling program, was investigated. The simulation results showed that, despite having near‐native backbone conformations among the starting decoys, the determinant of their refinement is side‐chain packing to a level that satisfies a minimum threshold of native contacts to allow efficient excursions toward the downhill refinement regime on the energy landscape. By repacking using SCWRL4 and by applying the RWplus statistical potential for structure identification, the SGLD/T‐ReX simulations achieved refinement to an average of 38% increase in the number of native contacts relative to the original I‐TASSER decoy sets and a 25% reduction in values of C_α root‐mean‐square deviation. The hybrid model succeeded in obtaining a sharper funnel to low‐energy states for a modeled target than the implicit solvent SGLD model; yet, structure identification remained roughly the same. Without meeting a threshold of near‐native packing of side chains, the T‐ReX simulations degrade the accuracy of the decoys, and subsequently, refinement becomes tantamount to the protein folding problem. Proteins 2013. 2012 Published by Wiley Periodicals, Inc.     

Contact order and ab initio protein structure prediction

Although much of the motivation for experimental studies of protein folding is to obtain insights for improving protein structure prediction, there has been relatively little connection between experimental protein folding studies and computational structural prediction work in recent years. In the present study, we show that the relationship between protein folding rates and the contact order (CO) of the native structure has implications for ab initio protein structure prediction. Rosetta ab initio folding simulations produce a dearth of high CO structures and an excess of low CO structures, as expected if the computer simulations mimic to some extent the actual folding process. Consistent with this, the majority of failures in ab initio prediction in the CASP4 (critical assessment of structure prediction) experiment involved high CO structures likely to fold much more slowly than the lower CO structures for which reasonable predictions were made. This bias against high CO structures can be partially alleviated by performing large numbers of additional simulations, selecting out the higher CO structures, and eliminating the very low CO structures; this leads to a modest improvement in prediction quality. More significant improvements in predictions for proteins with complex topologies may be possible following significant increases in high-performance computing power, which will be required for thoroughly sampling high CO conformations (high CO proteins can take six orders of magnitude longer to fold than low CO proteins). Importantly for such a strategy, simulations performed for high CO structures converge much less strongly than those for low CO structures, and hence, lack of simulation convergence can indicate the need for improved sampling of high CO conformations. The parallels between Rosetta simulations and folding in vivo may extend to misfolding: The very low CO structures that accumulate in Rosetta simulations consist primarily of local up-down beta-sheets that may resemble precursors to amyloid formation.     

Unification of the folding mechanisms of non-two-state and two-state proteins

We have collected the kinetic folding data for non-two-state and two-state globular proteins reported in the literature, and investigated the relationships between the folding kinetics and the native three-dimensional structure of these proteins. The rate constants of formation of both the intermediate and the native state of non-two-state folders were found to be significantly correlated with protein chain length and native backbone topology, which is represented by the absolute contact order and sequence-distant native pairs. The folding rate of two-state folders, which is known to be correlated with the native backbone topology, apparently does not correlate significantly with protein chain length. On the basis of a comparison of the folding rates of the non-two-state and two-state folders, it was found that they are similarly dependent on the parameters that reflect the native backbone topology. This suggests that the mechanisms behind non-two-state and two-state folding are essentially identical. The present results lead us to propose a unified mechanism of protein folding, in which folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the intermediate. Apparently, two-state folding is merely a simplified version of hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.     

C‐terminal deletion of leucine‐rich repeats from YopM reveals a heterogeneous distribution of stability in a cooperatively folded protein

Terminal deletions of units from α‐helical repeat proteins have provided insight into the physical origins of their cooperativity. To test if the same principles governing cooperativity apply to β‐sheet‐containing repeat proteins, we have created a series of C‐terminal deletion constructs from a large leucine‐rich repeat (LRR) protein, YopM. We have examined the structure and stability of the resulting deletion constructs by a combination of solution spectroscopy, equilibrium denaturation studies, and limited proteolysis. Surprisingly, a high degree of nonuniformity was found in the stability distribution of YopM. Unlike previously studied repeat proteins, we identified several key LRR that on deletion disrupt nearby structure, at distances as far away as up to three repeats, in YopM. This partial unfolding model is supported by limited proteolysis studies and by point substitution in repeats predicted to be disordered as a result of deletion of adjacent repeats. We show that key internal‐ and terminal‐caps must be present to maintain the structural integrity in adjacent regions (roughly four LRRs long) of decreased stability. The finding that full‐length YopM maintains a high level of cooperativity in equilibrium unfolding underscores the importance of interfacial interactions in stabilizing locally unstable regions of structure.     

Determination of the conformation of folding initiation sites in proteins by computer simulation

Experimental evidence and theoretical models both suggest that protein folding begins by specific short regions of the polypeptide chain intermittently assuming conformations close to their final ones. The independent folding properties and small size of these folding initiation sites make them suitable subjects for computational methods aimed at deriving structure from sequence. We have used a torsion space Monte Carlo procedure together with an all-atom free energy function to investigate the folding of a set of such sites. The free energy function is derived by a potential of mean force analysis of experimental protein structures. The most important contributions to the total free energy are the local main chain electrostatics, main chain hydrogen bonds, and the burial of nonpolar area. Six proposed independent folding units and four control peptides 11–14 residues long have been investigated. Thirty Monte Carlo simulations were performed on each peptide, starting from different random conformations. Five of the six folding units adopted conformations close to the experimental ones in some of the runs. None of the controls did so, as expected. The generated conformations which are close to the experimental ones have among the lowest free energies encountered, although some less native like low free energy conformations were also found. The effectiveness of the method on these peptides, which have a wide variety of experimental conformations, is encouraging in two ways: First, it provides independent evidence that these regions of the sequences are able to adopt native like conformations early in folding, and therefore are most probably key components of the folding pathways. Second, it demonstrates that available simulation methods and free energy functions are able to produce reasonably accurate structures. Extensions of the methods to the folding of larger portions of proteins are suggested. © 1995 Wiley-Liss, Inc.     

OPUS‐SSF: A side‐chain‐inclusive scoring function for ranking protein structural models

We introduce a side‐chain‐inclusive scoring function, named OPUS‐SSF, for ranking protein structural models. The method builds a scoring function based on the native distributions of the coordinate components of certain anchoring points in a local molecular system for peptide segments of 5, 7, 9, and 11 residues in length. Differing from our previous OPUS‐CSF [Xu et al., Protein Sci. 2018; 27: 286–292], which exclusively uses main chain information, OPUS‐SSF employs anchoring points on side chains so that the effect of side chains is taken into account. The performance of OPUS‐SSF was tested on 15 decoy sets containing totally 603 proteins, and 571 of them had their native structures recognized from their decoys. Similar to OPUS‐CSF, OPUS‐SSF does not employ the Boltzmann formula in constructing scoring functions. The results indicate that OPUS‐SSF has achieved a significant improvement on decoy recognition and it should be a very useful tool for protein structural prediction and modeling.     

The role of aromatic residues in the hydrophobic core of the villin headpiece subdomain

Small autonomously folding proteins are of interest as model systems to study protein folding, as the same molecule can be used for both experimental and computational approaches. The question remains as to how well these minimized peptide model systems represent larger native proteins. For example, is the core of a minimized protein tolerant to mutation like larger proteins are? Also, do minimized proteins use special strategies for specifying and stabilizing their folded structure? Here we examine these questions in the 35‐residue autonomously folding villin headpiece subdomain (VHP subdomain). Specifically, we focus on a cluster of three conserved phenylalanine (F) residues F47, F51, and F58, that form most of the hydrophobic core. These three residues are oriented such that they may provide stabilizing aromatic–aromatic interactions that could be critical for specifying the fold. Circular dichroism and 1D‐NMR spectroscopy show that point mutations that individually replace any of these three residues with leucine were destabilized, but retained the native VHP subdomain fold. In pair‐wise replacements, the double mutant that retains F58 can adopt the native fold, while the two double mutants that lack F58 cannot. The folding of the double mutant that retains F58 demonstrates that aromatic–aromatic interactions within the aromatic cluster are not essential for specifying the VHP subdomain fold. The ability of the VHP subdomain to tolerate mutations within its hydrophobic core indicates that the information specifying the three dimensional structure is distributed throughout the sequence, as observed in larger proteins. Thus, the VHP subdomain is a legitimate model for larger, native proteins.     

One short cysteine‐rich sequence pattern – two different disulfide‐bonded structures – a molecular dynamics simulation study

The nematocyst walls of Hydra are formed by proteins containing small cysteine‐rich domains (CRDs) of ~25 amino acids. The first CRD of nematocyst outer all antigen (NW1) and the C‐terminal CRD of minicollagen‐1 (Mcol1C) contain six cysteines at identical sequence positions, however adopt different disulfide bonded structures. NW1 shows the disulfide connectivities C2‐C14/C6‐C19/C10‐C18 and Mcol1C C2‐C18/C6‐C14/C10‐C19. To analyze if both show structural preferences in the open, non‐disulfide bonded form, which explain the formation of either disulfide connectivity pattern, molecular dynamics (MD) simulations at different temperatures were performed. NW1 maintained in the 100‐ns MD simulations at 283 K a rather compact fold that is stabilized by specific hydrogen bonds. The Mcol1C structure fluctuated overall more, however stayed most of the time also rather compact. The analysis of the backbone Φ/ψ angles indicated different turn propensities for NW1 and Mcol1C, which mostly can be explained based on published data about the influence of different amino acid side chains on the local backbone conformation. Whereas a folded precursor mechanism may be considered for NW1, Mcol1C may fold according to the quasi‐stochastic folding model involving disulfide bond reshuffling and conformational changes, locking the native disulfide conformations. The study further demonstrates the power of MD simulations to detect local structural preferences in rather dynamic systems such as the open, non‐disulfide bonded forms of NW1 and Mcol1C, which complement published information from NMR backbone residual dipolar couplings. Because the backbone structural preferences encoded by the amino acid sequence embedding the cysteines influence which disulfide connectivities are formed, the data are generally interesting for a better understanding of oxidative folding and the design of disulfide stabilized therapeutics. Copyright © 2015 European Peptide Society and John Wiley & Sons, Ltd.