Cyclic retro-inverso dipeptides with two aromatic side chains. II. Conformational analysis.

The conformations of cis and trans cyclic retro-inverso dipeptides--2-[(4-hydroxy)benzyl]-5-benzyl-4,6(1H,2H,3H,5H)-pyrimidinedi one (c[mTyr-gPhe]), and 2-benzyl-5-amino-5-[(4-hydroxy)benzyl]-4,6(1H,2H,3H,5H)-pyrimidinedione (c[mTyr-gPhe]), and 2-benzyl-5-amino-5-[(4-hydroxy)benzyl]-4,6(1H,2H,3H,5H)-pyrimidinedione (c[(alpha-amino)mTyr-gPhe])--and the parent cyclic dipeptides--c[tyrosyl-phenylalanine] (cis-c[L-Tyr-L-Phe]) and c[tyrosyl-D-phenylalanine] (trans-c[L-Tyr-D-Phe])--were studied by using 1H-nmr spectroscopy and semiempirical energy calculations. In the cis compounds of all the cyclic retro-inverso and parent dipeptides, the most stable conformer has both aromatic side chains sharing the space over the backbone ring in a "face-to-face" fashion. All the trans compounds predominantly assume a "sandwich" conformation in which the two aromatic rings are folded back over the backbone ring on opposite sides. However, different conformational preferences were observed for the backbones between the retro-inverso and parent cyclic dipeptides. The parent cyclic dipeptide trans-c[L-Tyr-D-Phe] adopts two types of boat structures with different side-chain orientations in almost equal amounts: one with the Tyr side chain in a pseudoaxial position and the Phe side chain in a pseudoequatorial position, the other with the Tyr side chain in a pseudoequatorial position and the Phe side chain in a pseudoaxial position. On the other hand, the cyclic retro-inverso dipeptides trans-c[mPhe-gTyr] and trans c[mTyr-gPhe] assume only one type of boat structure in which the malonyl side chain is in a pseudoequatorial and the gem-diamino side chain is in a pseudoaxial position. In addition to the preferred conformations, the conformational energies of the C alpha--C beta bonds in the malonyl and gem-diamino residues were estimated from the temperature variation of vicinal 1H--1H coupling constants for the H--C alpha--C beta--H groupings observed for the trans isomers of cyclic retro-inverso dipeptides. The energies were evaluated to be 1.1 and 1.8 kcal mol-1 for the malonyl and gem-diamino residues, respectively. Applying these energies to the parent cyclic dipeptide trans-c[L-Tyr-D-Phe], the observed fractions of three side-chain conformations are reasonably reproduced. The conformational energies as well as conformational properties of the molecules estimated in this investigation may be useful to refine force constants for both parent and retro-inverso peptides with aromatic side chains.     

Thermodynamic stability of beta-peptide helices and the role of cyclic residues

Beta-peptides are emerging as an attractive class of peptidomimetic molecules. In contrast to naturally occurring alpha-peptides, short oligomers of beta-amino acids (comprising just 4-6 monomers) exhibit stable secondary structures that make them amenable for quantitative, concerted experimental and theoretical studies of the effects of particular chemical interactions on structure. In this work, molecular simulations are used to study the thermodynamic stability of helical conformations formed by beta-peptides containing varying proportions of acyclic (beta(3)) and cyclic (ACH) residues. More specifically, several beta-peptides differing only in their content of cyclic residues are considered in this work. Previous computational studies of beta-peptides have relied mostly on energy minimization of molecular dynamics simulations. In contrast, our study relies on density-of-states based Monte Carlo simulations to calculate the free energy and examine the stability of various folded structures of these molecules along a well-defined order parameter. By resorting to an expanded-ensemble formalism, we are able to determine the free energy required to unfold specific molecules, a quantity that could be measured directly through single-molecule force spectroscopy. Simulations in both implicit and explicit solvents have permitted a systematic study of the role of cyclic residues and electrostatics on the stability of secondary structures. The molecules considered in this work are shown to exhibit stable H-14 helical conformations and, in some cases, relatively stable H-12 conformations, thereby suggesting that solvent quality may be used to manipulate the hydrogen-bonding patterns and structure of these peptides.     

Conformational studies on enkephalins using the MOLS technique

Conformational studies of two linear enkephalin molecules, Met-enkephalin and Leu-enkephalin, have been carried using the mutually orthogonal Latin squares (MOLS) technique with the ECEPP/3 force field. This technique was developed recently in our laboratory to perform an unbiased search of the conformational space of peptides and to locate low energy conformations. The present study identified all the folds predicted by other studies, and in addition picked up other energetically favorable structures. The results suggest that the peptide backbone exists as a mixture of folded and unfolded forms (approximately 50% each). The study also provides information on the distribution of the low energy conformations that we have classified on the basis of structural motifs, backbone hydrogen-bonding patterns, and root mean square deviations in atomic positions.     

Starting‐structure dependence of nanosecond timescale intersubstate transitions and reproducibility of MD‐derived order parameters

Factors affecting the accuracy of molecular dynamics (MD) simulations are investigated by comparing generalized order parameters for backbone NH vectors of the B3 immunoglobulin‐binding domain of streptococcal protein G (GB3) derived from simulations with values obtained from NMR spin relaxation (Yao L, Grishaev A, Cornilescu G, Bax A, J Am Chem Soc 2010;132:4295‐4309.). Choices for many parameters of the simulations, such as buffer volume, water model, or salt concentration, have only minor influences on the resulting order parameters. In contrast, seemingly minor conformational differences in starting structures, such as orientations of sidechain hydroxyl groups, resulting from applying different protonation algorithms to the same structure, have major effects on backbone dynamics. Some, but not all, of these effects are mitigated by increased sampling in simulations. Most discrepancies between simulated and experimental results occur for residues located at the ends of secondary structures and involve large amplitude nanosecond timescale transitions between distinct conformational substates. These transitions result in autocorrelation functions for bond vector reorientation that do not converge when calculated over individual simulation blocks, typically of length similar to the overall rotational diffusion time. A test for convergence before averaging the order parameters from different blocks results in better agreement between order parameters calculated from different sets of simulations and with NMR‐derived order parameters. Thus, MD‐derived order parameters are more strongly affected by transitions between conformational substates than by fluctuations within individual substates themselves, while conformational differences in the starting structures affect the frequency and scale of such transitions. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Conformational transitions in eosinophil cationic protein: a molecular dynamics study in aqueous environment

Extensive molecular dynamics simulations have been performed on eosinophil cationic protein (ECP). The two structures found in the crystallographic dimer (ECPA and ECPB) have been independently simulated. A small difference in the pattern of the sidechain hydrogen bonds in the starting structure has resulted in interesting differences in the conformations accessed during the simulations. In one simulation (ECPB), a stable equilibrium conformation was obtained and in the other (ECPA), conformational transitions at the level of sidechain interactions were observed. The conformational transitions exhibit the involvement of the solvent (water) molecules with a pore-like construct in the equilibrium conformation and an opening for a large number of water molecules during the transition phase. The details of these transitions are examined in terms of intra-protein hydrogen bonds, protein-water networks and the residence times of water molecules on the polar atoms of the protein. These properties show some significant differences in the region between the N-terminal helix and the loop before the C-terminal strand as a function of different conformations accessed during the simulations. However, the stable hydrogen bonds, the protein-water networks, and the hydration patterns in most part of the protein including the active site are very much similar in both the simulations, indicating the fact that these are intrinsic properties of proteins.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

Disaccharide conformational flexibility. II. Molecular dynamics simulations of sucrose

Molecular dynamics simulations have been used to study the motions in vacuum of the disaccharide sucrose. Ensembles of trajectories were calculated for each of the five local minimum energy conformations identified in the adiabatic conformational energy mapping of this molecule. The model sucrose molecules were found to exhibit a variety of motions, although the global minimum energy conformation was found to be dynamically stable, and no transitions away from this structure were observed to occur spontaneously. In all but one of these vacuum trajectories, the intramolecular hydrogen bond between residues was maintained, in accord with recent nmr studies of this molecule in aqueous solution. Considerable flexibility of the furanoid ring was found in the trajectories. No "flips" to the opposite puckering for this ring were found in the simulations starting from the global minimum, although such a transition was observed for a trajectory initiated with one of the higher local minimum energy conformations. Overall, the observed structural fluctuations were consistent with the experimental picture of sucrose as a relatively rigid molecule.     

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.     

Conformational analyses of sandostatin analogs containing stereochemical changes in positions 6 or 8

We report the conformational analysis by 1H nmr in DMSO and computer simulations involving distance geometry and molecular dynamics simulations of analogs of the cyclic octapeptide D-Phe1-c[Cys2-Phe3-D-Trp4-Lys5-Thr6-Cys 7]-Thr8-ol (sandostatin, octreotide). The analogs D-Phe1-c[Cys2-Phe3-D-Trp4-Lys5-Xaa6-Cys 7]-Xbb8-NH2 (Xaa = allo-Thr, D-allo-Thr, D-beta-Hyv, beta-Hyv, D-Thr, and Xbb = Thr or Xaa = Thr and Xbb = allo-Thr, D-allo-Thr, beta-Hyv, D-Thr) contain stereochemical changes in the Thr residues in positions 6 and 8, which allow us to investigate the influence of the stereochemistry within these residues on conformation and binding affinity. The molecular dynamics simulations provide insight into the conformational flexibility of these analogs. The compounds with (S)-configuration at the C(alpha) of residue 6 adopt beta-sheet structures containing a type II' beta-turn with D-Trp in the i+1 position, and these conformations are "folded" about residues 6 and 3. The structures are very similar to those observed for sandostatin, and the disulfide bridge results in a close proximity of the H(alpha) protons of residues 7 and 2, which confirms earlier observations that a disulfide bridge is a good mimic for a cis peptide bond. The compounds with (R)-configuration at the C(alpha) of residue 6 adopt considerably different backbone conformations. The structures observed for these analogs contain either a beta-turn about residue Lys and Xaa6 or a gamma-turn about the Xaa6 residue. These compounds do not exhibit significant binding to the somatostatin receptors, while the compounds with (S) configuration in position 6 bind potently to the sst2, 3, and 5 receptors. The nmr spectra of analogs with (R) or (S) configuration at the C(alpha) of residue 8 are strikingly similar to each other. We have demonstrated that the chemical shifts of protons of residues 3, 4, 5, and 6, which are part of the type II' beta-turn, and especially the effect on the Lys gamma-protons are considerably different in active molecules as compared to inactive analogs. Since the presence of a type II' beta-turn is crucial for the binding to the receptors, the chemical shifts, the amide temperature coefficients of the Thr residue and the medium strength NOE between LysNH and ThrNH can be extremely useful as an initial screening tool to separate the active molecules from inactive analogs.     

Phosphorylation-induced conformational changes in a mitogen-activated protein kinase substrate. Implications for tyrosine hydroxylase activation

Mitogen-activated protein (MAP) kinase-mediated phosphorylation of specific residues in tyrosine hydroxylase leads to an increase in enzyme activity. However, the mechanism whereby phosphorylation affects enzyme turnover is not well understood. We used a combination of fluorescence resonance energy transfer (FRET) measurements and molecular dynamics simulations to explore the conformational free energy landscape of a 10-residue MAP kinase substrate found near the N terminus of the enzyme. This region is believed to be part of an autoregulatory sequence that overlies the active site of the enzyme. FRET was used to measure the effect of phosphorylation on the ensemble of peptide conformations, and molecular dynamics simulations generated free energy profiles for both the unphosphorylated and phosphorylated peptides. We demonstrate how FRET transfer efficiencies can be calculated from molecular dynamics simulations. For both the unphosphorylated and phosphorylated peptides, the calculated FRET efficiencies are in excellent agreement with the experimentally determined values. Moreover, the FRET measurements and molecular simulations suggest that phosphorylation causes the peptide backbone to change direction and fold into a compact structure relative to the unphosphorylated state. These results are consistent with a model of enzyme activation where phosphorylation of the MAP kinase substrate causes the N-terminal region to adopt a compact structure away from the active site. The methods we employ provide a general framework for analyzing the accessible conformational states of peptides and small molecules. Therefore, they are expected to be applicable to a variety of different systems.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

Exploration of the binding of proton pump inhibitors to human P450 2C9 based on docking and molecular dynamics simulation

Human P450 protein CYP2C9 is one of the major drug-metabolizing isomers, contributing to the oxidation of 16% of the drugs currently in clinical use. To examine the interaction mechanisms between CYP2C9 and proton pump inhibitions (PPIs), we used molecular docking and molecular dynamics (MD) simulation methods to investigate the conformations and interactions around the binding sites of PPIs/CYPP2C9. Results from molecular docking and MD simulations demonstrate that nine PPIs adopt two different conformations (extended and U-bend structures) at the binding sites and position themselves far above the heme of 2C9. The presence of PPIs changes the secondary structures and residue flexibilities of 2C9. Interestingly, at the binding sites of all PPI–CYP2C9 complexes except for Lan/CYP2C9, there are hydrogen-bonding networks made of PPIs, water molecules, and some residues of 2C9. Moreover, there are strong hydrophobic interactions at all binding sites for PPIs/2C9, which indicate that electrostatic interactions and hydrophobic interactions appear to be important for stabilizing the binding sites of most PPIs/2C9. However, in the case of Lan/2C9, the hydrophobic interactions are more important than the electrostatic interactions for stabilizing the binding site. In addition, an interesting conformational conversion from extended to U-bend structures was observed for pantoprazole, which is attributed to an H-bond interaction in the binding pocket, an internal π–π stacking interaction, and an internal electrostatic interaction of pantoprazole.     

Protein flexibility: multiple molecular dynamics simulations of insulin chain B

Multiple molecular dynamics simulations totaling more than 100 ns were performed on chain B of insulin in explicit solvent at 300 K and 400 K. Despite some individual variations, a comparison of the protein dynamics of each simulation showed similar trends and most structures were consistent with NMR experimental values, even at the elevated temperature. The importance of packing interactions in determining the conformational transitions of the protein was observed, sometimes resulting in conformations induced by localized hydrophobic interactions. The high temperature simulation generated a more diverse range of structures with similar elements of secondary structure and populated conformations to the simulations at room temperature. A broad sampling of the conformational space of insulin chain B illustrated a wide range of conformational states with many transitions at room temperature in addition to the conformational states observed experimentally. The T-state conformation associated with insulin activity was consistently present and a possible mechanism of behavior was suggested.     

The sequence TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase displays structural ambivalence and interconverts between alpha-helical and beta-hairpin conformations mediated by collapsed conformational states.

The peptide TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase adopts a helical conformation in the crystal structure and is a site for two hydrated helical segments, which are thought to be helical folding intermediates. Overlapping sequences of four to five residues from the peptide, sample both helical and strand conformations in known protein structures, which are dissimilar to glyceraldehyde-3-phosphate dehydrogenase suggesting that the peptide may have a structural ambivalence. Molecular dynamics simulations of the peptide sequence performed for a total simulation time of 1.2 micros, starting from the various initial conformations using GROMOS96 force field under NVT conditions, show that the peptide samples a large number of conformational forms with transitions from alpha-helix to beta-hairpin and vice versa. The peptide, therefore, displays a structural ambivalence. The mechanism from alpha-helix to beta-hairpin transition and vice versa reveals that the compact bends and turns conformational forms mediate such conformational transitions. These compact structures including helices and hairpins have similar hydrophobic radius of gyration (Rgh) values suggesting that similar hydrophobic interactions govern these conformational forms. The distribution of conformational energies is Gaussian with helix sampling lowest energy followed by the hairpins and coil. The lowest potential energy of the full helix may enable the peptide to take up helical conformation in the crystal structure of the glyceraldehyde-3-phosphate dehydrogenase, even though the peptide has a preference for hairpin too. The relevance of folding and unfolding events observed in our simulations to hydrophobic collapse model of protein folding are discussed.     

Intrinsic α‐helical and β‐sheet conformational preferences: A computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle (-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (N–C_α–C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the “bridge” region of-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between -ψ conformations and τ.     

Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin.

Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used to sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.     

Molecular dynamics studies of a DNA-binding protein: 1. A comparison of the trp repressor and trp aporepressor aqueous simulations.

The results of two 30-ps molecular dynamics simulations of the trp repressor and trp aporepressor proteins are presented in this paper. The simulations were obtained using the AMBER molecular mechanical force field and in both simulations a 6-A shell of TIP3P waters surrounded the proteins. The trp repressor protein is a DNA-binding regulatory protein and it utilizes a helix-turn-helix (D helix-turn-E helix) motif to interact with DNA. The trp aporepressor, lacking two molecules of the L-tryptophan corepressor, cannot bind specifically to DNA. Our simulations show that the N- and C-termini and the residues in and near the helix-turn-helix motifs are the most mobile regions of the proteins, in agreement with the X-ray crystallographic studies. Our simulations also find increased mobility of the residues in the turn-D helix-turn regions of the proteins. We find the average distance separating the DNA-binding motifs to be larger in the repressor as compared to the aporepressor. In addition to examining the protein residue fluctuations and deviations with respect to X-ray structures, we have also focused on backbone dihedral angles and corepressor hydrogen-bonding patterns in this paper.     

Reversible scaling of dihedral angle barriers during molecular dynamics to improve structure prediction of cyclic peptides.

Peptide cyclization or chemical cross-linking has frequently been used to restrict the conformational freedom of a peptide, for example, to enhance its capacity for selective binding to a target receptor molecule. Structure prediction of cyclic peptides is important to evaluate possible conformations prior to synthesis. Because of the conformational constraints imposed by cyclization low energy conformations of cyclic peptides can be separated by large energy barriers. In order to improve the conformational search properties of molecular dynamics (MD) simulations a potential scaling method has been designed. The approach consists of several consecutive MD simulations with a specific lowering of dihedral energy barriers and reduced nonbonded interactions between atoms separated by three atoms followed by gradually scaling the potential until the original barriers are reached. Application to four cyclic penta- and hexa-peptide test cases and a protein loop of known structure indicates that the potential scaling method is more efficient and faster in locating low energy conformations than standard MD simulations. Combined with a generalized Born implicit solvation model the low energy cyclic peptide conformations and the loop structure are in good agreement with experiment. Applications in the presence of explicit water molecules during the simulations showed also improved convergence to structures close to experiment compared with regular MD.     

High resolution approach to the native state ensemble kinetics and thermodynamics

Many biologically interesting functions such as allosteric switching or protein-ligand binding are determined by the kinetics and mechanisms of transitions between various conformational substates of the native basin of globular proteins. To advance our understanding of these processes, we constructed a two-dimensional free energy surface (FES) of the native basin of a small globular protein, Trp-cage. The corresponding order parameters were defined using two native substructures of Trp-cage. These calculations were based on extensive explicit water all-atom molecular dynamics simulations. Using the obtained two-dimensional FES, we studied the transition kinetics between two Trp-cage conformations, finding that switching process shows a borderline behavior between diffusive and weakly-activated dynamics. The transition is well-characterized kinetically as a biexponential process. We also introduced a new one-dimensional reaction coordinate for the conformational transition, finding reasonable qualitative agreement with the two-dimensional kinetics results. We investigated the distribution of all the 38 native nuclear magnetic resonance structures on the obtained FES, analyzing interactions that stabilize specific low-energy conformations. Finally, we constructed a FES for the same system but with simple dielectric model of water instead of explicit water, finding that the results were surprisingly similar in a small region centered on the native conformations. The dissimilarities between the explicit and implicit model on the larger-scale point to the important role of water in mediating interactions between amino acid residues.     

Application of biasing‐potential replica‐exchange simulations for loop modeling and refinement of proteins in explicit solvent

Comparative or homology modeling of a target protein based on sequence similarity to a protein with known structure is widely used to provide structural models of proteins. Depending on the target‐template similarity these model structures may contain regions of limited structural accuracy. In principle, molecular dynamics (MD) simulations can be used to refine protein model structures and also to model loop regions that connect structurally conserved regions but it is limited by the currently accessible simulation time scales. A recently developed biasing potential replica exchange (BP‐REMD) method was used to refine loops and complete decoy protein structures at atomic resolution including explicit solvent. In standard REMD simulations several replicas of a system are run in parallel at different temperatures allowing exchanges at preset time intervals. In a BP‐REMD simulation replicas are controlled by various levels of a biasing potential to reduce the energy barriers associated with peptide backbone dihedral transitions. The method requires much fewer replicas for efficient sampling compared with T‐REMD. Application of the approach to several protein loops indicated improved conformational sampling of backbone dihedral angle of loop residues compared to conventional MD simulations. BP‐REMD refinement simulations on several test cases starting from decoy structures deviating significantly from the native structure resulted in final structures in much closer agreement with experiment compared to conventional MD simulations. Proteins 2010. © 2010 Wiley‐Liss, Inc.