Search tips
Search criteria

Results 1-25 (75)

Clipboard (0)

Select a Filter Below

Year of Publication
more »
Document Types
1.  Are accurate computations of the 13C′ shielding feasible at the DFT level of theory? 
The goal of this study is twofold. First, to investigate the relative influence of the main structural factors affecting the computation of the 13C′ shielding, namely, the conformation of the residue itself and the next nearest-neighbor effects. Second, to determine whether calculation of the 13C′ shielding at the DFT level of theory, with an accuracy similar to that of the 13Cα shielding, is feasible with the existing computational resources. The DFT calculations, carried out for a large number of possible conformations of the tripeptide Ac-GXY-NMe, with different combinations of X and Y residues, enable us to conclude that the accurate computation of the 13C′ shielding for a given residue X depends on the: (i) (φ,ψ) backbone torsional angles of X; (ii) side-chain conformation of X; (iii) (φ,ψ) torsional angles of Y; and (iv) identity of residue Y. Consequently, DFT-based quantum mechanical calculations of the 13C′ shielding, with all these factors taken into account, are two orders of magnitude more CPU demanding than the computation, with similar accuracy, of the 13Cα shielding. Despite not considering the effect of the possible hydrogen bond interaction of the carbonyl oxygen, this work contributes to our general understanding of the main structural factors affecting the accurate computation of the 13C′ shielding in proteins and may spur significant progress in effort to develop new validation methods for protein structures.
PMCID: PMC3902030  PMID: 24403017
2.  Improvement of the treatment of loop structures in the UNRES force field by inclusion of coupling between backbone- and side-chain-local conformational states 
Journal of chemical theory and computation  2013;9(10):10.1021/ct4004977.
The UNited RESidue (UNRES) coarse-grained model of polypeptide chains, developed in our laboratory, enables us to carry out millisecond-scale molecular-dynamics simulations of large proteins effectively. It performs well in ab initio predictions of protein structure, as demonstrated in the last Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP10). However, the resolution of the simulated structure is too coarse, especially in loop regions, which results from insufficient specificity of the model of local interactions. To improve the representation of local interactions, in this work we introduced new side-chain-backbone correlation potentials, derived from a statistical analysis of loop regions of 4585 proteins. To obtain sufficient statistics, we reduced the set of amino-acid-residue types to five groups, derived in our earlier work on structurally optimized reduced alphabets, based on a statistical analysis of the properties of amino-acid structures. The new correlation potentials are expressed as one-dimensional Fourier series in the virtual-bond-dihedral angles involving side-chain centroids. The weight of these new terms was determined by a trial-and-error method, in which Multiplexed Replica Exchange Molecular Dynamics (MREMD) simulations were run on selected test proteins. The best average root-mean-square deviations (RMSDs) of the calculated structures from the experimental structures below the folding-transition temperatures were obtained with the weight of the new side-chain-backbone correlation potentials equal to 0.57. The resulting conformational ensembles were analyzed in detail by using the Weighted Histogram Analysis Method (WHAM) and Ward's minimum-variance clustering. This analysis showed that the RMSDs from the experimental structures dropped by 0.5 Å on average, compared to simulations without the new terms, and the deviation of individual residues in the loop region of the computed structures from their counterparts in the experimental structures (after optimum superposition of the calculated and experimental structure) decreased by up to 8 Å. Consequently, the new terms improve the representation of local structure.
PMCID: PMC3836199  PMID: 24273465
proteins; local interactions; statistical potentials; physicochemical properties
3.  A unified coarse-grained model of biological macromolecules based on mean-field multipole–multipole interactions 
Journal of Molecular Modeling  2014;20(8):2306.
A unified coarse-grained model of three major classes of biological molecules—proteins, nucleic acids, and polysaccharides—has been developed. It is based on the observations that the repeated units of biopolymers (peptide groups, nucleic acid bases, sugar rings) are highly polar and their charge distributions can be represented crudely as point multipoles. The model is an extension of the united residue (UNRES) coarse-grained model of proteins developed previously in our laboratory. The respective force fields are defined as the potentials of mean force of biomacromolecules immersed in water, where all degrees of freedom not considered in the model have been averaged out. Reducing the representation to one center per polar interaction site leads to the representation of average site–site interactions as mean-field dipole–dipole interactions. Further expansion of the potentials of mean force of biopolymer chains into Kubo’s cluster-cumulant series leads to the appearance of mean-field dipole–dipole interactions, averaged in the context of local interactions within a biopolymer unit. These mean-field interactions account for the formation of regular structures encountered in biomacromolecules, e.g., α-helices and β-sheets in proteins, double helices in nucleic acids, and helicoidally packed structures in polysaccharides, which enables us to use a greatly reduced number of interacting sites without sacrificing the ability to reproduce the correct architecture. This reduction results in an extension of the simulation timescale by more than four orders of magnitude compared to the all-atom representation. Examples of the performance of the model are presented.
FigureComponents of the Unified Coarse Grained Model (UCGM) of biological macromolecules
PMCID: PMC4139597  PMID: 25024008
Coarse-graining; Mean-field approach; Multipole–multipole interactions; Proteins; Nucleic acids; Polysaccharides
4.  Local vs global motions in protein folding 
It is of interest to know whether local fluctuations in a polypeptide chain play any role in the mechanism by which the chain folds to the native structure of a protein. This question is addressed by analyzing folding and non-folding trajectories of a protein; as an example, the analysis is applied to the 37-residue triple β-strand WW domain from the Formin binding protein 28 (FBP28) (PDB ID: 1E0L). Molecular dynamics (MD) trajectories were generated with the coarse-grained united-residue force field, and one- and two-dimensional free-energy landscapes (FELs) along the backbone virtual-bond angle θ and backbone virtual-bond-dihedral angle γ of each residue, and principal components, respectively, were analyzed. The key residues involved in the folding of the FBP28 WW domain are elucidated by this analysis. The correlations between local and global motions are found. It is shown that most of the residues in the folding trajectories of the system studied here move in a concerted fashion, following the dynamics of the whole system. This demonstrates how the choice of a pathway has to involve concerted movements in order for this protein to fold. This finding also sheds light on the effectiveness of principal component analysis (PCA) for the description of the folding dynamics of the system studied. It is demonstrated that the FEL along the PCs, computed by considering only several critically-placed residues, can correctly describe the folding dynamics.
PMCID: PMC3727290  PMID: 23914144
5.  Extension of UNRES force field to treat polypeptide chains with D-amino-acid residues 
Coarse-grained force fields for protein simulations are usually designed and parameterized to treat proteins composed of natural L-amino-acid residues. However, D-amino-acid residues occur in bacterial, fungal (e.g., gramicidins), as well as human-designed proteins. For this reason, we have extended the UNRES coarse-grained force field developed in our laboratory to treat systems with D-amino-acid residues. We developed the respective virtual-bond-torsional and double-torsional potentials for rotation about the Cα · · · Cα virtual-bond axis and two consecutive Cα · · · Cα virtual-bond axes, respectively, as functions of virtual-bond-dihedral angles γ. In turn, these were calculated as potentials of mean force (PMFs) from the diabatic energy surfaces of terminally-blocked model compounds for glycine, alanine, and proline. The potential-energy surfaces were calculated by using the ab initio method of molecular quantum mechanics at the Møller-Plesset (MP2) level of theory and the 6-31G(d,p) basis set, with the rotation angles of the peptide groups about Ci-1α⋯Ciα(λ(1)) and Ciα⋯Ci+1α(λ(2)) used as variables, and the energy was minimized with respect to the remaining degrees of freedom. The PMFs were calculated by numerical integration for all pairs and triplets with all possible combinations of types (glycine, alanine, and proline) and chirality (D or L); however, symmetry relations reduce the number of non-equivalent torsional potentials to 13 and the number of double-torsional potentials to 63 for a given C-terminal blocking group. Subsequently, one- (for torsional) and two-dimensional (for double-torsional potentials) Fourier series were fitted to the PMFs to obtain analytical expressions. It was found that the torsional potentials of the x-Y and X-y types, where X and Y are Ala or Pro, respectively, and a lowercase letter denotes D-chirality, have global minima for small absolute values of γ, accounting for the double-helical structure of gramicidin A, which is a dimer of two chains, each possessing an alternating D-Tyr-L-Tyr sequence, and similar peptides. The side-chain and correlation potentials for D-amino-acid residues were obtained by applying the reflection about the Ci-1α⋯Ciα⋯Ci+1α plane to the respective potentials for the L-amino-acid residues.
PMCID: PMC3982868  PMID: 24729761
6.  CheShift-2 Resolves a Local Inconsistency Between Two X-ray Crystal Structures 
Journal of biomolecular NMR  2012;54(2):193-198.
Since chemical shifts provide important and relatively accessible information about protein structure in solution, a Web server, CheShift-2, was developed for structure interrogation, based on a quantum mechanics database of 13Cα chemical shifts. CheShift-2 to a local inconsistency between two X-ray crystal structures (PDB IDs 1IKN and INFI) of the complex between the p65/p50 heterodimer of NFκB and its inhibitor IκBα. The availability of NMR resonance assignments that included the region of the inconsistency provided an opportunity for independent validation of the CheShift-2 server. Application of the server showed that the 13Cψ chemical shift measured for the Gly270-Pro281 sequence close to the C-terminus of IκBα were unequivocally consistent with the backbone structure modeled in the 1IKN structure, and were inconsistent with the 1NFI structure. Previous NOE measurements had demonstrated that the position of a tryptophan ring in the region immediately N-terminal in this region was not consistent with either structure. Subsequent recalculation of the local structure in this region, based on the electron density of the deposited structure factors for 1IKN, confirmed that the local backbone structure was best modeled by 1IKN, but that the rotamer of Trp258 is consistent with the 1NFI structure, including the presence of a hydrogen bond between the ring NεH of Trp258 and the backbone carbonyl group of Gln278. The consensus between all of these measures suggests that the CheShift-2 server operates well under circumstances in which backbone chemical shifts are available but where local plasticity may render X-ray structural data ambiguous.
PMCID: PMC3471536  PMID: 22945426
7.  Relation between free energy landscapes of proteins and dynamics 
By examining the molecular dynamics (MD) of protein folding trajectories, generated with the coarse-grained UNRES force field, for the B-domain of staphylococcal protein A and the triple β-strand WW domain from the Formin binding protein 28 (FBP), by principal component analysis (PCA), it is demonstrated how different free energy landscapes (FELs) and folding pathways of trajectories can be, even though they appear to be very similar by visual inspection of the time-dependence of the root-mean-square deviation (rmsd). Approaches to determine the minimal dimensionality of FELs for a correct description of protein folding dynamics are discussed. The correlation between the amplitude of the fluctuations of proteins and the dimensionality of the FELs is shown. The advantage of internal coordinate PCA over Cartesian PCA for small proteins is also illustrated.
PMCID: PMC3633568  PMID: 23620713
8.  Mean-field interactions between nucleic-acid-base dipoles can drive the formation of the double helix 
Physical review letters  2013;110(9):098101.
A proposed coarse-grained model of nucleic acids demonstrates that average interactions between base dipoles, together with chain connectivity and excluded-volume interactions, are sufficient to form double-helical structures of DNA and RNA molecules. Additionally, local interactions determine helix handedness and direction of strand packing. This result, and earlier research on reduced protein models, suggest that mean-field multipole-multipole interactions are the principal factors responsible for the formation of regular structure of biomolecules.
PMCID: PMC3627500  PMID: 23496746
9.  Hidden protein folding pathways in free-energy landscapes uncovered by network analysis 
A network analysis is used to uncover hidden folding pathways in free-energy landscapes usually defined in terms of such arbitrary order parameters as root-mean-square deviation from the native structure, radius of gyration, etc. The analysis has been applied to molecular dynamics (MD) trajectories of the B-domain of staphylococcal protein A, generated with the coarse-grained united-residue (UNRES) force field in a broad range of temperatures (270K ≤ T ≤ 325K). Thousands of folding pathways have been identified at each temperature. Out of these many folding pathways, several most probable ones were selected for investigation of the conformational transitions during protein folding. Unlike other conformational space network (CSN) methods, a node in the CSN variant implemented in this work is defined according to the nativelikeness class of the structure, which defines the similarity of segments of the compared structures in terms of secondary-structure, contact-pattern, and local geometry, as well as the overall geometric similarity of the conformation under consideration to that of the reference (experimental) structure. Our previous findings, regarding the folding model and conformations found at the folding-transition temperature for protein A (Maisuradze et al., J. Am. Chem. Soc. 132, 9444, 2010), were confirmed by the conformational space network analysis. In the methodology and in the analysis of the results, the shortest path identified by using the shortest-path algorithm corresponds to the most probable folding pathway in the conformational space network.
PMCID: PMC3376395  PMID: 22715321
10.  Determination of effective potentials for the stretching of Cα ⋯ Cα virtual bonds in polypeptide chains for coarse-grained simulations of proteins from ab initio energy surfaces of N-methylacetamide and N-acetylpyrrolidine 
The potentials of mean force (PMF’s) for the deformation of the Cα ⋯ Cα virtual bonds in polypeptide chains were determined from the diabatic energy surfaces of N-methylacetamide (modeling regular peptide groups) and N-acetylpyrrolidine (modeling the peptide groups preceding proline), calculated at the Møller-Plesset (MP2) ab initio level of theory with the 6-31G(d,p) basis set. The energy surfaces were expressed in the Cα ⋯ Cα virtual-bond length (d) and the H-N-Cα ⋯ C′ improper dihedral angle (α) that describes the pyramidicity of the amide nitrogen, or in the Cα-C′(O)-N-Cα dihedral angle (ω) and the angle α. For each grid point, the potential energy was minimized with respect to all remaining degrees of freedom. The PMF’s obtained from the (d, α) energy surfaces produced realistic free-energy barriers to the trans-cis transition (10 kcal/mol and 13 kcal/mol for the regular and proline peptide groups, respectively, compared to 12.6 – 13.9 kcal/mol and 17.3 – 19.6 kcal/mol determined experimentally for glycylglycine and N-acylprolines, respectively), while those obtained from the (ω, α) energy maps produced either low-quality PMF curves when direct Boltzmann summation was implemented to compute the PMF’s or too-flat curves with too-low free-energy barriers to the trans-cis transition if harmonic extrapolation was used to estimate the contributions to the partition function. An analytical bimodal logarithmic-Gaussian expression was fitted to the PMF’s, and the potentials were implemented in the UNRES force field. Test Langevin-dynamics simulations were carried out for the Gly-Gly and Gly-Pro dipeptides, which showed a 106-fold increase of the simulated rate of the trans-cis isomerization with respect to that measured experimentally; effectively the same result was obtained with the analytical Kramers theory of reaction rate applied to the UNRES representation of the peptide groups. Application of Kramers’ theory to compute the rate constants from the all-atom ab initio energy surfaces of the model compounds studied resulted in isomerization rates close to the experimental values, which demonstrates that the increase of the isomerization rate in UNRES simulations results solely from averaging out the secondary degrees of freedom.
PMCID: PMC3475191  PMID: 23087598
protein folding; trans-cis isomerization of peptide bonds; UNRES force field; potential of mean force; ab initio molecular quantum mechanics
11.  Effects of mutation, truncation and temperature on the folding kinetics of a WW domain 
Journal of molecular biology  2012;420(4-5):350-365.
The purpose of this work is to show how mutation, truncation and change of temperature can influence the folding kinetics of a protein. This is accomplished by principal component analysis (PCA) of molecular dynamics (MD)-generated folding trajectories of the triple β-strand WW domain from the Formin binding protein 28 (FBP) [PDB: 1E0L] and its full-size, and singly- and doubly-truncated mutants at temperatures below and very close to the melting point. The reasons for biphasic folding kinetics [i.e., coexistence of slow (three-state) and fast (two-state) phases], including the involvement of a solvent-exposed hydrophobic cluster and another delocalized hydrophobic core in the folding kinetics, are discussed. New folding pathways are identified in free-energy landscapes determined in terms of principal components for full-size mutants. Three-state folding is found to be a main mechanism for folding FBP28 WW domain and most of the full-size and truncated mutants. The results from the theoretical analysis are compared to those from experiment. Agreements and discrepancies between the theoretical and experimental results are discussed. Because of its importance in understanding protein kinetics and function, the diffusive mechanism by which FBP28 WW domain and its full-size and truncated mutants explore their conformational space is examined in terms of the mean-square displacement, (MSD), and PCA eigenvalue spectrum analyses. Subdiffusive behavior is observed for all studied systems.
PMCID: PMC3586707  PMID: 22560992
12.  Identification of formation of initial native structure in onconase from an unfolded state 
Biochemistry  2011;51(1):521-532.
In the oxidative folding of onconase, the stabilization of intermediates early in the folding process gives rise to efficient formation of its biologically active form. To identify the residues responsible for initial formation of structured intermediates, the transition from an ensemble of unstructured three-disulfide species, 3SU, to a single structured three-disulfide intermediate species, des-[30–75] or 3SF, at pH 8.0, 25 °C was examined. This transition was first monitored by far-UV CD spectroscopy at pH 8.0, 25 °C, showing that it occurs with formation of secondary structure, presumably due to native interactions. The time-dependence of formation of native-like structure was then followed by NMR spectroscopy after arresting the transition at different times by lowering the pH to 3 and then acquiring 1H,15N – HSQC spectra at pH 3, 16 °C to identify amide hydrogens that become part of native-like structure. H/D exchange was utilized to reduce the intensity of resonances from backbone amide hydrogens not involved in structure, without allowing exchange of backbone amide hydrogens involved in initial structure. Six hydrogen-bonding residues, namely, Tyr38, Lys49, Ser82, Cys90, Glu91, and Ala94 were identified as involved in the earliest detectable native-like structure before complete formation of des-[30–75], and are further stabilized later in the formation of this intermediate through S-S/SH interchange. By observing the stabilization of the structures of these residues by their neighboring residues, the initial, native-like structural elements formed in this transition have been identified, providing details of the initial events in the oxidative folding of onconase.
PMCID: PMC3254794  PMID: 22142378
13.  Like-charged residues at the ends of oligoalanine sequences might induce a chain reversal 
Biopolymers  2011;97(4):240-249.
We have examined the effect of like-charged residues on the conformation of an oligoalanine sequence. This was facilitated by CD and NMR spectroscopic and differential scanning calorimetric (DSC) measurements, and molecular dynamics calculations, of the following three alanine-based peptides: Ac-K-(A)5-K-NH2 (KAK5), Ac-K-(A)4-K-NH2 (KAK4), Ac-K-(A)3-K-NH2 (KAK3), where A and K denote alanine and lysine residues, respectively. Our earlier studies suggested that the presence of like-charged residues at the end of a short polypeptide chain composed of nonpolar residues can induce a chain reversal. For all three peptides, canonical MD simulations with NMR-derived restraints demonstrate the presence of ensembles of structures with a tendency to form a chain reversal. The KAK3 peptide exhibits a bent shape with its ends close to each other, while KAK4 and KAK5 are more extended. In the KAK5 peptide, the lysine residues do not have any influence on each other and are very mobile. Nevertheless, the tendency to form a more or less pronounced chain reversal is observed and it seems to be stable in all three peptides. This chain reversal seems to be caused by screening of the nonpolar core from the solvent by the hydrated charged residues
PMCID: PMC3371584  PMID: 22161955
alanine-based peptides; NMR spectroscopy; molecular dynamics; conformational studies; bent-like structure
14.  A Study of the α-Helical Intermediate Preceding the Aggregation of the Amino-Terminal Fragment of the β Amyloid Peptide (Aβ1–28) 
The journal of physical chemistry. B  2011;115(44):12978-12983.
The β amyloid (Aβ) peptide aggregates to form β-rich structures that are known to trigger Alzheimer’s disease. Experiments suggest that an α-helical intermediate precedes the formation of these aggregates. However, a description at the molecular level of the α-to-β transition has not been obtained. Because it has been proposed that the transition might be initiated in the amino-terminal region of Aβ, we studied the aggregation of the 28-residue amino-terminal fragment of Aβ (Aβ1–28) using molecular dynamics and a coarse-grained force field. Simulations starting from extended and helical conformations showed that oligomerization is initiated by formation of intermolecular β -sheets between the residues in the N-terminal regions. In simulations starting from the α-helical conformation, forcing residues 17–21 to remain in the initial (helical) conformation prevents aggregation but allows for the formation of dimers, indicating that oligomerization, initiated along the non-helical N-terminal regions, cannot progress without the α-to-β transition propagating along the chains.
PMCID: PMC3236598  PMID: 21939202
Amyloids; Aβ peptide; UNRES force field; molecular dynamics; Alzheimer’s disease
15.  Coarse-grained force field; general folding theory 
Physical Chemistry Chemical Physics  2011;13(38):16890-16901.
We review the coarse-grained UNited RESidue (UNRES) force field for the simulations of protein structure and dynamics, which is being developed in our laboratory over the last several years. UNRES is a physics-based force field, the prototype of which is defined as a potential of mean force of polypeptide chains in water, where all the degrees of freedom except the coordinates of α-carbon atoms and side-chain centers have been integrated out. We describe the initial implementation of UNRES to protein-structure prediction formulated as a search for the global minimum of the potential-energy function and its subsequent molecular dynamics and extensions of molecular-dynamics implementation, which enabled us to study protein-folding pathways and thermodynamics, as well as to reformulate the protein-structure prediction problem as a search for the conformational ensemble with the lowest free energy at temperatures below the folding-transition temperature. Applications of UNRES to study biological problems are also described.
PMCID: PMC3362049  PMID: 21643583
16.  Influence of the Length of the Alanine Spacer on the Acidic–Basic Properties of the Ac–Lys–(Ala)n–Lys–NH2 Peptides (n = 0, 1, 2, …, 5) 
Journal of Solution Chemistry  2012;41(10):1738-1746.
By using the potentiometric titration method, we have determined the pKa values of the two terminal lysine groups in six alanine-based peptides differing in the length of the alanine chain: Ac–Lys–Lys–NH2 (KK), Ac–Lys–Ala–Lys–NH2 (KAK), Ac–Lys–Ala–Ala–Lys–NH2 (KAK2), Ac–Lys–Ala–Ala–Ala–Lys–NH2 (KAK3), Ac–Lys–Ala–Ala–Ala–Ala–Lys–NH2 (KAK4), and Ac–Lys–Ala–Ala–Ala–Ala–Ala–Lys–NH2 (KAK5) in aqueous solution. For each compound, the model of two stepwise acid–base equilibria was fitted to the potentiometric-titration data. As expected, the pKa values of the lysine groups increase with increasing length of the alanine spacer, which means that the influence of the electrostatic field between one charged group on the other decreases with increasing length of the alanine spacer. However, for KAK3, the pKa1 value (8.20) is unusually small and pKa2 (11.41) is remarkably greater than pKa1, suggesting that the two groups are close to each other and, in turn, that a chain-reversal conformation is present for this peptide. Starting with KAK3, the differences between pKa1 and pKa2 decrease; however, for the longest peptide (KAK5), the values of pKa1 and pKa2 still differ by about 1 unit, i.e., by more than the value of log10 (4) = 0.60 that is a limiting value for the pKa difference of dicarboxylic acids with increasing methylene-spacer length. Consequently, some interactions between the two charged groups are present and, in turn, a bent shape occurs even for the longest of the peptides studied.
PMCID: PMC3510421  PMID: 23204596
Alanine-based peptides; Potentiometric titration; Acid–base equilibria
17.  Assessing the Accuracy of Protein Structures by Quantum Mechanical Computations of 13Cα Chemical Shifts 
Accounts of Chemical Research  2009;42(10):1545-1553.
Two major techniques have been used to determine the three-dimensional structures of proteins: x-ray diffraction and NMR spectroscopy. In particular, the validation of NMR-derived protein structures is one of the most challenging problems in NMR spectroscopy. Therefore, researchers have proposed a plethora of methods to determine the accuracy and reliability of protein structures. Despite these proposals, there is a growing need for more sophisticated, physics-based structure validation methods. This approach will enable us to (a) characterize the “quality” of the NMR-derived ensemble as a whole by a single parameter, (b) unambiguously identify flaws in the sequence at a residue level, and (c) provide precise information, such as sets of backbone and side-chain torsional angles, that we can use to detect local flaws.
Rather than reviewing all of the existing validation methods, this Account describes the contributions of our research group toward a solution of the long-standing problem of both global and local structure validation of NMR-derived protein structures. We emphasize a recently introduced physics-based methodology that makes use of observed and computed 13Cα chemical shifts (at the DFT level of theory) for an accurate validation of protein structures in solution and in crystals. By assessing the ability of computed 13Cα chemical shifts to reproduce observed 13Cα chemical shifts of a single or ensemble of structures in solution and in crystals, we accomplish a global validation by using the conformationally-averaged root-mean-square-deviation, ca-rmsd, as a scoring function. In addition, the method enables us to provide local validation by identifying a set of individual amino acid conformations for which the computed and observed 13Cα chemical shifts do not agree within a certain error range and may represent a non-reliable fold of the protein model.
Although it is computationally intensive, our validation method has several advantages, which we illustrate through a series of applications. This method makes use of the 13Cα chemical shifts, not shielding, that are ubiquitous to proteins and can be computed precisely from the φ, ψ, and χ torsional angles. There is no need for a priori knowledge of the oligomeric state of the protein, and no knowledge-based information or additional NMR data are required. The primary limitation at this point is the computational cost of such calculations. However, we anticipate that enhancements both in the speed of calculating these chemical shifts coupled with ever increasing computational power should soon make this a standard method accessible to the general NMR community.
PMCID: PMC3396562  PMID: 19572703
18.  Simple physics-based analytical formulas for the potentials of mean force of the interaction of amino-acid side chains in water. VI. Oppositely-charged side chains 
The journal of physical chemistry. B  2011;115(19):6130-6137.
The two-site coarse-grained model for the interactions of charged side chains, to be used with our coarse-grained UNRES force field for protein simulations proposed in the accompanying paper, has been extended to pairs of oppositely-charged side chains. The potentials of mean force of four pairs of molecules modeling charged amino-acid side chains, i.e., propionate – n-pentylamine cation (for aspartic acid – lysine), butyrate…n-pentylamine cation (for glutamic acid – lysine), propionate –1-butylguanidine (for aspartic acid – arginine), and butyrate – 1-butylguanidine (for glutamic acid – arginine) pairs were determined by umbrella-sampling molecular dynamics simulations in explicit water as functions of distance and orientation, and the analytical expression was fitted to the potentials of mean force. Compared to pairs of like-charged side chains discussed in the accompanying paper, an average quadrupole-quadrupole interaction term had to be introduced to reproduce the Coulombic interactions, and a multi-state model of charge distribution had to be introduced to fit the potentials of mean force of all oppositely-charged pairs well. The model reproduces all salt-bridge minima and, consequently, is likely to improve the performance of the UNRES force field.
PMCID: PMC3093716  PMID: 21500791
charged side chains; new model of side-chain – side-chain interactions; potential of mean force; molecular dynamics; umbrella sampling; multipole expansion
19.  Simple physics-based analytical formulas for the potentials of mean force of the interaction of amino-acid side chains in water. V. Like-charged side chains 
The journal of physical chemistry. B  2011;115(19):6119-6129.
A new model of side-chain – side-chain interactions for charged side-chains of amino acids, to be used in the UNRES force-field, has been developed, in which a side chain consists of a nonpolar and a charged site. The interaction energy between the nonpolar sites is composed of a Gay-Berne and a cavity term; the interaction energy between the charged sites consists of a Lennard-Jones term, a Coulombic term, a Generalized-Born term, and a cavity term, while the interaction energy between the nonpolar and charged sites is composed of a Gay-Berne and a polarization term. We parameterized the energy function for the models of all six pairs of natural like-charged amino-acid side chains, namely propionate-propionate (for the aspartic acid-aspartic acid pair), butyrate-butyrate (for the glutamic acid-glutamic acid pair), propionate-butyrate (for the aspartic acid-glutamic acid pair), pentylamine cation-pentylamine cation (for the lysine-lysine pair), 1-butylguanidine cation-1-butylguanidine cation (for the arginine-arginine pair), and pentylamine cation-1-butylguanidine cation (for the lysine-arginine pair). By using umbrella-sampling molecular dynamics simulations in explicit TIP3P water, we determined the potentials of mean force of the above-mentioned pairs as functions of distance and orientation and fitted analytical expressions to them. The positions and depths of the contact minima and the positions and heights of the desolvation maxima, including their dependence on the orientation of the molecules were well represented by analytical expressions for all systems. The values of the parameters of all the energy components are physically reasonable, which justifies use of such potentials in coarse-grain protein-folding simulations.
PMCID: PMC3099398  PMID: 21500792
charged side chains; new model of side-chain – side-chain interactions; potential of mean force; molecular dynamics; umbrella sampling
20.  CheShift-2: graphic validation of protein structures 
Bioinformatics  2012;28(11):1538-1539.
Summary: The differences between observed and predicted 13Cα chemical shifts can be used as a sensitive probe with which to detect possible local flaws in protein structures. For this reason, we previously introduced CheShift, a Web server for protein structure validation. Now, we present CheShift-2 in which a graphical user interface is implemented to render such local flaws easily visible. A series of applications to 15 ensembles of conformations illustrate the ability of CheShift-2 to locate the main structural flaws rapidly and accurately on a per-residue basis. Since accuracy plays a central role in CheShift predictions, the treatment of histidine (His) is investigated here by exploring which form of His should be used in CheShift-2.
Availability: CheShift-2 is free of charge for academic use and can be accessed from
Supplementary information: Supplementary data are available at the Bioinformatics online.
PMCID: PMC3356844  PMID: 22495749
21.  PDZ binding to the BAR domain of PICK1 is elucidated by coarse-grained molecular dynamics 
Journal of molecular biology  2010;405(1):298-314.
A key regulator of AMPA (α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid) receptor traffic, PICK1 is also known to interact with over 40 other proteins, including receptors, transporters, and ionic channels, and to be active mostly as a homodimer. The current lack of a complete PICK1 structure determined at atomic resolution hinders the elucidation of its functional mechanisms. Here, we identify interactions between the component PDZ and BAR domains of PICK1 by calculating possible binding sites for the PDZ domain of PICK1, PICK1-PDZ, to the homology-modeled crescent-shaped dimer of the PICK1-BAR domain using multiplexed replica-exchange molecular dynamics (MREMD) and canonical molecular dynamics (MD) simulations with the coarse-grained UNRES force field. The MREMD results show that the preferred binding site for the single PDZ domain is the concave cavity of the BAR dimer. A second possible binding site is near the N-terminus of the BAR domain that is linked directly to the PDZ domain. Subsequent short MD simulations, used to determine how the PICK1-PDZ domain moves to the preferred binding site on the BAR domain of PICK1, revealed that initial hydrophobic interactions drive the progress of the simulated binding. Thus, the concave face of the BAR dimer accommodates the PDZ domain first by weak hydrophobic interactions, and then the PDZ domain slides to the center of the concave face, where more favorable hydrophobic interactions take over.
PMCID: PMC3008210  PMID: 21050858
PICK1; binding; hydrophobic interactions; UNRES force field; molecular dynamics
22.  Modification and optimization of the united-residue (UNRES) potential-energy function for canonical simulations. I. Temperature dependence of the effective energy function and tests of the optimization method with single training proteins 
We report the modification and parameterization of the united-residue (UNRES) force field for energy-based protein-structure prediction and protein-folding simulations. We tested the approach on three training proteins separately: 1E0L (β), 1GAB (α), and 1E0G (α + β). Heretofore, the UNRES force field had been designed and parameterized to locate native-like structures of proteins as global minima of their effective potential-energy surfaces, which largely neglected the conformational entropy because decoys composed of only lowest-energy conformations were used to optimize the force field. Recently, we developed a mesoscopic dynamics procedure for UNRES, and applied it with success to simulate protein folding pathways. How ever, the force field turned out to be largely biased towards α-helical structures in canonical simulations because the conformational entropy had been neglected in the parameterization. We applied the hierarchical optimization method developed in our earlier work to optimize the force field, in which the conformational space of a training protein is divided into levels each corresponding to a certain degree of native-likeness. The levels are ordered according to increasing native-likeness; level 0 corresponds to structures with no native-like elements and the highest level corresponds to the fully native-like structures. The aim of optimization is to achieve the order of the free energies of levels, decreasing as their native-likeness increases. The procedure is iterative, and decoys of the training protein(s) generated with the energy-function parameters of the preceding iteration are used to optimize the force field in a current iteration. We applied the multiplexing replica exchange molecular dynamics (MREMD) method, recently implemented in UNRES, to generate decoys; with this modification, conformational entropy is taken into account. Moreover, we optimized the free-energy gaps between levels at temperatures corresponding to a predominance of folded or unfolded structures, as well as to structures at the putative folding-transition temperature, changing the sign of the gaps at the transition temperature. This enabled us to obtain force fields characterized by a single peak in the heat capacity at the transition temperature. Furthermore, we introduced temperature dependence to the UNRES force field; this is consistent with the fact that it is a free-energy and not a potential-energy function.
PMCID: PMC3236617  PMID: 17201450
ab initio protein folding; folding transition; thermodynamic hypothesis; potential-function optimization; hierarchical energy landscapes; foldability
23.  Mechanism of fiber assembly; treatment of Aβ-peptide aggregation with a coarse-grained united-residue force field 
Journal of molecular biology  2010;404(3):537-552.
The mechanism of growth of fibrils of the β-amyloid peptide (Aβ) was studied by means of a physics-based coarse-grained united-residue (UNRES) model and molecular dynamics (MD) simulations. To identify the mechanism of monomer addition to an Aβ1–40 fibril, an unstructured monomer was placed at a 20 Å distance from a fibril template, and allowed to interact freely with it. The monomer was not biased towards the fibril conformation, by either the force field or the MD algorithm. By using a coarse-grained model with replica exchange MD, a longer time scale was accessible making it possible to observe how the monomers probe different binding modes during their search towards the fibril conformation. Although different assembly pathways were seen, they all follow a dock-lock mechanism, with two distinct locking stages, which is consistent with data from experiments on fibril elongation. Whereas these experiments have not been able to characterize the conformations populating the different stages, we have been able to describe these different stages explicitly by following free monomers as they dock onto a fibril template and adopt the fibril conformation; i.e., we describe fibril elongation step by step, at the molecular level. During the first stage of the assembly, “docking”, the monomer tries different conformations. After docking, the monomer is locked into the fibril through two different locking stages. In the first stage the monomer forms hydrogen bonds with the fibril template along one of the strands in a two-stranded β hairpin; in the second stage, hydrogen bonds are formed along the second strand, locking the monomer into the fibril structure. The data reveal a free-energy barrier separating the two locking stages. The importance of hydrophobic interactions and hydrogen bonds in the stability of the Aβ fibril structure was examined by carrying out additional canonical MD simulations of oligomers with different numbers of chains (4 to 16 chains) with the fibril structure as the initial conformation. The data confirm that the structures are stabilized largely by hydrophobic interactions and show that the intermolecular hydrogen bonds are highly stable and contribute to the stability of the oligomers as well.
PMCID: PMC2981693  PMID: 20888834
amyloids; Aβ peptide; Alzheimer’s disease; hydrophobic interactions; UNRES force field; molecular dynamics
24.  Towards crystal structure prediction of complex organic compounds – a report on the fifth blind test 
Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories – a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.
PMCID: PMC3222142  PMID: 22101543
25.  Towards crystal structure prediction of complex organic compounds – a report on the fifth blind test 
The results of the fifth blind test of crystal structure prediction, which show important success with more challenging large and flexible molecules, are presented and discussed.
Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories – a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.
PMCID: PMC3222142  PMID: 22101543
prediction; blind test; polymorph; crystal structure prediction

Results 1-25 (75)