Substrate inhibition imposes fitness penalty at high protein stability

Proteins are only moderately stable. It has long been debated whether this narrow range of stabilities is solely a result of neutral drift towards lower stability or purifying selection against excess stability is also at work - for which no experimental evidence was found so far. Here we show that mutations outside the active site in the essential E. coli enzyme adenylate kinase result in stability-dependent increase in substrate inhibition by AMP, thereby impairing overall enzyme activity at high stability. Such inhibition caused substantial fitness defects not only in the presence of excess substrate but also under physiological conditions. In the latter case, substrate inhibition caused differential accumulation of AMP in the stationary phase for the inhibition prone mutants. Further, we show that changes in flux through Adk could accurately describe the variation in fitness effects. Taken together, these data suggest that selection against substrate inhibition and hence excess stability may have resulted in a narrow range of optimal stability observed for modern proteins.Comment: 30 pages, 6 figures, 1 table, Supplementary figures and tables - 6 page

Stereochemical Criteria for Prediction of the Effects of Proline Mutations on Protein Stability

When incorporated into a polypeptide chain, proline (Pro) differs from all other naturally occurring amino acid residues in two important respects. The φ dihedral angle of Pro is constrained to values close to −65° and Pro lacks an amide hydrogen. Consequently, mutations which result in introduction of Pro can significantly affect protein stability. In the present work, we describe a procedure to accurately predict the effect of Pro introduction on protein thermodynamic stability. Seventy-seven of the 97 non-Pro amino acid residues in the model protein, CcdB, were individually mutated to Pro, and the in vivo activity of each mutant was characterized. A decision tree to classify the mutation as perturbing or nonperturbing was created by correlating stereochemical properties of mutants to activity data. The stereochemical properties including main chain dihedral angle φ and main chain amide H-bonds (hydrogen bonds) were determined from 3D models of the mutant proteins built using MODELLER. We assessed the performance of the decision tree on a large dataset of 163 single-site Pro mutations of T4 lysozyme, 74 nsSNPs, and 52 other Pro substitutions from the literature. The overall accuracy of this algorithm was found to be 81% in the case of CcdB, 77% in the case of lysozyme, 76% in the case of nsSNPs, and 71% in the case of other Pro substitution data. The accuracy of Pro scanning mutagenesis for secondary structure assignment was also assessed and found to be at best 69%. Our prediction procedure will be useful in annotating uncharacterized nsSNPs of disease-associated proteins and for protein engineering and design

Role of Water in the Enzymatic Catalysis: Study of ATP plus AMP -> 2ADP Conversion by Adenylate Kinase

The catalytic conversion ATP + AMP -> 2ADP by the enzyme adenylate kinase (ADK) involves the binding of one ATP. molecule to the LID domain and one AMP molecule to the NMP domain. The latter is followed by a. phosphate transfer and then the release of two ADP molecules. We have computed a novel two-dimensional configurational free energy surface (2DCFES), with one reaction coordinate each for the LID and the NMP domain motions, while considering explicit water interactions. Our computed 2DCFES clearly reveals the existence of a stable half-open half-closed (HOHC) intermediate stale of the enzyme. Cycling of the enzyme through the HOHC state reduces the conformational free energy barrier for. the reaction by about 20 kJ/mol. We find that the stability of the HOHC state (missed in all earlier studies with implicit solvent model) is largely because of the increase of specific interactions of the polar amino acid side chains with water, particularly with the arginine and the histidine residues. Free energy surface of the LID domain is rather rugged, which can conveniently slow down LID's conformational motion, thus facilitating a new substrate capture after the product release in the catalytic cycle

Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of adenosine triphosphate (ATP) and adenosine monophosphate (AMP) to adenosine diphosphate (ADP) by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of ATP binding domain (LID) and AMP binding domain (NMP) domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motion of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis, and multidimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two-dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed root mean square deviation variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder-type theoretical model to explain the observed dynamic coupling between the motion of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations. (C) 2011 American Institute of Physics. doi:10.1063/1.3516588

Supplementary Methods from Trade-offs between microbial growth phases lead to frequency-dependent and non-transitive selection

Document containing complete definition of the model, derivations of main results, and supplementary results (includes text and three figures)

Stereochemical criteria for prediction of the effects of proline mutations on protein stability.

When incorporated into a polypeptide chain, proline (Pro) differs from all other naturally occurring amino acid residues in two important respects. The phi dihedral angle of Pro is constrained to values close to -65 degrees and Pro lacks an amide hydrogen. Consequently, mutations which result in introduction of Pro can significantly affect protein stability. In the present work, we describe a procedure to accurately predict the effect of Pro introduction on protein thermodynamic stability. Seventy-seven of the 97 non-Pro amino acid residues in the model protein, CcdB, were individually mutated to Pro, and the in vivo activity of each mutant was characterized. A decision tree to classify the mutation as perturbing or nonperturbing was created by correlating stereochemical properties of mutants to activity data. The stereochemical properties including main chain dihedral angle phi and main chain amide H-bonds (hydrogen bonds) were determined from 3D models of the mutant proteins built using MODELLER. We assessed the performance of the decision tree on a large dataset of 163 single-site Pro mutations of T4 lysozyme, 74 nsSNPs, and 52 other Pro substitutions from the literature. The overall accuracy of this algorithm was found to be 81% in the case of CcdB, 77% in the case of lysozyme, 76% in the case of nsSNPs, and 71% in the case of other Pro substitution data. The accuracy of Pro scanning mutagenesis for secondary structure assignment was also assessed and found to be at best 69%. Our prediction procedure will be useful in annotating uncharacterized nsSNPs of disease-associated proteins and for protein engineering and design