Biomedical subjects

Bernard R Brooks

Publications and source records attributed to Bernard R Brooks.

7 recordsLinked to original sources

PubMed2003

An improved method for nonbonded list generation: rapid determination of near-neighbor pairs.

Generation of the list of near-neighbor pairs of atoms not bonded to each other is a key feature of many programs for calculating the energy and energy derivatives for large molecules. Because this step can take a significant amount of CPU time, more efficient nonbonded list generation can speed up the energy calculations. In this article, a novel nonbonded list generation algorithm, BYCC, is introduced. It combines certain features of other algorithms and achieves more rapid nonbonded list generation; a factor of approximately 2.5 for a molecule of 5000 atoms with a cutoff in the 10 A range is obtained on Hewlett-Packard (HP) and Alpha processors, without greatly increasing memory requirements.

Algorithms↗

PubMed2003

Annealing function of GroEL: structural and bioinformatic analysis.

The Escherichia coli chaperonin system, GroEL-GroES, facilitates folding of substrate proteins (SPs) that are otherwise destined to aggregate. The iterative annealing mechanism suggests that the allostery-driven GroEL transitions leading to changes in the microenvironment of the SP constitutes the annealing action of chaperonins. To describe the molecular basis for the changes in the nature of SP-GroEL interactions we use the crystal structures of GroEL (T state), GroEL-ATP (R state) and the GroEL-GroES-(ADP)(7) (R" state) complex to determine the residue-specific changes in the accessible surface area and the number of tertiary contacts as a result of the T-->R-->R" transitions. We find large changes in the accessible area in many residues in the apical domain, but relatively smaller changes are associated with residues in the equatorial domain. In the course of the T-->R transition the microenvironment of the SP changes which suggests that GroEL is an annealing machine even without GroES. This is reflected in the exposure of Glu386 which loses six contacts in the T-->R transition. We also evaluate the conservation of residues that participate in the various chaperonin functions. Multiple sequence alignments and chemical sequence entropy calculations reveal that, to a large extent, only the chemical identities and not the residues themselves important for the nominal functions (peptide binding, nucleotide binding, GroES and substrate protein release) are strongly conserved. Using chemical sequence entropy, which is computed by classifying aminoacids into four types (hydrophobic, polar, positively charged and negatively charged) we make several new predictions that are relevant for peptide binding and annealing function of GroEL. We identify a number of conserved peptide binding sites in the apical domain which coincide with those found in the 1.7 A crystal structure of 'mini-chaperone' complexed with the N-terminal tag. Correlated mutations in the HSP60 family, that might control allostery in GroEL, are also strongly conserved. Most importantly, we find that charged solvent-exposed residues in the T state (Lys 226, Glu 252 and Asp 253) are strongly conserved. This leads to the prediction that mutating these residues, that control the annealing function of the SP, can decrease the efficacy of the chaperonin function.

Algorithms↗

PubMed2003

A core-weighted fitting method for docking atomic structures into low-resolution maps: application to cryo-electron microscopy.

Cryo-electron microscopy of "single particles" is a powerful method to analyze structures of large macromolecular assemblies that are not amenable to investigation by traditional X-ray crystallographic methods. A key step in these studies is to obtain atomic interpretations of multiprotein complexes by fitting atomic structures of individual components into maps obtained from electron microscopic data. Here, we report the use of a "core-weighting" method, combined with a grid-threading Monte Carlo (GTMC) approach for this purpose. The "core" of an individual structure is defined to represent the part where the density distribution is least likely to be altered by other components that comprise the macromolecular assembly of interest. The performance of the method has been evaluated by its ability to determine the correct fit of (i) the alpha-chain of the T-cell receptor variable domain into a simulated map of the alphabeta complex at resolutions between 5 and 40 A, and (ii) the E2 catalytic domain of the pyruvate dehydrogenase into an experimentally determined map, at 14 A resolution, of the icosahedral complex formed by 60 copies of this enzyme. Using the X-ray structures of the two test cases as references, we demonstrate that, in contrast to more traditional methods, the combination of the core-weighting method and the grid-threading Monte Carlo approach can identify the correct fit reliably and rapidly from the low-resolution maps that are typical of structures determined with the use of single-particle electron microscopy.

Animals↗

PubMed2002

Molecular architecture and mechanism of an icosahedral pyruvate dehydrogenase complex: a multifunctional catalytic machine.

Electron cryo-microscopy of 'single particles' is a powerful method to determine the three-dimensional (3D) architectures of complex cellular assemblies. The pyruvate dehydrogenase multi-enzyme complex couples the activity of three component enzymes (E1, E2 and E3) in the oxidative decarboxylation of pyruvate to generate acetyl-CoA, linking glycolysis and the tricarboxylic acid cycle. We report here a 3D model for an 11 MDa, icosahedral pyruvate dehydrogenase sub-complex, obtained by combining a 28 A structure derived from electron cryo-microscopy with previously determined atomic coordinates of the individual E1 and E2 components. A key feature is that the E1 molecules are located on the periphery of the assembly in an orientation that allows each of the 60 mobile lipoyl domains tethered to the inner E2 core to access multiple E1 and E2 active sites from inside the icosahedral complex. This unexpected architecture provides a highly efficient mechanism for active site coupling and catalytic rate enhancement by the motion of the lipoyl domains in the restricted annular region between the inner core and outer shell of the complex.

Catalysis↗

PubMed2002

Hydrogen bonding in high-resolution protein structures: a new method to assess NMR protein geometry.

An analysis of backbone hydrogen bonds has been performed on nine high-resolution protein X-ray crystal structures. Backbone hydrogen-bond geometry is compared in the context of X-ray crystal structure resolution. A strong correlation between the hydrogen-bond distance, R(HO), and the hydrogen-bond angle, theta(NHO), is observed when the X-ray crystal structure resolution is <1.00 A. Ab initio calculations were performed to substantiate these results. The angle and distance limits found in our correlation for the backbone hydrogen-bond geometry can be used to evaluate the quality of protein structures and for further NMR structure refinement.

Crystallography, X-Ray↗

PubMed2002

Direct observation of the folding and unfolding of a beta-hairpin in explicit water through computer simulation.

The cooperative folding and unfolding of a beta-hairpin structure are observed in explicit water at native folding conditions through self-guided molecular dynamics simulation. The folded structure agrees excellently with the NMR NOE data. After going through a fully hydrated state, the peptide folds into a beta-hairpin structure in a highly cooperative process. During the folding process it is observed that side chain interaction occurs first, while intrapeptide hydrogen bonds only form at the final stage. On the contrary, the unfolding process starts with the breaking of interstrand hydrogen bonds. Energetic analysis indicates that the driving force of the folding is the intrapeptide interaction, while the solvent interaction opposes the folding.

Amino Acid Sequence↗

PubMed2002

Simulations of membranes and other interfacial systems using P2(1) and Pc periodic boundary conditions.

We demonstrate the ease and utility of simulating heterogeneous interfacial systems with P2(1) and Pc periodic boundary conditions which allow, for example, lipids in a membrane to switch leaflets. In preliminary tests, P2(1) was shown to yield equivalent results to P1 in simulations of bulk water, a water/vacuum interface, and pure DPPC bilayers with an equal number of lipids per leaflet; equivalence of Pc and P1 was also demonstrated for the former two systems. P2(1) was further tested in simulations involving the spreading of an octane film on water, and equilibration of a DPPC bilayer from an initial condition containing different numbers of lipids in the two leaflets. Lastly, a simulation in P2(1) of a DOPC/melittin membrane showed significant passage of lipids to the melittin-containing leaflet from the initial distribution, and lends insight into the condensation of lipids by melittin.

1,2-Dipalmitoylphosphatidylcholine↗