Molecular Modeling of Proteins: A Strategy for Energy Minimization by Molecular Mechanics in the AMBER Force Field

Abstract

Energy minimization is an important step in molecular modeling of proteins. In this study, we sought to develop a minimization strategy which would give the best final structures with the shortest computer time in the AMBER force field. In the all-atom model, we performed energy minimization of the melittin (mostly α-helical) and cardiotoxin (mostly β-sheet and β-turns) crystal structures by both constrained and unconstrained pathways. In the constrained path, which has been recommended in the energy minimization of proteins, hydrogens were relaxed first, followed by the side chains of amino acid residues, and finally the whole molecule. Despite the logic of this approach, however, the structures minimized by the unconstrained path fit the experimental structures better than those minimized by constrained paths. Moreover, the unconstrained path saved considerable computer time. We also compared the effects of the steepest descents and conjugate gradients algorithms in energy minimization. Previously, steepest descents has been used in the initial stages of minimization and conjugate gradients in the final stages of minimization. We therefore studied the effect on the final structure of performing an initial minimization by steepest descents. The structures minimized by conjugate gradients alone resembled the structures minimized initially by the steepest descents and subsequently by the conjugate gradients algorithms. Thus an initial minimization using steepest descents is wasteful and unnecessary, especially when starting from the crystal structure. Based on these results, we propose the use of an unconstrained path and conjugate gradients for energy minimization of proteins. This procedure results in low energy structures closer to the experimental structures, and saves about 70–80% of computer time. This procedure was applied in building models of lysozyme mutants. The crystal structure of native T₄ lysozyme was mutated to three different mutants and the structures were minimized. The minimized structures closely fit the crystal structures of the respective mutants (< 0.3 Å root-mean-square, RMS, deviation in the position of all heavy atoms). These results confirm the efficiency of the proposed minimization strategy in modeling closely related homologs. To determine the reliability of the united atom approximation, we also performed all of the above minimizations with united atom models. This approximation gave structures with similar but slightly higher RMS deviations than the all-atom model, but gave further savings of60-70% in computer time. However, we feel further investigation is essential to determine the reliability of this approximation. Finally, to determine the limitation of the procedure, we built the melittin molecule interactively in an α-helical conformation and this model showed an RMS deviation greater than 2.8 Å when compared to the melittin crystal structure. This model was minimized by various strategies. None of the minimized structures converged towards the crystal structure. Thus, although the proposed method seems to give valid structures starting from closely related crystal structures, it cannot predict the native structure when the starting structure is far from the native structure. From these results, we recommend the use of the proposed strategy of minimizing by an unconstrained path using the conjugate gradients algorithm, but only for modeling of closely related structural homologs of proteins.     

Critical Parameters for Efficient Sonication and Improved Chromatin Immunoprecipitation of High Molecular Weight Proteins

Solubilization of cross-linked cells followed by chromatin shearing is essential for successful chromatin immunoprecipitation (ChIP). However, this task, typically accomplished by ultrasound treatment, may often become a pitfall of the process, due to inconsistent results obtained between different experiments under seemingly identical conditions. To address this issue we systematically studied ultrasound-mediated cell lysis and chromatin shearing, identified critical parameters of the process and formulated a generic strategy for rational optimization of ultrasound treatment. We also demonstrated that whereas ultrasound treatment required to shear chromatin to within a range of 100–400 bp typically degrades large proteins, a combination of brief sonication and benzonase digestion allows for the generation of similarly sized chromatin fragments while preserving the integrity of associated proteins. This approach should drastically improve ChIP efficiency for this class of proteins.     

Q-FADD: A Mechanistic Approach for Modeling the Accumulation of Proteins at Sites of DNA Damage

The repair of DNA damage requires the ordered recruitment of many different proteins that are responsible for signaling and subsequent repair. A powerful and widely used tool for studying the orchestrated accumulation of these proteins at damage sites is laser microirradiation in live cells, followed by monitoring the accumulation of the fluorescently labeled protein in question. Despite the widespread use of this approach, there exists no rigorous method for characterizing the recruitment process quantitatively. Here, we introduce a diffusion model that explicitly accounts for the unique sizes and shapes of individual nuclei and uses two variables: D_eff, the effective coefficient of diffusion, and F, the fraction of mobile protein that accumulates at sites of DNA damage. Our model quantitatively describes the accumulation of three test proteins, poly-ADP-ribose polymerases 1 and 2 (PARP1/2) and histone PARylation factor 1. D_eff for PARP1, as derived by our approach, is 6× greater than for PARP2 and in agreement with previous literature reports using fluorescence correlation spectroscopy and fluorescence recovery after photobleaching. Our data indicate that histone PARylation factor 1 arrives at sites of DNA damage independently of either PARP. Importantly, our model, which can be applied to existing data, allows for the direct comparison of the coefficient of diffusion for any DNA repair protein between different cell types, obtained in different laboratories and by different methods, and also allows for the interrogation of cell-to-cell variability.     

Identification of Fibrinogen-Binding Proteins of Aspergillus fumigatus Using Proteomic Approach

Aspergillus fumigatus, the main etiological agent for various forms of human aspergillosis, gets access to the respiratory system of human host by inhalation of airborne conidia. These conidia possibly adhere to extracellular matrix (ECM) proteins. Among the ECM proteins involved in adherence, fibrinogen is thought to be crucial. Here, we studied whether A. fumigatus three-week culture filtrate (3wcf) proteins promote binding of A. fumigatus to ECM proteins and promote fungal growth. We observed that incubation of ECM with 3wcf proteins led to dose- and time-dependent increase in adherence of conidia to the ECM. In order to identify the catalogue of fibrinogen-binding A. fumigatus proteins, we carried out fibrinogen affinity blotting using two-dimensional gel electrophoresed 3wcf proteins. A total of 15 fibrinogen-binding protein spots corresponding to 7 unique proteins were identified in 3wcf using matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF-TOF). Among these, 4 proteins, namely, beta-glucosidase, alpha-mannosidase, pectate lyase A and oryzin precursor were predicted to have cell wall or extracellular localization, whereas amidase family protein and two hypothetical proteins did not display the signal sequence. This study reports seven novel fibrinogen-binding proteins of A. fumigatus, some of which could be further explored for targeting the adhesion phenomenon as antifungal strategy.     

Stability Curve Prediction of Homologous Proteins Using Temperature-Dependent Statistical Potentials

The unraveling and control of protein stability at different temperatures is a fundamental problem in biophysics that is substantially far from being quantitatively and accurately solved, as it requires a precise knowledge of the temperature dependence of amino acid interactions. In this paper we attempt to gain insight into the thermal stability of proteins by designing a tool to predict the full stability curve as a function of the temperature for a set of 45 proteins belonging to 11 homologous families, given their sequence and structure, as well as the melting temperature () and the change in heat capacity () of proteins belonging to the same family. Stability curves constitute a fundamental instrument to analyze in detail the thermal stability and its relation to the thermodynamic stability, and to estimate the enthalpic and entropic contributions to the folding free energy. In summary, our approach for predicting the protein stability curves relies on temperature-dependent statistical potentials derived from three datasets of protein structures with targeted thermal stability properties. Using these potentials, the folding free energies () at three different temperatures were computed for each protein. The Gibbs-Helmholtz equation was then used to predict the protein''s stability curve as the curve that best fits these three points. The results are quite encouraging: the standard deviations between the experimental and predicted ''s, ''s and folding free energies at room temperature () are equal to 13 , 1.3 ) and 4.1 , respectively, in cross-validation. The main sources of error and some further improvements and perspectives are briefly discussed.     

Dissecting the Molecular Properties of Prokaryotic Flotillins

Flotillins are universally conserved proteins that are present in all kingdoms of life. Recently it was demonstrated that the B. subtilis flotillin YuaG (FloT) has a direct influence on membrane domain formation by orchestrating lipid domains. Thereby it allocates a proper environment for diverse cellular machineries. YuaG creates platforms for signal transduction, processes crucial for biofilm formation, sporulation, competence, secretion, and others. Even though, flotillins are an emerging topic of research in the field of microbiology little is known about the molecular architecture of prokaryotic flotillins. All flotillins share common structural elements and are tethered to the membrane N’- terminally, followed by a so called PHB domain and a flotillin domain. We show here that prokaryotic flotillins are, similarly to eukaryotic flotillins, tethered to the membrane via a hairpin loop. Further it is demonstrated by sedimentation assays that B. subtilis flotillins do not bind to the membrane via their PHB domain contrary to eukaryotic flotillins. Size exclusion chromatography experiments, blue native PAGE and cross linking experiments revealed that B. subtilis YuaG can oligomerize into large clusters via the PHB domain. This illustrates an important difference in the setup of prokaryotic flotillins compared to the organization of eukaryotic flotillins.     

Pharmacophore Modeling for Anti-Chagas Drug Design Using the Fragment Molecular Orbital Method

BackgroundChagas disease, caused by the parasite Trypanosoma cruzi, is a neglected tropical disease that causes severe human health problems. To develop a new chemotherapeutic agent for the treatment of Chagas disease, we predicted a pharmacophore model for T. cruzi dihydroorotate dehydrogenase (TcDHODH) by fragment molecular orbital (FMO) calculation for orotate, oxonate, and 43 orotate derivatives.Conclusions/SignificanceFMO-based interaction energy analyses revealed a pharmacophore model for TcDHODH inhibitor. Hydrogen bond acceptor pharmacophores correspond to Lys43 and Lys214, hydrogen bond donor and acceptor pharmacophores correspond to Asn67 and Asn194, and the aromatic ring pharmacophore corresponds to FMN, which shows important characteristics of compounds that inhibit TcDHODH. In addition, the Lys214 residue is not conserved between TcDHODH and human DHODH. Our analysis suggests that these orotate derivatives should preferentially bind to TcDHODH, increasing their selectivity. Our results obtained by pharmacophore modeling provides insight into the structural requirements for the design of TcDHODH inhibitors and their development as new anti-Chagas drugs.     

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.     

Exploring the Linkage Dependence of Polyubiquitin Conformations Using Molecular Modeling

Posttranslational modification of proteins by covalent attachment of a small protein ubiquitin (Ub) or a polymeric chain of Ub molecules (called polyubiquitin) is involved in controlling a vast variety of processes in eukaryotic cells. The question of how different polyubiquitin signals are recognized is central to understanding the specificity of various types of polyubiquitination. In polyubiquitin, monomers are linked to each other via an isopeptide bond between the C-terminal glycine of one Ub and a lysine of the other. The functional outcome of polyubiquitination depends on the particular lysine involved in chain formation and appears to rely on linkage-dependent conformation of polyubiquitin. Thus, K48-linked chains, a universal signal for proteasomal degradation, under physiological conditions adopt a closed conformation where functionally important residues L8, I44, and V70 are sequestered at the interface between two adjacent Ub monomers. By contrast, K63-linked chains, which act as a nonproteolytic regulatory signal, adopt an extended conformation that lacks hydrophobic interubiquitin contact. Little is known about the functional roles of the so-called “noncanonical” chains (linked via K6, K11, K27, K29, or K33, or linked head-to-tail), and no structural information on these chains is available, except for information on the crystal structure of the head-to-tail-linked diubiquitin (Ub₂). In this study, we use molecular modeling to examine whether any of the noncanonical chains can adopt a closed conformation similar to that in K48-linked polyubiquitin. Our results show that the eight possible Ub₂ chains can be divided into two groups: chains linked via K6, K11, K27, or K48 are predicted to form a closed conformation, whereas chains linked via K29, K33, or K63, or linked head-to-tail are unable to form such a contact due to steric occlusion. These predictions are validated by the known structures of K48-, K63-, and head-to-tail-linked chains. Our study also predicts structural models for Ub₂ chains linked via K6, K11, or K27. The implications of these findings for linkage-selective recognition of noncanonical polyubiquitin signals by various receptors are discussed.     

Isolation of Metallothionein Genes and in silico Structural Characterization of Their Proteins Using Molecular Modeling from Yak (Bos grunniens)

Yak metallothioneins (BgMTs) are cysteine-rich metal-chelating proteins with highly conserved cysteine residues in their amino acid sequences. The 3D structures of the Cd(7)-BgMTs reconstructed by molecular modeling included two domains: the β-domain with M(3)(S(cys))(9) metal-thiolate clusters and the α-domain with M(4)(S(cys))(11) metal-thiolate clusters. An unusual variant was found at position 30 (Cys30→Ser30) in BgMT-III, which is usually conserved in the mammalian MT-I/-II (Cys29) and MT-III (Cys30). The variant residue of BgMT-III may play a key role in yak genetic evolution, metal-binding activity, dynamic conformation, and heavy metal metabolism. BgMT-III contained a Thr insertion at position 5 (T(5)), which may loosen the structure of the β-domain of BgMT-III, and a conserved C(6)PCP(9) motif, which may provide an interacting surface for protein-protein interactions. There is also an acidic hexapeptide insertion (E(55)GAEAE(60)) that could regulate the particular interdomain interactions and lead to the conformational change in the β-domain.     

Structural and Thermodynamic Analysis of HIV-1 Fusion Inhibition Using Small gp41 Mimetic Proteins

Development of effective inhibitors of the fusion between HIV-1 and the host cell membrane mediated by gp41 continues to be a grand challenge due to an incomplete understanding of the molecular and mechanistic details of the fusion process. We previously developed single-chain, chimeric proteins (named covNHR) that accurately mimic the N-heptad repeat (NHR) region of gp41 in a highly stable coiled-coil conformation. These molecules bind strongly to peptides derived from the gp41 C-heptad repeat (CHR) and are potent and broad HIV-1 inhibitors. Here, we investigated two covNHR variants differing in two mutations, V10E and Q123R (equivalent to V38E and Q40R in gp41 sequence) that reproduce the effect of HIV-1 mutations associated with resistance to fusion inhibitors, such as T20 (enfuvirtide). A detailed calorimetric analysis of the binding between the covNHR proteins and CHR peptides (C34 and T20) reveals drastic changes in affinity due to the mutations as a result of local changes in interactions at the site of T20 resistance. The crystallographic structure of the covNHR:C34 complex shows a virtually identical CHR–NHR binding interface to that of the post-fusion structure of gp41 and underlines an important role of buried interfacial water molecules in binding affinity and in development of resistance against CHR peptides. Despite the great difference in affinity, both covNHR variants demonstrate strong inhibitory activity for a wide variety of HIV-1 strains. These properties support the high potential of these covNHR proteins as new potent HIV-1 inhibitors. Our results may guide future inhibition approaches.     

脂肪乳剂理化性质和稳定性的影响因素浅析

静脉注射脂肪乳剂是由两种互不相溶的液体(油相和水相)混合而成,是水不溶性药物的良好栽体,可以提高易水解药物的稳定性并且减少各种毒副作用.脂肪乳荆属于热力学不稳定体系,其会表现出各种不稳定的现象如聚集,絮凝,聚结等,最终可能出现相分离.本文就处方、工艺及包装等影响脂肪乳剂理化性质及稳定性的主要因素进行综述.     

Modeling the Effect of Amino Acids and Copper on Monoclonal Antibody Productivity and Glycosylation: A Modular Approach

In manufacturing monoclonal antibodies (mAbs), it is crucial to be able to predict how process conditions and supplements affect productivity and quality attributes, especially glycosylation. Supplemental inputs, such as amino acids and trace metals in the media, are reported to affect cell metabolism and glycosylation; quantifying their effects is essential for effective process development. We aim to present and validate, through a commercially relevant cell culture process, a technique for modeling such effects efficiently. While existing models can predict mAb production or glycosylation dynamics under specific process configurations, adapting them to new processes remains challenging, because it involves modifying the model structure and often requires some mechanistic understanding. Here, a modular modeling technique for adapting an existing model for a fed-batch Chinese hamster ovary (CHO) cell culture process without structural modifications or mechanistic insight is presented. Instead, data is used, obtained from designed experimental perturbations in media supplementation, to train and validate a supplemental input effect model, which is used to “patch” the existing model. The combined model can be used for model-based process development to improve productivity and to meet product quality targets more efficiently. The methodology and analysis are generally applicable to other CHO cell lines and cell types.