Helices in biomolecules.
Identical objects, regularly assembled, form a helix, which is the principal motif of nucleic acids, proteins, and viral capsids.
Biopolymers↗
Biomedical subjects
Kevin Cahill
Publications and source records attributed to Kevin Cahill.
Rydberg-London potential for diatomic molecules and unbonded atom pairs.
We propose and test a pair potential that is accurate at all relevant distances and simple enough for use in large-scale computer simulations. A combination of the Rydberg potential from spectroscopy and the London inverse-sixth-power energy, the proposed form fits spectroscopically determined potentials better than the Morse, Varnshi, and Hulburt-Hirschfelder potentials and much better than the Lennard-Jones and harmonic potentials. At long distances, it goes smoothly to the London force appropriate for gases and preserves van der Waals's "continuity of the gas and liquid states," which is routinely violated by coefficients assigned to the Lennard-Jones 6-12 form.
Journal Article↗
Alternative splicing and genomic stability.
Alternative splicing allows an organism to make different proteins in different cells at different times, all from the same gene. In a cell that uses alternative splicing, the total length of all the exons is much shorter than in a cell that encodes the same set of proteins without alternative splicing. This economical use of exons makes genes more stable during reproduction and development because a genome with a shorter exon length is more resistant to harmful mutations. Genomic stability may be the reason why higher vertebrates splice alternatively. For a broad class of alternatively spliced genes, a formula is given for the increase in their stability.
Alternative Splicing↗
On the kinematics of protein folding.
We offer simple solutions to three kinematic problems that occur in the folding of proteins. We show how to construct suitably local elementary Monte Carlo moves, how to close a loop, and how to fold a loop without breaking the bond that closes it.
Algorithms↗
Proteins wriggle.
We propose an algorithmic strategy for improving the efficiency of Monte Carlo searches for the low-energy states of proteins. Our strategy is motivated by a model of how proteins alter their shapes. In our model, when proteins fold under physiological conditions, their backbone dihedral angles change synchronously in groups of four or more to avoid steric clashes and respect the kinematic conservation laws. They wriggle; they do not thrash. We describe a simple algorithm that can be used to incorporate wriggling in Monte Carlo simulations of protein folding. We have tested this wriggling algorithm against a code in which the dihedral angles are varied independently (thrashing). Our standard of success is the average root-mean-square distance (rmsd) between the alpha-carbons of the folding protein and those of its native structure. After 100,000 Monte Carlo sweeps, the relative decrease in the mean rmsd, as one switches from thrashing to wriggling, rises from 11% for the protein 3LZM with 164 amino acids (aa) to 40% for the protein 1A1S with 313 aa and 47% for the protein 16PK with 415 aa. These results suggest that wriggling is useful and that its utility increases with the size of the protein. One may implement wriggling on a parallel computer or a computer farm.
Algorithms↗