Search

The stability of folded proteins is critical to their biological function and for the efficacy of protein therapeutics. Predicting the energetic effects of protein mutations can improve our fundamental understanding of structural biology, the molecular basis of diseases, and possible routes to addressing those diseases with biological drugs. Identifying the effect of single amino acid point mutations on the thermodynamic equilibrium between the folded and unfolded states of a protein can pinpoint residues of critical importance that should be avoided in the process of improving other properties (affinity, solubility, viscosity, etc.) and suggest changes at other positions for increasing stability in protein engineering. Multiple computational tools have been developed for in silico predictions of protein stability in recent years, ranging from sequence-based empirical approaches to rigorous physics-based free energy methods. In this work, we show that FEP+, which is a free energy perturbation method based on all-atom molecular dynamics simulations, can provide accurate thermal stability predictions for a wide range of biologically relevant systems. Significantly, the FEP+ approach, while originally developed for relative binding free energies of small molecules to proteins and not specifically fitted for protein stability calculations, performs well compared to other methods that were fitted specifically to predict protein stability. Here, we present the broadest validation of a rigorous free energy-based approach applied to protein stability reported to date: 700+ single-point mutations spanning 10 different protein targets. Across the entire data set, we correctly classify the mutations as stabilizing or destabilizing in 84% of the cases, and obtain statistically significant predictions as compared with experiment [average error of ~1.6kcal/mol and coefficient of determination (R 2 ) of 0.40]. This study demonstrates, for the first time in a large-scale validation, that rigorous free energy calculations can be used to predict changes in protein stability from point mutations without parameterization or system-specific customization, although further improvements should be possible with additional sampling and a better representation of the unfolded state of the protein. Here, we describe the FEP+ method as applied to protein stability calculations, summarize the large-scale retrospective validation results, and discuss limitations of the method, along with future directions for further improvements., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

A new formulation of mometasone furoate (MF) for administration by dry powder inhaler (DPI) was evaluated for the treatment of asthma. A 12-week, double-blind, placebo-controlled dose-ranging study compared the efficacy and safety of three doses of MF DPI (100, 200 and 400 mcg b.i.d) with beclomethasone dipropionate (BDP) 168 mcg b.i.d. administered by metered dose inhaler in 365 adult or adolescent patients being treated with inhaled glucocorticoids. The mean change from baseline to endpoint (last treatment visit) for forced expiratory volume in 1 sec (FEV1) was the primary efficacy variable. Secondary efficacy variables included other objective measures of pulmonary function [forced vital capacity (FVC), forced expiratory flow 25-75% (FEV25-75%.) and peak expiratory flow rate (PEFR)] as well as subjective measures of therapeutic response (patients' daily evaluation of asthma symptoms and physicians' evaluation). At endpoint, all four active treatments were significantly more effective than placebo (P < 0.01) in improving FEV1 (MF DPI 5 to 7%, BDP 3%, placebo -6.6%) and all other measures of pulmonary function (FVC: MF DPI 4 to 5%, BDP 2%, placebo -4.7%; FEF25-75%: MF DPI 6 to 18%, BDP 7.5%, placebo -9.5%; PEFR (AM): MF DPI 5 to 10%, BDP 5.7%, placebo -7%). A consistent trend was observed for better improvement in patients treated with MF DPI 200 mcg b.i.d. than with MF DPI 100 mcg b.i.d., with no apparent additional benefit of MF DPI 400 mcg b.i.d. Results for the MF DPI 100 mcg b.i.d. and BDP 168 mcg b.i.d. treatment groups were similar. Patients' and physicians' subjective evaluations of symptoms found similar improvement in the MF DPI 200 and 400 mcg b.i.d. treatment groups, which were slightly better than that in the MF DPI 100 mcg b.i.d. group. Symptoms tended to worsen in the placebo group. MF DPI was well tolerated at all dose levels and the most frequently reported treatment-related adverse effects were headache, pharyngitis and oral candidiasis. No evidence of HPA-axis suppression was detected in any treatment group. In summary, all doses of MF DPI were well tolerated and significantly improved lung function and MF DPI 400 mcg (200 mcg b.i.d.) was the optimal dose in this study of patients with moderate persistent asthma.

We use the thermodynamic integration technique to calculate the free energy associated with the Tyr82-->Phe82 mutation (Y82F) in the protein FKBP-12, both free and bound to known inhibitor FK506 (tacrolimis). We find that the net difference in free energy for the two changes is 0.85 kcal/mol, with the binding of FK506 relatively more favorable for the native protein than the mutant. This net energy compares very favorably with the experimentally measured value of 0.60 kcal/mol. The results indicate that the relatively better binding of FK506 to the native protein is driven by the favorable entropy associated with the release of water molecules from the protein when the ligand binds. For a variety of reasons, modest size of the system, smallness of the change being examined, rapid convergence of the ensemble that needs to be determined and use of statistical estimates to control sampling, we have been able to carry out atypically reliable and reproducible free energy calculations for this protein system. Free energy changes for the two simulations (Y82F FKBP-12/FK506 and Y82F FKBP-12) have been calculated a total of eight times each, to compare a variety of different methodological choices and to ensure that the results are statistically significant. Detailed analysis of the free energy results has been carried out, and indicates that even when applicable, deconvolution of the total free energy into components can be very difficult, that the statistical error estimates can give a reasonable bound on the error in a simulation, and that one must be careful to use the same simulation protocol in all simulations being compared.

Two-dimensional nuclear Overhauser effect (2D NOE) spectra have been used as the experimental basis for determining the solution structure of the duplex [d(GTATATAC)]2 employing restrained molecular dynamics (rMD) simulations. The MARDIGRAS algorithm has been employed to construct a set of 233 interproton distance constraints via iterative complete relaxation matrix analysis utilizing the peak intensities from the 2D NOE spectra obtained for different mixing times and model structures. The upper and lower bounds for each of the constraints, defining size of a flat-well potential function term used in the rMD simulations, were conservatively chosen as the largest or smallest value calculated by MARDIGRAS. Three different starting models were utilized in several rMD calculations: energy-minimized A-DNA, B-DNA, and a structure containing wrinkled D-DNA in the interior. Considerable effort was made to define the appropriate force constants to be employed with the NOE terms in the AMBER force field, using as criteria the average constraints deviation, the constraints violation energy and the total energy. Of the 233 constraints, one was generated indirectly, but proved to be crucial in defining the structure: the cross-strand A5-H2 A5-H2 distance. As those two protons resonate isochronously for the self-complementary duplex, the distance cannot be determined directly. However, the general pattern of 2D NOE peak intensities, spin-lattice relaxation time (T1) values, and 31P nuclear magnetic resonance spectra lead to use of the A3-H2 A7-H2 distance for A5-H2 A5-H2 as well. Five rMD runs, with different random number seeds, were made for each of the three starting structures with the full distance constraint set. The average structure from all 15 runs and the five-structure averages from each starting structure were all quite similar. Two rMD runs for each starting structure were made with the A5-H2 A5-H2 constraint missing. The average of these six rMD runs revealed differences in structure, compared to that with the full set of constraints, primarily for the middle two base-pairs involving the missing cross-strand constraint but global deviations also were found. Conformational analysis of the resulting structures revealed that the inner four to six base-pairs differed in structure from the termini. Furthermore, an alternating structure was suggested with features alternating for the A-T and T-A steps.

A recently suggested method for refinement of structural data obtained from two-dimensional nuclear magnetic resonance experiments using molecular dynamics (MD) is explored. In this method, the time-averaged values of the appropriate internal co-ordinates of the molecule, calculated from the MD trajectory, are driven by restraints towards the experimental target values. This contrasts with most refinement procedures currently in use, where restraints are applied based on the instantaneous values of the appropriate co-ordinates. Both refinement methods are applied to the EcoRI restriction site DNA hexamer d(GAATTC)2, using target nuclear Overhauser enhancement distances derived from a one nanosecond unrestrained MD simulation of this structure. The resulting refined structures are compared to the results of the unrestrained MD trajectory, which serves as our "experimental" data. We show that although both methods can yield an average structure with the correct gross morphology, the new method allows both a much more realistic picture of inherent flexibility, and reproduces fine conformational detail better, such as sequence dependency. We also analyze the very long MD trajectory generated here (longer than any previously reported for a DNA oligomer), and find that significantly shorter simulations, typical of those frequently performed, may not yield acceptably reliable values for certain structural parameters.

We have derived a complete set of atomic charges for DNA from very high resolution, low temperature, single-crystal X-ray diffraction data, collected for a variety of nucleosides and nucleotides: cytidine; deoxycytidine 5'-monophosphate; deoxythymidine; guanosine 5'-monophosphate; deoxyadenosine; adenosine. This set of charges represents the first experimentally based parameterization of an important term in the energy function used in most modeling of DNA. The resulting charges are in good agreement with chemical intuition and experimental observations. They also agree qualitatively with the theoretically derived values now commonly used, but numerous and significant quantitative differences are observed. Possible reasons for the quantitative disagreement are discussed. An averaged set of charges (derived from the experimental results), which can be used in DNA modeling calculations, is presented.