A Novel Domain Assembly Routine for Creating Full-Length Models of Membrane Proteins from Known Domain Structures

Key aspects of the past 30 years of protein design

Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.

Structural model of tissue factor (TF) and TF-factor VIIa complex in a lipid membrane: A combined experimental and computational study

Tissue factor (TF) is a membrane protein involved in blood coagulation. TF initiates a cascade of proteolytic reactions, ultimately leading to the formation of a blood clot. The first reaction consists of the binding of the coagulation factor VII and its conversion to the activated form, FVIIa. Here, we combined experimental, i.e. quartz crystal microbalance with dissipation monitoring and neutron reflectometry, and computational, i.e. molecular dynamics (MD) simulation, methods to derive a complete structural model of TF and TF/FVIIa complex in a lipid bilayer. This model shows that the TF transmembrane domain (TMD), and the flexible linker connecting the TMD to the extracellular domain (ECD), define the location of the ECD on the membrane surface. The average orientation of the ECD relative to the bilayer surface is slightly tilted towards the lipid headgroups, a conformation that we suggest is promoted by phosphatidylserine lipids, and favours the binding of FVIIa. On the other hand, the formation of the TF/FVIIa complex induces minor changes in the TF structure, and reduces the conformational freedom of both TF and FVIIA. Altogether we describe the protein-protein and protein-lipid interactions favouring blood coagulation, but also instrumental to the development of new drugs.

Ensuring scientific reproducibility in bio-macromolecular modeling via extensive, automated benchmarks

Each year vast international resources are wasted on irreproducible research. The scientific community has been slow to adopt standard software engineering practices, despite the increases in high-dimensional data, complexities of workflows, and computational environments. Here we show how scientific software applications can be created in a reproducible manner when simple design goals for reproducibility are met. We describe the implementation of a test server framework and 40 scientific benchmarks, covering numerous applications in Rosetta bio-macromolecular modeling. High performance computing cluster integration allows these benchmarks to run continuously and automatically. Detailed protocol captures are useful for developers and users of Rosetta and other macromolecular modeling tools. The framework and design concepts presented here are valuable for developers and users of any type of scientific software and for the scientific community to create reproducible methods. Specific examples highlight the utility of this framework, and the comprehensive documentation illustrates the ease of adding new tests in a matter of hours.

Co-occurring gain-of-function mutations in HER2 and HER3 modulate HER2/HER3 activation, oncogenesis, and HER2 inhibitor sensitivity

Activating mutations in HER2 (ERBB2) drive the growth of a subset of breast and other cancers and tend to co-occur with HER3 (ERBB3) missense mutations. The HER2 tyrosine kinase inhibitor neratinib has shown clinical activity against HER2-mutant tumors. To characterize the role of HER3 mutations in HER2-mutant tumors, we integrate computational structural modeling with biochemical and cell biological analyses. Computational modeling predicts that the frequent HER3E928G kinase domain mutation enhances the affinity of HER2/HER3 and reduces binding of HER2 to its inhibitor neratinib. Co-expression of mutant HER2/HER3 enhances HER2/HER3 co-immunoprecipitation and ligand-independent activation of HER2/HER3 and PI3K/AKT, resulting in enhanced growth, invasiveness, and resistance to HER2-targeted therapies, which can be reversed by combined treatment with PI3Kα inhibitors. Our results provide a mechanistic rationale for the evolutionary selection of co-occurring HER2/HER3 mutations and the recent clinical observations that HER3 mutations are associated with a poor response to neratinib in HER2-mutant cancers.

Order and disorder-An integrative structure of the full-length human growth hormone receptor

Because of its small size (70 kilodalton) and large content of structural disorder (>50%), the human growth hormone receptor (hGHR) falls between the cracks of conventional high-resolution structural biology methods. Here, we study the structure of the full-length hGHR in nanodiscs with small-angle x-ray scattering (SAXS) as the foundation. We develop an approach that combines SAXS, x-ray diffraction, and NMR spectroscopy data obtained on individual domains and integrate these through molecular dynamics simulations to interpret SAXS data on the full-length hGHR in nanodiscs. The hGHR domains reorient freely, resulting in a broad structural ensemble, emphasizing the need to take an ensemble view on signaling of relevance to disease states. The structure provides the first experimental model of any full-length cytokine receptor in a lipid membrane and exemplifies how integrating experimental data from several techniques computationally may access structures of membrane proteins with long, disordered regions, a widespread phenomenon in biology.

Domain-Swap Dimerization of Acanthamoeba castellanii CYP51 and a Unique Mechanism of Inactivation by Isavuconazole

Cytochromes P450 (P450, CYP) metabolize a wide variety of endogenous and exogenous lipophilic molecules, including most drugs. Sterol 14α-demethylase (CYP51) is a target for antifungal drugs known as conazoles. Using X-ray crystallography, we have discovered a domain-swap homodimerization mode in CYP51 from a human pathogen, Acanthamoeba castellanii CYP51 (AcCYP51). Recombinant AcCYP51 with a truncated transmembrane helix was purified as a heterogeneous mixture corresponding to the dimer and monomer units. Spectral analyses of these two populations have shown that the CO-bound ferrous form of the dimeric protein absorbed at 448 nm (catalytically competent form), whereas the monomeric form absorbed at 420 nm (catalytically incompetent form). AcCYP51 dimerized head-to-head via N-termini swapping, resulting in formation of a nonplanar protein-protein interface exceeding 2000 Ų with a total solvation energy gain of -35.4 kcal/mol. In the dimer, the protomers faced each other through the F and G α-helices, thus blocking the substrate access channel. In the presence of the drugs clotrimazole and isavuconazole, the AcCYP51 drug complexes crystallized as monomers. Although clotrimazole-bound AcCYP51 adopted a typical CYP monomer structure, isavuconazole-bound AcCYP51 failed to refold 74 N-terminal residues. The failure of AcCYP51 to fully refold upon inhibitor binding in vivo would cause an irreversible loss of a structurally aberrant enzyme through proteolytic degradation. This assumption explains the superior potency of isavuconazole against A. castellanii The dimerization mode observed in this work is compatible with membrane association and may be relevant to other members of the CYP family of biologic, medical, and pharmacological importance. SIGNIFICANCE STATEMENT: We investigated the mechanism of action of antifungal drugs in the human pathogen Acanthamoeba castellanii. We discovered that the enzyme target [Acanthamoeba castellanii sterol 14α-demethylase (AcCYP51)] formed a dimer via an N-termini swap, whereas drug-bound AcCYP51 was monomeric. In the AcCYP51-isavuconazole complex, the protein target failed to refold 74 N-terminal residues, suggesting a fundamentally different mechanism of AcCYP51 inactivation than only blocking the active site. Proteolytic degradation of a structurally aberrant enzyme would explain the superior potency of isavuconazole against A. castellanii.

Macromolecular modeling and design in Rosetta: recent methods and frameworks

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.