Emanuele Paci

CHARMM: The biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model‐building capabilities. The CHARMM program is applicable to problems involving a much broader class of many‐particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical‐molecular mechanical force fields, to all‐atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

The major determinant of exendin‐4/glucagon‐like peptide 1 differential affinity at the rat glucagon‐like peptide 1 receptor N‐terminal domain is a hydrogen bond from SER‐32 of exendin‐4

Exendin‐4 (exenatide, Ex4) is a high‐affinity peptide agonist at the glucagon‐like peptide‐1 receptor (GLP‐1R), which has been approved as a treatment for type 2 diabetes. Part of the drug/hormone binding site was described in the crystal structures of both GLP‐1 and Ex4 bound to the isolated N‐terminal domain (NTD) of GLP‐1R. However, these structures do not account for the large difference in affinity between GLP‐1 and Ex4 at this isolated domain, or for the published role of the C‐terminal extension of Ex4. Our aim was to clarify the pharmacology of GLP‐1R in the context of these new structural data.

Three key residues form a critical contact network in a protein folding transition state

Determining how a protein folds is a central problem in structural biology. The rate of folding of many proteins is determined by the transition state, so that a knowledge of its structure is essential for understanding the protein folding reaction. Here we use mutation measurements—which determine the role of individual residues in stabilizing the transition state—as restraints in a Monte Carlo sampling procedure to determine the ensemble of structures that make up the transition state. We apply this approach to the experimental data for the 98-residue protein acylphosphatase, and obtain a transition-state ensemble with the native-state topology and an average root-mean-square deviation of 6 Å from the native structure. Although about 20 residues with small positional fluctuations form the structural core of this transition state, the native-like contact network of only three of these residues is sufficient to determine the overall fold of the protein. This result reveals how a nucleation mechanism involving a small number of key residues can lead to folding of a polypeptide chain to its unique native-state structure.

Pulling geometry defines the mechanical resistance of a β-sheet protein

Proteins show diverse responses when placed under mechanical stress. The molecular origins of their differing mechanical resistance are still unclear, although the orientation of secondary structural elements relative to the applied force vector is thought to have an important function. Here, by using a method of protein immobilization that allows force to be applied to the same all-β protein, E2lip3, in two different directions, we show that the energy landscape for mechanical unfolding is markedly anisotropic. These results, in combination with molecular dynamics (MD) simulations, reveal that the unfolding pathway depends on the pulling geometry and is associated with unfolding forces that differ by an order of magnitude. Thus, the mechanical resistance of a protein is not dictated solely by amino acid sequence, topology or unfolding rate constant, but depends critically on the direction of the applied extension.

Mechanical Unfolding of a Titin Ig Domain: Structure of Unfolding Intermediate Revealed by Combining AFM, Molecular Dynamics Simulations, NMR and Protein Engineering

The mechanical unfolding of an immunoglobulin domain from the human muscle protein titin (TI I27) has been shown to proceed via a metastable intermediate in which the A-strand is detached. The structure and properties of this intermediate are characterised in this study. A conservative destabilising mutation in the A-strand has no effect on the unfolding force, nor the dependence of the unfolding force on the pulling speed, indicating that the unfolding forces measured in an AFM experiment are those required for the unfolding of the intermediate and not the native state. A mutant of TI I27 with the A-strand deleted (TI I27-A) is studied by NMR and standard biophysical techniques, combined with protein engineering. Molecular dynamics simulations show TI I27-A to be a good model for the intermediate. It has a structure very similar to the native state, and is surprisingly stable. Comparison with a Phi-value analysis of the unfolding pathway clearly shows that the protein unfolds by a different pathway under an applied force than on addition of denaturant.

ORAC: A Molecular dynamics program to simulate complex molecular systems with realistic electrostatic interactions

In this study, we present a new molecular dynamics program for simulation of complex molecular systems. The program, named ORAC, combines state‐of‐the‐art molecular dynamics (MD) algorithms with flexibility in handling different types and sizes of molecules. ORAC is intended for simulations of molecular systems and is specifically designed to treat biomolecules efficiently and effectively in solution or in a crystalline environment. Among its unique features are: (i) implementation of reversible and symplectic multiple time step algorithms (or r‐RESPA, reversible reference system propagation algorithm) specifically designed and tuned for biological systems with periodic boundary conditions; (ii) availability for simulations with multiple or single time steps of standard Ewald or smooth particle mesh Ewald (SPME) for computation of electrostatic interactions; and (iii) possibility of simulating molecular systems in a variety of thermodynamic ensembles. We believe that the combination of these algorithms makes ORAC more advanced than other MD programs using standard simulation algorithms. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1848–1862, 1997

Transition states for protein folding have native topologies despite high structural variability

We present a structural analysis of the folding transition states of three SH3 domains. Our results reveal that the secondary structure is not yet fully formed at this stage of folding and that the solvent is only partially excluded from the interior of the protein. Comparison of the members of the transition state ensemble with a database of native folds shows that, despite substantial local variability, the transition state structures can all be classified as having the topology characteristic of an SH3 domain. Our results suggest a mechanism for folding in which the formation of a network of interactions among a subset of hydrophobic residues ensures that the native topology is generated. Such a mechanism enables high fidelity in folding while minimizing the need to establish a large number of specific interactions in the conformational search.

Structures and relative free energies of partially folded states of proteins

The ability of proteins to fold to well defined compact structures is one of the most remarkable examples of the effect of natural selection on biological molecules. To understand their properties, including the stability, the mechanism of folding, and the possibilities of misfolding and association, it is necessary to know the protein free energy landscape. We use NMR data as restraints in a Monte Carlo sampling procedure to determine the ensemble of structures populated by human α-lactalbumin in the presence of increasing concentrations of urea. The ensembles of structures that represent the partially folded states of the protein show that two structural cores, corresponding to portions of the α and β domains of the native protein, are preserved even when the native-like interactions that define their existence are substantially weakened. Analysis of the network of residual contacts reveals the presence of a complex interface region between the two structural cores and indicates that the development of specific interactions within this interface is the key step in achieving the native structure. The relative probabilities of the conformations determined from the NMR data are used to construct a coarse-grained free energy landscape for α-lactalbumin in the absence of urea. The form of the landscape, together with the existence of distinct cores, supports the concept that robustness and modularity are the properties that make possible the folding of complex proteins.

Native and non‐native interactions along protein folding and unfolding pathways

Gō‐type models, which include only native contact interactions in the energy function, are being used increasingly to describe the protein folding reaction. To investigate the validity of such models, we determine the role of native and non‐native interactions along folding and unfolding pathways. For this purpose, we use a molecular mechanics energy function with an implicit solvation model (an effective energy function or potential of mean force) that can be expressed in a pairwise decomposable form. We find that for the native state and a wide range of other configurations, the contact energy is an accurate description, in part due to the cancellation of non‐zero contributions from more distant residues. However, significant errors in the energy are introduced for non‐native structures if the energy is calculated from the native contacts alone. Non‐native contacts tend to make a significant contribution, particularly for molten globules and collapsed states along the unfolding pathways. The implication of these results for the use of Gō‐type models in studies of protein folding are discussed. Proteins 2002;47:379–392.

Unraveling the molecular basis of subunit specificity in P pilus assembly by mass spectrometry

P pili are multisubunit fibers essential for the attachment of uropathogenic Escherichia coli to the kidney. These fibers are formed by the noncovalent assembly of six different homologous subunit types in an array that is strictly defined in terms of both the number and order of each subunit type. Assembly occurs through a mechanism termed “donor-strand exchange (DSE)” in which an N-terminal extension (Nte) of one subunit donates a β-strand to an adjacent subunit, completing its Ig fold. Despite structural determination of the different subunits, the mechanism determining specificity of subunit ordering in pilus assembly remained unclear. Here, we have used noncovalent mass spectrometry to monitor DSE between all 30 possible pairs of P pilus subunits and their Ntes. We demonstrate a striking correlation between the natural order of subunits in pili and their ability to undergo DSE in vitro. The results reveal insights into the molecular mechanism by which subunit ordering during the assembly of this complex is achieved.