Rafael C. Bernardi
During the PhD work, Rafael worked mainly with molecular dynamics studies of anesthetics and their effect in both biological membranes and the TREK-1 potassium channel. This work awarded him the Best Thesis Award in Biophysics and Biotechnology in 2010 by the Secretariat of Strategic Affair of the Brazilian Presidency. He is also co-worker in several studies simulating biological membranes, and he has been responsible for the first QM/MM dynamics of lipid membranes in which some lipids were studied in the DFT level. For the last couple of years, he had been working with molecular modeling of cellulose, cellulases and cellulosomes, first at INMETRO (Brazil), in a group that aims to develop a new second-generation biofuel, and more recently, at TCBG.
Home Department: Beckman Institute for Advanced Science and Technology.
Office Address: Beckman Institute, Room 3157.
Office Phone: (217) 244-0177
- PhD in Biophysics (2010) - Federal University of Rio de Janeiro - Brazil
- Masters in Physics (2007) - Brazilian Center for Physical Research - Brazil
- BS Physics (2005) - State University of Londrina - Brazil
Advisor: Dr. Pedro Geraldo Pascutti and Dr. Werner Treptow
Advisor: Dr. Carlton Anthony Taft
Advisor: Dr. André Tsutomu Ota
- Second-generation Biofuel
- Cellulose, Cellulases and Cellulosomes
- Local Anesthetics
- Biological Membranes
- Protein Folding
QwikMD - Gateway to Easy Simulation
To assist experimentalists and any novice to MD to overcome the initial learning currve barrier of MD simulation software, we developed QwikMD, a user interface that connects the widely employed and user-friendly molecular graphics program VMD to the widely adopted MD program NAMD.
By enabling easy control of every step, QwikMD, meets even the needs of experts in the field, increasing the efficiency and quality of their work.
Easy Setup of MD Simulations
Changes in Protonation
Protocols for MD and SMD
Live View Simulations
Integrated Basic Analysis
Info Buttons to Guide Novices
Available on Amazon Cloud
Full log Capabilities
Biofuels are one of the most studied alternative sources of fossil-derived fuels, especially due to its smaller effect on greenhouse gas accumulation. The interest in the so-called second-generation biofuels has increased in the last few years, since they can be produced using the waste of agricultural production and thus do not raise ethical issues. Enzymatic hydrolysis to release carbohydrates units from plant cell wall polysaccharides has been one of the most promising strategies to produce the aforementioned new renewable biofuels. However the natural resistance of the plant biomass, usually referred to as biomass recalcitrance, makes the industrial process not cost-effective. To address this problem we are performing a joint computational and experimental study to unravel the mechanism of biomass degradation by cellulases, the enzymes found in some microbes that can digest plant fibers.
In another joint study we expect to elucidate the cooperative mechanism of cellulases complexes known as cellulosomes. The plant cell wall polysaccharides, cellulose and xylan represent the most abundant reservoirs of fixed carbon in nature and thus have considerable potential as a renewable energy source. These polysaccharides are highly recalcitrant to degradation; however fungi and bacteria have evolved complex enzymatic strategies for the degradation and fermentation of these polysaccharides. One of the important paradigms for enzymatic degradation of plant cell walls is the cellulosomal system, elaborated by certain strains of bacteria. These extracellular proteinaceous assemblages consist of a scaffoldin backbone onto which cohesin-containing catalytic modules and carbohydrate binding modules are appended. Analogous to a Swiss army knife, these cellulosomes contain a plethora of different catalytic and substrate binding activities that facilitate the degradation of plant cell wall material. It was reported that Clostridium thermocellum, the most studied cellulosomal organism, exhibits one of the highest rates of cellulose utilization known in nature, and the cellulosomal system of this bacterium is reported to display a specific activity against crystalline cellulose that is 50-fold higher than the corresponding non-cellusomal fungal system in Trichoderma reesei. To understand the efficient mechanism that leads cellulosomal machineries in organisms to degrade biomass, we are studying the interaction of the subunits of the cellulosome to comprehend their synergistic mechanism.
One of the Strongest Biomolecular Interactions
Studying membrane lipids in the DFT level using a QM/MM approach
From very simple models, such as dielectric interfaces in the 1980s and 1990s, to all-atoms simulations with thousands of lipids, going through several methods and models developed, biological membranes as well as their proteins have become one of the most studied systems by molecular dynamics simulations, especially because of the huge advances in the computer capability. Membranes are especially important to drug development, since several membrane proteins are drug targets. Moreover, the membrane is a remarkable barrier to a majority of the drugs that act inside the cell.
Numerous studies attempt to understand how could we better use computational methods to study such regions, for instance, using more accurate methods of molecular modeling, such as quantum mechanical (QM) calculations. Meanwhile, because this interface has electrostatic properties that, at first, could not be studied by classical molecular modeling, the application of QM methods seems to be a more accurate approach. Furthermore, the behavior of the charge distribution over a drug surface, especially those containing aromatic groups, will adapt to the media when crossing a membrane interface, opposite of the classical MD approach, where it remains constant.
To study biological systems, especially enzymatic reactions, a hybrid QM/MM approach in molecular dynamics have been largely used. However, despite its importance in several biochemical processes, to the best of our knowledge, the very same approach has not been used to study the lipid membranes interface using a quantum mechanical treatment for the lipids. Understand how amphiphilic drugs interact, at the level of molecular orbitals, with the lipids in this region is essential for the complete knowledge of such drugs mechanism.
We choose to study how the local anesthetic benzocaine behaves in such region, once its mechanism of action may be related to its stabilization on the membrane/water interface. The QM/MM simulations were carried out taking advantage of an integration of the Car-Parrinello method, to treat the benzocaine as well as two of the membrane lipids and the GROMOS force field to simulate the remaining system. Our calculations showed an overlap of the atomic orbitals of the anesthetic and the lipid, suggesting a strong hydrogen bond between the amine terminal of the drug and the palmitate group of the lipid.
One of the main paradigms of molecular biology tells us that the three-dimensional structure of proteins defines its function and dynamics. Such three-dimensional structures, in turn, derive from the amino acid sequence itself, through the folding process. A protein's structure determines its activity and properties, thus knowing such conformation on an atomic level is essential for both basic and applied studies of protein function and dynamics. However, the acquisition of such structures by experimental methods is slow and expensive, and current computational methods mostly depend on previously known structures to determine new ones. We developed a new software called GSAFold that applies the Generalized Simulated Annealing (GSA) algorithm on ab-initio protein structure prediction. The GSA is a stochastic search algorithm employed in energy minimization and used in global optimization problems, especially those that depend on long-range interactions, such as gravity models and conformation optimization of small molecules. This new implementation applies, for the first time in ab-initio protein structure prediction, an analytical inverse for the Visitation function of GSA. It also employs the broadly used NAMD Molecular Dynamics package to carry out energy calculations, allowing the user to select different force fields and parameterizations. Moreover, the software also allows the execution of several simulations simultaneously. Applications that depend on protein structures include rational drug design and structure-based protein function prediction.
- Energy Biosciences Institute - University of Illinois and University of California Berkeley
- Gaub Lab - Ludwig-Maximilians Universitat Munchen - LMU - Germany
- Bayer Group - Weizmann Institute of Science - Israel
- Laboratory for Molecular Modeling and Dynamics - Federal University of Rio de Janeiro - UFRJ - Brazil
- Laboratory of Biotechnology - National Institute of Metrology, Quality and Technology - INMETRO - Brazil
24. QwikMD Integrative Molecular Dynamics Toolkit for Novices and Experts; JV Ribeiro, RC Bernardi, T Rudack, JE Stone ,JC Phillips,PL Freddolino, K Schulten; Scientific Reports, 2016
23. Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes; BC Goh, JA Hadden, RC Bernardi, A Singharoy, R McGreevy, T Rudack, CK Cassidy, K Schulten; Annual Review of Biophysics, 2016
22. Easy and Fast Setup of Molecular Dynamics Simulations: Combining VMD and NAMD for Experimentalists; JV Ribeiro, RC Bernardi, T Rudack, K Schulten; Biophysical Journal, 2016
21. Cellulose degradation in the human gut: Ruminococcus champanellensis expands the cellulosome paradigm; I Cann, RC Bernardi, RI Mackie; Environmental Microbiology, 2016
20. Mapping mechanical force propagation through biomolecular complexes; C Schoeler, RC Bernardi, KH Malinowska, E Durner, W Ott, EA Bayer, K Schulten, MA Nash, HE Gaub ; Nano Letters, 2015
19. Enhanced sampling techniques in molecular dynamics simulations of biological systems; RC Bernardi, MCR Melo, K Schulten; Biochimica et Biophysica Acta (BBA), 2015
18. Molecular dynamics simulations of large macromolecular complexes; JR Perilla, BC Goh, CK Cassidy, B Liu, RC Bernardi, T Rudack, H Yu, Z Wu, K Schulten; Current Opinion in Structural Biology , 2015
17. Mechanism of Inhibition of Glycoside Hydrolases Investigated by Molecular Dynamics Simulations; RC Bernardi, I Cann, E Imsand, D Clark, K Schulten; Biophysical Journal, 2015
16. Ultrastable cellulosome-adhesion complex tightens under load; C Schoeler, KH Malinowska, RC Bernardi, LF Milles, MA Jobst, E Durner, W Ott, DB Fried, EA Bayer, K Schulten, HE Gaub, MA Nash; Nature Communications, 2014
15. Molecular dynamics study of enhanced Man5B enzymatic activity; RC Bernardi, I Cann, K Schulten; Biotechnology for Biofuels, 2014
14. Large Scale Structure Sampling for Protein Fold Prediction using the Generalized Simulated Annealing; MCR Melo, RC Bernardi, PG Pascutti; Biophysical Journal, 2013
13. Molecular Dynamics Studies of Buckminsterfullerene Derivatives as Drug Carriers; RC Bernardi, K Schulten, PG Pascutti; Biophysical Journal, 2013
12. The Structural Dynamics of the Flavivirus Fusion Peptide-Membrane Interaction; YS Mendes, NS Alves, TLF Souza, IP Sousa, ML Bianconi, RC Bernardi, PG, et. al.; PloS one, 2012
11. GSAFold: A new application of GSA to protein structure prediction; MCR Melo, RC Bernardi, TVA Fernandes, PG Pascutti; Proteins: Structure, Function, and Bioinformatics, 2012
10. Hybrid QM/MM Molecular Dynamics Study of Benzocaine in a Membrane Environment: How Does a Quantum Mechanical Treatment of Both Anesthetic and Lipids Affect Their Interaction; RC Bernardi, PG Pascutti; Journal of Chemical Theory and Computation, 2012
9. The role of helices 5 and 6 on the human beta1-adrenoceptor activation mechanism; LVB Hoelz, AAST Ribeiro, RC Bernardi, BAC Horta, MG Albuquerque, et. al.; Molecular Simulation, 2012
8. QM/MM Molecular Dynamics Methods Applied to Investigate Cellulose Fibers Hydration; RC Bernardi, MCR Melo, PG Pascutti, Biophysical Journal, 2012
7. New Developments on Generalized Simulated Annealing Applied to ab-initio Protein Structure Prediction; MC Melo, TV Fernandes, RC Bernardi, PG Pascutti; Biophysical Journal, 2012
6. Dynamical behaviour of the human beta1-adrenoceptor under agonist binding; LVB Hoelz, RC Bernardi, BAC Horta, JQ Araújo, MG Albuquerque, et. al.; Molecular Simulation, 2011
5. The Ionic Lock Activation Mechanism of the Human B1-adrenoceptor; LVB Hoelz, RC Bernardi, MG Albuquerque, JFM Silva; Drugs of the Future, 2010
4. Molecular dynamics study of biomembrane/local anesthetics interactions; RC Bernardi, DEB Gomes, R Gobato, CA Taft, AT Ota, PG Pascutti; Molecular Physics 2009
3. Density functional and molecular dynamics simulations of local anesthetics in 0.9% NaCl solution; RC Bernardi, DEB Gomes, AS Ito, AT Ota, PG Pascutti, C Taft; Molecular Simulation, 2007
2. Water solvent and local anesthetics: A computational study; RC Bernardi, DEB Gomes, PG Pascutti, AS Ito, CA Taft, AT Ota; International Journal of Quantum Chemistry, 2007
1. Theoretical studies on water-tetracaine interaction; RC Bernardi, DEB Gomes, PG Pascutti, AS Ito, AT Ota; International journal of quantum chemistry, 2006