Protein and Drug Design Algorithms Using Improved Biophysical Modeling

This thesis focuses on the development of algorithms that will allow protein design calculations to incorporate more realistic modeling assumptions. Protein design algorithms search large sequence spaces for protein sequences that are biologically and medically useful. Better modeling could improve the chance of success in designs and expand the range of problems to which these algorithms are applied. I have developed algorithms to improve modeling of backbone flexibility (DEEPer) and of more extensive continuous flexibility in general (EPIC and LUTE). I’ve also developed algorithms to perform multistate designs, which account for effects like specificity, with provable guarantees of accuracy (COMETS), and to accommodate a wider range of energy functions in design (EPIC and LUTE).

Using a cartographic modeling language to manipulate spectral satellite imagery

Land related information about the Earth's surface is commonIJ found in two forms: (1) map infornlation and (2) satellite image da ta. Satellite imagery provides a good visual picture of what is on the ground but complex image processing is required to interpret features in an image scene. Increasingly, methods are being sought to integrate the knowledge embodied in mop information into the interpretation task, or, alternatively, to bypass interpretation and perform biophysical modeling directly on derived data sources. A cartographic modeling language, as a generic map analysis package, is suggested as a means to integrate geographical knowledge and imagery in a process-oriented view of the Earth. Specialized cartographic models may be developed by users, which incorporate mapping information in performing land classification. In addition, a cartographic modeling language may be enhanced with operators suited to processing remotely sensed imagery. We demonstrate the usefulness of a cartographic modeling language for pre-processing satellite imagery, and define two nerv cartographic operators that evaluate image neighborhoods as post-processing operations to interpret thematic map values. The language and operators are demonstrated with an example image classification task.

A model of the dynamic relationship between blood flow and volume changes during brain activation

The temporal relationship between changes in cerebral blood flow (CBF) and cerebral blood volume (CBV) is important in the biophysical modeling and interpretation of the hemodynamic response to activation, particularly in the context of magnetic resonance imaging and the blood oxygen level-dependent signal. Grubb et al. (1974) measured the steady state relationship between changes in CBV and CBF after hypercapnic challenge. The relationship CBV proportional to CBFPhi has been used extensively in the literature. Two similar models, the Balloon (Buxton et al., 1998) and the Windkessel (Mandeville et al., 1999), have been proposed to describe the temporal dynamics of changes in CBV with respect to changes in CBF. In this study, a dynamic model extending the Windkessel model by incorporating delayed compliance is presented. The extended model is better able to capture the dynamics of CBV changes after changes in CBF, particularly in the return-to-baseline stages of the response.

Using Helix-coil Models to Study Protein Unfolded States

An abstract of a thesis devoted to using helix-coil models to study unfolded states.\\

Research on polypeptide unfolded states has received much more attention in the last decade or so than it has in the past. Unfolded states are thought to be implicated in various

misfolding diseases and likely play crucial roles in protein folding equilibria and folding rates. Structural characterization of unfolded states has proven to be

much more difficult than the now well established practice of determining the structures of folded proteins. This is largely because many core assumptions underlying

folded structure determination methods are invalid for unfolded states. This has led to a dearth of knowledge concerning the nature of unfolded state conformational

distributions. While many aspects of unfolded state structure are not well known, there does exist a significant body of work stretching back half a century that

has been focused on structural characterization of marginally stable polypeptide systems. This body of work represents an extensive collection of experimental

data and biophysical models associated with describing helix-coil equilibria in polypeptide systems. Much of the work on unfolded states in the last decade has not been devoted

specifically to the improvement of our understanding of helix-coil equilibria, which arguably is the most well characterized of the various conformational equilibria

that likely contribute to unfolded state conformational distributions. This thesis seeks to provide a deeper investigation of helix-coil equilibria using modern

statistical data analysis and biophysical modeling techniques. The studies contained within seek to provide deeper insights and new perspectives on what we presumably

know very well about protein unfolded states. \\

Chapter 1 gives an overview of recent and historical work on studying protein unfolded states. The study of helix-coil equilibria is placed in the context

of the general field of unfolded state research and the basics of helix-coil models are introduced.\\

Chapter 2 introduces the newest incarnation of a sophisticated helix-coil model. State of the art modern statistical techniques are employed to estimate the energies

of various physical interactions that serve to influence helix-coil equilibria. A new Bayesian model selection approach is utilized to test many long-standing

hypotheses concerning the physical nature of the helix-coil transition. Some assumptions made in previous models are shown to be invalid and the new model

exhibits greatly improved predictive performance relative to its predecessor. \\

Chapter 3 introduces a new statistical model that can be used to interpret amide exchange measurements. As amide exchange can serve as a probe for residue-specific

properties of helix-coil ensembles, the new model provides a novel and robust method to use these types of measurements to characterize helix-coil ensembles experimentally

and test the position-specific predictions of helix-coil models. The statistical model is shown to perform exceedingly better than the most commonly used

method for interpreting amide exchange data. The estimates of the model obtained from amide exchange measurements on an example helical peptide

also show a remarkable consistency with the predictions of the helix-coil model. \\

Chapter 4 involves a study of helix-coil ensembles through the enumeration of helix-coil configurations. Aside from providing new insights into helix-coil ensembles,

this chapter also introduces a new method by which helix-coil models can be extended to calculate new types of observables. Future work on this approach could potentially

allow helix-coil models to move into use domains that were previously inaccessible and reserved for other types of unfolded state models that were introduced in chapter 1.

Modeling the biophysical impacts of global change in mountain biosphere reserves

Mountains and mountain societies provide a wide range of goods and services to humanity, but they are particularly sensitive to the effects of global environmental change. Thus, the definition of appropriate management regimes that maintain the multiple functions of mountain regions in a time of greatly changing climatic, economic, and societal drivers constitutes a significant challenge. Management decisions must be based on a sound understanding of the future dynamics of these systems. The present article reviews the elements required for an integrated effort to project the impacts of global change on mountain regions, and recommends tools that can be used at 3 scientific levels (essential, improved, and optimum). The proposed strategy is evaluated with respect to UNESCO's network of Mountain Biosphere Reserves (MBRs), with the intention of implementing it in other mountain regions as well. First, methods for generating scenarios of key drivers of global change are reviewed, including land use/land cover and climate change. This is followed by a brief review of the models available for projecting the impacts of these scenarios on (1) cryospheric systems, (2) ecosystem structure and diversity, and (3) ecosystem functions such as carbon and water relations. Finally, the cross-cutting role of remote sensing techniques is evaluated with respect to both monitoring and modeling efforts. We conclude that a broad range of techniques is available for both scenario generation and impact assessments, many of which can be implemented without much capacity building across many or even most MBRs. However, to foster implementation of the proposed strategy, further efforts are required to establish partnerships between scientists and resource managers in mountain areas.

Biophysical 3D modeling approach of the Canary Current System: a ROMS-PISCES coupling study

Máster Universitario en Oceanografía

Computational modeling of prefrontal cortex circuits

Dissertation presented to obtain the Ph.D degree in Biology

Vascular disease modeling using induced pluripotent stem cells: Focus in Hutchinson-Gilford Progeria Syndrome

Transparency document related to this article can be found online at http://dx.doi.org/10.1016/j.bbrc.2015.10.014

Comment on recent modeling studies of astrocyte-neuron metabolic interactions.

Recent years have seen a surge in mathematical modeling of the various aspects of neuron-astrocyte interactions, and the field of brain energy metabolism is no exception in that regard. Despite the advent of biophysical models in the field, the long-lasting debate on the role of lactate in brain energy metabolism is still unresolved. Quite the contrary, it has been ported to the world of differential equations. Here, we summarize the present state of this discussion from the modeler's point of view and bring some crucial points to the attention of the non-mathematically proficient reader.

Molecular approaches for structural characterization of Bothrops L-amino acid oxidases with antiprotozoal activity: cDNA cloning, comparative sequence analysis, and molecular modeling

Two L-amino acid oxidases (LAAOs) were identified by random sequencing of cDNA libraries from the venom glands of Bothrops moojeni (BmooLAAO) and Bothrops jararacussu (Bjussu LAAO). Phylogenetic analysis involving other SV-LAAOs showed sequence identities within the range 83-87% being closely related to those from Agkistrodon and Trimeresurus. Molecular modeling experiments indicated the FAD-binding, substrate-binding, and helical domains of Bmoo and Bjussu LAAOs. The RMS deviations obtained by the superposition of those domains and that from Calloselasma rhodostoma LAAO crystal structure confirm the high degree of structural similarity between these enzymes. Purified BjussuLAAO-I and BmooLAAO-I exhibited antiprotozoal activities which were demonstrated to be hydrogen-peroxide mediated. This is the first report on the isolation and identification of cDNAs encoding LAAOs from Bothrops venom. The findings here reported contribute to the overall structural elucidation of SV-LAAOs and will advance the understanding on their mode of action. (c) 2006 Elsevier B.V. All rights reserved.

Modeling the spatial and temporal heterogeneity of deforestation-driven carbon emissions: the INPE-EM framework applied to the Brazilian Amazon

We present a generic spatially explicit modeling framework to estimate carbon emissions from deforestation (INPE-EM). The framework incorporates the temporal dynamics related to the deforestation process and accounts for the biophysical and socioeconomic heterogeneity of the region under study. We build an emission model for the Brazilian Amazon combining annual maps of new clearings, four maps of biomass, and a set of alternative parameters based on the recent literature. The most important results are as follows: (a) Using different biomass maps leads to large differences in estimates of emission; for the entire region of the Brazilian Amazon in the last decade, emission estimates of primary forest deforestation range from 0.21 to 0.26 similar to Pg similar to C similar to yr-1. (b) Secondary vegetation growth presents a small impact on emission balance because of the short duration of secondary vegetation. In average, the balance is only 5% smaller than the primary forest deforestation emissions. (c) Deforestation rates decreased significantly in the Brazilian Amazon in recent years, from 27 similar to Mkm2 in 2004 to 7 similar to Mkm2 in 2010. INPE-EM process-based estimates reflect this decrease even though the agricultural frontier is moving to areas of higher biomass. The decrease is slower than a non-process instantaneous model would estimate as it considers residual emissions (slash, wood products, and secondary vegetation). The average balance, considering all biomass, decreases from 0.28 in 2004 to 0.15 similar to Pg similar to C similar to yr-1 in 2009; the non-process model estimates a decrease from 0.33 to 0.10 similar to Pg similar to C similar to yr-1. We conclude that the INPE-EM is a powerful tool for representing deforestation-driven carbon emissions. Biomass estimates are still the largest source of uncertainty in the effective use of this type of model for informing mechanisms such as REDD+. The results also indicate that efforts to reduce emissions should focus not only on controlling primary forest deforestation but also on creating incentives for the restoration of secondary forests.

Molecular modeling and small angle X-ray scattering studies of Hoplosternum littorale cathodic haemoglobin

Considerable interest is currently focused on fish haemoglobins in order to identify the structural basis for their diversity of functional behavior. Hoplosternum littorale is a catfish that presents bimodal gill (water)/gut (air) -breathing, which allows this species to survive in waters with low oxygen content. The hemolysate of this fish showed the presence of two main haemoglobins, cathodic and anodic. This work describes structural features analyzed here by integration of molecular modeling with small angle X-ray scattering. Here is described a molecular model for the cathodic haemoglobin in the unliganded and liganded states. The models were determined by molecular modeling based on the high-resolution crystal structure of fish haemoglobins. The structural models for both forms of H. littorale haemoglobin were compared to human haemoglobin. (C) 2004 Elsevier B.V. All rights reserved.

Parmodel: a web server for automated comparative modeling of proteins

Parmodel is a web server for automated comparative modeling and evaluation of protein structures. The aim of this tool is to help inexperienced users to perform modeling, assessment, visualization, and optimization of protein models as well as crystallographers to evaluate structures solved experimentally. It is subdivided in four modules: Parmodel Modeling, Parmodel Assessment, Parmodel Visualization, and Parmodel Optimization. The main module is the Parmodel Modeling that allows the building of several models ford a same protein in a reduced time, through the distribution of modeling processes on a Beowulf cluster. Parmodel automates and integrates the main softwares used in comparative modeling as MODELLER, Whatcheck, Procheck, Raster3D, Molscript, and Gromacs. This web server is freely accessible at http://www.biocristalografia.df.ibilce.unesp.br/tools/parmodel. (C) 2004 Elsevier B.V. All rights reserved.