Physical scoring function based on AMBER force field and Poisson-Boltzmann implicit solvent for protein structure prediction

A well-behaved physics-based all-atom scoring function for protein structure prediction is analyzed with several widely used all-atom decoy sets. The scoring function, termed AMBER/Poisson-Boltzmann (PB), is based on a refined AMBER force field for intramolecular interactions and an efficient PB model for solvation interactions. Testing on the chosen decoy sets shows that the scoring function, which is designed to consider detailed chemical environments, is able to consistently discriminate all 62 native crystal structures after considering the heteroatom groups, disulfide bonds, and crystal packing effects that are not included in the decoy structures. When NMR structures are considered in the testing, the scoring function is able to discriminate 8 out of 10 targets. In the more challenging test of selecting near-native structures, the scoring function also performs very well: for the majority of the targets studied, the scoring function is able to select decoys that are close to the corresponding native structures as evaluated by ranking numbers and backbone Calpha root mean square deviations. Various important components of the scoring function are also studied to understand their discriminative contributions toward the rankings of native and near-native structures. It is found that neither the nonpolar solvation energy as modeled by the surface area model nor a higher protein dielectric constant improves its discriminative power. The terms remaining to be improved are related to 1-4 interactions. The most troublesome term is found to be the large and highly fluctuating 1-4 electrostatics term, not the dihedral-angle term. These data support ongoing efforts in the community to develop protein structure prediction methods with physics-based potentials that are competitive with knowledge-based potentials.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

Molecular dynamics simulation of mammalian 15S-lipoxygenase with AMBER force field

A molecular dynamics simulation study of mononuclear iron 15S-lipoxygenase (15S-LOX) from rabbit reticulocytes was performed to investigate its structure and dynamics; newly developed AMBER force field parameters were employed for the first coordination sphere of the catalytic iron (II). The results obtained from this study demonstrate that the structural features of the catalytic iron coordination site are in good agreement with available data obtained from experiments. The motional flexibility of the N-terminal β-barrel domain is greater than the C-terminal catalytic domain; flexibility was assessed in terms of B-factors and secondary structure calculations. The significant features obtained for the relative motional flexibility of these two domains of 15S-LOX in solution as well as the isolated C-terminal domain were analyzed in terms of radius of gyration and maximum diameter, which correlated well with the structural flexibility of 15-lipoxygenase-1 in solution as probed by small-angle X-ray scattering. The motional flexibility indicates interdomain motion between the N-terminal β-barrel and the C-terminal catalytic domain; this was further verified by the evaluation of central bending in the solvated LOX molecule, which identified an unstructured stretch of amino acids as the interdomain linker. The average bending angle confirmed significant central bending between these two domains, which was linked to the high degree of motional freedom of the N-terminal β-barrel domain in aqueous solutions. This can be considered to have biological relevance for membrane binding as well as for regulating the catalytic domain.     

Gel filtration of myelin using 2,2,2-trichloroethanol as the solvent

The use of 2, 2, 2-trichloroethanol as a solvent for myelin from both the central and peripheral nervous systems is described. Concentrated, optically clear solutions of lyophilized myelin in this solvent are stable for weeks. The preparation of highly concentrated myelin proteins by gel filtration in trichloroethanol is described.     

Gel filtration of myelin using 2,2,2-trichloroethanol as the solvent.

The use of 2,2,2-trichloroethanol as a solvent for myelin from both the central and peripheral nervous systems is described. Concentrated, optically clear solutions of lyophilized myelin in this solvent are stable for weeks. The preparation of highly concentrated myelin proteins by gel filtration in trichloroethanol is described.     

Reverse bimodal action of 2,2,2-trifluoroethanol on Rumex crispus seed germination

Light-stimulated germination of R. crispus L. seed is inhibited by pretreatment with 5 m M 2,2,2-trifluoroethanol, whereas at higher concentrations (0.1 M ) germination is greatly enhanced. The action is the reverse of the usual promotion-inhibition order of bimodal response to applied substances. In darkness, 0.1 M trifluoroethanol promotes germination. Germination of eight other species (including R. obtusifolius L. and R. acetosella L.) is inhibited by levels of trifluoroethanol that promote germination of R. crispus . Ethanol inhibition of light-stimulated germination of R. crispus can be partially offset by simultaneous treatment with 0.1 M trifluoroethanol. Similarly, inhibition of the light-stimulated germination of R. crispus seeds by 5 m M trifluoroethanol can be partially offset by simultaneous treatment with 2-propanol. The inhibitory and promotive actions of trifluoroethanol may be caused by a specific enzyme inhibition and an anaesthetic effect, respectively.     

Molecular dynamics simulations of a DMSO/water mixture using the AMBER force field

Due to its protective properties of biological samples at low temperatures and under desiccation, dimethyl sulfoxide (DMSO) in aqueous solutions has been studied widely by many experimental approaches and molecular dynamics (MD) simulations. In the case of the latter, AMBER is among the most commonly used force fields for simulations of biomolecular systems; however, the parameters for DMSO published by Fox and Kollman in 1998 have only been tested for pure liquid DMSO. We have conducted an MD simulation study of DMSO in a water mixture and computed several structural and dynamical properties such as of the mean density, self-diffusion coefficient, hydrogen bonding and DMSO and water ordering. The AMBER force field of DMSO is seen to reproduce well most of the experimental properties of DMSO in water, with the mixture displaying strong and specific water ordering, as observed in experiments and multiple other MD simulations with other non-polarizable force fields.

Open image in new window

Graphical abstract Hydration structure within hydrogen-bonding distance around a DMSOmolecule

     

Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers

We present here the parmbsc0 force field, a refinement of the AMBER parm99 force field, where emphasis has been made on the correct representation of the alpha/gamma concerted rotation in nucleic acids (NAs). The modified force field corrects overpopulations of the alpha/gamma = (g+,t) backbone that were seen in long (more than 10 ns) simulations with previous AMBER parameter sets (parm94-99). The force field has been derived by fitting to high-level quantum mechanical data and verified by comparison with very high-level quantum mechanical calculations and by a very extensive comparison between simulations and experimental data. The set of validation simulations includes two of the longest trajectories published to date for the DNA duplex (200 ns each) and the largest variety of NA structures studied to date (15 different NA families and 97 individual structures). The total simulation time used to validate the force field includes near 1 mus of state-of-the-art molecular dynamics simulations in aqueous solution.     

Comparison of the effects of 2,2,2-trifluoroethanol on peptide and protein structure and function

The co-solvent 2,2,2-trifluoroethanol (TFE) has been often used to aid formation of secondary structure in solution peptides or alternately as a denaturant within protein folding studies. Hen egg white lysozyme (HEWL) and a synthetic model peptide defining HEWL helix-4 were used as comparative model systems to systematically investigate the effect of increasing TFE concentrations on the structure of proteins and peptides. HEWL was analyzed using NMR, far-UV CD and fluorescence spectroscopy; with correlation of these results towards changes in enzymatic activity and the helix-4 peptide was analysed using NMR. Data illustrates two conflicting modes of interaction: Low TFE concentrations stabilize tertiary structure, observed from an increase in the number of NMR NOE contacts. Higher TFE concentrations denatured HEWL with the loss of lysozyme tertiary structure. The effects of TFE upon secondary structural elements within HEWL are distinct from those observed for the helix-4 peptide. This illustrates a dissimilar interaction of TFE towards both protein and peptide at equivalent TFE concentrations. The concentration that TFE promotes stabilization over denaturation is likely to be protein dependent although the structural action can be extrapolated to other protein systems with implications for the use of TFE in structural stability studies.     

Analysis and purification of alcohol-sensitive chiral compounds using 2,2,2-trifluoroethanol as a modifier in supercritical fluid chromatography

2,2,2-Trifluoroethanol (TFE) is evaluated as an alternative modifier for the analysis and purification of alcohol-sensitive chiral compounds using supercritical fluid chromatography (SFC). Four chiral compounds, selected for their sensitivity to alcohols, in addition to a variety of standard chiral compounds were analyzed by SFC using TFE with polysaccharide and Pirkle-type chiral stationary phases (CSPs) to produce selectivities (alpha) and resolutions (Rs) as high as 1.4 and 7.2. A preparative isolation of 2-phenylglutaric anhydride was achieved using TFE as the mobile phase modifier to produce clean enantiomers.     

Assessing the accuracy of simple field based root strength measurements

Background and aims

Root reinforcement of slopes is a key control on landslide triggering, yet remains difficult to measure. Dozens of studies have utilised a wide range of testing methods to understand the tensile strength properties to estimate root reinforcement. We present a systematic attempt to evaluate the simple and efficient field spring scale method.

Methods

This study compared different testing methods to assess the strength of Picea sitchensis roots. We tested roots in the field using a spring scale and with two different pre-treatments in the laboratory using a universal testing machine. Root strengths were assessed under different testing conditions in laboratory and field experiments using different pre-testing treatments, stress concentration at clamps, and spacer types.

Results

Tensile strengths measured in the laboratory and field, with different spacers and clamping stress concentrations were not significantly different in our testing. Roots that were incompletely rehydrated after being dried were significantly stronger than wetter roots, suggesting moisture is one of the dominant controls on root strength.

Conclusions

Our results suggest that there is a large range of natural variability in root tensile strength. The field-based, spring scale method produces results that are indistinguishable from those of more precise universal testing machines.     

Molecular modeling of proteins: a strategy for energy minimization by molecular mechanics in the AMBER force field.

Energy minimization is an important step in molecular modeling of proteins. In this study, we sought to develop a minimization strategy which would give the best final structures with the shortest computer time in the AMBER force field. In the all-atom model, we performed energy minimization of the melittin (mostly alpha-helical) and cardiotoxin (mostly beta-sheet and beta-turns) crystal structures by both constrained and unconstrained pathways. In the constrained path, which has been recommended in the energy minimization of proteins, hydrogens were relaxed first, followed by the side chains of amino acid residues, and finally the whole molecule. Despite the logic of this approach, however, the structures minimized by the unconstrained path fit the experimental structures better than those minimized by constrained paths. Moreover, the unconstrained path saved considerable computer time. We also compared the effects of the steepest descents and conjugate gradients algorithms in energy minimization. Previously, steepest descents has been used in the initial stages of minimization and conjugate gradients in the final stages of minimization. We therefore studied the effect on the final structure of performing an initial minimization by steepest descents. The structures minimized by conjugate gradients alone resembled the structures minimized initially by the steepest descents and subsequently by the conjugate gradients algorithms. Thus an initial minimization using steepest descents is wasteful and unnecessary, especially when starting from the crystal structure. Based on these results, we propose the use of an unconstrained path and conjugate gradients for energy minimization of proteins. This procedure results in low energy structures closer to the experimental structures, and saves about 70-80% of computer time. This procedure was applied in building models of lysozyme mutants. The crystal structure of native T4 lysozyme was mutated to three different mutants and the structures were minimized. The minimized structures closely fit the crystal structures of the respective mutants (less than 0.3 A root-mean-square, RMS, deviation in the position of all heavy atoms). These results confirm the efficiency of the proposed minimization strategy in modeling closely related homologs. To determine the reliability of the united atom approximation, we also performed all of the above minimizations with united atom models. This approximation gave structures with similar but slightly higher RMS deviations than the all-atom model, but gave further savings of 60-70% in computer time. However, we feel further investigation is essential to determine the reliability of this approximation.(ABSTRACT TRUNCATED AT 400 WORDS)     

Reactivation and refolding of rabbit muscle creatine kinase denatured in 2,2,2-trifluoroethanol solutions

The unfolding and refolding of creatine kinase (ATP:creatine N-phosphotransferase (CK), EC 2.7.3.2) during denaturation and reactivation by trifluoroethanol (TFE) have been studied. Significant aggregation was observed when CK was denatured at TFE concentrations between 10% and 40% (v/v). 50% TFE (v/v) was used to study the denaturation and unfolding of CK. The activity loss of CK was a very quick process, as was the marked conformational changes during denaturation followed by fluorescence emission spectra and far-ultraviolet CD spectra. DTNB modification and size exclusion chromatography were used to find that CK dissociated and was in its monomer state after denaturation with 50% TFE. Reactivation and refolding were observed after 80-fold dilution of the denatured CK into 0.05 M Tris-HCl buffer, pH 8.0. The denatured CK recovered about 38% activity following a two phase course (k(1)=4.82+/-0.41x10(-3) s(-1), k(2)=0.60+/-0.01x10(-3) s(-1)). Intrinsic fluorescence maximum intensity changes showed that the refolding process also followed biphasic kinetics (k(1)=4.34+/-0.27x10(-3) s(-1), k(2)=0.76+/-0.02x10(-3) s(-1)) after dilution into the proper solutions. The far-ultraviolet CD spectra ellipticity changes at 222 nm during the refolding process also showed a two phase course (k(1)=4.50+/-0.07x10(-3) s(-1), k(2)=1.13+/-0.05x10(-3) s(-1)). Our results suggest that TFE can be used as a reversible denaturant like urea and GuHCl. The 50% TFE induced CK denaturation state, which was referred to as the 'TFE state', and the partially refolded CK are compared with the molten globule state. The aggregation caused by TFE during denaturation is also discussed in this paper.     

Conformational changes and inactivation of bovine carbonic anhydrase II in 2,2,2-trifluoroethanol solutions

Changes in unfolding and enzymatic activity of bovine carbonic anhydrase II (BCA II) in different concentrations of 2,2,2-trifluoroethanol (TFE) were investigated by 1-anilino-8-naphthalenesulfonate (ANS) fluorescence emission spectra, far-UV CD spectra, and enzyme activity. The results showed that the activity and conformation of BCA II changed according to the concentration of TFE. Significant aggregation was observed when BCA II was denatured at TFE concentrations between 10 and 35% (v/v). When the concentration of TFE exceeded 40%, the aggregation of BCA II was not very obvious. The activity of BCA II decreased almost to zero as the TFE concentration reached 26%. The ANS fluorescence spectra indicated the tertiary conformations of BCA II were more stable in solutions with TFE concentrations lower than 15% (v/v) and higher than 40% (v/v). Far-UV CD spectra showed that high concentrations (higher than 25%) of TFE could induce BCA II to form more alpha-helix structures and caused these structures to be in relatively stable states. The native conformation of BCA II being destroyed after its inactivity indicated that the active sites of BCA II is situated in a limited region and has more flexibility than the whole enzyme molecule.     

Spectroscopic evidence for backbone desolvation of helical peptides by 2,2,2-trifluoroethanol: an isotope-edited FTIR study

2,2,2-Trifluoroethanol (TFE) is widely used to induce helix formation in peptides and proteins, but the mechanism behind this effect is still poorly understood. Several recent papers have proposed that TFE acts by selectively desolvating the peptide backbone groups of the helix state. Infrared (IR) spectroscopy of the amide I band of polypeptides can be used to probe both secondary structure and backbone solvation, making this technique well suited for addressing the effect of TFE on polypeptide conformation. In this paper, we report the IR spectra as a function of TFE concentration for an alanine-rich peptide based on the repeat (AAKAA)(n)(). The IR spectra confirm that TFE desolvates the helical state of the peptide to a greater extent than the random coil state. Moreover, using a series of specifically (13)C-labeled peptides, the precise residues desolvated in the presence of TFE were identified. The residues most desolvated by TFE are the alanines located at position i - 4 in the sequence, where i is a lysine residue. This pattern of desolvation is consistent with molecular dynamics simulations which predict strong interactions between the lysine side chain at position n and the backbone carbonyl of the alanine at position i - 4. This is the first direct spectroscopic evidence of specific desolvation of helix backbone atoms in model alanine-rich peptides.     

Distinguish protein decoys by using a scoring function based on a new AMBER force field, short molecular dynamics simulations, and the generalized born solvent model

Recent works have shown the ability of physics-based potentials (e.g., CHARMM and OPLS-AA) and energy minimization to differentiate the native protein structures from large ensemble of non-native structures. In this study, we extended previous work by other authors and developed an energy scoring function using a new set of AMBER parameters (also recently developed in our laboratory) in conjunction with molecular dynamics and the Generalized Born solvent model. We evaluated the performance of our new scoring function by examining its ability to distinguish between the native and decoy protein structures. Here we present a systematic comparison of our results with those obtained with use of other physics-based potentials by previous authors. A total of 7 decoy sets, 117 protein sequences, and more than 41,000 structures were evaluated. The results of our study showed that our new scoring function represents a significant improvement over previously published physics-based scoring functions.     

Interaction of 2,2,2-trifluoroethanol with proteins: calorimetric, densimetric and surface tension approach

The thermal denaturation of hen egg-white lysozyme was studied in the presence of 2,2,2-trifluoroethanol (TFE) at various pH values using micro differential scanning calorimetry. Quantitative thermodynamic parameters accompanying the thermal transitions were evaluated. It is observed that thermal unfolding of lysozyme in the presence of TFE upto a concentration of 4.0 mol dm(-3) follows a two-state denaturation mechanism as indicated by the equality of van't Hoff and calorimetric enthalpies. The finer details of interaction were studied by measuring the partial molar volume of some constituent amino acids and glycine peptides from water to aqueous TFE at 298.15 K. The physico-chemical properties of aqueous TFE: apparent molar heat capacities, apparent molar volumes and surface tension were measured to understand the intrinsic properties of the cosolvent as well. From the correlation among the thermal unfolding data on lysozyme in aqueous TFE, calculated preferential interaction parameters, physico chemical properties of aqueous TFE and partial molar volumes of transfer, it is concluded that both solvent mediated effect and direct interaction constitute the mechanism of TFE-protein interactions.     

Destabilisation of native tertiary structural interactions is linked to helix-induction by 2,2,2-trifluoroethanol in proteins

The effect of 2,2,2-trifluoroethanol (TFE) on the structure of an all β-sheet protein, cardiotoxin analogue 111 (CTX III) from the Taiwan cobra (Naja naja atra) is studied. It is found that high concentrations ( > 80% v/v) of TFE induced a β-sheet to -helix structural transition. It is found that in denatured and reduced CTX III (rCTX III) helical conformation is induced even upon addition of low concentrations ( > 10% v/v) of TFE. Using three other proteins, namely, ribonuclease A (RNase A), lysozyme and -lactalbumin, it is been observed that helix-induction by TFE is intricately linked to drastic destabilization of native tertiary structural interactions in the proteins.     

Equilibrium and folding simulations of NS4B H2 in pure water and water/2,2,2-trifluoroethanol mixed solvent: examination of solvation models

The structural stability and preference of a protein are highly sensitive to the environment accommodating it. In this work, the solvation effect on the structure and folding dynamics of a small peptide, NS4B H2, was studied by computer simulation. The native structure of NS4B H2 was solved previously in 50 % v/v water/2,2,2-trifluoroethanol (TFE) mixed solvent. In this work, both pure water and water/TFE cosolvent were utilized. The force field parameters for water were taken from the TIP3P water model, and those for TFE were generated following the routine of the general AMBER force field (GAFF). The simulated structure of NS4B H2 in the mixed solvent is quite in line with experimental data, while in pure water it undergoes a large structural deformation. The generalized Born (GB) model was also investigated by tuning the dielectric constant to match experimental measurements. However, the results show that its performance was less satisfactory. Two independent direct folding simulations of NS4B H2 in explicit water/TFE cosolvent were carried out, both of which resulted in successful folding. Investigation of the distribution of solvent molecules around the peptide indicates that folding is triggered by the aggregation of TFE on the peptide surface.