Climate–ecosystem modelling made easy: The Land Sites Platform

Dynamic Global Vegetation Models (DGVMs) provide a state-of-the-art process-based approach to study the complex interplay between vegetation and its physical environment. For example, they help to predict how terrestrial plants interact with climate, soils, disturbance and competition for resources. We argue that there is untapped potential for the use of DGVMs in ecological and ecophysiological research. One fundamental barrier to realize this potential is that many researchers with relevant expertize (ecology, plant physiology, soil science, etc.) lack access to the technical resources or awareness of the research potential of DGVMs. Here we present the Land Sites Platform (LSP): new software that facilitates single-site simulations with the Functionally Assembled Terrestrial Ecosystem Simulator, an advanced DGVM coupled with the Community Land Model. The LSP includes a Graphical User Interface and an Application Programming Interface, which improve the user experience and lower the technical thresholds for installing these model architectures and setting up model experiments. The software is distributed via version-controlled containers; researchers and students can run simulations directly on their personal computers or servers, with relatively low hardware requirements, and on different operating systems. Version 1.0 of the LSP supports site-level simulations. We provide input data for 20 established geo-ecological observation sites in Norway and workflows to add generic sites from public global datasets. The LSP makes standard model experiments with default data easily achievable (e.g., for educational or introductory purposes) while retaining flexibility for more advanced scientific uses. We further provide tools to visualize the model input and output, including simple examples to relate predictions to local observations. The LSP improves access to land surface and DGVM modelling as a building block of community cyberinfrastructure that may inspire new avenues for mechanistic ecosystem research across disciplines.     

Transient Exposure of Pulmonary Surfactant to Hyaluronan Promotes Structural and Compositional Transformations into a Highly Active State

Pulmonary surfactant is a lipid-protein complex that lowers surface tension at the respiratory air-liquid interface, stabilizing the lungs against physical forces tending to collapse alveoli. Dysfunction of surfactant is associated with respiratory pathologies such as acute respiratory distress syndrome or meconium aspiration syndrome where naturally occurring surfactant-inhibitory agents such as serum, meconium, or cholesterol reach the lung. We analyzed the effect of hyaluronan (HA) on the structure and surface behavior of pulmonary surfactant to understand the mechanism for HA-promoted surfactant protection in the presence of inhibitory agents. In particular, we found that HA affects structural properties such as the aggregation state of surfactant membranes and the size, distribution, and order/packing of phase-segregated lipid domains. These effects do not require a direct interaction between surfactant complexes and HA and are accompanied by a compositional reorganization of large surfactant complexes that become enriched with saturated phospholipid species. HA-exposed surfactant reaches very high efficiency in terms of rapid and spontaneous adsorption of surfactant phospholipids at the air-liquid interface and shows significantly improved resistance to inactivation by serum or cholesterol. We propose that physical effects pertaining to the formation of a meshwork of interpenetrating HA polymer chains are responsible for the changes in surfactant structure and composition that enhance surfactant function and, thus, resistance to inactivation. The higher resistance of HA-exposed surfactant to inactivation persists even after removal of the polymer, suggesting that transient exposure of surfactant to polymers like HA could be a promising strategy for the production of more efficient therapeutic surfactant preparations.     

Development of a novel MATLAB-based framework for implementing mechanical joint stability constraints within OpenSim musculoskeletal models

The Static Optimization (SO) solver in OpenSim estimates muscle activations and forces that only equilibrate applied moments. In this study, SO was enhanced through an open-access MATLAB interface, where calculated muscle activations can additionally satisfy crucial mechanical stability requirements. This Stability-Constrained SO (SCSO) is applicable to many OpenSim models and can potentially produce more biofidelic results than SO alone, especially when antagonistic muscle co-contraction is required to stabilize body joints. This hypothesis was tested using existing models and experimental data in the literature. Muscle activations were calculated by SO and SCSO for a spine model during two series of static trials (i.e. simulation 1 and 2), and also for a lower limb model (supplementary material 2). In simulation 1, symmetric and asymmetric flexion postures were compared, while in simulation 2, various external load heights were compared, where increases in load height did not change the external lumbar flexion moment, but necessitated higher EMG activations. During the tasks in simulation 1, the predicted muscle activations by SCSO demonstrated less average deviation from the EMG data (6.8% −7.5%) compared to those from SO (10.2%). In simulation 2, SO predicts constant muscle activations and forces, while SCSO predicts increases in the average activations of back and abdominal muscles that better match experimental data. Although the SCSO results are sensitive to some parameters (e.g. musculotendon stiffness), when considering the strategy of the central nervous system in distributing muscle forces and in activating antagonistic muscles, the assigned activations by SCSO are more biofidelic than SO.     

Partner-specific prediction of RNA-binding residues in proteins: A critical assessment

RNA-protein interactions play essential roles in regulating gene expression. While some RNA-protein interactions are “specific”, that is, the RNA-binding proteins preferentially bind to particular RNA sequence or structural motifs, others are “non-RNA specific.” Deciphering the protein-RNA recognition code is essential for comprehending the functional implications of these interactions and for developing new therapies for many diseases. Because of the high cost of experimental determination of protein-RNA interfaces, there is a need for computational methods to identify RNA-binding residues in proteins. While most of the existing computational methods for predicting RNA-binding residues in RNA-binding proteins are oblivious to the characteristics of the partner RNA, there is growing interest in methods for partner-specific prediction of RNA binding sites in proteins. In this work, we assess the performance of two recently published partner-specific protein-RNA interface prediction tools, PS-PRIP, and PRIdictor, along with our own new tools. Specifically, we introduce a novel metric, RNA-specificity metric (RSM), for quantifying the RNA-specificity of the RNA binding residues predicted by such tools. Our results show that the RNA-binding residues predicted by previously published methods are oblivious to the characteristics of the putative RNA binding partner. Moreover, when evaluated using partner-agnostic metrics, RNA partner-specific methods are outperformed by the state-of-the-art partner-agnostic methods. We conjecture that either (a) the protein-RNA complexes in PDB are not representative of the protein-RNA interactions in nature, or (b) the current methods for partner-specific prediction of RNA-binding residues in proteins fail to account for the differences in RNA partner-specific versus partner-agnostic protein-RNA interactions, or both.     

Highly specific inactivation of triosephosphate isomerase from Trypanosoma cruzi

We searched for molecules that selectively inactivate homodimeric triosephosphate isomerase from Trypanosoma cruzi (TcTIM), the parasite that causes Chagas' disease. We found that some benzothiazoles inactivate the enzyme. The most potent were 3-(2-benzothiazolylthio)-propanesulfonic acid, 2-(p-aminophenyl)-6-methylbenzothiazole-7-sulfonic acid, and 2-(2-4(4-aminophenyl)benzothiazole-6-methylbenzothiazole-7-sulfonic acid. Half-maximal inactivation by these compounds was attained with 33, 56, and 8 microM, respectively; in human TIM, half-maximal inactivation required 422 microM, 3.3 mM, and 1.6 mM. In TcTIM, the effect of the benzothiazoles decreased as the concentration of the enzyme was increased. TcTIM has a cysteine (Cys 15) at the dimer interface, whereas human TIM has methionine in that position. In M15C human TIM, the benzothiazole concentrations that caused half-maximal inactivation were much lower than in the wild type. The overall findings suggest that the benzothiazoles perturb the interactions between the two subunits of TcTIM through a process in which the interface cysteine is central in their deleterious action.     

A Mutation within the Extended X Loop Abolished Substrate-induced ATPase Activity of the Human Liver ATP-binding Cassette (ABC) Transporter MDR3

The human multidrug resistance protein 3 (MDR3/ABCB4) belongs to the ubiquitous family of ATP-binding cassette (ABC) transporters and is located in the canalicular membrane of hepatocytes. There it flops the phospholipids of the phosphatidylcholine (PC) family from the inner to the outer leaflet. Here, we report the characterization of wild type MDR3 and the Q1174E mutant, which was identified previously in a patient with progressive familial intrahepatic cholestasis type 3 (PFIC-3). We expressed different variants of MDR3 in the yeast Pichia pastoris, purified the proteins via tandem affinity chromatography, and determined MDR3-specific ATPase activity in the presence or absence of phospholipids. The ATPase activity of wild type MDR3 was stimulated 2-fold by liver PC or 1,2-dioleoyl-sn-glycero-3-phosphatidylethanolamine lipids. Furthermore, the cross-linking of MDR3 with a thiol-reactive fluorophore blocked ATP hydrolysis and exhibited no PC stimulation. Similarly, phosphatidylethanolamine, phosphatidylserine, and sphingomyelin lipids did not induce an increase of wild type MDR3 ATPase activity. The phosphate analogues beryllium fluoride and aluminum fluoride led to complete inhibition of ATPase activity, whereas orthovanadate inhibited exclusively the PC-stimulated ATPase activity of MDR3. The Q1174E mutation is located in the nucleotide-binding domain in direct proximity of the leucine of the ABC signature motif and extended the X loop, which is found in ABC exporters. Our data on the Q1174E mutant demonstrated basal ATPase activity, but PC lipids were incapable of stimulating ATPase activity highlighting the role of the extended X loop in the cross-talk of the nucleotide-binding domain and the transmembrane domain.     

A three-dimensional vertex dynamics cell model of space-filling polyhedra simulating cell behavior in a cell aggregate

We developed a three-dimensional (3D) cell model of a multicellular aggregate consisting of several polyhedral cells to investigate the deformation and rearrangement of cells under the influence of external forces. The polyhedral cells fill the space in the aggregate without gaps or overlaps, consist of contracting interfaces and maintain their volumes. The interfaces and volumes were expressed by 3D vertex coordinates. Vertex movements obey equations of motion that rearrange the cells to minimize total free energy, and undergo an elementary process that exchanges vertex pair connections when vertices approach each other. The total free energy includes the interface energy of cells and the compression or expansion energy of cells. Computer simulations provided the following results: An aggregate of cells becomes spherical to minimize individual cell surface areas; Polygonal interfaces of cells remain flat; Cells within the 3D cell aggregate can move and rearrange despite the absence of free space. We examined cell rearrangement to elucidate the viscoelastic properties of the aggregate, e.g. when an external force flattens a cell aggregate (e.g. under centrifugation) its component cells quickly flatten. Under a continuous external force, the cells slowly rearrange to recover their original shape although the cell aggregate remains flat. The deformation and rearrangement of individual cells is a two-step process with a time lag. Our results showed that morphological and viscoelastic properties of the cell aggregate with long relaxation time are based on component cells where minimization of interfacial energy of cells provides a motive force for cell movement.     

Enhancing Antibody Fc Heterodimer Formation through Electrostatic Steering Effects: APPLICATIONS TO BISPECIFIC MOLECULES AND MONOVALENT IgG

Naturally occurring IgG antibodies are bivalent and monospecific. Bispecific antibodies having binding specificities for two different antigens can be produced using recombinant technologies and are projected to have broad clinical applications. However, co-expression of multiple light and heavy chains often leads to contaminants and pose purification challenges. In this work, we have modified the CH3 domain interface of the antibody Fc region with selected mutations so that the engineered Fc proteins preferentially form heterodimers. These novel mutations create altered charge polarity across the Fc dimer interface such that coexpression of electrostatically matched Fc chains support favorable attractive interactions thereby promoting desired Fc heterodimer formation, whereas unfavorable repulsive charge interactions suppress unwanted Fc homodimer formation. This new Fc heterodimer format was used to produce bispecific single chain antibody fusions and monovalent IgGs with minimal homodimer contaminants. The strategy proposed here demonstrates the feasibility of robust production of novel Fc-based heterodimeric molecules and hence broadens the scope of bispecific molecules for therapeutic applications.     

Characterization of Protein–Protein Interfaces

We analyze the characteristics of protein–protein interfaces using the largest datasets available from the Protein Data Bank (PDB). We start with a comparison of interfaces with protein cores and non-interface surfaces. The results show that interfaces differ from protein cores and non-interface surfaces in residue composition, sequence entropy, and secondary structure. Since interfaces, protein cores, and non-interface surfaces have different solvent accessibilities, it is important to investigate whether the observed differences are due to the differences in solvent accessibility or differences in functionality. We separate out the effect of solvent accessibility by comparing interfaces with a set of residues having the same solvent accessibility as the interfaces. This strategy reveals residue distribution propensities that are not observable by comparing interfaces with protein cores and non-interface surfaces. Our conclusions are that there are larger numbers of hydrophobic residues, particularly aromatic residues, in interfaces, and the interactions apparently favored in interfaces include the opposite charge pairs and hydrophobic pairs. Surprisingly, Pro-Trp pairs are over represented in interfaces, presumably because of favorable geometries. The analysis is repeated using three datasets having different constraints on sequence similarity and structure quality. Consistent results are obtained across these datasets. We have also investigated separately the characteristics of heteromeric interfaces and homomeric interfaces.