Accurate and Rigorous Prediction of the Changes in Protein Free Energies in a Large-Scale Mutation Scan.

The prediction of mutation-induced free-energy changes in protein thermostability or protein-protein binding is of particular interest in the fields of protein design, biotechnology, and bioengineering. Herein, we achieve remarkable accuracy in a scan of 762 mutations estimating changes in protein thermostability based on the first principles of statistical mechanics. The remaining error in the free-energy estimates appears to be due to three sources in approximately equal parts, namely sampling, force-field inaccuracies, and experimental uncertainty. We propose a consensus force-field approach, which, together with an increased sampling time, leads to a free-energy prediction accuracy that matches those reached in experiments. This versatile approach enables accurate free-energy estimates for diverse proteins, including the prediction of changes in the melting temperature of the membrane protein neurotensin receptor 1.

Binding affinities controlled by shifting conformational equilibria: opportunities and limitations.

Conformational selection is an established mechanism in molecular recognition. Despite its power to explain binding events, it is hardly used in protein/ligand design to modulate molecular recognition. Here, we explore the opportunities and limitations of design by conformational selection. Using appropriate thermodynamic cycles, our approach predicts the effects of a conformational shift on binding affinity and also allows one to disentangle the effects induced by a conformational shift from other effects influencing the binding affinity. The method is assessed and applied to explain the contribution of a conformational shift on the binding affinity of six ubiquitin mutants showing different conformational shifts in six different complexes.

pmx: Automated protein structure and topology generation for alchemical perturbations.

Computational protein design requires methods to accurately estimate free energy changes in protein stability or binding upon an amino acid mutation. From the different approaches available, molecular dynamics-based alchemical free energy calculations are unique in their accuracy and solid theoretical basis. The challenge in using these methods lies in the need to generate hybrid structures and topologies representing two physical states of a system. A custom made hybrid topology may prove useful for a particular mutation of interest, however, a high throughput mutation analysis calls for a more general approach. In this work, we present an automated procedure to generate hybrid structures and topologies for the amino acid mutations in all commonly used force fields. The described software is compatible with the Gromacs simulation package. The mutation libraries are readily supported for five force fields, namely Amber99SB, Amber99SB*-ILDN, OPLS-AA/L, Charmm22*, and Charmm36.

A designed conformational shift to control protein binding specificity.

In a conformational selection scenario, manipulating the populations of binding-competent states should be expected to affect protein binding. We demonstrate how in silico designed point mutations within the core of ubiquitin, remote from the binding interface, change the binding specificity by shifting the conformational equilibrium of the ground-state ensemble between open and closed substates that have a similar population in the wild-type protein. Binding affinities determined by NMR titration experiments agree with the predictions, thereby showing that, indeed, a shift in the conformational equilibrium enables us to alter ubiquitin's binding specificity and hence its function. Thus, we present a novel route towards designing specific binding by a conformational shift through exploiting the fact that conformational selection depends on the concentration of binding-competent substates.

Influence of the nucleobase and anchimeric assistance of the carboxyl acid groups in the hydrolysis of amino acid nucleoside phosphoramidates.

Nucleoside phosphoramidates (NPs) are a class of nucleotide analogues that has been developed as potential antiviral/antitumor prodrugs. Recently, we have shown that some amino acid nucleoside phosphoramidates (aaNPs) can act as substrates for viral polymerases like HIV-1 RT. Herein, we report the synthesis and hydrolysis of a series of new aaNPs, containing either natural or modified nucleobases to define the basis for their differential reactivity. Aqueous stability, kinetics, and hydrolysis pathways were studied by NMR spectroscopy at different solution pD values (5-7) and temperatures. It was observed that the kinetics and mechanism (P-N and/or P-O bond cleavage) of the hydrolysis reaction largely depend on the nature of the nucleobase and amino acid moieties. Aspartyl NPs were found to be more reactive than Gly or ß-Ala NPs. For aspartyl NPs, the order of reactivity of the nucleobase was 1-deazaadenine>7-deazaadenine>adenine>thymine≥3-deazaadenine. Notably, neutral aqueous solutions of Asp-1-deaza-dAMP degraded spontaneously even at 4 °C through exclusive P-O bond hydrolysis (a 50-fold reactivity difference for Asp-1-deaza-dAMP vs. Asp-3-deaza-dAMP at pD 5 and 70 °C). Conformational studies by NMR spectroscopy and molecular modeling suggest the involvement of the protonated N3 atom in adenine and 1- and 7-deazaadenine in the intramolecular catalysis of the hydrolysis reaction through the rare syn conformation.

Reactivity of amino acid nucleoside phosphoramidates: a mechanistic quantum chemical study.

Recent experimental evidence (Maiti et al. Chem.-Eur. J., submitted) indicates that hydrolysis of nucleoside phosphoramidates is subjected to anchimeric influence by carboxyl moieties in the leaving group but also by the base in the nucleotide. A quantum chemical analysis of these findings is presented. First the intrinsic hydrolysis mechanism is investigated for simplified model compounds, and then both amino acid and nucleoside substituents are included. It is found that hydrolysis is assisted by the α-carboxyl group via formation of a five-membered intermediate and that the barrier for the reaction of this intermediate toward the product state can be influenced by the nucleobase. The adenine base protonated on N3 interacts with the transition state and considerably lowers the barrier for hydrolysis. The influence of several base modifications is explained by calculating the pK(a) for protonation on N3.

Normal mode analysis of Trp RNA binding attenuation protein: structure and collective motions.

The Trp RNA-binding protein (TRAP) has a toroidal topology and a perfect 11-fold symmetry, which makes it an excellent candidate for a vibrational study of elastic properties. Normal mode analysis in combination with correlation matrix calculations was used to detect collective low-frequency motions in TRAP. The results reveal the presence of highly correlated modes at the lower end of the spectrum, which directly reflect the annular and toroidal topology. The integral of the correlations over the low-frequency torsional part of the vibrational spectrum further demonstrates the relative rigidity of the 11 monomer building blocks of TRAP. The internal flexibility of each monomer and the effects of Trp-binding were also examined. The study clearly shows the determining influence of symmetry and topology on the elastic properties and also offers a detailed view on the Trp affinity of TRAP.

Calculation of binding free energies.

Molecular dynamics simulations enable access to free energy differences governing the driving force underlying all biological processes. In the current chapter we describe alchemical methods allowing the calculation of relative free energy differences. We concentrate on the binding free energies that can be obtained using non-equilibrium approaches based on the Crooks Fluctuation Theorem. Together with the theoretical background, the chapter covers practical aspects of hybrid topology generation, simulation setup, and free energy estimation. An important aspect of the validation of a simulation setup is illustrated by means of calculating free energy differences along a full thermodynamic cycle. We provide a number of examples, including protein-ligand and protein-protein binding as well as ligand solvation free energy calculations.

Molecular Dynamic Indicators of the Photoswitching Properties of Green Fluorescent Proteins.

Reversibly photoswitchable fluorescent proteins (RSFPs) are highly useful probes for a range of applications including diffraction-unlimited fluorescence microscopy. It was previously shown that reversible photoswitching not only involves cis-trans isomerization and protonation-deprotonation of the chromophore but also results in a marked difference in ß-barrel flexibility. In this work, we performed flexibility profiling and functional mode analysis (FMA) using molecular dynamics calculations to study how the flexibility of the RSFP ß-barrel influences the photoswitching properties of several fluorescent proteins. We also used Partial Least-Squared (PLS) FMA to detect promising mutation sites for the modulation of photoswitching properties of RSFPs. Our results show that the flexibility of RSFP does depend on its state with a systematically higher flexibility in the dark state compared to the bright state. In particular our method highlights the importance of Val157 in Dronpa, which upon mutation yields a striking difference in the collective motions of the two mutants. Overall, we show that PLS-FMA yields information, complementary to static structures, that can guide the rational design of fluorescent proteins.

Improving the identification rate of data independent label-free quantitative proteomics experiments on non-model crops: a case study on apple fruit.

Complex peptide extracts from non-model crops are troublesome for proper identification and quantification. To increase the identification rate of label free DIA experiments of Braeburn apple a new workflow was developed where a DDA database was constructed and linked to the DIA data. At a first level, parent masses found in DIA were searched in the DDA database based on their mass to charge ratio and retention time; at a second level, masses of fragmentation ions were compared for each of the linked spectrum. Following this workflow, a tenfold increase of peptides was identified from a single DIA run. As proof of principle, the designed workflow was applied to determine the changes during a storage experiment, achieving a two-fold identification increase in the number of significant peptides. The corresponding protein families were divided into nine clusters, representing different time profiles of changes in abundances during storage. Up-regulated protein families already show a glimpse of important pathways affecting aging during long-term storage, such as ethylene synthesis, and responses to abiotic stresses and their influence on the central metabolism. BIOLOGICAL SIGNIFICANCE: Proteomics research on non-model crops causes additional difficulties in identifying the peptides present in, often complex, samples. This work proposes a new workflow to retrieve more identifications from a set of quantitative data, based on linking DIA and DDA data at two consecutive levels. As proof of principle, a storage experiment on Braeburn apple resulted in twice as much identified storage related peptides. Important proteins involved in central metabolism and stress are significantly up-regulated after long term storage. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

The harmonic analysis of cylindrically symmetric proteins: a comparison of Dronpa and a DNA sliding clamp.

The harmonic analysis of two types of proteins with cylindrical symmetry is performed by the Standard Force Field Normal Mode Analysis and by the elastic network model. For both proteins the global elastic modes are assigned to their characteristic topologies. Dronpa is a rigid ß-barrel structure, presenting the twisting, bending and breathing motion of a cylindrical rod. The ß sliding clamp of Escherichia coli is a hexagonal ß-wheel, consisting of rigid segments. In its spectrum four classes of vibrations are identified which are characteristic of an elastic torus. Correlation diagrams and RMSF analysis are compared. The results provide not only a comprehensive validation of the use of both methods to describe the elastic behavior according to the low-frequency normal modes, but also depict the correlated motions of ß-barrel and ß-wheel proteins. The harmonic flexibility of the Dronpa protein is compared to the principal components of molecular dynamics (MD) simulation. A functionally important localized cleft opening mode is found, which is not detected by harmonic analysis.

Molecular dynamics in principal component space.

A molecular dynamics algorithm in principal component space is presented. It is demonstrated that sampling can be improved without changing the ensemble by assigning masses to the principal components proportional to the inverse square root of the eigenvalues. The setup of the simulation requires no prior knowledge of the system; a short initial MD simulation to extract the eigenvectors and eigenvalues suffices. Independent measures indicated a 6-7 times faster sampling compared to a regular molecular dynamics simulation.

Improved Replica Exchange Method for Native-State Protein Sampling.

We present a new replica exchange method, designed for optimal native state protein sampling in explicit solvent, called replica exchange with flexible tempering (REFT). The method was built upon the recently introduced replica exchange with solute tempering (REST). The potential function is adapted to direct the conformational search toward interdomain movements and the flexible portions of the protein. We demonstrate the improved sampling efficiency of REFT compared to the original REST for the bacteriophage T4 lysozyme.

Structural basis for the role of LYS220 as proton donor for nucleotidyl transfer in HIV-1 reverse transcriptase.

Biochemical studies by Castro et al. have recently revealed a crucial role for a general acid in the catalysis of nucleic acid transfer in distinct classes of polymerases. For HIV-RT LYS220 was identified as proton donor. This was unanticipated from a structural point of view, since in all ternary crystal structures of HIV-RT LYS220 are too distant from the active site to fulfill this role. In this work molecular dynamics simulations were used to reveal the dynamics of HIV-RT and to provide structural evidence for the role of LYS220. During a 1µs molecular dynamics simulation LYS220 migrates toward the active site and occupies several positions enabling direct and water mediated proton transfer towards pyrophosphate. A combination of quantum mechanical and molecular mechanics methods was used to validate the different modes of interaction.

Hydrolysis of aspartic acid phosphoramidate nucleotides: a comparative quantum chemical study.

L-Aspartic acid has recently been found to be a good leaving group during HIV reverse transcriptase catalyzed incorporation of deoxyadenosine monophosphate (dAMP) in DNA. This showed that L-Asp is a good mimic for the pyrophosphate moiety of deoxyadenosine triphosphate. The present work explores the thermochemistry and mechanism for hydrolysis of several models for L-aspartic-dAMP using B3LYP/DGDZVP, MP2/6-311++G** and G3MP2 level of theory. The effect of the new compound is gradually investigated: starting from a simple methyl amine leaving group up to the aspartic acid leaving group. The enzymatic environment was mimicked by involving two Mg(2+) ions and some important active site residues in the reaction. All reactions are compared to the corresponding O-coupled leaving group, which is methanol for methyl amine and malic acid for aspartic acid. With methyl amine as a leaving group a tautomeric associative or tautomeric dissociative mechanism is preferred and the barrier is lower than the comparable mechanism with methanol as a leaving group. The calculations on the aspartic acid in the enzymatic environment show that qualitatively the mechanism is the same as for triphosphate but the barrier for hydrolysis by the associative mechanism is higher for L-aspartic-dAMP than for L-malic-dAMP and pyrophosphate.

How Is cis-trans Isomerization Controlled in Dronpa Mutants? A Replica Exchange Molecular Dynamics Study.

The reversibly photoactivatable green fluorescent protein analog Dronpa holds great promise as a marker for various new cellular imaging applications. Using a replica exchange method which combines both Hamiltonian and temperature exchanges, the ground-state dynamics of Dronpa and two mutants with increased switching kinetics, Val157Gly and Met159Thr, were compared. The dominant chromophore state was found to be the cis isomer in all three proteins. The simulation data suggest that both mutations strongly increase the chromophore flexibility and cis-trans isomerization rate. We identify three key amino acids, Val157, Met159, and Phe173, which are able to impede the bottom hula-twist transition path, depending on their position and rotameric state. We believe our insights will help to understand the switching process and provide useful information for the design of new variants with improved fluorescence properties.