Single-site mutation and secondary structure stability: an isodesmic reaction approach. The case of unnatural amino acid mutagenesis Ala-->Lac

A method is described to evaluate backbone interactions in proteins via computational unnatural amino acid mutagenesis. Several N-acetyl polyalanyl amides (AcA(n)NH(2)) were optimized in the representative helical (3(10)-, 4(13)-, and a "hybrid" kappa-helix, n = 7, 9, 10, 14) and hairpin (two- and three-stranded antiparallel beta-sheets with type I turns betaalphaalphaepsilon, n = 6, 9, 10) conformations, and extended conformers of N-acetyl polyalanyl methylamides (n = 2, 3) were used to derive multistranded beta-sheet fragments. Subsequently, each residue of every model structure was substituted, one at a time, with l-lactic acid. The resulting mutant structures were again optimized, and group-transfer energies DeltaE(GT) were obtained as heats of the isodesmic reactions: AcA(n)NHR + AcOMe --> AcA(x)LacA(y)NHR + AcNHMe (R = H, CH(3)). These group-transfer energies correlate with the degree of charge polarization of the substituted peptide linkages as measured by the difference Deltae in H and O Mulliken populations in HN-C=O and with the H-bond distances in the "wild-type" structures. A good correlation obtains for the HF/3-21G and B3LYP/6-31G* group-transfer energies. The destabilization effects are interpreted in terms of loss of interstrand and intrastrand H-bonds, decrease in Lewis basicity of the C=O group, and O...O repulsion. On the basis of several comparisons of Ala --> Lac DeltaE(GT)'s with heats of the NH --> CH(2) substitutions, the latter contribution is estimated (B3LYP/6-31G*) to range between 1.5 and 2.4 kcal mol(-1), a figure close to the recent experimental DeltaDeltaG(o) value of 2.6 kcal mol(-1) (McComas, C. C.; Crowley, B. M.; Boger, D. L. J. Am.Chem. Soc. 2003, 125, 9314). The partitioning yields the following maximum values of the electronic association energy of H-bonds in the examined sample of model structures (B3LYP/6-31G* estimates): 3(10)-helix D(e) = -1.7 kcal mol(-1), alpha-helix D(e) = -3.8 kcal mol(-1), beta-sheet D(e) = -6.1 kcal mol(-1). The premise of experimental evaluations of the backbone-backbone H-bonding that Ala --> Lac substitution in proteins is isosteric (e.g., Koh, J. T.; Cornish, V. W.; Schultz, P. G. Biochemistry 1997, 36, 11314) is often but not always corroborated. Examination of the integrity of H-bonding pattern and phi(i), psi(i) distribution identified several mutants with significant distortions of the "wild-type" structure resulting inter alia from the transitions between i, i + 3 and i, i + 4 H-bonding in helices, observed previously in the crystallographic studies of depsipeptides (Ohyama, T.; Oku, H.; Hiroki, A.; Maekawa, Y.; Yoshida, M.; Katakai, R. Biopolymers 2000, 54, 375; Karle, I. L.; Das, C.; Balaram, P. Biopolymers 2001, 59, 276). Thus, the isodesmic reaction approach provides a simple way to gauge how conformation of the polypeptide chain and dimensions of the H-bonding network affect the strength of backbone-backbone C=O...HN bonds. The results indicate that the stabilization provided by such interactions increases on going from 3(10)-helix to alpha-helix to beta-sheet.     

Coordination number of percolating clusters

The coordination number of percolating clusters and their backbones, on the square and simple cubic lattices, is investigated numerically. At the percolation threshold, the average coordination number of the clusters is 2.239 ± 0.012 and 2.079 ± 0.016 for d = 2 and d = 3, respectively. The corresponding numbers for the backbones are 2.402 ± 0.016 and 2.246 ± 0.034.     

Boundary conditions at a fluid-solid interface

We study the boundary conditions at a fluid-solid interface using molecular dynamics simulations covering a broad range of fluid-solid interactions and fluid densities and both simple and chain-molecule fluids. The slip length is shown to be independent of the type of flow, but rather is related to the fluid organization near the solid, as governed by the fluid-solid molecular interactions.     

Hidden valley structure of Ising spin glasses

The hidden valley structure of the phase space of short-range spin glasses is investigated. In the spirit of a similar investigation for liquids by Stillinger and Weber, valleys are defined as the set of states which by steepest descent lead to the same metastable state. The lifetime of valleys and the size of spin clusters involved in intervalley transitions on thermal Monte Carlo trajectories are calculated as a function of temperature.     

Mechanical unfolding of ubiquitin molecules

Mechanical stretching of ubiquitin and of its several repeats are studied through molecular-dynamics simulations. A Go-type model [H. Abe and N. Go, Biopolymers 20, 1013 (1981)] with a realistic contact map and with Lennard-Jones contact interactions is used. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on polyubiquitins stretched by various linkages. The terminal-to-terminal stretching of polyubiquitin results in peak forces similar to those measured for titin-based polyproteins and of a magnitude that matches measurements. Consistent with the experimental measurements, the simulated peak forces depend on the pulling speed logarithmically when thermal fluctuations are explicitly introduced. These results validate the application of topology-based models in the study of the mechanical stretching of proteins.     

Tightening of knots in proteins

We perform theoretical studies of stretching of 20 proteins with knots within a coarse-grained model. The knot's ends are found to jump to well defined sequential locations that are associated with sharp turns, whereas in homopolymers they diffuse around and eventually slide off. The waiting times of the jumps are increasingly stochastic as the temperature is raised. Knots typically do not return to their native locations when a protein is released after stretching.     

Coarse-grained molecular dynamics simulations of nanopatterning with multivalent inks

A coarse-grained molecular dynamics (MD) model is developed to study the multivalent, or multisite, binding of small functionalized dendrimer molecules to beta-cyclodextrin-terminated self-assembled monolayers, the so-called "molecular printboards" used to print "ink" molecules on surfaces with a high degree of positional control and specificity. Some current and future bionanotechnology applications are in the creation of nanoparticle assemblies, directed protein assembly, platforms for biosensing, and cell:surface attachment. The coarse-grained model allows us to probe up to microsecond timescales and model ink diffusion, crucial for the application of the printboard in, for example, medical diagnostics. Recent all-atom MD simulations identified and quantified the molecular strain limiting the stability of nanopatterns created with small dendrimer inks, and explained the different patterns obtained experimentally with different dendrimer inks. In the present work, the all-atom simulations are "scaled up" to longer timescales via coarse graining, without incurring significant additional computational expense, and, crucially, without significant loss in atom-scale detail, the coarse-grained MD simulations yielding properties similar to those obtained from the all-atom simulations. The anchoring of the ink molecules to the monolayer is of multivalent nature and the degree of multivalency shows a sharp dependence on temperature, control of temperature thus providing a further operational "switch" for directed molecular assembly. The computational protocol developed can, in principle, be extended to model any multivalent assembly, for example, virus-cell complexation.