Next Article in Journal
Synthesis and Magnetic Properties of Stable Radical Derivatives Carrying a Phenylacetylene Unit
Next Article in Special Issue
An Investigation of Atomic Structures Derived from X-ray Crystallography and Cryo-Electron Microscopy Using Distal Blocks of Side-Chains
Previous Article in Journal
Nectandra grandiflora By-Products Obtained by Alternative Extraction Methods as a Source of Phytochemicals with Antioxidant and Antifungal Properties
Previous Article in Special Issue
Exploring Protein Cavities through Rigidity Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Segmenting Proteins into Tripeptides to Enhance Conformational Sampling with Monte Carlo Methods

LAAS-CNRS, Université de Toulouse, CNRS, 31400 Toulouse, France
*
Author to whom correspondence should be addressed.
Current address: Computer Science Department, Princeton University, Princeton, NJ 08540, USA.
Molecules 2018, 23(2), 373; https://doi.org/10.3390/molecules23020373
Submission received: 13 December 2017 / Revised: 22 January 2018 / Accepted: 1 February 2018 / Published: 9 February 2018

Abstract

:
This paper presents an approach to enhance conformational sampling of proteins employing stochastic algorithms such as Monte Carlo (MC) methods. The approach is based on a mechanistic representation of proteins and on the application of methods originating from robotics. We outline the general ideas of our approach and detail how it can be applied to construct several MC move classes, all operating on a shared representation of the molecule and using a single mathematical solver. We showcase these sampling techniques on several types of proteins. Results show that combining several move classes, which can be easily implemented thanks to the proposed approach, significantly improves sampling efficiency.

Graphical Abstract

1. Introduction

The study of proteins is essential to understanding living organisms, since they are responsible for most processes inside cells. Proteins are also of interest for other domains. For instance, proteins are pharmaceutical targets of drugs, their catalytic properties are exploited in biotechnology, and they are used as components of nano-devices in the rising field of bionanotechnology. A better understanding of the relationship between structural and dynamic features of proteins and their function is fundamental in all of these domains.
Protein investigations are very challenging because of the large size and flexibility of these biological molecules. X-ray crystallography [1] is the most widely-used experimental method to obtain atomic-scale data. However, it only provides a static picture. In this respect, nuclear magnetic resonance (NMR) spectroscopy is a more complete technique, since it is able to provide relevant structural and dynamic data of proteins in solution [2,3,4]. These data are nevertheless limited and averaged. Therefore, computational methods, which have been proposed over the last decades, are used as a complement to experimental methods to model proteins and to simulate their behavior. They are mostly based on molecular dynamics simulations or Monte Carlo methods [5,6], although alternative methods have been proposed in recent years [7,8], some of which are originating from robotics [9,10,11].
This paper presents an approach to enhance conformational exploration methods (note that a preliminary version of this work was presented at the ASME IDETC/CIE 2012 conference [12]). The approach is based on a mechanistic view of proteins [13]. The idea is to cut the protein into small fragments of three amino acid residues, which we refer to as tripeptides. Each fragment can be represented as a kinematic chain, similar to a robotic manipulator. Such a representation enables performing local deformations of the protein backbone model that preserve bond geometry, using closed-form inverse kinematics (IK) solvers [14]. This work focuses on the application of this approach for devising Monte Carlo move classes. Nevertheless, the tripeptide-based protein representation introduced below can be exploited within other types of methods, as further discussed in the conclusion.
The basic Monte Carlo (MC) method [6,15] explores the conformational space through a random walk. At each iteration, the protein conformation is randomly perturbed, and the trial move is accepted or rejected with a probability that depends on the potential energies of the old and the new states. The main difficulty involving the application of the MC method to proteins lies in devising suitable trial move classes for complex chain-like molecules. An effective move class would yield a good acceptance rate (therefore avoiding futile, expensive energy evaluations), while enabling the exploration of large regions of the conformational space. Along more than three decades, significant efforts have been devoted to the development of trial move classes to enhance the efficiency of MC methods applied to proteins and other chain-like polymers [16,17,18,19,20,21,22]. A comprehensive review was published by Vitalis and Pappu a few years ago [23]. Building on these important contributions, the approach presented in this paper enables devising and combining different types of move classes that can be easily implemented using a unique molecular representation and a single solver.
The paper presents the general aspects of the mechanistic protein representation using the tripeptide-based decomposition, and explains how to implement several move classes based on this representation. Note however that the application of this approach to structure prediction or other related problems is out of the scope of this paper. Our aim is to show the interest of the proposed tripeptide-based model to enhance protein backbone sampling.

2. Results and Discussion

2.1. Implemented Move Classes and Parameter Settings

The methods presented in this paper have been implemented in a software package called MoMA (for Molecular Motion Algorithms), which includes modeling tools and algorithms to sample conformations and transition paths of biomolecules. The open-source code of MoMA is not available yet. Nevertheless, binaries can be provided upon request.
We have implemented four move classes. Three of them, producing fixed-end moves, are based on the proposed approach. The other one is simply based on the variation of a single bond torsion at each iteration. More precisely:
-
The simplest class of trial moves, largely applied to sample the conformation of chain-like molecules, consists of perturbing a randomly selected bond torsion and then propagating the motion toward the end of the chain. Such moves, usually called pivot moves, are named here OneTorsion moves. They are illustrated in Figure 1a.
-
The second move class is named ConRot, since it is inspired by the concerted rotations proposed by Dodd et al. [17]. It has been implemented using the proposed tripeptide-based model as follows: an amino-acid residue is randomly selected and one of its bond torsions ( ϕ or ψ ) is randomly perturbed; the backbone conformation of the next three residues (the next tripeptide) is modified by inverse kinematics in order to maintain fixed ends (see Section 3.2 for details). The move class is illustrated in Figure 1b.
-
The third move class, called OneParticle moves, corresponds to the simplest move class involving tripeptide-based particle perturbations, as described in Section 3.2, and illustrated in Figure 1c.
-
The last move class, called Hinge moves, corresponds to the rigid-body block moves described in Section 3.2, and illustrated in Figure 1d. The number of consecutive particles affected by the move is randomly sampled at each iteration between 3 and 10 (i.e., moves involve between 9 and 30 residues).
These four move classes have been applied within a basic MC method, using the Metropolis criterion [15] to accept or to reject trial moves. At each iteration, the algorithm randomly chooses between performing either a backbone move or a side-chain move. A side-chain move consists of randomly selecting a side-chain and perturbing all of its dihedral angles χ i . We have tested the four backbone move classes individually as well as a simple combination of all of them. The Mixed move class selects one of the four backbone move classes with equal probability at each iteration.
At each iteration of the MC method, the conformational parameters involved in the applied move class (bond torsions or oriented particle poses) are perturbed by adding a random value to their original value, sampled in the interval [ δ , δ ] . Thus, the parameter δ defines the maximum perturbation step-size. For a meaningful comparative analysis, the values used for the perturbation of bond torsions ( δ b ), particle translations ( δ p t ) and particle rotations ( δ p r ) in each move class have to be chosen in such a way that they will produce average atom displacements of similar length. In general, this also implies that MC acceptance rates will be similar for all move classes. The values used in this work for the two test systems introduced below are presented in Table 1.
For energy evaluation, we use an in-house implementation of the AMBER parm96 force-field [24] with an implicit representation of the solvent using the Generalized Born (GB) approximation. A geometric filter is applied before energy evaluation with the aim of improving computational efficiency (note that our implementation of the energy force-field is not optimal in terms of computing time.). After applying each trial move, the model is checked for atom overlaps. A trial move is rejected if the distance between two non-bonded atoms is less than 70% of the van der Waals equilibrium distance [25]. In addition, for ConRot, OneParticle and Hinge move classes, a trial move is rejected if the IK solver fails to find a solution. All the tests have been performed at a temperature of 300 K.

2.2. Test Systems

We have chosen two different types of proteins to evaluate the performance of the move classes. The first one is a SH3 domain [26]. More precisely, we have used the SH3 domain of obscurin, represented in Figure 2a. This is a small globular protein composed of 68 amino-acid residues. It presents a relatively rigid beta-barrel-like core and two flexible loops. Its crystal structure is available in the Protein Data Bank (PDB ID: 1V1C). The second test system is an intrinsically disordered protein (IDP). More precisely, we have considered a disordered region of the Sic1 protein [27], containing 77 residues. A model of this protein region is shown in Figure 2b. The model was generated using the Flexible-Meccano method [28] for sampling a statistically probable backbone conformation, and SCWRL4 [29] for placing the side-chains. Note that, although the application of MC methods to sample the extremely vast conformational space of IDPs may appear unrealistic, it can be interesting for more specific problems such as the analysis of the energy landscape around a biologically functional conformation of the protein (e.g., a conformation of the IDP while interacting with another molecule).
The two structures were locally energy-minimized before running the tests. Starting from these initial conformations, the MC method was iterated until the generation of 2 × 10 6 accepted samples for each molecule, using the four move classes individually or combined together. Each test was repeated three times. All tests using the same settings provided similar results.
The results presented in this paper are not aimed at providing new insights into these biological systems, but to serve as a proof of concept and to show the interest of the proposed approach for different types of proteins.

2.3. Computational Performance

Table 2 contains results on the computational performance of the methods, averaged over the 3 runs. It shows the MC acceptance rate for backbone moves, the total number of iterations (considering backbone and side-chain moves) and the overall CPU time for each system and move class. Tests were run on a single Intel® Xeon® E5-1650 processor at 3.2 GHz. These results show that, for a similar MC acceptance rate, generating a given number of samples ( 2 × 10 6 in this case) requires more computational resources using the OneTorsion move class compared to the other move classes, which apply fixed-end motions. The reason is that, although the OneTorsion move class does not require solving inverse kinematics, a significant amount of computing time is needed for propagating atom motions along the chain by forward kinematics and to recompute interaction energies between atom pairs. In other words, the computing time needed by the ConRot, OneParticle and Hinge move classes to solve inverse kinematics is largely compensated since only the positions of a small number of atoms need to be updated after each move, which also decreases the cost of energy recalculation. Nevertheless, computing time is probably not the most important indicator of the performance of the different move classes. Indeed, the numbers in Table 2 have to be analyzed together with other data related to the quality of the sampling strategy in terms of energy and distance distributions, as discussed below.

2.4. Distribution of Sampled States

Figure 3 shows plots aimed at comparing the performance of the different move classes in terms of conformational space coverage for the SH3 domain and the Sic1 protein. They represent the projection of the sampled states on two dimensions: the distance with respect to the initial structure and the potential energy. The distance is measured as the root mean square deviation (RMSD) of the bond torsion angles from the initial conformation. For clarity reasons, only one in every one hundred sampled conformations has been plotted (i.e., 20,000 samples from each run). Each plot includes the results of the three runs.
Several conclusions can be extracted from the analysis of these plots: (1) the Mixed move class shows the best performance. The combination of move-classes provides samples that have lower energies, and conformational coverage (distances from a reference conformation) is comparable or better than for the individual move classes; (2) move classes involving a larger number of atoms, i.e., OneTorsion and Hinge, show better performance on the disordered Sic1 protein, which presents higher conformational variability than the globular SH3 domain. For the Sic1 protein, these move classes reach lower-energy regions compared to ConRot and OneParticle; (3) as it could be expected, the performance of ConRot and OneParticle moves is very similar for both test systems; and (4) the Hinge move class shows slightly poorer exploration capabilities compared to the other move classes. The maximum distances of the samples to the reference conformation are smaller when using this move class alone. The difference is more significant for the SH3 domain.

2.5. Exploration Efficiency Analysis

The distance vs. energy plots permit comparing the different move classes in terms of coverage, but they do not allow us to quantitatively evaluate their exploration efficiency. Furthermore, these plots illustrate the performance at the end of a long MC run, but do not show the short-term exploration capabilities of the different move classes, which are of interest for applications where MC methods are run for a short time, such as structural relaxation. This section presents additional tests allowing such an analysis.

2.5.1. Time Dependent RMSD Function

The time dependent RMSD function, rmsd ( τ ) , is aimed at showing the rapidity of the exploration process. It represents the average distance between conformations separated by τ MC steps:
rmsd ( τ ) = RMSD ( q k , q k + τ ) k ,
where q k is the configuration of the system after the kth trial in the MC simulation, and X ( k ) k denotes the average of the observable X over all k.
Figure 4 shows the time dependent RMSD function for each move class and for the two molecules, using the results from the aforementioned MC simulations ( 2 × 10 6 accepted samples), for τ values up to 30,000 MC steps. It provides interesting additional information with respect to the analysis in the previous section, and allows for highlighting more significant differences between the move classes: (1) it confirms that the Mixed move class performs better than the other individual move classes; (2) although the OneParticle and ConRot move classes perform very similarly, it appears that the OneParticle move class explores a bit faster than ConRot. In the case of the SH3 domain, OneParticle can even compete with the Mixed move class on small time scales; and (3) the Hinge move class performs almost two times worse than the Mixed move class for both proteins, while the OneTorsion move class performs almost four times worse. This last result, together with results in the previous section, show that, although OneTorsion moves may provide good performance in terms of coverage (particularly for the disordered Sic1 protein), convergence can be slow.

2.5.2. Autocorrelation Time

Autocorrelation is a statistical tool that computes the correlation of a time series with a lagged version of itself. In the context of a MC simulation, it is used to characterize how new accepted states gradually become independent from the previous ones. Formally, the autocorrelation of an observable O with a lag τ is defined by:
acf ( τ ) = E [ ( O k O ) ( O k + τ O ) ] Var ( O ) ,
where E [ X ] denotes the expected value of X, Var ( O ) denotes the variance of O. O k and O k + τ represent two series of observations of O during the MC run with a lag of τ steps. This formula only holds true when the simulation is stationary. In the context of MC, this requirement implies that the simulation is long enough for the autocorrelation value to converge. Since this would require very long computing time for large systems such as the SH3 domain or the Sic1 protein, we have used a smaller system for this test. Inspired by related work [22,30], we have used a 14-alanine molecule with a structural constraint: starting from a low energy stable state, corresponding to an α -helical conformation, the two end residues were blocked in position and orientation so that the system could only fluctuate around this state. This constraint, associated with the small size of the system, restricts the applicable move classes to ConRot, OneParticle, and a Mixed move class involving both of them.
Starting from the helical conformation, the MC method was iterated until 10 8 samples were accepted for the two move classes individually or combined together (no side-chain moves were performed). Figure 5 shows the average autocorrelation function over the 20 central dihedral angles for each move class. We can observe that the Mixed move class yields the fastest independence time, whereas the autocorrelation of ConRot seems to stagnate when approaching zero.
For a more precise evaluation, we computed the characteristic autocorrelation time, also called the integral autocorrelation time. It is defined by:
τ i n t = 1 2 + τ = 1 acf ( τ ) .
In practice, to reduce computational cost and noisy values, the summation is usually truncated. We have summed for all τ 5 τ i n t . The autocorrelation times of the 20 central dihedral angles are plotted in Figure 6, and the statistics for each move class are summarized in Table 3. These results allow us to quantify the improvement of autocorrelation time between the different move classes. The OneParticle move class has an autocorrelation time about 1.8 times lower on average than the ConRot move class, while the Mixed move class consistently outperforms the OneParticle move class by a factor of 2. These results confirm that the Mixed move class is much more efficient than the individual move classes.

2.6. Additional Results for Ubiquitin

We performed additional tests on human ubiquitin, a well studied system with experimental and computational methods [31]. The structure of this protein (PDB ID: 1UBQ) is shown in Figure 7. Our main goal was to compare, qualitatively, the results obtained with the methods proposed in this paper to those presented by Bottaro et al. [22]. The authors of this related work propose an extension of the concerted rotation move class, called CRISP (Concerted Rotations Involving Self-consistent Proposals), and use ubiquitin to demonstrate the ability of their enhanced MC method to reproduce structural variations around the native conformation, with a performance comparable to molecular dynamics simulations. Note that the CRISP move class is compatible with the approach presented in this paper, and could be implemented using the proposed framework.
For the comparative analysis, we applied the ConRot and OneParticle move classes, as well as a Mixed move class involving both of them. Three MC simulations of 15 × 10 6 steps each were performed for each move class, starting from a relaxed state of the protein. The plots presented below display average values of the three runs. Each MC simulation randomly alternates between the compared move class (probability of 20%), a side-chain perturbation (probability of 75%), or a OneTorsion move (probability of 5%). We used the same step sizes than for the SH3 domain (see Table 1).
Figure 8 shows the root mean squared fluctuations (RMSF) of the C α atoms. The RMSF is calculated as the RMSD from the mean (equilibrium) position. As for the previous experiments, the figure shows that the ConRot and OneParticle move classes perform similarly, and that the Mixed move class explores more widely than the individual move classes. In this case, ConRot slightly outperforms OneParticle. This can be explained by the compactness of the ubiquitin structure with respect to the SH3 domain and the Sic1 protein, which favours local moves involving few atoms. Interestingly, although we are using a different force-field and a less sophisticated sampling method for the side-chains, the C α RMSF profiles are very similar to those presented in the referred publication by Bottaro et al. [22].
The efficiency of the different move classes can be roughly evaluated and compared by looking at the evolution of the RMSF with the number of iterations. The RMSF of a particular atom is expected to grow with the simulation time, and eventually converge to a maximum value when the energy basing around the equilibrium state has been completely explored. Figure 9 shows the evolution of the cumulative RMSF (i.e., the sum of the RMSF values of all C α atoms) with the length of the simulation. The plot confirms a faster exploration of the Mixed move class, which allows reaching a cumulative RMSF value around 50 Å after 15 × 10 6 MC steps, whereas the individual ConRot and OneParticle move classes only attain 40 Å. These values are also very similar to the ones obtained by Bottaro et al. [22]. In fact, a direct comparison of Figure 9 in this paper and Figure 4 in the referred publication suggests a faster convergence of our Mixed move class with respect to the CRISP move class, which requires around 10 times more iterations to reach 50 Å ( 5 × 10 7 MC steps instead of 5 × 10 6 MC steps). However, such a direct comparison should be put into perspective because the models and the implementations are different.

3. Materials and Methods

3.1. Protein Model

3.1.1. Mechanistic Model

The conformation (i.e., spatial arrangement) of a protein can be defined by the Cartesian coordinates of all its constituent atoms, or by a vector of internal coordinates that represent the relative position of bonded atoms. These internal coordinates correspond to the bond lengths, bond angles and bond torsions. Since the bond lengths and bond angles vary very slightly at room temperature, they are often considered to be constant parameters in molecular simulations [32]. Under such assumptions, the bond torsions are the only degrees of freedom of the molecule. An additional simplification of molecular models is to consider that double bonds, such as peptide bonds in proteins, are rigid connections (i.e., the dihedral angle ω associated with the peptide bond torsion is constant). In summary, the variable parameters that define the conformation of a protein backbone are the pairs of dihedral angles, ϕ and ψ , of all its constituent amino-acid residues. The conformation of the side-chains is determined by a variable number of dihedral angles χ i for each residue.
In this work, we consider the aforementioned rigid geometry assumption. Nevertheless, small variations of bond lengths, bond angles and peptide bond torsions, which can be important from a structural point of view [33], can be considered within our approach, as discussed below.
Using the internal coordinate representation described above, proteins can be modeled as articulated mechanisms. The bodies of the mechanism correspond to rigidly-bonded atom groups, and the joints are the bond torsions. The kinematic chains corresponding to the protein backbone and side-chains can then be modeled using standard conventions usually applied in robotics. In this work, we have used the modified Denavit–Hartenberg (mDH) convention [34]. Following this convention, a Cartesian coordinate system F i is attached to each rigid atom group. The relative location of consecutive frames in a kinematic chain can be then defined by a homogeneous transformation matrix of the form:
i 1 T i = cos θ i sin θ i 0 a i 1 sin θ i cos α i 1 cos θ i cos α i 1 sin α i 1 d i sin α i 1 sin θ i sin α i 1 cos θ i sin α i 1 cos α i 1 d i cos α i 1 0 0 0 1 .
The elements of i 1 T i depend on the bond geometry. Assuming constant bond lengths and bond angles, the parameters d i and α i 1 (which can be directly computed from them) are constant. The bond torsion angle θ i is the only variable parameter. Figure 10 illustrates the method to assign the frames and to obtain the mDH parameters in a general case, and when peptide bond torsion angles ( ω ) are considered to have an arbitrary constant value. Note that, in the general case, when all the bond torsion angles are variable, the parameters a i = 0 .

3.1.2. Decomposition into Tripeptides

The main idea proposed in this paper is to segment the protein chain into fragments of three amino acid residues, which we refer to as tripeptides. The reason for choosing such a subdivision is that each tripeptide backbone involves six degrees of freedom (three pairs of ϕ , ψ angles), corresponding to the shortest fragment with full instantaneous mobility of the end-frame relatively to the base-frame (i.e., the relative position and orientation of the two frames can change in any direction). Figure 11 illustrates this idea. Figure 11a shows a protein model with a ribbon representing the backbone embedded in the model of the protein surface. Figure 11b represents the protein backbone trace with the frames corresponding to the connections between consecutive tripeptides. Figure 11c,d represent the chemical and the mechanistic models of the backbone of a single tripeptide, respectively. As depicted in the figure, the tripeptide backbone can be seen as a robotic manipulator with six revolute joints. The base of the manipulator corresponds to the first body of the tripeptide backbone (i.e., the first rigid atom group in the backbone of the first amino-acid residue). Since tripeptides are linked through rigid peptide bonds, the location of the end effector of tripeptide i can be determined from the base-frame of tripeptide i + 1 by a constant transformation. Given the location of the base-frame and the end-frame, the conformation of a tripeptide backbone can be determined by efficient analytical inverse kinematics (IK). The IK problem formulation and the particular solver applied in this work are described in the next section. Consequently, the conformation of the whole protein backbone can be determined from the pose of a single reference frame for each tripeptide, the one attached to the first body of each tripeptide backbone. In the following, we will refer to these reference frames as (oriented) particles. The last affirmation is true for all the protein backbone except two short fragments at the N-terminal and C-terminal ends of the chain. Since the choice of the first residue for the decomposition into tripeptides is arbitrary (and may change during the conformational exploration process), the polypeptide chain model involves two terminal fragments, containing up to three residues, which require a particular treatment. The conformations of these terminal fragments are directly defined by their internal bond torsions.

3.2. Devising Move Classes

This section presents a unified approach for devising different move classes that can be employed in the context of a Monte Carlo framework. These classes all utilize the tripeptide-based representations described above. The principle consists of perturbing the pose (position and orientation) of a set of particles, and then adapting the conformation of several tripeptides in order to keep the integrity of the molecular chain while maintaining the local geometry of the bonds (i.e., constant bond lengths and bond angles). Several strategies can be considered for perturbing the pose of particles. The number of particles selected for perturbation and the correlation/uncorrelation between the motion direction of several particles leads to different move classes, more or less local, and more or less collective.
Next, we explain how to implement simple general-purpose move classes. Here, we only deal with unbiased random moves. However, our approach is also suitable for devising biased moves. In such a case, the selection of the particles to be perturbed and the motion directions would be determined depending on the specific context. For instance, moves could be devised to deform proteins while simulating the interaction with other molecules, or for applications such as all-atom model fitting into lower-resolution electron density maps.
The move classes presented below, as well as other ones such as pivot moves applied to a single bond torsion, can be combined within a higher-level sampling protocol that selects a move class at each iteration. In addition, the constraints on bond lengths, bond angles, and peptide bond torsions imposed by the tripeptide-based model can be relaxed. This can be done by performing separate MC moves on these parameters, or by slightly perturbing the geometry of some or all the bonds in a tripeptide before applying the IK solver. Note that a similar approach to include variable bond geometry has been applied in related work [35].
Side-chain moves, which are usually performed separately from backbone moves, are not treated in this section. Side-chain conformations can be sampled by simple perturbations of the bond torsion angles χ i or following more sophisticated approaches [19,36]. Moreover, different strategies can be adopted to combine backbone and side-chain trial moves in a suitable manner [37].

3.2.1. Perturbing Particles

Different move classes can be easily implemented from the proposed tripeptide-based model. As an example, three move classes involving perturbations of one or several particles as explained next. As explained in Section 2.1, we have implemented several of them as well as other frequently-used moves for an empirical evaluation.
One-particle moves: Local moves can be implemented by perturbing the pose of a single particle, as depicted in Figure 1c. Such perturbation implies that the two tripeptides linked through this particle (i.e., with end-frame or base-frame defined from it) change their backbone conformation. In other words, 12 consecutive bond torsions are modified, while the rest of the protein conformation remains unchanged. This type of moves has a similar effect to other local, fixed-end move classes [17,18,38] proposed since the seminal work of Gō and Scheraga [16].
Flexible fragment moves: The previous move class can be extended to larger fragments by applying perturbations to a set of n consecutive particles. If the particles are perturbed independently (i.e., in different random directions), the backbone of n + 1 tripeptides is affected by the move. This move class has a similar effect to moves based on the cyclic coordinate descent method (CCD) [39]. Such moves consist of breaking the chain at an arbitrary point, performing a random perturbation of several bond torsions at one of the sides, and applying CCD to close the chain again (but with a different, perturbed conformation). An interesting advantage of the proposed move class with respect to CCD-based moves is that the deformation of sub-fragments can be easily modulated (e.g., larger perturbations for the middle particles in the fragment and smaller perturbations near the ends).
Rigid-body block moves: A simple variant of the previous move class may produce a very different effect, as illustrated in Figure 1d. In this case, n consecutive particles are also perturbed, but the perturbations are correlated in such a way that the particles do not move with respect to each other. Indeed, the perturbation is applied to a fictitious rigid body formed by the set of n particles. In principle, a random translation and rotation around an arbitrary axis could be tried. Nevertheless, it may be more interesting to apply moves that simulate hinge motions. Note that only two “hinge” tripeptides (i.e., the tripeptides preceding and following the selected particle sequence) change their conformation. Hinge-like moves, called Closed Rigid-body Rotation Under Bond-angle Restraints (CRRUBAR) moves [20], have been shown to be particularly efficient for sampling conformations of proteins. Although the proposed method involves more complex algebraic operations than CRRUBAR, it presents the advantage that bond angles do not need to be distorted.

3.2.2. Solving Inverse Kinematics for a Tripeptide

Once the location of the particles is set, obtaining the conformation of each tripeptide requires solving an inverse kinematics (IK) problem for the chain corresponding to its backbone. As explained above, the model of a tripeptide backbone is similar to a six-revolute (6R) serial manipulator with general geometry. Its conformation is defined by the vector of six angles { θ 1 , θ 2 , θ 3 , θ 4 , θ 5 , θ 6 } , corresponding to the three pairs of ϕ , ψ bond torsion angles. The IK problem is to find values for this vector of angles such that they solve the following matrix equation: 0 T 1 1 T 2 2 T 3 3 T 4 4 T 5 5 T 6 = T r e l , where T r e l is the homogeneous transformation matrix representing the relative pose of two consecutive particles, and each i 1 T i is the transformation matrix between consecutive atom groups in the chain, which is a function of θ i .
The method applied in this work for solving the IK problem for a general 6R serial kinematic chain has been adapted from the solver proposed by Renaud [40,41]. This solver is based on algebraic elimination theory, and develops an ad hoc resultant formulation inspired by the work of Lee and Liang [42,43]. Starting from the matrix equation representing the IK problem, the elimination procedure leads to an 8-by-8 quadratic polynomial matrix in one variable. The problem can then be treated as a generalized eigenvalue problem, as was previously proposed by Manocha and Canny [44], for which efficient and robust solutions are available [45]. Our implementation applies the Schur factorization from the linear algebra package LAPACK [46]. Further details on the applied IK solver are provided in a technical report by Renaud [41].
This solver has been successfully applied in previous works on protein and polymer modeling [47,48]. The advantage of this semi-analytical method with respect to numerical (optimization-based) methods, such as CCD, is that it provides the exact solution in a single iteration, not suffering from slow convergence issues. The solver is very computationally efficient, requiring about 0.2 milliseconds on a single processor. Note, however, that our approach is not dependent on this solver, so that other IK methods [17,44,49,50] could be applied.
In general, the IK problem for a 6R serial kinematic chain has a finite number of solutions (up to 16 in the most general case). All the solutions correspond to geometrically valid conformations of the tripeptide backbone with fixed ends defined by the pose of the particles. Depending on the type of application, several strategies can be adopted to select one of the solutions. The simplest strategy within a MC method is to select one at random. However, if detailed balance needs to be satisfied for a correct sampling of equilibrium fluctuations in the canonical ensemble, some works [17,19,30,49] recommend to take one of the solutions with a probability that depends on the Boltzmann factor and another term, called the Jacobian, which attempts to correct for the non-uniformity in the distribution of the torsion angles introduced by the fixed-end moves. Alternatively, when the moves are performed for finding energy minima, all the conformations can be evaluated in order to keep the best one (in terms of the Boltzmann factor, for instance). Otherwise, when the goal is to simulate continuous motions, the closest conformation to the one prior to the perturbation is selected in order to minimize the jumps in conformational space (if none of the solutions remains within a distance threshold that depends on the perturbation step-size, the local move is rejected). Note that this strategy obeys detained balance if a “reverse proximity criterion” is satisfied [51].

4. Conclusions

We have presented a unified approach to devise efficient Monte Carlo move classes for chain-like molecules. Based on a subdivision of the protein into tripeptides and on the application of efficient 6R inverse kinematics, many of the move classes that have been proposed in the last decades to enhance protein backbone sampling can be easily implemented. The approach is general, being applicable to proteins on any size, shape and topology. First results for proteins from different structural classes (globular and disordered) have been presented as a proof of concept. The overall conclusion of these results is that mixing move classes provides better results than using a single move class, as also suggested by related work on MC methods. Combining different move classes is a straightforward task following the proposed approach.
In this paper, we have only presented results of simple, unbiased MC sampling. However, the approach could be implemented within more advanced methods for conformational sampling [52] or global optimization [53]. Future work could also involve the implementation of a more sophisticated sampling strategy performing different types of trial moves selected among a variety of parameterized move classes.
We also envisage investigating the application of our approach to enhance other types of conformational exploration methods. More precisely, we plan to introduce the tripeptide-based representation within path-planning-based methods to compute motions associated with protein–ligand interactions [54]. Other research directions would concern the development of methods to refine predictions of protein–protein docking methods or to fit structures into electron density maps from cryo-electron microscopy (cryo-EM) using flexible protein models. For these applications, the proposed tripeptide-based representation is particularly interesting since it enables a very simple implementation of biased moves that deform regions of the protein with respect to interactions with other molecules or to locally improve structure fitting. Computational protein design is another potential application area for the proposed approach. It has been shown that introducing local backbone perturbation may significantly improve the results provided by protein design methods [55]. Using our approach, different classes of local moves around mutated residues can be easily implemented inside stochastic algorithms for protein design.

Acknowledgments

The authors would like to thank Pau Bernadó for providing us with the model of the intrinsically disordered protein Sic1. This work has been partially supported by the French National Research Agency (ANR) under project ProtiCAD (project number ANR-12-MONU-0015).

Author Contributions

T.S. and J.C. conceived the methods; L.D., I.A.-B., M.V. and J.C. implemented the methods; L.D. and J.C. conceived and designed the experiments; L.D. performed the experiments and analyzed the data; L.D. and J.C. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Woolfson, M.M. An Introduction to X-ray Crystallography; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  2. Kay, L.E. NMR studies of protein structure and dynamics. J. Magn. Reson. 2005, 173, 193–207. [Google Scholar] [CrossRef] [PubMed]
  3. Lange, O.F.; Lakomek, N.A.; Farès, C.; Schröder, G.F.; Walter, K.F.A.; Becker, S.; Meiler, J.; Grubmüller, H.; Griesinger, C.; de Groot, B.L. Recognition Dynamics Up to Microseconds Revealed from an RDC-Derived Ubiquitin Ensemble in Solution. Science 2008, 320, 1471–1475. [Google Scholar] [CrossRef] [PubMed]
  4. Fenwick, R.B.; van den Bedem, H.; Fraser, J.S.; Wright, P.E. Integrated description of protein dynamics from room-temperature X-ray crystallography and NMR. Proc. Natl. Acad. Sci. USA 2014, 111, E445–E454. [Google Scholar] [CrossRef] [PubMed]
  5. Leach, A.R. Molecular Modelling: Principles and Applications; Pearson Education: Harlow, UK, 2001. [Google Scholar]
  6. Frenkel, D.; Smit, B. Understanding Molecular Simulations: From Algorihtms to Applications; Academic Press: London, UK, 2002. [Google Scholar]
  7. Wales, D. Energy Landscapes: Applications to Clusters, Biomolecules and Glasses; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  8. Donald, B.R. Algorithms in Structural Molecular Biology; The MIT Press: Cambridge, MA, USA, 2011. [Google Scholar]
  9. Al-Bluwi, I.; Siméon, T.; Cortés, J. Motion Planning Algorithms for Molecular Simulations: A Survey. Comput. Sci. Rev. 2012, 6, 125–143. [Google Scholar] [CrossRef]
  10. Gipson, B.; Hsu, D.; Kavraki, L.; Latombe, J.C. Computational models of protein kinematics and dynamics: Beyond simulation. Ann. Rev. Anal. Chem. 2012, 5, 273–291. [Google Scholar] [CrossRef] [PubMed]
  11. Shehu, A.; Plaku, E. A Survey of Computational Treatments of Biomolecules by Robotics-Inspired Methods Modeling Equilibrium Structure and Dynamic. J. Artif. Intell. Res. 2016, 57, 509–572. [Google Scholar]
  12. Cortés, J.; Al-Bluwi, I. A Robotics Approach To Enhance Conformational Sampling Of Proteins. In Proceedings of the ASME 2012 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference, Chicago, IL, USA, 12–15 August 2012; Volume 4, pp. 1177–1186. [Google Scholar]
  13. Parsons, D.; Canny, J. Geometric Problems in Molecular Biology and Robotics. Proc. Int. Conf. Intell. Syst. Mol. Biol. 1994, 2, 322–330. [Google Scholar] [PubMed]
  14. Siciliano, B.; Khatib, O. (Eds.) Springer Handbook of Robotics; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  15. Metropolis, N.; Rosenbluth, A.; Rosenbluth, M.; Teller, A.; Teller, E. Equation of State Calculations by Fast Computing Machines. J. Chem. Phys. 1953, 21, 1087–1092. [Google Scholar] [CrossRef]
  16. Gō, N.; Scheraga, H. Ring Closure and Local Conformational Deformations of Chain Molecules. Macromolecules 1970, 3, 178–187. [Google Scholar] [CrossRef]
  17. Dodd, L.R.; Boone, T.D.; Theodorou, D.N. A Concerted Rotation Algorithm for Atomistic Monte Carlo Simulation of Polymer Melts and Glasses. Mol. Phys. 1993, 78, 961–996. [Google Scholar] [CrossRef]
  18. Leontidis, E.; de Pablo, J.J.; Laso, M.; Suter, U.W. A Critical Evaluation of Novel Algorithms for the Off-Lattice Monte Carlo Simulation of Condensed Polymer Phases. Adv. Polym. Sci. 1994, 116, 283–318. [Google Scholar]
  19. Wu, N.G.; Deem, M.W. Efficient Monte Carlo Methods for Cyclic Peptides. Mol. Phys. 1999, 94, 559–580. [Google Scholar] [CrossRef]
  20. Betancourt, M.R. Efficient Monte Carlo Trial Moves for Polypeptide Simulations. J. Chem. Phys. 2005, 123, 174905. [Google Scholar] [CrossRef] [PubMed]
  21. Davis, I.W.; Arendall, W.B., III; Richardson, D.C.; Richardson, J.S. The Backrub Motion: How Protein Backbone Shrugs When a Sidechain Dances. Structure 2006, 14, 265–274. [Google Scholar] [CrossRef] [PubMed]
  22. Bottaro, S.; Boomsma, W.; Johansson, K.E.; Andreetta, C.; Hamelryck, T.; Ferkinghoff-Borg, J. Subtle Monte Carlo Updates in Dense Molecular Systems. J. Chem. Theory Comput. 2012, 8, 695–702. [Google Scholar] [CrossRef] [PubMed]
  23. Vitalis, A.; Pappu, R.V. Methods for Monte Carlo Simulations of Biomacromolecules. Annu. Rep. Comput. Chem. 2009, 5, 49–76. [Google Scholar] [PubMed]
  24. Kollman, P.; Dixon, R.; Cornell, W.; Fox, T.; Chipot, C.; Pohorille, A. The development/application of a “minimalist” organic/biochemical molecular mechanic force field using a combination of ab initio calculations and experimental data. In Computer Simulation of Biomolecular Systems; van Gunsteren, W., Weiner, P., Wilkinson, A., Eds.; Springer: Dordrecht, The Netherlands, 1997; Volume 3, pp. 83–96. [Google Scholar]
  25. Bondi, A. Van der Waals Volumes and Radii. J. Phys. Chem. 1964, 68, 441–451. [Google Scholar] [CrossRef]
  26. Kurochkina, N.; Guha, U. SH3 domains: Modules of protein–protein interactions. Biophys. Rev. 2013, 5, 29–39. [Google Scholar] [CrossRef] [PubMed]
  27. Brocca, S.; Samalíkovà, M.; Uversky, V.N.; Lotti, M.; Vanoni, M.; Alberghina, L.; Grandori, R. Order propensity of an intrinsically disordered protein, the cyclin-dependent-kinase inhibitor Sic1. Proteins 2009, 76, 731–746. [Google Scholar] [CrossRef] [PubMed]
  28. Bernadó, P.; Blanchard, L.; Timmins, P.; Marion, D.; Ruigrok, R.; Blackledge, M. A structural model for unfolded proteins from residual dipolar couplings and small-angle X-ray scattering. Proc. Natl. Acad. Sci. USA 2005, 102, 17002–17005. [Google Scholar] [CrossRef] [PubMed]
  29. Krivov, G.G.; Shapovalov, M.V.; Dunbrack, R.L., Jr. Improved Prediction of Protein Side-chain Conformations with SCWRL4. Proteins Struct. Funct. Bioinf. 2009, 77, 778–795. [Google Scholar] [CrossRef] [PubMed]
  30. Ulmschneider, J.P.; Jorgensen, W.L. Monte Carlo backbone sampling for polypeptides with variable bond angles and dihedral angles using concerted rotations and a Gaussian bias. J. Chem. Phys. 2003, 118, 4261–4271. [Google Scholar] [CrossRef]
  31. Pickart, C.M.; Eddins, M.J. Ubiquitin: Structures, functions, mechanisms. Biochim. Biophys. Acta (BBA) Mol. Cell Res. 2004, 1695, 55–72. [Google Scholar] [CrossRef] [PubMed]
  32. Scott, R.; Scheraga, H. Conformational Analysis of Macromolecules. II. The Rotational Isomeric States of the Normal Hydrocarbons. J. Chem. Phys. 1966, 44, 3054–3069. [Google Scholar] [CrossRef]
  33. Ulmschneider, J.P.; Jorgensen, W.L. Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation. J. Am. Chem. Soc. 2004, 126, 1849–1857. [Google Scholar] [CrossRef] [PubMed]
  34. Craig, J.J. Introduction to Robotics; Addison-Wesley: Reading, MA, USA, 1989. [Google Scholar]
  35. Nilmeier, J.; Hua, L.; Coutsias, E.A.; Jacobson, M.P. Assessing Protein Loop Flexibility by Hierarchical Monte Carlo Sampling. J. Chem. Theory Comput. 2011, 7, 1564–1574. [Google Scholar] [CrossRef] [PubMed]
  36. Nilmeier, J.; Jacobson, M.P. Multiscale Monte Carlo Sampling of Protein Sidechains: Application to Binding Pocket Flexibility. J. Chem. Theory Comput. 2008, 4, 835–846. [Google Scholar] [CrossRef] [PubMed]
  37. Nilmeier, J.; Jacobson, M.P. Monte Carlo Sampling with Hierarchical Move Sets: POSH Monte Carlo. J. Chem. Theory Comput. 2009, 5, 1968–1984. [Google Scholar] [CrossRef] [PubMed]
  38. Wu, N.G.; Deem, M.W. Analytical Rebridging Monte Carlo: Application to cis/trans Isomerization in Proline-Containing, Cyclic Peptides. J. Chem. Phys. 1999, 111, 6625–6632. [Google Scholar] [CrossRef]
  39. Canutescu, A.A.; Dunbrack, R.L., Jr. Cyclic Coordinate Descent: A Robotics Algorithm for Protein Loop Closure. Protein Sci. 2003, 12, 963–972. [Google Scholar] [CrossRef] [PubMed]
  40. Renaud, M. A simplified inverse kinematic model calculation method for all 6R type manipulators. In Current Advances in Mechanical Design and Production VII; Hassan, M.F., Megahed, S.M., Eds.; Pergamon: New York, NY, USA, 2000; pp. 57–66. [Google Scholar]
  41. Renaud, M. Calcul des Modèles Géométriques Inverses des Robots Manipulateurs 6R; Rapport LAAS 06332; LAAS: Toulouse, France, 2006. [Google Scholar]
  42. Lee, H.Y.; Liang, C.G. A New Vector Theory for the Analysis of Spatial Mechanisms. Mech. Mach. Theory 1988, 23, 209–217. [Google Scholar] [CrossRef]
  43. Lee, H.Y.; Liang, C.G. Displacement Analysis of the General Spatial 7-Link 7R Mechanisms. Mech. Mach. Theory 1988, 23, 219–226. [Google Scholar]
  44. Manocha, D.; Canny, J.F. Efficient Inverse Kinematics for General 6R Manipulators. IEEE Trans. Robot. Autom. 1994, 10, 648–657. [Google Scholar] [CrossRef]
  45. Golub, G.H.; Van Loan, C.F. Matrix Computations, 3rd ed.; Johns Hopkins University Press: Baltimore, MD, USA, 1996. [Google Scholar]
  46. Anderson, E.; Bai, Z.; Bischof, C.; Blackford, S.; Demmel, J.; Dongarra, J.; Du Croz, J.; Greenbaum, A.; Hammarling, S.; McKenney, A.; et al. LAPACK Users’ Guide, 3rd ed.; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1999. [Google Scholar]
  47. Cortés, J.; Siméon, T.; Remaud-Siméon, M.; Tran, V. Geometric Algorithms for the Conformational Analysis of Long Protein Loops. J. Comput. Chem. 2004, 25, 956–967. [Google Scholar] [CrossRef] [PubMed]
  48. Cortés, J.; Carrión, S.; Curco, D.; Renaud, M.; Aleman, C. Relaxation of Amorphous Multichain Polymer Systems using Inverse Kinematics. Polymer 2010, 51, 4008–4014. [Google Scholar] [CrossRef]
  49. Dinner, A.R. Local deformations of polymers with nonplanar rigid main-chain internal coordinates. J. Comput. Chem. 2000, 21, 1132–1144. [Google Scholar] [CrossRef]
  50. Coutsias, E.A.; Seok, C.; Jacobson, M.P.; Dill, K.A. A Kinematic View of Loop Closure. J. Comput. Chem. 2004, 25, 510–528. [Google Scholar] [CrossRef] [PubMed]
  51. Mezei, M. Efficient Monte Carlo Sampling of Long Molecular Chains Using Local Moves, Tested on a Solvated Lipid Bilayer. J. Chem. Phys. 2003, 118, 3874–3879. [Google Scholar] [CrossRef]
  52. Okamoto, Y. Generalized-ensemble algorithms: Enhanced sampling techniques for Monte Carlo and molecular dynamics simulations. J. Mol. Graph. Model. 2004, 22, 425–439. [Google Scholar] [CrossRef] [PubMed]
  53. Carr, J.M.; Wales, D.J. Global optimization and folding pathways of selected alpha-helical proteins. J. Chem. Phys. 2005, 123, 234901. [Google Scholar] [CrossRef] [PubMed]
  54. Cortés, J.; Le, D.; Iehl, R.; Siméon, T. Simulating ligand-induced conformational changes in proteins using a mechanical disassembly method. Phys. Chem. Chem. Phys. 2010, 12, 8268–8276. [Google Scholar] [CrossRef] [PubMed]
  55. Mandell, D.J.; Kortemme, T. Backbone flexibility in computational protein design. Curr. Opin. Biotechnol. 2009, 20, 420–428. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Not Available.
Figure 1. Illustration of the implemented move classes. (a) pivot move, which is the most frequently used move class within MC methods applied to chain-like molecules; (b) concerted rotation, as implemented using the proposed tripeptide-based model; (c) local move produced by the perturbation of the particle associated with one tripeptide; and (d) hinge move produced by a correlated perturbation of a set of particles. The protein backbone trace before and after the move are drawn in black and blue, respectively. Tripeptides that change their conformation are represented with dashed lines.
Figure 1. Illustration of the implemented move classes. (a) pivot move, which is the most frequently used move class within MC methods applied to chain-like molecules; (b) concerted rotation, as implemented using the proposed tripeptide-based model; (c) local move produced by the perturbation of the particle associated with one tripeptide; and (d) hinge move produced by a correlated perturbation of a set of particles. The protein backbone trace before and after the move are drawn in black and blue, respectively. Tripeptides that change their conformation are represented with dashed lines.
Molecules 23 00373 g001
Figure 2. Proteins used for the analysis of the methods: (a) SH3 domain; (b) Sic1 protein.
Figure 2. Proteins used for the analysis of the methods: (a) SH3 domain; (b) Sic1 protein.
Molecules 23 00373 g002
Figure 3. Projection of sampled states on distance vs. energy plots for the SH3 domain and the Sic1 protein.
Figure 3. Projection of sampled states on distance vs. energy plots for the SH3 domain and the Sic1 protein.
Molecules 23 00373 g003
Figure 4. Plots of the function rmsd ( τ ) for each move class during MC runs for the SH3 domain (left) and the Sic1 protein (right).
Figure 4. Plots of the function rmsd ( τ ) for each move class during MC runs for the SH3 domain (left) and the Sic1 protein (right).
Molecules 23 00373 g004
Figure 5. Plots of the average autocorrelation function over the 20 central dihedral angles of 14-alanine for the ConRot, OneParticle, and Mixed move classes.
Figure 5. Plots of the average autocorrelation function over the 20 central dihedral angles of 14-alanine for the ConRot, OneParticle, and Mixed move classes.
Molecules 23 00373 g005
Figure 6. Autocorrelation times τ i n t (in MC steps) of the 20 central dihedral angles of 14-alanine for the ConRot, OneParticle, and Mixed move classes.
Figure 6. Autocorrelation times τ i n t (in MC steps) of the 20 central dihedral angles of 14-alanine for the ConRot, OneParticle, and Mixed move classes.
Molecules 23 00373 g006
Figure 7. Structure of ubiquitin.
Figure 7. Structure of ubiquitin.
Molecules 23 00373 g007
Figure 8. Ubiquitin C α root mean squared fluctuations (RMSF) from 15 × 10 6 steps of MC simulations using ConRot, OneParticle, and Mixed move classes.
Figure 8. Ubiquitin C α root mean squared fluctuations (RMSF) from 15 × 10 6 steps of MC simulations using ConRot, OneParticle, and Mixed move classes.
Molecules 23 00373 g008
Figure 9. Evolution of the cumulative C α root mean squared fluctuations (RMSF) of ubiquitin with the number of MC steps using ConRot, OneParticle, and Mixed move classes.
Figure 9. Evolution of the cumulative C α root mean squared fluctuations (RMSF) of ubiquitin with the number of MC steps using ConRot, OneParticle, and Mixed move classes.
Molecules 23 00373 g009
Figure 10. Geometric model of the protein backbone: (a) in a general case; (b) in the particular case of a rigid peptide bond. In the general case, the mDH parameters involved in the homogeneous transformation matrices between consecutive atom groups are obtained directly from bond lengths and bond angles. In the case of a rigid peptide bond, this requires some (simple) geometric operations. Note that the origin of a frame attached to an atom group does not necessarily correspond to an atom position.
Figure 10. Geometric model of the protein backbone: (a) in a general case; (b) in the particular case of a rigid peptide bond. In the general case, the mDH parameters involved in the homogeneous transformation matrices between consecutive atom groups are obtained directly from bond lengths and bond angles. In the case of a rigid peptide bond, this requires some (simple) geometric operations. Note that the origin of a frame attached to an atom group does not necessarily correspond to an atom position.
Molecules 23 00373 g010
Figure 11. Illustration of the protein subdivision approach. Fragments of three amino-acid residues are treated as kinematic chains, similar to robotic manipulators.
Figure 11. Illustration of the protein subdivision approach. Fragments of three amino-acid residues are treated as kinematic chains, similar to robotic manipulators.
Molecules 23 00373 g011
Table 1. Perturbation step-sizes.
Table 1. Perturbation step-sizes.
OneTorsionConRotOneParticleHinge
δ b δ b δ pt δ pr δ pr
SH3 domain0.01 rad0.025 rad0.05 Å0.003 rad0.01 rad
Sic1 protein0.02 rad0.025 rad0.05 Å0.003 rad0.02 rad
Table 2. Computational performance.
Table 2. Computational performance.
Move ClassAcc. Rate# IterationsT CPU
SH3 domainOneTorsion0.68 3.28 × 10 6 63 h
ConRot0.56 3.55 × 10 6 51 h
OneParticle0.42 4.06 × 10 6 56 h
Hinge0.59 3.46 × 10 6 57 h
Mixed0.56 3.58 × 10 6 57 h
Sic1 proteinOneTorsion0.56 3.54 × 10 6 89 h
ConRot0.65 3.25 × 10 6 63 h
OneParticle0.52 3.66 × 10 6 69 h
Hinge0.53 3.60 × 10 6 75 h
Mixed0.57 3.49 × 10 6 74 h
Table 3. Autocorrelation times τ i n t (in MC steps) of the ConRot, OneTorsion, and Mixed move classes.
Table 3. Autocorrelation times τ i n t (in MC steps) of the ConRot, OneTorsion, and Mixed move classes.
Move ClassAverageMinMedianMax
ConRot6016512531514,070
OneParticle3426343121511,281
Mixed15181579853706

Share and Cite

MDPI and ACS Style

Denarie, L.; Al-Bluwi, I.; Vaisset, M.; Siméon, T.; Cortés, J. Segmenting Proteins into Tripeptides to Enhance Conformational Sampling with Monte Carlo Methods. Molecules 2018, 23, 373. https://doi.org/10.3390/molecules23020373

AMA Style

Denarie L, Al-Bluwi I, Vaisset M, Siméon T, Cortés J. Segmenting Proteins into Tripeptides to Enhance Conformational Sampling with Monte Carlo Methods. Molecules. 2018; 23(2):373. https://doi.org/10.3390/molecules23020373

Chicago/Turabian Style

Denarie, Laurent, Ibrahim Al-Bluwi, Marc Vaisset, Thierry Siméon, and Juan Cortés. 2018. "Segmenting Proteins into Tripeptides to Enhance Conformational Sampling with Monte Carlo Methods" Molecules 23, no. 2: 373. https://doi.org/10.3390/molecules23020373

APA Style

Denarie, L., Al-Bluwi, I., Vaisset, M., Siméon, T., & Cortés, J. (2018). Segmenting Proteins into Tripeptides to Enhance Conformational Sampling with Monte Carlo Methods. Molecules, 23(2), 373. https://doi.org/10.3390/molecules23020373

Article Metrics

Back to TopTop