Biomedical subjects

Marcos R Betancourt

Publications and source records attributed to Marcos R Betancourt.

3 recordsLinked to original sources

PubMed2005

Efficient Monte Carlo trial moves for polypeptide simulations.

A new move set for the Monte Carlo simulations of polypeptide chains is introduced. It consists of a rigid rotation along the (C(alpha)) ends of an arbitrary long segment of the backbone in such a way that the atoms outside this segment remain fixed. This fixed end move, or FEM, alters only the backbone dihedral angles phi and psi and the C(alpha) bond angles of the segment ends. Rotations are restricted to those who keep the alpha bond angles within their maximum natural range of approximately +/-10 degrees. The equations for the angular intervals (tau) of the allowed rigid rotations and the equations required for satisfying the detailed balance condition are presented in detail. One appealing property of the FEM is that the required number of calculations is minimal, as it is evident from the simplicity of the equations. In addition, the moving backbone atoms undergo considerable but limited displacements of up to 3 A. These properties, combined with the small number of backbone angles changed, lead to high acceptance rates for the new conformations and make the algorithm very efficient for sampling the conformational space. The FEMs, combined with pivot moves, are used in a test to fold a group of coarse-grained proteins with lengths of up to 200 residues.

Algorithms↗

PubMed2004

Local propensities and statistical potentials of backbone dihedral angles in proteins.

The following three issues concerning the backbone dihedral angles of protein structures are presented. (1) How do the dihedral angles of the 20 amino acids depend on the identity and conformation of their nearest residues? (2) To what extent are the native dihedral angles determined by local (dihedral) potentials? (3) How to build a knowledge-based potential for a residue's dihedral angles, considering the identity and conformation of its nearest residues? We find that the dihedral angle distribution for a residue can significantly depend on the identity and conformation of its adjacent residues. These correlations are in sharp contrast to the Flory isolated-pair hypothesis. Statistical potentials are built for all combinations of residue triplets and depend on the dihedral angles between consecutive residues. First, a low-resolution potential is obtained, which only differentiates between the main populated basins in the dihedral angle density plots. Minimization of the dihedral potential for 125 test proteins reveals that most native alpha-helical residues (89%) and a large fraction of native beta-sheet residues (47%) adopt conformations close to their native one. For native loop residues, the percentage is 48%. It is also found that this fraction is higher for residues away from the ends of alpha or beta secondary structure elements. In addition, a higher resolution potential is built as a function of dihedral angles by a smoothing procedure and continuous functions interpolations. Monte Carlo energy minimization with this potential results in a lower fraction for native beta-sheet residues. Nevertheless, because of the higher flexibility and entropy of beta structures, they could be preferred under the influence of non-local interactions. In general, most alpha-helices and many beta-sheets are strongly determined by the local potential, while the conformations in loops and near the end of beta-sheets are more influenced by non-local interactions.

Amino Acids↗

PubMed2003

A reduced protein model with accurate native-structure identification ability.

A protein model that is simple enough to be used in protein-folding simulations but accurate enough to identify a protein native fold is described. Its geometry consists of describing the residues by one, two, or three pseudoatoms, depending on the residue size. Its energy is given by a pairwise, knowledge-based potential obtained for all the pseudoatoms as a function of their relative distance. The pseudoatomic potential is also a function of the primary chain separation and residue order. The model is tested by gapless threading on a large, representative set of known protein and decoy structures obtained from the "Decoys 'R' Us" database. It is also tested by threading on gapped decoys generated for proteins with many homologs. The gapless threading tests show near 98% native-structure recognition as the lowest energy structure and almost 100% as one of the three lowest energy structures for over 2200 test proteins. In decoy threading tests, the model recognized the majority of the native structures. It is also able to recognize native structures among gapped decoys, in spite of close structural similarities. The results indicate that the pseudoatomic model has native recognition ability similar to comparable atomic-based models but much better than equivalent residue-based models.

Algorithms↗