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.     

Mimicking the action of folding chaperones by Hamiltonian replica-exchange molecular dynamics simulations: application in the refinement of de novo models

The efficiency of using a variant of Hamiltonian replica‐exchange molecular dynamics (Chaperone H‐replica‐exchange molecular dynamics [CH‐REMD]) for the refinement of protein structural models generated de novo is investigated. In CH‐REMD, the interaction between the protein and its environment, specifically, the electrostatic interaction between the protein and the solvating water, is varied leading to cycles of partial unfolding and refolding mimicking some aspects of folding chaperones. In 10 of the 15 cases examined, the CH‐REMD approach sampled structures in which the root‐mean‐square deviation (RMSD) of secondary structure elements (SSE‐RMSD) with respect to the experimental structure was more than 1.0 Å lower than the initial de novo model. In 14 of the 15 cases, the improvement was more than 0.5 Å. The ability of three different statistical potentials to identify near‐native conformations was also examined. Little correlation between the SSE‐RMSD of the sampled structures with respect to the experimental structure and any of the scoring functions tested was found. The most effective scoring function tested was the DFIRE potential. Using the DFIRE potential, the SSE‐RMSD of the best scoring structures was on average 0.3 Å lower than the initial model. Overall the work demonstrates that targeted enhanced‐sampling techniques such as CH‐REMD can lead to the systematic refinement of protein structural models generated de novo but that improved potentials for the identification of near‐native structures are still needed. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Stabilization provided by neighboring strands is critical for the mechanical stability of proteins

Single-molecule force spectroscopy studies and steered molecular dynamics simulations have revealed that protein topology and pulling geometry play important roles in determining the mechanical stability of proteins. Most studies have focused on local interactions that are associated with the force-bearing β-strands. Interactions mediated by neighboring strands are often overlooked. Here we use Top7 and barstar as model systems to illustrate the critical importance of the stabilization effect provided by neighboring β-strands on the mechanical stability. Using single-molecule atomic force microscopy, we showed that Top7 and barstar, which have similar topology in their force-bearing region, exhibit vastly different mechanical-stability characteristics. Top7 is mechanically stable and unfolds at ∼150 pN, whereas barstar is mechanically labile and unfolds largely below 50 pN. Steered molecular dynamics simulations revealed that stretching force peels one force-bearing strand away from barstar to trigger unfolding, whereas Top7 unfolds via a substructure-sliding mechanism. This previously overlooked stabilization effect from neighboring β-strands is likely to be a general mechanism in protein mechanics and can serve as a guideline for the de novo design of proteins with significant mechanical stability and novel protein topology.     

CASP10–BCL::Fold efficiently samples topologies of large proteins

During CASP10 in summer 2012, we tested BCL::Fold for prediction of free modeling (FM) and template‐based modeling (TBM) targets. BCL::Fold assembles the tertiary structure of a protein from predicted secondary structure elements (SSEs) omitting more flexible loop regions early on. This approach enables the sampling of conformational space for larger proteins with more complex topologies. In preparation of CASP11, we analyzed the quality of CASP10 models throughout the prediction pipeline to understand BCL::Fold's ability to sample the native topology, identify native‐like models by scoring and/or clustering approaches, and our ability to add loop regions and side chains to initial SSE‐only models. The standout observation is that BCL::Fold sampled topologies with a GDT_TS score > 33% for 12 of 18 and with a topology score > 0.8 for 11 of 18 test cases de novo. Despite the sampling success of BCL::Fold, significant challenges still exist in clustering and loop generation stages of the pipeline. The clustering approach employed for model selection often failed to identify the most native‐like assembly of SSEs for further refinement and submission. It was also observed that for some β‐strand proteins model refinement failed as β‐strands were not properly aligned to form hydrogen bonds removing otherwise accurate models from the pool. Further, BCL::Fold samples frequently non‐natural topologies that require loop regions to pass through the center of the protein. Proteins 2015; 83:547–563. © 2015 Wiley Periodicals, Inc.     

Fourier transform infrared spectroscopic analysis of protein secondary structures

Infrared spectroscopy is one of the oldest and well established experimental techniques for the analysis of secondary structure of polypeptides and proteins. It is convenient, non-destructive, requires less sample preparation, and can be used under a wide variety of conditions. This review introduces the recent developments in Fourier transform infrared （FTIR） spectroscopy technique and its applications to protein structural studies. The experimental skills, data analysis, and correlations between the FTIR spectroscopic bands and protein secondary structure components are discussed. The applications of FTIR to the second- ary structure analysis, conformational changes, structural dynamics and stability studies of proteins are also discussed.     

Role of active site rigidity in activity: MD simulation and fluorescence study on a lipase mutant

Relationship between stability and activity of enzymes is maintained by underlying conformational flexibility. In thermophilic enzymes, a decrease in flexibility causes low enzyme activity while in less stable proteins such as mesophiles and psychrophiles, an increase in flexibility is associated with enhanced enzyme activity. Recently, we identified a mutant of a lipase whose stability and activity were enhanced simultaneously. In this work, we probed the conformational dynamics of the mutant and the wild type lipase, particularly flexibility of their active site using molecular dynamic simulations and time-resolved fluorescence techniques. In contrast to the earlier observations, our data show that active site of the mutant is more rigid than wild type enzyme. Further investigation suggests that this lipase needs minimal reorganization/flexibility of active site residues during its catalytic cycle. Molecular dynamic simulations suggest that catalytically competent active site geometry of the mutant is relatively more preserved than wild type lipase, which might have led to its higher enzyme activity. Our study implies that widely accepted positive correlation between conformation flexibility and enzyme activity need not be stringent and draws attention to the possibility that high enzyme activity can still be accomplished in a rigid active site and stable protein structures. This finding has a significant implication towards better understanding of involvement of dynamic motions in enzyme catalysis and enzyme engineering through mutations in active site.     

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1–2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug‐specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans‐membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α‐helical MPs as well as β‐barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge‐based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade. Proteins 2015; 83:1–24. © 2014 Wiley Periodicals, Inc.     

Recent developments in structural proteomics for protein structure determination

The major challenges in structural proteomics include identifying all the proteins on the genome-wide scale, determining their structure-function relationships, and outlining the precise three-dimensional structures of the proteins. Protein structures are typically determined by experimental approaches such as X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. However, the knowledge of three-dimensional space by these techniques is still limited. Thus, computational methods such as comparative and de novo approaches and molecular dynamic simulations are intensively used as alternative tools to predict the three-dimensional structures and dynamic behavior of proteins. This review summarizes recent developments in structural proteomics for protein structure determination; including instrumental methods such as X-ray crystallography and NMR spectroscopy, and computational methods such as comparative and de novo structure prediction and molecular dynamics simulations.     

Chaperone discovery

Molecular chaperones assist de novo protein folding and facilitate the refolding of stress‐denatured proteins. The molecular chaperone concept was coined nearly 35 years ago, and since then, tremendous strides have been made in understanding how these factors support protein folding. Here, we focus on how various chaperone proteins were first identified to play roles in protein folding. Examples are used to illustrate traditional routes of chaperone discovery and point out their advantages and limitations. Recent advances, including the development of folding biosensors and promising methods for the stabilization of proteins in vivo, provide new routes for chaperone discovery.     

Designing cooperativity into the designed protein Top7

The topology of the designed protein Top7 is not found in natural proteins. Top7 shows signatures of non‐cooperative folding in both experimental studies and computer simulations. In particular, molecular dynamics of coarse‐grained structure‐based models of Top7 show a well‐populated C‐terminal folding‐intermediate. Since most similarly sized globular proteins are cooperative folders, the non‐natural topology of Top7 has been suggested as a reason for its non‐cooperative folding. Here, we computationally examine the folding of Top7 with the intent of making it cooperative. We find that its folding cooperativity can be increased in two ways: (a) Optimization of packing interactions in the N‐terminal half of the protein enables further folding of the C‐terminal intermediate. (b) Reduction in the packing density of the C‐terminal region destabilizes the intermediate. In practice, these strategies are implemented in our Top7 model through modifications to the contact‐map. These modifications do not alter the topology of Top7 but result in cooperative folding. Amino‐acid mutations that mimic these modifications also lead to a significant increase in folding cooperativity. Finally, we devise a method to randomize the sizes of amino‐acids within the same topology, and confirm that the structure of Top7 makes its folding sensitive to packing interactions. In contrast, the ribosomal protein S6, which has secondary structure similar to Top7, has folding which is much less sensitive to packing perturbations. Thus, it should be possible to make a sequence fold cooperatively to the structure of Top7, but to do so its side‐chain packing needs to be carefully designed. Proteins 2014; 82:364–374. © 2013 Wiley Periodicals, Inc.     

Native Topology of the Designed Protein Top7 is Not Conducive to Cooperative Folding

Many single-domain proteins with <100 residues fold cooperatively; but the recently designed 92-residue Top7 protein exhibits clearly non-two-state behaviors. In apparent agreement with experiment, we found that coarse-grained, native-centric chain models, including potentials with and without elementary desolvation barriers, predicted that Top7 has a stable intermediate state in which the C-terminal fragment is folded while the rest of the chain remains disordered. We observed noncooperative folding in Top7 models that incorporated nonnative hydrophobic interactions as well. In contrast, free energy profiles deduced from models with desolvation barriers for a set of thirteen natural proteins with similar chain lengths and secondary structure elements suggested that they fold much more cooperatively than Top7. Buttressed by related studies on smaller natural proteins with chain lengths of ∼40 residues, our findings argue that the de novo native topology of Top7 likely imposed a significant restriction on the cooperativity achievable by any design for this target structure.     

iATTRACT: Simultaneous global and local interface optimization for protein–protein docking refinement

Protein‐protein interactions are abundant in the cell but to date structural data for a large number of complexes is lacking. Computational docking methods can complement experiments by providing structural models of complexes based on structures of the individual partners. A major caveat for docking success is accounting for protein flexibility. Especially, interface residues undergo significant conformational changes upon binding. This limits the performance of docking methods that keep partner structures rigid or allow limited flexibility. A new docking refinement approach, iATTRACT, has been developed which combines simultaneous full interface flexibility and rigid body optimizations during docking energy minimization. It employs an atomistic molecular mechanics force field for intermolecular interface interactions and a structure‐based force field for intramolecular contributions. The approach was systematically evaluated on a large protein‐protein docking benchmark, starting from an enriched decoy set of rigidly docked protein–protein complexes deviating by up to 15 Å from the native structure at the interface. Large improvements in sampling and slight but significant improvements in scoring/discrimination of near native docking solutions were observed. Complexes with initial deviations at the interface of up to 5.5 Å were refined to significantly better agreement with the native structure. Improvements in the fraction of native contacts were especially favorable, yielding increases of up to 70%. Proteins 2015; 83:248–258. © 2014 Wiley Periodicals, Inc.     

Conformational Effects on the Circular Dichroism of Human Carbonic Anhydrase II: A Multilevel Computational Study

Circular Dichroism (CD) spectroscopy is a powerful method for investigating conformational changes in proteins and therefore has numerous applications in structural and molecular biology. Here a computational investigation of the CD spectrum of the Human Carbonic Anhydrase II (HCAII), with main focus on the near-UV CD spectra of the wild-type enzyme and it seven tryptophan mutant forms, is presented and compared to experimental studies. Multilevel computational methods (Molecular Dynamics, Semiempirical Quantum Mechanics, Time-Dependent Density Functional Theory) were applied in order to gain insight into the mechanisms of interaction between the aromatic chromophores within the protein environment and understand how the conformational flexibility of the protein influences these mechanisms. The analysis suggests that combining CD semi empirical calculations, crystal structures and molecular dynamics (MD) could help in achieving a better agreement between the computed and experimental protein spectra and provide some unique insight into the dynamic nature of the mechanisms of chromophore interactions.     

Understanding peptide competitive inhibition of botulinum neurotoxin a binding to SV2 protein via molecular dynamics simulations

Botulinum neurotoxins (BoNTs) are known as the most toxic natural substances. Synaptic vesicle protein 2 (SV2) has been proposed to be a protein receptor for BoNT/A. Recently, two short peptides (BoNT/A‐A2 and SV2C‐A3) were designed to inhibit complex formation between the BoNT/A receptor‐binding domain (BoNT/A‐RBD) and the synaptic vesicle protein 2C luminal domain (SV2C‐LD). In this article, the two peptide complex systems are studied by molecular dynamics (MD) simulations. The structural stability analysis indicates that BoNT/A‐A2 system is more stable than SV2C‐A3 system. The conformational analysis implies that the β‐sheet in BoNT/A‐A2 system maintains its secondary structure but the two β‐strands in SV2C‐A3 system have remarkable conformational changes. Based on the calculation of hydrogen bonds, hydrophobic interactions and cation‐π interactions, it is found that the internal hydrogen bonds play crucial roles in the structural stability of the peptides. Because of the stable secondary structure, the β‐sheet in BoNT/A‐A2 system establishes effective interactions at the interface and inhibits BoNT/A‐RBD binding to SV2C‐LD. In contrast, without other β‐strands forming internal hydrogen bonds, the two isolated β‐strands in SV2C‐A3 system become the random coil. This conformational change breaks important hydrogen bonds and weakens cation‐π interaction in the interface, so the complex formation is only partially inhibited by the two β‐strands. These results are consistent with experimental studies and may be helpful in understanding the inhibition mechanisms of peptide inhibitors. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 597–608, 2015.     

Dynameomics: large-scale assessment of native protein flexibility

Structure is only the first step in understanding the interactions and functions of proteins. In this paper, we explore the flexibility of proteins across a broad database of over 250 solvated protein molecular dynamics simulations in water for an aggregate simulation time of approximately 6 micros. These simulations are from our Dynameomics project, and these proteins represent approximately 75% of all known protein structures. We employ principal component analysis of the atomic coordinates over time to determine the primary axis and magnitude of the flexibility of each atom in a simulation. This technique gives us both a database of flexibility for many protein fold families and a compact visual representation of a particular protein's native-state conformational space, neither of which are available using experimental methods alone. These tools allow us to better understand the nature of protein motion and to describe its relationship to other structural and dynamical characteristics. In addition to reporting general properties of protein flexibility and detailing many dynamic motifs, we characterize the relationship between protein native-state flexibility and early events in thermal unfolding and show that flexibility predicts how a protein will begin to unfold. We provide evidence that fold families have conserved flexibility patterns, and family members who deviate from the conserved patterns have very low sequence identity. Finally, we examine novel aspects of highly inflexible loops that are as important to structural integrity as conventional secondary structure. These loops, which are difficult if not impossible to locate without dynamic data, may constitute new structural motifs.     

Rotameric preferences of a protein spin label at edge‐strand β‐sheet sites

Protein spin labeling to yield the nitroxide‐based R1 side chain is a powerful method to measure protein dynamics and structure by electron spin resonance. However, R1 measurements are complicated by the flexibility of the side chain. While analysis approaches for solvent‐exposed α‐helical environment have been developed to partially account for flexibility, similar work in β‐sheets is lacking. The goal of this study is to provide the first essential steps for understanding the conformational preferences of R1 within edge β‐strands using X‐ray crystallography and double electron electron resonance (DEER) distance measurements. Crystal structures yielded seven rotamers for a non‐hydrogen‐bonded site and three rotamers for a hydrogen‐bonded site. The observed rotamers indicate contextual differences in R1 conformational preferences compared to other solvent‐exposed environments. For the DEER measurements, each strand site was paired with the same α‐helical site elsewhere on the protein. The most probable distance observed by DEER is rationalized based on the rotamers observed in the crystal structure. Additionally, the appropriateness of common molecular modeling methods that account for R1 conformational preferences are assessed for the β‐sheet environment. These results show that interpretation of R1 behavior in β‐sheets is difficult and indicate further development is needed for these computational methods to correctly relate DEER distances to protein structure at edge β‐strand sites.     

A molecular dynamics study of the C-terminal fragment of the L7/L12 ribosomal protein

The crystallographic dimer of the C-terminal fragment (CTF) of the L7/L12 ribosomal protein has been subjected to molecular dynamics (MD) simulations. A 90 picosecond (ps) trajectory for the protein dimer, 19 water molecules and two counter ions has been calculated at constant temperature. Effects of intermolecular interactions on the structure and dynamics have been studied. The exact crystallographic symmetry is lost and the atomic fluctuations differ from one monomer to the other. The average MD structure is more stable than the X-ray one, as judged by accessible surface area and energy calculations. Crystal (non-dimeric) interactions have been simulated in another 40 ps trajectory by using harmonic restraints to represent intermolecular hydrogen bonds. The conformational changes with respect ot the X-ray structure are then virtually suppressed.The unrestrained dimer trajectory has been scanned for cooperative motions involving secondary structure elements. The intrinsic collective motions of the monomer are transmitted via intermolecular contacts to the dimer structure.The existence of a stable dimeric form of CTF, resembling the crystallographic one, has been documented. At the cost of fairly small energy expenditure the dimer has considerable conformational flexibility. This flexibility may endow the dimer with some functional potential as an energy transducer.     

The dynamic duo: Combining NMR and small angle scattering in structural biology

Structural biology provides essential information for elucidating molecular mechanisms that underlie biological function. Advances in hardware, sample preparation, experimental methods, and computational approaches now enable structural analysis of protein complexes with increasing complexity that more closely represent biologically entities in the cellular environment. Integrated multidisciplinary approaches are required to overcome limitations of individual methods and take advantage of complementary aspects provided by different structural biology techniques. Although X‐ray crystallography remains the method of choice for structural analysis of large complexes, crystallization of flexible systems is often difficult and does typically not provide insights into conformational dynamics present in solution. Nuclear magnetic resonance spectroscopy (NMR) is well‐suited to study dynamics at picosecond to second time scales, and to map binding interfaces even of large systems at residue resolution but suffers from poor sensitivity with increasing molecular weight. Small angle scattering (SAS) methods provide low resolution information in solution and can characterize dynamics and conformational equilibria complementary to crystallography and NMR. The combination of NMR, crystallography, and SAS is, thus, very useful for analysis of the structure and conformational dynamics of (large) protein complexes in solution. In high molecular weight systems, where NMR data are often sparse, SAS provides additional structural information and can differentiate between NMR‐derived models. Scattering data can also validate the solution conformation of a crystal structure and indicate the presence of conformational equilibria. Here, we review current state‐of‐the‐art approaches for combining NMR, crystallography, and SAS data to characterize protein complexes in solution.     

TROSY and CRINEPT: NMR with large molecular and supramolecular structures in solution

TROSY and CRINEPT are new techniques for solution NMR studies of molecular and supramolecular structures. They allow the collection of high-resolution spectra of structures with molecular weights >100 kDa, significantly extending the range of macromolecular systems that can be studied by NMR in solution. TROSY has already been used to map protein-protein interfaces, to conduct structural studies on membrane proteins and to study nucleic acid conformations in multimolecular assemblies. These techniques will help us to investigate the conformational states of individual macromolecular components and will support de novo protein structure determination in large supramolecular structures.