Search

We employ all-atom structure-based models with a force field with multiple energetic basins for the C-terminal (residues 166-226) of the mammalian prion protein. One basin represents the known alpha-helical (αH) structure while the other represents the same residues in a left-handed beta-helical (LHBH) conformation. The LHBH structure has been proposed to help describe one class of in vitro grown fibrils, as well as possibly self-templating the conversion of normal cellular prion protein to the infectious form. Yet, it is unclear how the protein may make this global rearrangement. Our results demonstrate that the conformation changes are not strongly limited by large-scale geometry modification and that there may exist an overall preference for the LHBH conformation. Furthermore, our model presents novel intermediate trapping conformations with twisted LHBH structure.

Uridine 5′-diphosphate-N-acetylglucosamine (UDP-GlcNAc) acyltransferase (LpxA) catalyzes a reversible reaction for adding anO-acyl group to the GlcNAc in UDP-GlcNAc in the first step of lipid A biosynthesis. Lipid A constitutes a major component of lipopolysaccharides, also referred to as endotoxins, which form the outer monolayer of the outer membrane of Gram-negative bacteria. Ligand-free and UDP-GlcNAc-bound crystal structures of LpxA fromBacteroides fragilisNCTC 9343, the most common pathogenic bacteria found in abdominal abscesses, have been determined and are presented here. The enzyme crystallizes in a cubic space group, with the crystallographic threefold axis generating the biological functional homotrimer and with each monomer forming a nine-rung left-handed β-helical (LβH) fold in the N-terminus followed by an α-helical motif in the C-terminus. The structure is highly similar to LpxA from other organisms. Yet, despite sharing a similar LβH structure with LpxAs fromEscherichia coliand others, previously unseen calcium ions are observed on the threefold axis inB. fragilisLpxA to help stabilize the trimeric assembly.

The PI led theoretical studies of correlated hybridization in transition metal complexes, compounds, and molecules, and of electron transport in DNA associated with nanoelectronic conformations attached to gold electrodes and in the presence of DNA repair proteins.

The left-handed beta helix is an intriguing structural motif in several known proteins. We attempt to elucidate the factors that contribute to its stability with a theoretical-computational approach. Combining a novel form of coarse-grained molecular dynamics with parallel tempering affords access to the regime of equilibration in short peptides. For example, this method has been used to reproduce important features of the helix-coil transition in polyalanine.Three-layered beta-helical fragments are formed from various sequences that are superposed upon two backbone templates taken from sections of ideal type-I and type-II beta helices. Native, native-like, and various homogeneous sequences are simulated and their stability analyzed relative to the initial structure.An all-atom potential energy and associated parameters drawn from the Assisted Model-Building and Energy Refinement package with slight modifications. The effects of aqueous solvent are treated with the generalized Born model and a recently-proposed hydrophobic potential of mean force. Larger conformational changes can be explored by treating non-bonded forces according to Brownian dynamics while simultaneously maintaining molecular geometry with a separate algorithm.

We propose models for in vitro grown mammalian prion protein fibrils based upon left handed beta helices formed both from the N-terminal and C-terminal regions of the proteinase resistant infectious prion core. The C-terminal threading onto a beta-helical structure is almost uniquely determined by fixing the cysteine disulfide bond on a helix corner. In comparison to known left handed helical peptides, the resulting model structures have similar stability attributes including relatively low root mean square deviations in all atom molecular dynamics, substantial side-chain-to-side-chain hydrogen bonding, good volume packing fraction, and low hydrophilic/hydrophobic frustration. For the N-terminus, we propose a new threading of slightly more than two turns, which improves upon the above characteristics relative to existing three turn beta-helical models. The N-terminal and C-terminal beta helices can be assembled into eight candidate models for the fibril repeat units, held together by large hinge (order 30 residues) domain swapping, with three amenable to fibril promoting domain swapping via a small (five residue) hinge on the N-terminal side. Small concentrations of the metastable C-terminal beta helix in vivo might play a significant role in templating the infectious conformation and in enhancing conversion kinetics for inherited forms of the disease and explain resistance (for canines) involving hypothesized coupling to the methionine 129 sulfur known to play a role in human disease.

We combine molecular dynamics simulations and density functional theory to analyze the electrical structure and transmission probability in four different DNA sequences under physiological conditions. The conductance in these sequences is primarily controlled by interstrand and intrastrand coupling between low-energy guanine orbitals. Insertion of adenine-thymine base pairs between the guanine-cytosine rich domains acts as a tunneling barrier. Our theory explains recent length dependent conductance data for individual DNA molecules in water.

The mechanism of site-specific recognition of DNA by proteins has been a long-standing issue. The DNA glycosylase MutY, for instance, must find the rare 8-oxoguanine-adenine mismatches among the large number of basepairs in the DNA. This protein has a [4Fe-4S] cluster, which is highly conserved in species as diverse as Escherichia Coli and Homo sapiens. The mixed-valent nature of this cluster suggests that charge transfer may play a role in MutY's function. We have studied the energetics of the charge transfer in Bacillus stearothermophilus MutY-DNA complex using multiscale calculation including density functional theory and molecular dynamics. The [4Fe-4S] cluster in MutY is found to undergo 2+ to 3+ oxidation when coupling to DNA through hole transfer, especially when MutY is near an oxoguanine modified base (oxoG). Employing the Marcus theory for electron transfer, we find near optimal Frank-Condon factors for electron transfer from MutY to oxoguanine modified base. MutY has modest selectivity for oxoguanine over guanine due to the difference in oxidation potential. The tunneling matrix element is significantly reduced with the mutation R149W, whereas the mutation L154F reduces the tunneling matrix element as well as the Frank-Condon factor. Both L154F and R149W mutations are known to dramatically reduce or eliminate repair efficiency. We suggest a scenario where the charge transfer leads to a stabilization of the specific binding conformation, which is likely the recognition mode, thus enabling it to find the damaged site efficiently.