Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates

Determination of Trends Underlying Aspartic Acid Isomerization in Intact Proteins Reveals Unusually Rapid Isomerization of Tau

Spontaneous chemical modifications in long-lived proteins can potentially change protein structure in ways that impact proteostasis and cellular health. For example, isomerization of aspartic acid interferes with protein turnover and is anticorrelated with cognitive acuity in Alzheimer's disease. However, few isomerization rates have been determined for Asp residues in intact proteins. To remedy this deficiency, we used protein extracts from SH-SY5Y neuroblastoma cells as a source of a complex, brain-relevant proteome with no baseline isomerization. Cell lysates were aged in vitro to generate isomers, and extracted proteins were analyzed by data-independent acquisition (DIA) liquid chromatography-mass spectrometry (LC-MS). Although no Asp isomers were detected at day 0, isomerization increased over time and was quantifiable for 105 proteins by day 50. Data analysis revealed that the isomerization rate is influenced by both primary sequence and secondary structure, suggesting that steric hindrance and backbone rigidity modulate isomerization. Additionally, we examined lysates extracted under gentle conditions to preserve protein complexes and found that protein-protein interactions often slow isomerization. Base catalysis was explored as a means to accelerate Asp isomerization due to findings of accelerated asparagine deamidation. However, no substantial rate enhancement was found for isomerization, suggesting fundamental differences in acid-base chemistry. With an enhanced understanding of Asp isomerization in proteins in general, we next sought to better understand Asp isomerization in tau. In vitro aging of monomeric and aggregated recombinant tau revealed that tau isomerizes significantly faster than any similar protein within our data set, which is likely related to its correlation with cognition in Alzheimer's disease.

Unusually Rapid Isomerization of Aspartic Acid in Tau

Spontaneous chemical modifications in long-lived proteins can potentially change protein structure in ways that impact proteostasis and cellular health. For example, isomerization of aspartic acid interferes with protein turnover and is anticorrelated with cognitive acuity in Alzheimer's disease. However, few isomerization rates have been determined for Asp residues in intact proteins. To remedy this deficiency, we used protein extracts from SH-SY5Y neuroblastoma cells as a source of a complex, brain-relevant proteome with no baseline isomerization. Cell lysates were aged in vitro to generate isomers, and extracted proteins were analyzed by data-independent acquisition (DIA) liquid chromatography-mass spectrometry (LC-MS). Although no Asp isomers were detected at Day 0, isomerization increased across time and was quantifiable for 105 proteins by Day 50. Data analysis revealed that isomerization rate is influenced by both primary sequence and secondary structure, suggesting that steric hindrance and backbone rigidity modulate isomerization. Additionally, we examined lysates extracted under gentle conditions to preserve protein complexes and found that protein-protein interactions often slow isomerization. Base catalysis was explored as a means to accelerate Asp isomerization due to findings of accelerated asparagine deamidation. However, no substantial rate enhancement was found for isomerization, suggesting fundamental differences in acid-base chemistry. With an enhanced understanding of Asp isomerization in proteins in general, we next sought to better understand Asp isomerization in tau. In vitro aging of monomeric and aggregated recombinant tau revealed that tau isomerizes significantly faster than any similar protein within our dataset, which is likely related to its correlation with cognition in Alzheimer's disease.

Insights into the Main Protease of SARS-CoV-2: Thermodynamic Analysis, Structural Characterization, and the Impact of Inhibitors

The main protease (Mpro) of SARS-CoV-2 is an essential enzyme for coronaviral maturation and is the target of Paxlovid, which is currently the standard-of-care treatment for COVID-19. There remains a need to identify new inhibitors of Mpro as viral resistance to Paxlovid emerges. Here, we report the use of native mass spectrometry coupled with 193 nm ultraviolet photodissociation (UVPD) and integrated with other biophysical tools to structurally characterize Mpro and its interactions with potential covalent inhibitors. The overall energy landscape was obtained using variable temperature nanoelectrospray ionization (vT-nESI), thus providing quantitative evaluation of inhibitor binding on the stability of Mpro. Thermodynamic parameters extracted from van't Hoff plots revealed that the dimeric complexes containing each inhibitor showed enhanced stability through increased melting temperatures as well as overall lower average charge states, giving insight into the basis for inhibition mechanisms.

Unraveling the peculiarities and development of novel inhibitors of leishmanial arginyl-tRNA synthetase

Aminoacylation by tRNA synthetase is a crucial part of protein synthesis and is widely recognized as a therapeutic target for drug development. Unlike the arginyl-tRNA synthetases (ArgRSs) reported previously, here, we report an ArgRS of Leishmania donovani (LdArgRS) that can follow the canonical two-step aminoacylation process. Since a previously uncharacterized insertion region is present within its catalytic domain, we implemented the splicing by overlap extension PCR (SOE-PCR) method to create a deletion mutant (ΔIns-LdArgRS) devoid of this region to investigate its function. Notably, the purified LdArgRS and ΔIns-LdArgRS exhibited different oligomeric states along with variations in their enzymatic activity. The full-length protein showed better catalytic efficiency than ΔIns-LdArgRS, and the insertion region was identified as the tRNA binding domain. In addition, a benzothiazolo-coumarin derivative (Comp-7j) possessing high pharmacokinetic properties was recognized as a competitive and more specific inhibitor of LdArgRS than its human counterpart. Removal of the insertion region altered the mode of inhibition for ΔIns-LdArgRS and caused a reduction in the inhibitor's binding affinity. Both purified proteins depicted variances in the secondary structural content upon ligand binding and thus, thermostability. Apart from the trypanosomatid-specific insertion and Rossmann fold motif, LdArgRS revealed typical structural characteristics of ArgRSs, and Comp-7j was found to bind within the ATP binding pocket. Furthermore, the placement of tRNAArg near the insertion region enhanced the stability and compactness of LdArgRS compared to other ligands. This study thus reports a unique ArgRS with respect to catalytic as well as structural properties, which can be considered a plausible drug target for the derivation of novel anti-leishmanial agents.

Mapping HDX-MS Data to Protein Conformations through Training Ensemble-Based Models

An original approach that adopts machine learning inference to predict protein structural information using hydrogen-deuterium exchange mass spectrometry (HDX-MS) is described. The method exploits an in-house optimization program that increases the resolution of HDX-MS data from peptides to amino acids. A system is trained using Gradient Tree Boosting as a type of machine learning ensemble technique to assign a protein secondary structure. Using limited training data we generate a discriminative model that uses optimized HDX-MS data to predict protein secondary structure with an accuracy of 75%. This research could form the basis for new methods exploiting artificial intelligence to model protein conformations by HDX-MS.

Ordered-domain unfolding of thermophilic isolated β subunit ATP synthase

The flexibility of the ATP synthase's β subunit promotes its role in the ATP synthase rotational mechanism, but its domains stability remains unknown. A reversible thermal unfolding of the isolated β subunit (Tβ) of the ATP synthase from Bacillus thermophilus PS3, tracked through circular dichroism and molecular dynamics, indicated that Tβ shape transits from an ellipsoid to a molten globule through an ordered unfolding of its domains, preserving the β-sheet residual structure at high temperature. We determined that part of the stability origin of Tβ is due to a transversal hydrophobic array that crosses the β-barrel formed at the N-terminal domain and the Rossman fold of the nucleotide-binding domain (NBD), while the helix bundle of the C-terminal domain is the less stable due to the lack of hydrophobic residues, and thus the more flexible to trigger the rotational mechanism of the ATP synthase.

pH-Dependent Changes in Structural Stabilities of Bt Cry1Ac Toxin and Contrasting Model Proteins following Adsorption on Montmorillonite

The environmental fate of insecticidal Cry proteins, including time-dependent conservation of biological properties, results from their structural stability in soils. The complex cascade of reactions involved in biological action requires Cry proteins to be in solution. However, the pH-dependent changes in conformational stability and the adsorption-desorption mechanisms of Cry protein on soil minerals remain unclear. We used Derjaguin-Landau-Verwey-Overbeek (DLVO) calculation and differential scanning calorimetry to interpret the driving forces and structural stabilities of Cry1Ac and two contrasting model proteins adsorbed by montmorillonite. The structural stability of Cry1Ac is closer to that of the "hard" protein, α-chymotrypsin, than that of the "soft" bovine serum albumin (BSA). The pH-dependent adsorption of Cry1Ac and α-chymotrypsin could be explained by DLVO theory, whereas the BSA adsorption deviated from it. Patch-controlled electrostatic attraction, hydrophobic effects, and entropy changes following protein unfolding on a mineral surface could contribute to Cry1Ac adsorption. Cry1Ac, like chymotrypsin, was partly denatured on montmorillonite, and its structural stability decreased with an increase in pH. Moreover, small changes in the conformational heterogeneity of both Cry1Ac and chymotrypsin were observed following adsorption. Conversely, adsorbed BSA was completely denatured regardless of the solution pH. The moderate conformational rearrangement of adsorbed Cry1Ac may partially explain why the insecticidal activity of Bt toxin appears to be conserved in soils, albeit for a relatively short time period.

Long-Term Biocatalyst Performance: Mechanistic Prediction and Continuous Non-Isothermal Testing

The assessment of the operational stability of biocatalysts by conventional direct determination of the total turnover number (TTN), a useful indicator of lifetime biocatalyst productivity, via continuous isothermal experiments tends to be time-consuming, material-intensive, and prone to disturbances, especially in case of rather stable catalysts. In the present work, we present and validate two alternative methods for estimating the TTN of a biocatalyst for any desired operating temperature. The first method is a mechanistic approach, built upon mathematical derivation of enzyme deactivation models derived from first principles, in which TTN can be calculated from two straightforward isothermal biochemical batch measurements. The second method relies on a few non-isothermal, continuous-mode experiments in conjunction with mathematical modeling to determine the intrinsic deactivation parameters of the biocatalyst. We verify both methods on the test case of TEM-1 β-lactamase-catalyzed penicillin G (Pen G) hydrolysis. Both alternative methods provide estimates of TTN which are typically within a factor of two to five or less of the values measured directly via lengthy, costly, and error-prone conventional isothermal aging tests. Therefore, both the mechanistic approach and the non-isothermal continuous approach are extremely valuable tools to enable calculation of catalyst cost contribution in continuous processing and to eliminate underperforming candidates in search of the most stable biocatalyst.

Wine Thermosensitive Proteins Adsorb First and Better on Bentonite during Fining: Practical Implications and Proposition of Alternative Heat Tests

Bentonite fining is the most popular treatment used to remove proteins in white and rosé wines. The usual heat test used to adjust the bentonite dose consists of heating the wine during 30 min at 80 °C. At this temperature, all of the proteins are unfolded, and this can lead to an overestimation of the dose. We have shown that proteins adsorb on bentonite in a specific order and, more importantly, that the proteins responsible for haze formation adsorb first. Fluorescence spectroscopy showed that this is due to the structural properties of proteins, which can be classified as hard and soft proteins. Alternative heat tests were performed at a lower temperature (40 °C) and showed a better correlation with accelerated aging. These tests were also less dependent upon the wine pH.

On the Thermal Stability of O6-Methylguanine-DNA Methyltransferase from Archaeon Pyrococcus kodakaraensis by Molecular Dynamics Simulations

We have employed molecular dynamics simulations to analyze the thermal stability of the O⁶-methylguanine-DNA methyltransferase (MGMT) protein, both hyperthermophilic archaeon Pyrococcus kodakaraensis (Pk-MGMT) and its mesophilic homologue pair, obtained from enterobacterium Escherichia coli (AdaC). This theoretical study was done at three different temperatures: 302, 371, and 450 K. The molecular dynamics has been performed in explicit aqueous solvent during a period of time of 95 ns, including periodic boundary conditions and constant pressure. The same procedure has been used for both proteins, and each simulation has been carried out by triplicate. Hence, we performed 18 simulations. In this way, we have done different analyses to explore the factors that may affect the thermal stability of Pk-MGMT. The structural behavior was analyzed using indicators such as root-mean-square deviation, radius of gyration, solvent-accessible surface area, hydrogen bonds, native contacts, secondary structure, and salt bridge formation. The results showed that when the temperature increases, the global atomic fluctuations increase too, which suggests that both proteins lose thermal stability, but as expected, this fact is highlighted in AdaC. Moreover, the contacts of the native state in AdaC are considerably lower than those found in Pk-MGMT at 450 K. Also, the structural studies showed that conserved and nonconserved salt bridges kept close contacts with the Pk-MGMT protein at high temperatures. These interaction types act as molecular staples and are mainly responsible to provide thermostability to the hyperthermophilic protein.

Comparison of α-Helix and β-Sheet Structure Adaptation to a Quantum Dot Geometry: Toward the Identification of an Optimal Motif for a Protein Nanoparticle Cover

While quantum dots (QDs) are useful as fluorescent labels, their application in biosciences is limited due to the stability and hydrophobicity of their surface. In this study, we tested two types of proteins for use as a cover for spherical QDs, composed of cadmium selenide. Pumilio homology domain (Puf), which is mostly α-helical, and leucine-rich repeat (LRR) domain, which is rich in β-sheets, were selected to determine if there is a preference for one of these secondary structure types for nanoparticle covers. The protein sequences were optimized to improve their interaction with the surface of QDs. The solubilization of the apoproteins and their assembly with nanoparticles required the application of a detergent, which was removed in subsequent steps. Finally, only the Puf-based cover was successful enough as a QD hydrophilic cover. We showed that a single polypeptide dimer of Puf, PufPuf, can form a cover. We characterized the size and fluorescent properties of the obtained QD:protein assemblies. We showed that the secondary structure of the Puf proteins was not destroyed upon contact with the QDs. We demonstrated that these assemblies do not promote the formation of reactive oxygen species during illumination of the nanoparticles. The data represent advances in the effort to obtain a stable biocompatible cover for QDs.

Dihedral angle preferences of amino acid residues forming various non-local interactions in proteins

In theory, a polypeptide chain can adopt a vast number of conformations, each corresponding to a set of backbone rotation angles. Many of these conformations are excluded due to steric overlaps. Ramachandran and coworkers were the first to look into this problem by plotting backbone dihedral angles in a two-dimensional plot. The conformational space in the Ramachandran map is further refined by considering the energetic contributions of various non-bonded interactions. Alternatively, the conformation adopted by a polypeptide chain may also be examined by investigating interactions between the residues. Since the Ramachandran map essentially focuses on local interactions (residues closer in sequence), out of interest, we have analyzed the dihedral angle preferences of residues that make non-local interactions (residues far away in sequence and closer in space) in the folded structures of proteins. The non-local interactions have been grouped into different types such as hydrogen bond, van der Waals interactions between hydrophobic groups, ion pairs (salt bridges), and ππ-stacking interactions. The results show the propensity of amino acid residues in proteins forming local and non-local interactions. Our results point to the vital role of different types of non-local interactions and their effect on dihedral angles in forming secondary and tertiary structural elements to adopt their native fold.

Temperature dependent dynamics of DegP-trimer: A molecular dynamics study

DegP is a heat shock protein from high temperature requirement protease A family, which reacts to the environmental stress conditions in an ATP independent way. The objective of the present analysis emerged from the temperature dependent functional diversity of DegP between chaperonic and protease activities at temperatures below and above 28 °C, respectively. DegP is a multimeric protein and the minimal functional unit, DegP-trimer, is of great importance in understanding the DegP pathway. The structural aspects of DegP-trimer with respect to temperature variation have been studied using molecular dynamics simulations (for 100 ns) and principal component analysis to highlight the temperature dependent dynamics facilitating its functional diversity. The DegP-trimer revealed a pronounced dynamics at both 280 and 320 K, when compared to the dynamics observed at 300 K. The LA loop is identified as the highly flexible region during dynamics and at extreme temperatures, the residues 46–80 of LA loop express a flip towards right (at 280) and left ( at 320 K) with respect to the fixed β-sheet connecting the LA loop of protease for which Phe46 acts as one of the key residues. Such dynamics of LA loop facilitates inter-monomeric interaction with the PDZ1 domain of the neighbouring monomer and explains its active participation when DegP exists as trimer. Hence, the LA loop mediated dynamics of DegP-trimer is expected to provide further insight into the temperature dependent dynamics of DegP towards the understanding of its assembly and functional diversity in the presence of substrate.

Molecular dynamics perspective on the protein thermal stability: a case study using SAICAR synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) Cα, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Deactivation of TEM-1 β-Lactamase Investigated by Isothermal Batch and Non-Isothermal Continuous Enzyme Membrane Reactor Methods

The thermal deactivation of TEM-1 β-lactamase was examined using two experimental techniques: a series of isothermal batch assays and a single, continuous, non-isothermal assay in an enzyme membrane reactor (EMR). The isothermal batch-mode technique was coupled with the three-state "Equilibrium Model" of enzyme deactivation, while the results of the EMR experiment were fitted to a four-state "molten globule model". The two methods both led to the conclusions that the thermal deactivation of TEM-1 β-lactamase does not follow the Lumry-Eyring model and that the T(eq) of the enzyme (the point at which active and inactive states are present in equal amounts due to thermodynamic equilibrium) is at least 10 °C from the T(m) (melting temperature), contrary to the idea that the true temperature optimum of a biocatalyst is necessarily close to the melting temperature.

High temperature unfolding of Bacillus anthracis amidase-03 by molecular dynamics simulations

The stability of amidase-03 structure (a cell wall hydrolase protein) from Bacillus anthracis was studied using classical molecular dynamics (MD) simulation. This protein (GenBank accession number: NP_844822) contains an amidase-03 domain which is known to exhibit the catalytic activity of N-acetylmuramoyl-L-alanine amidase (digesting MurNAc-Lalanine linkage of bacterial cell wall). The amidase-03 enzyme has stability at high temperature due to the core formed by the combination of several secondary structure elements made of beta-sheets. We used root-mean-square-displacement (RMSD) of the simulated structure from its initial state to demonstrate the unfolding of the enzyme using its secondary structural elements. Results show that amidase-03 unfolds in transition state ensemble (TSE). The data suggests that alpha-helices unfold before beta-sheets from the core during simulation.

Dynamics of lysozyme structure network: probing the process of unfolding

Recently we showed that the three-dimensional structure of proteins can be investigated from a network perspective, where the amino acid residues represent the nodes in the network and the noncovalent interactions between them are considered for the edge formation. In this study, the dynamical behavior of such networks is examined by considering the example of T4 lysozyme. The equilibrium dynamics and the process of unfolding are followed by simulating the protein at 300 K and at higher temperatures (400 K and 500 K), respectively. The snapshots of the protein structure from the simulations are represented as protein structure networks in which the strength of the noncovalent interactions is considered an important criterion in the construction of edges. The profiles of the network parameters, such as the degree distribution and the size of the largest cluster (giant component), were examined as a function of interaction strength at different temperatures. Similar profiles are seen at all the temperatures. However, the critical strength of interaction (Icritical) and the size of the largest cluster at all interaction strengths shift to lower values at 500 K. Further, the folding/unfolding transition is correlated with contacts evaluated at Icritical and with the composition of the top large clusters obtained at interaction strengths greater than Icritical. Finally, the results are compared with experiments, and predictions are made about the residues, which are important for stability and folding. To summarize, the network analysis presented in this work provides insights into the details of the changes occurring in the protein tertiary structure at the level of amino acid side-chain interactions, in both the equilibrium and the unfolding simulations. The method can also be employed as a valuable tool in the analysis of molecular dynamics simulation data, since it captures the details at a global level, which may elude conventional pairwise interaction analysis.

On the temperature and pressure dependence of a range of properties of a type of water model commonly used in high-temperature protein unfolding simulations

Molecular dynamics simulations of protein folding and unfolding are often carried out at temperatures (400-600 K) that are much higher than physiological or room temperature to speed up the (un)folding process. Use of such high temperatures changes both the protein and solvent properties considerably, compared to physiological or room temperature. Water models designed for use in conjunction with biomolecules, such as the simple point charge (SPC) model, have generally been calibrated at room temperature and pressure. To determine the distortive effect of high simulation temperatures on the behavior of such "room temperature" water models, the structural, dynamic, and thermodynamic properties of the much-used SPC water model are investigated in the temperature range from 300 to 500 K. Both constant pressure and constant volume conditions, as used in protein simulations, were analyzed. We found that all properties analyzed change markedly with increasing temperature, but no phase transition in this temperature range was observed.

Effective energy function for proteins in solution

A Gaussian solvent-exclusion model for the solvation free energy is developed. It is based on theoretical considerations and parametrized with experimental data. When combined with the CHARMM 19 polar hydrogen energy function, it provides an effective energy function (EEF1) for proteins in solution. The solvation model assumes that the solvation free energy of a protein molecule is a sum of group contributions, which are determined from values for small model compounds. For charged groups, the self-energy contribution is accounted for primarily by the exclusion model. Ionic side-chains are neutralized, and a distance-dependent dielectric constant is used to approximate the charge-charge interactions in solution. The resulting EEF1 is subjected to a number of tests. Molecular dynamics simulations at room temperature of several proteins in their native conformation are performed, and stable trajectories are obtained. The deviations from the experimental structures are similar to those observed in explicit water simulations. The calculated enthalpy of unfolding of a polyalanine helix is found to be in good agreement with experimental data. Results reported elsewhere show that EEF1 clearly distinguishes correctly from incorrectly folded proteins, both in static energy evaluations and in molecular dynamics simulations and that unfolding pathways obtained by high-temperature molecular dynamics simulations agree with those obtained by explicit water simulations. Thus, this energy function appears to provide a realistic first approximation to the effective energy hypersurface of proteins.

Tertiary contacts in alpha-lactalbumin at pH 7 and pH 2: a molecular dynamics study

Molecular dynamics simulations of alpha-lactalbumin were performed under conditions of neutral pH and low pH in order to study the acid-induced molten globule state. Through the use of experimental techniques such as NMR and CD spectroscopy, molten globules have been characterized as being compact intermediates with secondary structure similar to that of the native protein but with tertiary structure that is disordered. The detailed structure of the molten globule state is unknown, however. Through the use of computer simulations we can study the structural changes which occur upon lowering pH. The simulations presented here differ from previous unfolding simulations in two important ways: the electrostatic interactions are treated more accurately than ever before, and artificially high temperatures are not used to force the protein to unfold. Simulations of 880 psec each were run at pH 7 (control simulation) and pH 2. We concentrate on the interesting changes in the tertiary interactions within the protein with lowering of pH. In particular, there is a loss of native tertiary contacts in the beta domain and interdomain region, and a large decrease in interdomain hydrogen bonds.

Identification of the binding surface on beta-lactamase for GroEL by limited proteolysis and MALDI-mass spectrometry

Escherichia coli beta-lactamase, alone or as a complex with GroEL at 48 degreesC, was partially digested with trypsin, endoproteinase Glu-C, or thermolysin. Peptides were analyzed by matrix-assisted laser desorption and ionization mass spectrometry and aligned with the known sequence. From the protease cleavage sites which become protected upon binding and those which become newly accessible, a model of the complex is proposed in which the carboxy-terminal helix has melted, two loops form the binding interface and the large beta-sheet become partially uncovered by the slight dislocation of other structural elements. This explains how hydrophobic surface on the substrate protein can become accessible while scarcely disrupting the hydrogen bond network of the native structure. An analysis of the GroEL-bound peptides bound after digestion of the beta-lactamase showed no obvious sequence motifs, indicating that binding is provided by hydrophobic patches in the three-dimensional structure.

Calculations on folding of segment B1 of streptococcal protein G

We present an investigation of the folding thermodynamics and mechanism of segment B1 of streptococcal protein G. Molecular dynamics simulations of the fully solvated protein are used to probe thermodynamically significant states at different stages of folding. We performed several unfolding simulations to generate a database of initial conditions. The database is analyzed and clustered. The cluster centers extracted from this database were then used as starting points for umbrella sampling of the folding free energy landscape under folding conditions. The resulting sampling was combined with the weighted histogram analysis method. One and two-dimensional free energy surfaces were constructed along several order parameters and used to analyze the folding process. Our findings indicate that an initial collapse precedes the formation of significant native structure. Elements of local structure originate in the regions of the protein shown to have higher H/2H exchange protection factors in early stages of folding. A non-native contact, observed experimentally at the N terminus of the alpha-helix in a peptide excised from the protein, is seen to pre-organize the chain in early stages of folding. Collapse and early structure formation yields a compact globule with a significant number of water molecules present. Desolvation of the protein core is coincident with the final stages of folding from the compact state.

Molecular dynamics simulations of hydrophobic collapse of ubiquitin

Nine nonnative conformations of ubiquitin, generated during two different thermal denaturation trajectories, were simulated under nearly native conditions (62 degrees C). The simulations included all protein and solvent atoms explicitly, and simulation times ranged from 1-2.4 ns. The starting structures had alpha-carbon root-mean-square deviations (RMSDs) from the crystal structure of 4-12 A and radii of gyration as high as 1.3 times that of the native state. In all but one case, the protein collapsed when the temperature was lowered and sampled conformations as compact as those reached in a control simulation beginning from the crystal structure. In contrast, the protein did not collapse when simulated in a 60% methanol:water mixture. The behavior of the protein depended on the starting structure: during simulation of the most native-like starting structures (<5 A RMSD to the crystal structure) the RMSD decreased, the number of native hydrogen bonds increased, and the secondary and tertiary structure increased. Intermediate starting structures (5-10 A RMSD) collapsed to the radius of gyration of the control simulation, hydrophobic residues were preferentially buried, and the protein acquired some native contacts. However, the protein did not refold. The least native starting structures (10-12 A RMSD) did not collapse as completely as the more native-like structures; instead, they experienced large fluctuations in radius of gyration and went through cycles of expansion and collapse, with improved burial of hydrophobic residues in successive collapsed states.

Simulations of the structural and dynamical properties of denatured proteins: the "molten coil" state of bovine pancreatic trypsin inhibitor

The dynamic nature of denatured, unfolded proteins makes it difficult to characterize their structures experimentally. To complement experiment and to obtain more detailed information about the structure and dynamic behavior of the denatured state, we have performed eleven 2.5 ns molecular dynamics simulations of reduced bovine pancreatic trypsin inhibitor (BPTI) at high temperature in water and a control simulation at 298 K, for a total of 30 ns of simulation time. In a neutral pH environment (acidic residues ionized), the unfolded protein structures were compact with an average radius of gyration 9% greater than the native state. The compact conformations resulted from the transient formation of non-native hydrophobic clusters, turns and salt bridges. However, when the acidic residues were protonated, the protein periodically expanded to a radius of gyration of 18 to 20 A. The early steps in unfolding were similar in the different simulations until passing through the major transition state of unfolding. Afterwards, unfolding proceeded through one of two general pathways with respect to secondary structure: loss of the C-terminal helix followed by loss of beta-structure or the opposite. To determine whether the protein preferentially sampled particular conformational substates in the denatured state, pairwise Calpha root-mean-square deviations were measured between all structures, but similar structures were found between only two trajectories. Yet, similar composite properties (secondary structure content, side-chain and water contacts, solvent accessible surface area, etc.) were observed for the structures that unfolded through different pathways. Somewhat surprisingly, the unfolded structures are in agreement with both past experiments suggesting that reduced BPTI is a random coil and more recent experiments providing evidence for non-random structure, demonstrating how ensembles of fluctuating structures can give rise to experimental observables that are seemingly at odds.

Molecular dynamics simulations of the unfolding of barnase in water and 8 M aqueous urea

Molecular dynamics simulations of barnase have been conducted both in water and in 8 M urea solution for 500 ps at 25 degrees C and for 2000 ps at 85 degrees C. The final structure of the aqueous simulation at room temperature matches closely the structure obtained by NMR and the experimentally observed protections from isotopic exchange. The comparison of the structures generated by the aqueous simulation at 85 degrees C reveals a trajectory composed of groups of geometrically related structures separated by narrow regions of rapid change in structure. The first of these regions displays changes in backbone rmsd to the crystal structure and solvent-accessible area suggestive of a transition state, while the properties observed during the final 300 ps of the simulation are consistent with a stable intermediate. These assignments were confirmed by calculation of the "progress along the reaction coordinate" phi-values using an empirical equation based on a linear response method. The pathway of unfolding defined in this fashion agrees well with the experimental results of site-directed mutagenesis in terms of secondary structure content of the transition state and the intermediate and reproduces the relative stability of the different elements of secondary structure. The results of the simulations in urea suggest a mechanism at the molecular level for its well-known enhancement of the denaturation of proteins. The analysis of radial distribution functions shows that the first solvation shell of the protein is enriched in urea relative to the bulk solvent. The displacement of water molecules allows greater exposure of hydrophobic side chains, as witnessed particularly in the analysis of solvent-accessible surface areas at the higher temperature. Almost all urea molecules in the first shell form at least one hydrogen bond with the protein. They provide a more favorable environment for accommodation of the remaining water molecules, and they facilitate the separation of secondary structure elements by acting as a bridge between groups previously forming intraprotein hydrogen bonds.

Dynamics of ribonuclease A and ribonuclease S: computational and experimental studies

RNase S is a complex consisting of two proteolytic fragments of RNase A: the S peptide (residues 1-20) and S protein (residues 21-124). RNase S and RNase A have very similar X-ray structures and enzymatic activities. Previous experiments have shown increased rates of hydrogen exchange and greater sensitivity to tryptic cleavage for RNase S relative to RNase A. It has therefore been asserted that the RNase S complex is considerably more dynamically flexible than RNase A. In the present study we examine the differences in the dynamics of RNase S and RNase A computationally, by MD simulations, and experimentally, using trypsin cleavage as a probe of dynamics. The fluctuations around the average solution structure during the simulation were analyzed by measuring the RMS deviation in coordinates. No significant differences between RNase S and RNase A dynamics were observed in the simulations. We were able to account for the apparent discrepancy between simulation and experiment by a simple model. According to this model, the experimentally observed differences in dynamics can be quantitatively explained by the small amounts of free S peptide and S protein that are present in equilibrium with the RNase S complex. Thus, folded RNase A and the RNase S complex have identical dynamic behavior, despite the presence of a break in polypeptide chain between residues 20 and 21 in the latter molecule. This is in contrast to what has been widely believed for over 30 years about this important fragment complementation system.

The role of protein-solvent interactions in protein unfolding

Protein unfolding occurs when the balance of forces between the protein's interaction with itself and the protein's interaction with its environment is disrupted. The disruption of this balance of forces may be as simple as a perturbance of the normal water structure around the protein. A decrease in the normal water-water interaction will result in an increase in the relative interaction of water with the protein. An increase in the number of interactions between water and the protein may initiate a protein's unfolding. This model for protein unfolding is supported by a range of recent experimental and computational data.

Molecular dynamics simulations of adipocyte lipid-binding protein: effect of electrostatics and acyl chain unsaturation

Molecular dynamics (MD) simulations have been performed on adipocyte lipid-binding protein, using the apo and holo forms, bound with stearic and oleic acid. The contribution of electrostatics to protein dynamics and ligand stabilization was assayed by perturbing the electrostatic charge of Arg106 and Arg126 (positive-->neutral) and the fatty acid (132H) headgroup (negative-->neutral). MD simulations for charged holo forms demonstrated significantly greater electrostatic binding energy and a more stabilized hydrogen bond network than simulations performed using neutral forms. Electrostatics, however, appeared to have little effect on fatty acid behavior, e.g., fluctuation of the dihedral head group; number of dihedral transitions within the acyl chain; and change in the end-to-end distance for fatty acid. Instead, fatty acid behavior appeared to be dictated by the presence or absence of an unsaturated bond within the acyl chain. A significantly greater number of transitions were observed during MD simulations in oleic than stearic acid. In addition, significantly greater fluctuation was observed for oleic acid, within the C2 headgroup and C9 and C11 dihedrals (which lie adjacent to the olefin bond of oleic acid). The dynamic behavior of the acyl chain may thereby be more a property of van der Waals contact, and the degree of acyl chain unsaturation, than a function of electrostatics. In the absence of fatty acid, an increase in distance between guanidino carbon centered atoms of Arg126 and Arg106 was observed during MD simulations of the charged apo form. This effect not observed with the neutral apo form or in any of the holo complexes and, presumbably, was a result of repulsion between the negatively charged arginine sidechains. Conserved waters reflected substantially lower mean-square displacement (msd) in all simulations, except the neutral apo form. This suggests that the presence of either charged amino acids or lipid provides increased order for water within the binding pocket. These results provide a dynamic perspective of the interactive nature within the FABP binding pocket regulated in a complex manner by the electrostatics within the binding cavity, acyl chain structure and behavior, and water energetics.

Theoretical studies of protein folding and unfolding

The mechanism of protein folding is being investigated theoretically by the use of both simplified and all-atom models of the polypeptide chain. Lattice heteropolymer simulations of the folding process have led to proposals for the folding mechanism and for the resolution of the Levinthal paradox. Both stability and rapid folding have been shown in model studies to result from the presence of a pronounced global energy minimum corresponding to the native state. Concomitantly, molecular dynamics simulations with detailed atomic models have been used to analyze the initial stages of protein unfolding. Results concerning possible folding intermediates and the role of water in the unfolding process have been obtained. The two types of theoretical approaches are providing information essential for an understanding of the mechanism of protein folding and are useful for the design of experiments to study the mechanism in different proteins.