Adjacent gauche stabilization in linear alkanes: implications for polymer models and conformational analysis

High-level ab initio quantum mechanical calculations are used to study various gauche conformational energies of n-pentane to n-decane. The destabilizing "pentane effect" (adjacent gauche states of opposite sign) for alkanes is confirmed, but the energies were found to depend slightly on chain length. In contrast, introducing an adjacent gauche of the same sign requires only 0.22-0.37 kcal/mol, approximately half of the single gauche state energy. This adjacent gauche stabilization should be taken into account when formulating or analyzing rotational isomeric models, carrying out conformational analysis, and developing force fields for alkanes, lipids, and related polymers.

A computational investigation of the geometrical structure and protodeboronation of boroglycine, H2N-CH2-B(OH)2

In this article the geometrical structure of the simple, achiral, alpha-amino boronic acid boroglycine, H2N-CH2-B(OH)2, was investigated using density functional theory (DFT), second-order Møller-Plesset (MP2) perturbation theory, and coupled cluster methodology with single- and double-excitations (CCSD); the effects of an aqueous environment were incorporated into the results by using a few explicit water molecules and/or self-consistent reaction field (SCRF) calculations with the IEF polarizable continuum model (PCM). Neutral reaction mechanisms were investigated for the direct protodeboronation (hydrolysis) of boroglycine (H2O+H2N-CH2-B(OH)2-->B(OH)3+H2N-CH3), for which DeltaH degrees 298 was -21.9 kcal/mol at the MP2(FC)/aug-cc-pVDZ level, and for the 1,2-carbon-to-nitrogen shift of the -B(OH)2 moiety (H2N-CH2-B(OH)2-->H3C-NH-B(OH)2), for which the corresponding value of DeltaH degrees 298 was -18.2 kcal/mol. A boron-oxygen double-bonded intermediate was found to play an important role in the 1,2-rearrangement mechanism.

Interfacing Q-Chem and CHARMM to perform QM/MM reaction path calculations

A hybrid quantum mechanical/molecular mechanical (QM/MM) potential energy function with Hartree-Fock, density functional theory (DFT), and post-HF (RIMP2, MP2, CCSD) capability has been implemented in the CHARMM and Q-Chem software packages. In addition, we have modified CHARMM and Q-Chem to take advantage of the newly introduced replica path and the nudged elastic band methods, which are powerful techniques for studying reaction pathways in a highly parallel (i.e., parallel/parallel) fashion, with each pathway point being distributed to a different node of a large cluster. To test our implementation, a series of systems were studied and comparisons were made to both full QM calculations and previous QM/MM studies and experiments. For instance, the differences between HF, DFT, MP2, and CCSD QM/MM calculations of H2O...H2O, H2O...Na+, and H2O...Cl- complexes have been explored. Furthermore, the recently implemented polarizable Drude water model was used to make comparisons to the popular TIP3P and TIP4P water models for doing QM/MM calculations. We have also computed the energetic profile of the chorismate mutase catalyzed Claisen rearrangement at various QM/MM levels of theory and have compared the results with previous studies. Our best estimate for the activation energy is 8.20 kcal/mol and for the reaction energy is -23.1 kcal/mol, both calculated at the MP2/6-31+G(d)//MP2/6-31+G(d)/C22 level of theory.

Allosteric transitions in the chaperonin GroEL are captured by a dominant normal mode that is most robust to sequence variations

The Escherichia coli chaperonin GroEL, which helps proteins to fold, consists of two heptameric rings stacked back-to-back. During the reaction cycle GroEL undergoes a series of allosteric transitions triggered by ligand (substrate protein, ATP, and the cochaperonin GroES) binding. Based on an elastic network model of the bullet-shaped double-ring chaperonin GroEL-(ADP)(7)-GroES structure (R''T state), we perform a normal mode analysis to explore the energetically favorable collective motions encoded in the R''T structure. By comparing each normal mode with the observed conformational changes in the R''T --> TR'' transition, a single dominant normal mode provides a simple description of this highly intricate allosteric transition. A detailed analysis of this relatively high-frequency mode describes the structural and dynamic changes that underlie the positive intra-ring and negative inter-ring cooperativity. The dynamics embedded in the dominant mode entails highly concerted structural motions with approximate preservation of sevenfold symmetry within each ring and negatively correlated ones between the two rings. The dominant normal mode (in comparison with the other modes) is robust to parametric perturbations caused by sequence variations, which validates its functional importance. Response of the dominant mode to local changes that mimic mutations using the structural perturbation method technique leads to a wiring diagram that identifies a network of key residues that regulate the allosteric transitions. Many of these residues are located in intersubunit interfaces, and may therefore play a critical role in transmitting allosteric signals between subunits.

Exploring SCC-DFTB paths for mapping QM/MM reaction mechanisms

A new first-order procedure for locating transition structures (TS) that employs hybrid quantum mechanical/molecular mechanical (QM/MM) potentials has been developed. This new technique (RPATh+RESD) combines the replica path method (RPATh) and standard reaction coordinate driving (RCD) techniques in an approach that both efficiently determines reaction barriers and successfully eliminates two key weaknesses of RCD calculations (i.e., hysteresis/discontinuities in the path and the sequential nature of the RCD procedure). In addition, we have extended CHARMM's QM/MM reaction pathway methods, the RPATh and nudged elastic band (NEB) methods, to incorporate SCC-DFTB wave functions. This newly added functionality has been applied to the chorismate mutase-catalyzed interconversion of chorismate to prephenate, which is a key step in the shikimate pathway of bacteria, fungi, and other higher plants. The RPATh+RESD barrier height (DeltaE=5.7 kcal/mol) is in good agreement with previous results from full-energy surface mapping studies (Zhang, X.; Zhang, X.; Bruice, T. C. Biochemistry 2005, 44, 10443-10448). Full reaction paths were independently mapped with RPATh and NEB methods and showed good agreement with the final transition state from the RPATh+RESD "gold standard" and previous high-level QM/MM transition states (Woodcock, H. L.; Hodoscek, M.; Gilbert, T. B.; Gill, P. M. W.; Schaefer, H. F.; Brooks, B. R. J. Comput. Chem. 2007, 28, 1485-1502). The SCC-DFTB TS geometry most closely approximates the MP2/6-31+G(d) QM/MM result. However, the barrier height is underestimated and possibly points to an area for improvement in SCC-DFTB parametrization. In addition, the steepest descents (SD) minimizer for the NEB method was modified to uncouple the in-path and off-path degrees of freedom during the minimization, which significantly improved performance. The convergence behavior of the RPATh and NEB was examined for SCC-DFTB wave functions, and it was determined that, in general, both methods converge at about the same rate, although the techniques used for convergence may be different. For instance, RPATh can effectively use the adopted basis Newton-Raphson (ABNR) minimizer, where NEB seems to require a combination of SD and ABNR.

Toward the mechanism of dynamical couplings and translocation in hepatitis C virus NS3 helicase using elastic network model

Hepatitis C virus NS3 helicase is an enzyme that unwinds double-stranded polynucleotides in an ATP-dependent reaction. It provides a promising target for small molecule therapeutic agents against hepatitis C. Design of such drugs requires a thorough understanding of the dynamical nature of the mechanochemical functioning of the helicase. Despite recent progress, the detailed mechanism of the coupling between ATPase activity and helicase activity remains unclear. Based on an elastic network model (ENM), we apply two computational analysis tools to probe the dynamical mechanism underlying the allosteric coupling between ATP binding and polynucleotide binding in this enzyme. The correlation analysis identifies a network of hot-spot residues that dynamically couple the ATP-binding site and the polynucleotide-binding site. Several of these key residues have been found by mutational experiments as functionally important, while our analysis also reveals previously unexplored hot-spot residues that are potential targets for future mutational studies. The conformational changes between different crystal structures of NS3 helicase are found to be dominated by the lowest frequency mode solved from the ENM. This mode corresponds to a hinge motion of the highly flexible domain 2. This motion simultaneously modulates the opening/closing of the domains 1-2 cleft where ATP binds, and the domains 2-3 cleft where the polynucleotide binds. Additionally, a small twisting motion of domain 1, observed in both mode 1 and the computed ATP binding induced conformational change, fine-tunes the binding affinity of the domains 1-3 interface for the polynucleotide. The combination of these motions facilitates the translocation of a single-stranded polynucleotide in an inchworm-like manner.

An ab initio study on the torsional surface of alkanes and its effect on molecular simulations of alkanes and a DPPC bilayer

Energies of 119 conformations of normal alkanes from butane to heptane were calculated at approximately the CCSD(T)/cc-pVQZ level. Energies of gauche (g) conformers relative to trans (t) decrease as chain length increases. In what is termed the "positive pentane effect", adjacent gauche conformers of the same sign are stabilized compared to nonadjacent conformers; e.g., for hexane the energies of tgt, tgg, and gtg are 0.600, 0.930, and 1.18 kcal/mol, respectively. Torsional terms in the CHARMM27 (C27) force field were fit to the calculated QM energies to yield a revised potential, C27r. Molecular dynamics simulations of normal alkanes (heptane, decane, tridecane, and pentadecane) with C27r yield higher populations of gauche states, increased transition rates, and improved agreement with experiment as compared to C27. In addition, C27r simulations of a hydrated DPPC lipid bilayer yield improved agreement with the experimental NMR deuterium order parameters for the aliphatic chain ends.

Coupling between allosteric transitions in GroEL and assisted folding of a substrate protein

Escherichia coli chaperonin, GroEL, helps proteins fold under nonpermissive conditions. During the reaction cycle, GroEL undergoes allosteric transitions in response to binding of a substrate protein (SP), ATP, and the cochaperonin GroES. Using coarse-grained representations of the GroEL and GroES structures, we explore the link between allosteric transitions and the folding of a model SP, a de novo-designed four-helix bundle protein, with low spontaneous yield. The ensemble of GroEL-bound SP is less structured than the bulk misfolded structures. Upon binding, which kinetically occurs in two stages, the SP loses not only native tertiary contacts but also experiences a decrease in helical content. During multivalent binding and the subsequent ATP-driven transition of GroEL the SP undergoes force-induced stretching. Upon encapsulation, which occurs upon GroES binding, the SP finds itself in a "hydrophilic" cavity in which it can reach the folded conformation. Surprisingly, we find that the yield of the native state in the expanded GroEL cavity is relatively small even after it remains in it for twice the spontaneous folding time. Thus, in accord with the iterative annealing mechanism, multiple rounds of binding, partial unfolding, and release of the SP are required to enhance the yield of the folded SP.

Long-range Lennard-Jones and electrostatic interactions in interfaces: application of the isotropic periodic sum method

Molecular dynamics (MD) simulations of heptane/vapor, hexadecane/vapor, water/vapor, hexadecane/water, and dipalmitoylphosphatidylcholine (DPPC) bilayers and monolayers are analyzed to determine the accuracy of treating long-range interactions in interfaces with the isotropic periodic sum (IPS) method. The method and cutoff (rc) dependences of surface tensions, density profiles, water dipole orientation, and electrostatic potential profiles are used as metrics. The water/vapor, heptane/vapor, and hexadecane/vapor interfaces are accurately and efficiently calculated with 2D IPS (rc=10 A). It is demonstrated that 3D IPS is not practical for any of the interfacial systems studied. However, the hybrid method PME/IPS [Particle Mesh Ewald for electrostatics and 3D IPS for Lennard-Jones (LJ) interactions] provides an efficient way to include both types of long-range forces in simulations of large liquid/vacuum and all liquid/liquid interfaces, including lipid monolayers and bilayers. A previously published pressure-based long-range LJ correction yields results similar to those of PME/IPS for liquid/liquid interfaces. The contributions to surface tension of LJ terms arising from interactions beyond 10 A range from 13 dyn/cm for the hexadecane/vapor interface to approximately 3 dyn/cm for hexadecane/water and DPPC bilayers and monolayers. Surface tensions of alkane/vapor, hexadecane/water, and DPPC monolayers based on the CHARMM lipid force fields agree very well with experiment, whereas surface tensions of the TIP3P and TIP4P-Ew water models underestimate experiment by 16 and 11 dyn/cm, respectively. Dipole potential drops (DeltaPsi) are less sensitive to long-range LJ interactions than surface tensions. However, DeltaPsi for the DPPC bilayer (845+/-3 mV proceeding from water to lipid) and water (547+/-2 mV for TIP4P-Ew and 521+/-3 mV for TIP3P) overestimate experiment by factors of 3 and 5, respectively, and represent expected deficiencies in nonpolarizable force fields.

Carbene Stabilization by Aryl Substituents. Is Bigger Better?

The geometries and relative stabilities of the singlet and triplet states of phenyl- (Cs), diphenyl- (C2), 1-naphthyl- (Cs), di(1-naphthyl)- (C2), and 9-anthryl-substituted (Cs) carbenes were investigated at the B3LYP/6-311+G(d,p) + ZPVE level of density functional theory. The singlet-triplet energy separations (DeltaEST), 2.7, 2.9, 3.4, 3.7, and 5.7 kcal/mol, respectively, after including an empirical correction (2.8 kcal/mol) based on the error in the computed singlet-triplet gap for methylene versus experiment, are in good agreement with available experimental values. Consistent with literature reports, triplet di(9-anthryl)carbene has a linear, D2d symmetrical, allene structure with 1.336 A C=C bond lengths and considerable biradical character. B3LYP favors such cumulene biradical structures and triplet spin states and predicts a large (>15 kcal/mol) "di(9-anthryl)carbene" singlet-triplet (biradical) energy gap. The resonance stabilization of both singlet and triplet carbenes increases modestly with the size of the arene substituent and overall, (di)arylcarbenes, both singlet and triplet, are better stabilized by bigger substituents. For example, methylene is stabilized more by a naphthyl than a phenyl group (singlets, 26.6 versus 24.4; and triplets, 20.9 versus 18.1 kcal/mol, respectively). The carbene geometries are affected by both steric effects and arene-carbene orbital interactions (sigma-p and p-pi). For instance, the central angles at the carbene are widened by a second arene group, which leads to increased s-character and shorter carbene bond lengths (i.e., C-C, C-H). In general, the aromaticity of the substituted rings in triplet carbenes is most affected by the presence of the unpaired electrons.

Sugar binding in lactose permease: anomeric state of a disaccharide influences binding structure

Lactose permease in Escherichia coli (LacY) transports both anomeric states of disaccharides but has greater affinity for alpha-sugars. Molecular dynamics (MD) simulations are used to probe the protein-sugar interactions, binding structures, and global protein motions in response to sugar binding by investigating LacY (the experimental mutant and wild-type) embedded in a fully hydrated lipid bilayer. A total of 12 MD simulations of 20-25 ns each with beta(alpha)-d-galactopyranosyl-(1,1)-beta-d-galactopyranoside (betabeta-(Galp)(2)) and alphabeta-(Galp)(2) result in binding conformational families that depend on the anomeric state of the sugar. Both sugars strongly interact with Glu126 and alphabeta-(Galp)(2) has a greater affinity to this residue. Binding conformations are also seen that involve protein residues not observed in the crystal structure, as well as those involved in the proton translocation (Phe118, Asn119, Asn240, His322, Glu325, and Tyr350). Common to nearly all protein-sugar structures, water acts as a hydrogen bond bridge between the disaccharide and protein. The average binding energy is more attractive for alphabeta-(Galp)(2) than betabeta-(Galp)(2), i.e. -10.7(+/-0.7) and -3.1(+/-1.0) kcal/mol, respectively. Of the 12 helices in LacY, helix-IV is the least stable with betabeta-(Galp)(2) binding resulting in larger distortion than alphabeta-(Galp)(2).

Engineering Metal Complexes of Chiral Pentaazacrowns as Privileged Reverse-turn Scaffolds

Reverse turns are common structural motifs and recognition sites in protein/protein interactions. The design of peptidomimetics is often based on replacing the amide backbone of peptides by a non-peptidic scaffold while retaining the biologic mode of action. This study evaluates the potential of metal complexes of chiral pentaazacrowns conceptually derived by reduction of cyclic pentapeptides as reverse-turn mimetics. The possible conformations of metal complexes of chiral pentaazacrown scaffolds have been probed by analysis of 28 crystal structures complexed with six different metals (Mn, Fe, Co, Ni, Cu, and Zn). The solvated structures as well as the impact of complexation with different metals/oxidation states have been examined with density functional theory (DFT) calculation as explicitly represented by interactions with a single water molecule. The results suggest that most reverse-turn motifs seen in proteins could be mimicked effectively with a subset of metal complexes of chiral pentaazacrown scaffolds with an RMSD of approximately 0.3 A. Due to the relatively fixed orientation of the pendant chiral side groups in these metal complexes, one can potentially elicit information about the receptor-bound conformation of the parent peptide from their binding affinities. The presence of 20 H-atoms on the pentaazacrown ring that could be functionalized as well as the conformational perturbations available from complexation with different metals offer a desirable diversity to probe receptors for reverse-turn recognition.

Dynamical motions of lipids and a finite size effect in simulations of bilayers

Molecular dynamics (MD) simulations of dipalmitoylphosphatidylcholine bilayers composed of 72 and 288 lipids are used to examine system size dependence on dynamical properties associated with the particle mesh Ewald (PME) treatment of electrostatic interactions. The lateral diffusion constant Dl is 2.92 x 10(-7) and 0.95 x 10(-7) cm2/s for 72 and 288 lipids, respectively. This dramatic finite size effect originates from the correlation length of lipid diffusion, which extends to next-nearest neighbors in the 288 lipid system. Consequently, diffusional events in smaller systems can propagate across the boundaries of the periodic box. The internal dynamics of lipids calculated from the PME simulations are independent of the system size. Specifically, reorientational correlation functions for the slowly relaxing phosphorus-glycerol hydrogen, phosphorus-nitrogen vectors, and more rapidly relaxing CH vectors in the aliphatic chains are equivalent for the 72 and 288 lipid simulations. A third MD simulation of a bilayer with 72 lipids using spherical force-shift electrostatic cutoffs resulted in interdigitated chains, thereby rendering this cutoff method inappropriate.

Structure of the Boronic Acid Dimer and the Relative Stabilities of Its Conformers

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order Møller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.

Isotropic periodic sum: a method for the calculation of long-range interactions

This work presents an accurate and efficient approach to the calculation of long-range interactions for molecular modeling and simulation. This method defines a local region for each particle and describes the remaining region as images of the local region statistically distributed in an isotropic and periodic way, which we call isotropic periodic images. Different from lattice sum methods that sum over discrete lattice images generated by periodic boundary conditions, this method sums over the isotropic periodic images to calculate long-range interactions, and is referred to as the isotropic periodic sum (IPS) method. The IPS method is not a lattice sum method and eliminates the need for a reciprocal space sum. Several analytic solutions of IPS for commonly used potentials are presented. It is demonstrated that the IPS method produces results very similar to that of Ewald summation, but with three major advantages, (1) it eliminates unwanted symmetry artifacts raised from periodic boundary conditions, (2) it can be applied to potentials of any functional form and to fully and partially homogenous systems as well as finite systems, and (3) it is more computationally efficient and can be easily parallelized for multiprocessor computers. Therefore, this method provides a general approach to an efficient calculation of long-range interactions for various kinds of molecular systems.

Low-frequency normal modes that describe allosteric transitions in biological nanomachines are robust to sequence variations

By representing the high-resolution crystal structures of a number of enzymes using the elastic network model, it has been shown that only a few low-frequency normal modes are needed to describe the large-scale domain movements that are triggered by ligand binding. Here we explore a link between the nearly invariant nature of the modes that describe functional dynamics at the mesoscopic level and the large evolutionary sequence variations at the residue level. By using a structural perturbation method (SPM), which probes the residue-specific response to perturbations (or mutations), we identify a sparse network of strongly conserved residues that transmit allosteric signals in three structurally unrelated biological nanomachines, namely, DNA polymerase, myosin motor, and the Escherichia coli chaperonin. Based on the response of every mode to perturbations, which are generated by interchanging specific sequence pairs in a multiple sequence alignment, we show that the functionally relevant low-frequency modes are most robust to sequence variations. Our work shows that robustness of dynamical modes at the mesoscopic level is encoded in the structure through a sparse network of residues that transmit allosteric signals.

Development of small molecules designed to modulate protein-protein interactions

Protein-protein interactions are ubiquitous, essential to almost all known biological processes, and offer attractive opportunities for therapeutic intervention. Developing small molecules that modulate protein-protein interactions is challenging, owing to the large size of protein-complex interface, the lack of well-defined binding pockets, etc. We describe a general approach based on the "privileged-structure hypothesis" [Che, Ph.D. Thesis, Washington University, 2003] - that any organic templates capable of mimicking surfaces of protein-recognition motifs are potential privileged scaffolds as protein-complex antagonists--to address the challenges inherent in the discovery of small-molecule inhibitors of protein-protein interactions.

Simulation-based methods for interpreting x-ray data from lipid bilayers

The fully hydrated liquid crystalline phase of the dimyristoylphosphatidycholine lipid bilayer at 30 degrees C was simulated using molecular dynamics with the CHARMM potential for five surface areas per lipid (A) in the range 55-65 A(2) that brackets the previously determined experimental area 60.6 A(2). The results of these simulations are used to develop a new hybrid zero-baseline structural model, denoted H2, for the electron density profile, rho(z), for the purpose of interpreting x-ray diffraction data. H2 and also the older hybrid baseline model were tested by fitting to partial information from the simulation and various constraints, both of which correspond to those available experimentally. The A, rho(z), and F(q) obtained from the models agree with those calculated directly from simulation at each of the five areas, thereby validating this use of the models. The new H2 was then applied to experimental dimyristoylphosphatidycholine data; it yields A = 60.6 +/- 0.5 A(2), in agreement with the earlier estimate obtained using the hybrid baseline model. The electron density profiles also compare well, despite considerable differences in the functional forms of the two models. Overall, the simulated rho(z) at A = 60.7 A(2) agrees well with experiment, demonstrating the accuracy of the CHARMM lipid force field; small discrepancies indicate targets for improvements. Lastly, a simulation-based model-free approach for obtaining A is proposed. It is based on interpolating the area that minimizes the difference between the experimental F(q) and simulated F(q) evaluated for a range of surface areas. This approach is independent of structural models and could be used to determine structural properties of bilayers with different lipids, cholesterol, and peptides.

Subcortical reorganization in amyotrophic lateral sclerosis

The cerebral cortex reorganizes in response to central or peripheral lesions. Although basal ganglia and cerebellum are key components of the network dedicated to movement control, their role in motor reorganization remains elusive. We therefore tested if slowly progressive neurodegenerative motor disease alters the subcortical functional anatomy of the basal ganglia-thalamo-cerebellar circuitry. Ten patients with amyotrophic lateral sclerosis (ALS) and ten healthy controls underwent functional magnetic resonance imaging (fMRI), while executing a simple finger flexion task. Cued by an acoustic trigger, they squeezed a handgrip force transducer with their right hand at 10% of their maximum voluntary contraction force. Movement frequency, amplitude, and force were controlled. Statistical parametric mapping of task-related BOLD-response revealed increased activation in ALS patients as compared to healthy controls. The main activation increases were found in the supplementary motor area, basal ganglia, brainstem, and cerebellum. These findings suggest that degeneration of cortical and spinal motor neurons in ALS leads to a recruitment of subcortical motor structures. These subcortical activation patterns strongly resemble functional activation in motor learning and might therefore represent adaptations of cortico-subcortical motor loops as a - albeit finally ineffective - mechanism to compensate for the ongoing loss of motor neurons in ALS.

Modeling protein conformational changes by iterative fitting of distance constraints using reoriented normal modes

Recently we have developed a normal-modes-based algorithm that predicts the direction of protein conformational changes given the initial state crystal structure together with a small number of pairwise distance constraints for the end state. Here we significantly extend this method to accurately model both the direction and amplitude of protein conformational changes. The new protocol implements a multisteps search in the conformational space that is driven by iteratively minimizing the error of fitting the given distance constraints and simultaneously enforcing the restraint of low elastic energy. At each step, an incremental structural displacement is computed as a linear combination of the lowest 10 normal modes derived from an elastic network model, whose eigenvectors are reorientated to correct for the distortions caused by the structural displacements in the previous steps. We test this method on a list of 16 pairs of protein structures for which relatively large conformational changes are observed (root mean square deviation >3 angstroms), using up to 10 pairwise distance constraints selected by a fluctuation analysis of the initial state structures. This method has achieved a near-optimal performance in almost all cases, and in many cases the final structural models lie within root mean square deviation of 1 approximately 2 angstroms from the native end state structures.

Residues in substrate proteins that interact with GroEL in the capture process are buried in the native state

We have used a bioinformatic approach to predict the natural substrate proteins for the Escherichia coli chaperonin GroEL based on two simple criteria. Natural substrate proteins should contain binding motifs similar in sequence to the mobile loop peptide of GroES that displaces the binding motif during the chaperonin cycle. Secondly, each substrate protein should contain multiple copies of the binding motif so that the chaperonin can perform "work" on the substrate protein. To validate these criteria, we have used a database of 252 proteins that have been experimentally shown to interact with the chaperonin machinery in vivo. More than 80% are identified by these criteria. The binding motifs of all 79 proteins in the database with a known three-dimensional structure are buried (<50% solvent-accessible surface area) in the native state. Our results show that the binding motifs are inaccessible in the native state but become solvent-exposed in unfolded state, thus enabling GroEL to distinguish between unfolded and native states. The structures of the binding motif in the native states of the substrate proteins include alpha-helices, beta-strands, and random coils. The diversity of secondary structures implies that there are large and varied conformational transitions in the recognition motifs after their displacement by the mobile loops of GroES.

Molecular structure of a 9-MDa icosahedral pyruvate dehydrogenase subcomplex containing the E2 and E3 enzymes using cryoelectron microscopy

The pyruvate dehydrogenase multienzyme complexes are among the largest multifunctional catalytic machines in cells, catalyzing the production of acetyl CoA from pyruvate. We have previously reported the molecular architecture of an 11-MDa subcomplex comprising the 60-mer icosahedral dihydrolipoyl acetyltransferase (E2) decorated with 60 copies of the heterotetrameric (alpha(2)beta(2)) 153-kDa pyruvate decarboxylase (E1) from Bacillus stearothermophilus (Milne, J. L. S., Shi, D., Rosenthal, P. B., Sunshine, J. S., Domingo, G. J., Wu, X., Brooks, B. R., Perham, R. N., Henderson, R., and Subramaniam, S. (2002) EMBO J. 21, 5587-5598). An annular gap of approximately 90 A separates the acetyltransferase catalytic domains of the E2 from an outer shell formed of E1 tetramers. Using cryoelectron microscopy, we present here a three-dimensional reconstruction of the E2 core decorated with 60 copies of the homodimeric 100-kDa dihydrolipoyl dehydrogenase (E3). The E2E3 complex has a similar annular gap of approximately 75 A between the inner icosahedral assembly of acetyltransferase domains and the outer shell of E3 homodimers. Automated fitting of the E3 coordinates into the map suggests excellent correspondence between the density of the outer shell map and the positions of the two best fitting orientations of E3. As in the case of E1 in the E1E2 complex, the central 2-fold axis of the E3 homodimer is roughly oriented along the periphery of the shell, making the active sites of the enzyme accessible from the annular gap between the E2 core and the outer shell. The similarities in architecture of the E1E2 and E2E3 complexes indicate fundamental similarities in the mechanism of active site coupling involved in the two key stages requiring motion of the swinging lipoyl domain across the annular gap, namely the synthesis of acetyl CoA and regeneration of the dithiolane ring of the lipoyl domain.