Servaas Michielssens

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.

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.

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

UNLABELLED 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.

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

Abstract 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.

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.

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.

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.

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.

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.