Carol Beth Post

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.

Energetics of Quasiequivalence: Computational Analysis of Protein-Protein Interactions in Icosahedral Viruses

Quaternary structure polymorphism found in quasiequivalent virus capsids provides a static framework for studying the dynamics of protein interactions. The same protein subunits are found in different structural environments within these particles, and in some cases, the molecular switching required for the polymorphic quaternary interactions is obvious from high-resolution crystallographic studies. Employing atomic resolution structures, molecular mechanics, and continuum electrostatic methods, we have computed association energies for unique subunit interfaces of three icosahedral viruses, black beetle virus, southern bean virus, and human rhinovirus 14. To quantify the chemical determinants of quasiequivalence, the energetic contributions of individual residues forming quasiequivalent interfaces were calculated and compared. The potential significance of the differences in stabilities at quasiequivalent interfaces was then explored with the combinatorial assembly approach. The analysis shows that the unique association energies computed for each virus serve as a sensitive basis set that may determine distinct intermediates and pathways of virus capsid assembly. The pathways for the quasiequivalent viruses displayed isoenergetic oligomers at specific points, suggesting that these may determine the quaternary structure polymorphism required for the assembly of a quasiequivalent particle.

Flavivirus Capsid Is a Dimeric Alpha-Helical Protein

ABSTRACT The capsid proteins of two flaviviruses, yellow fever virus and dengue virus, were expressed in Escherichia coli and purified to near homogeneity suitable for biochemical characterization and structure determination by nuclear magnetic resonance. The oligomeric properties of the capsid protein in solution were investigated. In the absence of nucleic acid, both proteins were predominately dimeric in solution. Further analysis of both proteins with far-UV circular dichroism spectroscopy indicated that they were largely alpha-helical. The secondary structure elements of the dengue virus capsid were determined by chemical shift indexing of the sequence-specific backbone resonance assignments. The dengue virus capsid protein devoid of its C-terminal signal sequence was found to be composed of four alpha helices. The longest alpha helix, 20 residues, is located at the C terminus and has an amphipathic character. In contrast, the N terminus was found to be unstructured and could be removed without disrupting the structural integrity of the protein.

Structure and self-association of the Rous sarcoma virus capsid protein.

BACKGROUND The capsid protein (CA) of retroviruses, such as Rous sarcoma virus (RSV), consists of two independently folded domains. CA functions as part of a polyprotein during particle assembly and budding and, in addition, forms a shell encapsidating the genomic RNA in the mature, infectious virus. RESULTS The structures of the N- and C-terminal domains of RSV CA have been determined by X-ray crystallography and solution nuclear magnetic resonance (NMR) spectroscopy, respectively. The N-terminal domain comprises seven alpha helices and a short beta hairpin at the N terminus. The N-terminal domain associates through a small, tightly packed, twofold symmetric interface within the crystal, different from those previously described for other retroviral CAs. The C-terminal domain is a compact bundle of four alpha helices, although the last few residues are disordered. In dilute solution, RSV CA is predominantly monomeric. We show, however, using electron microscopy, that intact RSV CA can assemble in vitro to form both tubular structures constructed from toroidal oligomers and planar monolayers. Both modes of assembly occur under similar solution conditions, and both sheets and tubes exhibit long-range order. CONCLUSIONS The tertiary structure of CA is conserved across the major retroviral genera, yet sequence variations are sufficient to cause change in associative behavior. CA forms the exterior shell of the viral core in all mature retroviruses. However, the core morphology differs between viruses. Consistent with this observation, we find that the capsid proteins of RSV and human immunodeficiency virus type 1 exhibit different associative behavior in dilute solution and assemble in vitro into different structures.

Molecular Basis for a Direct Interaction between the Syk Protein-tyrosine Kinase and Phosphoinositide 3-Kinase

After engagement of the B cell receptor for antigen, the Syk protein-tyrosine kinase becomes phosphorylated on multiple tyrosines, some of which serve as docking sites for downstream effectors with SH2 or other phosphotyrosine binding domains. The most frequently identified binding partner for catalytically active Syk identified in a yeast two-hybrid screen was the p85 regulatory subunit of phosphoinositide 3-kinase. The C-terminal SH2 domain of p85 was sufficient for mediating an interaction with tyrosine-phosphorylated Syk. Interestingly, this domain interacted with Syk at phosphotyrosine 317, a site phosphorylated in trans by the Src family kinase, Lyn, and identified previously as a binding site for c-Cbl. This site interacted preferentially with the p85 C-terminal SH2 domain compared with the c-Cbl tyrosine kinase binding domain. Molecular modeling studies showed a good fit between the p85 SH2 domain and a peptide containing phosphotyrosine 317. Tyr-317 was found to be essential for Syk to support phagocytosis mediated by FcγRIIA receptors expressed in a heterologous system. These studies establish a new type of p85 binding site that can exist on proteins that serve as substrates for Src family kinases and provide a molecular explanation for observations on direct interactions between Syk and phosphoinositide 3-kinase.

Molecular dynamics simulations of native and substrate-bound lysozyme: A study of the average structures and atomic fluctuations

Molecular dynamics simulations of hen egg-white lysozyme in the free and substrate-bound states are reported and the nature of the average structures and atomic fluctuations are analyzed. Crystallographic water molecules of structural importance, as determined by hydrogen-bonding, were included in the simulations. Comparisons are made between the dynamics and the X-ray results for the atomic positions, the main-chain and side-chain dihedral angles, and the hydrogen-bonding geometry. Improvements over earlier simulations in the potential energy function and methodology resulted in stable trajectories with the C alpha co-ordinates within 1.5 A of the starting X-ray structure. Structural features analyzed in the simulations agreed well with the X-ray results except for some surface residues. The Asx chi 2 dihedral distribution and the geometry of hydrogen bonding at reverse turns show differences; possible causes are discussed. The relation between the magnitudes and time-scales of the residue fluctuations and secondary structural features, such as helices beta-sheets and coiled loops, is examined. Significant differences in the residue mobilities between the simulations of the free and substrate-bound states were found in a region of the enzyme that is in direct contact with the substrate and in a region that is distant from the active-site cleft. The dynamic behavior of the structural water molecules is analyzed by examining the correlation between the fluctuations of the water oxygens and the lysozyme heavy-atoms to which they are hydrogen-bonded.

Design, Synthesis, and Preclinical Evaluation of Prostate-Specific Membrane Antigen Targeted 99mTc-Radioimaging Agents

The high mortality and financial burden associated with prostate cancer can be partly attributed to a lack of sensitive screening methods for detection and staging of the disease. Guided by in silico docking studies using the crystal structure of PSMA, we designed and synthesized a series of PSMA-targeted (99m)Tc-chelate complexes for imaging PSMA-expressing human prostate cancer cells (LNCaP cell line). Of the six targeted radioimaging agents synthesized, three were found to bind LNCaP cells with low nanomolar affinity. Moreover, the same three PSMA-targeted imaging agents were shown to localize primarily to LNCaP tumor xenografts in nu/nu mice, with an average of 9.8 +/- 2.4% injected dose/g tissue accumulating in the tumor and only 0.11% injected dose/g tissue retained in the muscle at 4 h postinjection. Collectively, these high affinity, PSMA-specific radioimaging agents demonstrate significant potential for use in localizing prostate cancer masses, monitoring response to therapy, detecting prostate cancer recurrence following surgery, and selecting patients for subsequent PSMA-targeted chemotherapy.

Design, Synthesis, and Biological Evaluation of Antiviral Agents Targeting Flavivirus Envelope Proteins

Flavivirus envelope proteins (E proteins) have been shown to play a pivotal role in virus assembly, morphogenesis, and infection of host cells. Inhibition of flavivirus infection of a host cell by means of a small molecule envelope protein antagonist is an attractive strategy for the development of antiviral agents. Virtual screening of the NCI chemical database using the dengue virus envelope protein structure revealed several hypothetical hit compounds. Bioassay results identified a class of thiazole compounds with antiviral potency in cell-based assays. Modification of these lead compounds led to a series of analogues with improved antiviral activity and decreased cytotoxicity. The most active compounds 11 and 36 were effective in the low micromolar concentration range in a cellular assay system.

Antiviral Compounds Discovered by Virtual Screening of Small–Molecule Libraries against Dengue Virus E Protein

Infection by the mosquito-borne dengue virus causes dengue fever and the sometimes fatal dengue hemorrhagic fever. The increasing number of dengue infections per year suggests that the virus is becoming more virulent and its transmission is expanding. Nevertheless, no effective treatment for dengue infection currently exists. In a search for antiviral agents effective against dengue virus, we investigated the potential of targeting a structural protein site rather than an enzymatic one. Using this approach, we now report the discovery of a small molecule ligand that inhibits viral growth. Our results also provide the first evidence that the binding site, a pocket located at the hinge between domains 1 and 2 of the envelope protein (E protein) on the virus surface, is a valid target for antiviral therapy. Ligand candidates were identified from libraries of approximately 142,000 compounds using a computational high-throughput screening protocol targeting this pocket of the E protein. Cell-based assays were conducted on 23 top-ranked compounds. Among four with good antiviral activity profiles, the compound P02 was found to inhibit viral reproduction at micromolar concentrations. Using saturation transfer difference NMR spectroscopy, we also show that the compound binds virus and competes for binding E protein with the known ligand N-octyl-beta-D-glucoside. Together, the results are consistent with an inhibition mechanism against maturation or host-cell entry mediated by ligand binding to the E-protein pocket. P02 is a promising lead compound for future development of an effective treatment against dengue virus and related flaviviruses.

Src kinase activation: A switched electrostatic network

Src tyrosine kinases are essential in numerous cell signaling pathways, and improper functioning of these enzymes has been implicated in many diseases. The activity of Src kinases is regulated by conformational activation, which involves several structural changes within the catalytic domain (CD): the orientation of two lobes of CD; rearrangement of the activation loop (A‐loop); and movement of an α‐helix (αC), which is located at the interface between the two lobes, into or away from the catalytic cleft. Conformational activation was investigated using biased molecular dynamics to explore the transition pathway between the active and the down‐regulated conformation of CD for the Src‐kinase family member Lyn kinase, and to gain insight into the interdependence of these changes. Lobe opening is observed to be a facile motion, whereas movement of the A‐loop motion is more complex requiring secondary structure changes as well as communication with αC. A key result is that the conformational transition involves a switch in an electrostatic network of six polar residues between the active and the down‐regulated conformations. The exchange between interactions links the three main motions of the CD. Kinetic experiments that would demonstrate the contribution of the switched electrostatic network to the enzyme mechanism are proposed. Possible implications for regulation conferred by interdomain interactions are also discussed.